From a8d591738909c2156b3f2cb9d81fdfaa4e5fa86c Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 11 Mar 2015 18:06:04 +0000
Subject: [PATCH] TomiyamaLift: Updated coefficient to make the model
continuous at EoH = 10.7 Correction provided by Juho Peltola
---
.../liftModels/TomiyamaLift/TomiyamaLift.C | 5 ++++-
.../liftModels/TomiyamaLift/TomiyamaLift.H | 6 +++++-
2 files changed, 9 insertions(+), 2 deletions(-)
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.C
index 310f2f30cd2..d4a4dfd17e5 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::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)
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.H b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/TomiyamaLift/TomiyamaLift.H
index 881d661bd3b..c179558dde9 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
@@ -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
--
GitLab