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

Correction provided by Juho Peltola

applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C

@@ -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)

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
@@ -35,6 +35,10 @@ Description
         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