diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C index 310f2f30..d4a4dfd1 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C @@ -63,7 +63,10 @@ Foam::tmp 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) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.H b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.H index f6787a62..c179558d 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.H +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.H @@ -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