research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

TomiyamaLift: Updated coefficient to make the model continuous at EoH = 10.7

Browse source 
Correction provided by Juho Peltola
This commit is contained in:
Henry 2015-03-11 18:09:12 +00:00
parent 883438d463
commit 96b1bb0b27
2 changed files with 11 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C
View file

@ -63,7 +63,10 @@ Foam::tmp<Foam::volScalarField> Foam::liftModels::TomiyamaLift::Cl() const
{
volScalarField EoH(pair_.EoH2());
volScalarField f(0.00105*pow3(EoH) - 0.0159*sqr(EoH) - 0.0204*EoH + 0.474);
volScalarField f
(
0.0010422*pow3(EoH) - 0.0159*sqr(EoH) - 0.0204*EoH + 0.474
);
return
neg(EoH - scalar(4))*min(0.288*tanh(0.121*pair_.Re()), f)

10
applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.H
View file

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2014 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2014-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -30,11 +30,15 @@ Description
Reference:
\verbatim
"Transverse migration of single bubbles in simple shear flows"
A Tomiyama, H Tamai, I Zun, S Hosokawa
Tomiyama, A., Tamai, H., Zun, I., Hosokawa, S.,
Chemical Engineering Science
Volume 57, Issue 11, Pages 1849-1858, June 2002
Volume 57, Issue 11, June 2002, pp. 1849-1858
\endverbatim
The coefficient for pow3(EoH) proposed by Tomiyama (2002) has been modified
to make the model continuous at EoH = 10.7 while maintaining the
lift coefficient proposed by Tomiyama (2002) when EoH >= 10.7.
SourceFiles
TomiyamaLift.C