From a80d9c7407b58af2bfc3ac0a5056e265996f5d77 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Tue, 10 Mar 2015 19:57:42 +0000
Subject: [PATCH] TomiyamaLift: Corrected coefficients to correspond exactly to
the original paper cited Resolved bug-report
http://www.openfoam.org/mantisbt/view.php?id=1564
---
.../liftModels/TomiyamaLift/TomiyamaLift.C | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C
index a789fd20791..310f2f30cd2 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C
@@ -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
@@ -67,8 +67,8 @@ Foam::tmp<Foam::volScalarField> Foam::liftModels::TomiyamaLift::Cl() const
return
neg(EoH - scalar(4))*min(0.288*tanh(0.121*pair_.Re()), f)
- + pos(EoH - scalar(4))*neg(EoH - scalar(10))*f
- + pos(EoH - scalar(10))*(-0.29);
+ + pos(EoH - scalar(4))*neg(EoH - scalar(10.7))*f
+ + pos(EoH - scalar(10.7))*(-0.288);
}
--
GitLab