From acb63f5141d269247b3267229395c630d5e53b14 Mon Sep 17 00:00:00 2001
From: william <william>
Date: Tue, 18 Mar 2014 15:03:41 +0000
Subject: [PATCH] BUG: twoPhaseEulerFoam: corrected the definintion of K in
turbulence models and in the segregated drag model
---
.../interfacialModels/dragModels/segregated/segregated.C | 6 +-----
.../phaseIncompressible/LES/Niceno/NicenoKEqn.C | 2 +-
.../phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C | 3 ++-
.../RAS/mixtureKEpsilon/mixtureKEpsilon.C | 8 ++++----
4 files changed, 8 insertions(+), 11 deletions(-)
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/segregated/segregated.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/segregated/segregated.C
index 92b82e537d2..7ae380da6dc 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/segregated/segregated.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/segregated/segregated.C
@@ -138,11 +138,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
(
pair_.rho()
*pair_.magUr()
- /(
- magGradI
- *max(alpha1*alpha2, sqr(residualAlpha_))
- *muI
- )
+ /(magGradI*muI)
);
volScalarField lambda(m_*ReI + n_*muAlphaI/muI);
diff --git a/src/TurbulenceModels/phaseIncompressible/LES/Niceno/NicenoKEqn.C b/src/TurbulenceModels/phaseIncompressible/LES/Niceno/NicenoKEqn.C
index 07002a6a1ec..3035699024b 100644
--- a/src/TurbulenceModels/phaseIncompressible/LES/Niceno/NicenoKEqn.C
+++ b/src/TurbulenceModels/phaseIncompressible/LES/Niceno/NicenoKEqn.C
@@ -182,7 +182,7 @@ tmp<volScalarField> NicenoKEqn<BasicTurbulenceModel>::bubbleG() const
tmp<volScalarField> bubbleG
(
- Cp_*gas*sqr(magUr)*fluid.drag(gas).K()/liquid.rho()
+ Cp_*sqr(magUr)*fluid.drag(gas).K()/liquid.rho()
);
return bubbleG;
diff --git a/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C b/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
index 8454103e3b8..83edf7de157 100644
--- a/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
+++ b/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
@@ -194,9 +194,10 @@ tmp<volScalarField> LaheyKEpsilon<BasicTurbulenceModel>::bubbleG() const
tmp<volScalarField> bubbleG
(
Cp_
+ *liquid*liquid.rho()
*(
pow3(magUr)
- + pow(fluid.drag(gas).K()*gas.d()/liquid.rho(), 4.0/3.0)
+ + pow(fluid.drag(gas).CdRe()*liquid.nu()/gas.d(), 4.0/3.0)
*pow(magUr, 5.0/3.0)
)
*gas
diff --git a/src/TurbulenceModels/phaseIncompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.C b/src/TurbulenceModels/phaseIncompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.C
index 23d69508bda..52cc2ecadfe 100644
--- a/src/TurbulenceModels/phaseIncompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.C
+++ b/src/TurbulenceModels/phaseIncompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.C
@@ -380,7 +380,7 @@ tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::Ct2() const
volScalarField beta
(
(6*this->Cmu_/(4*sqrt(3.0/2.0)))
- *alphag*fluid.drag(gas).K()/liquid.rho()
+ *fluid.drag(gas).K()/liquid.rho()
*(liquidTurbulence.k_/liquidTurbulence.epsilon_)
);
volScalarField Ct0((3 + beta)/(1 + beta + 2*gas.rho()/liquid.rho()));
@@ -488,10 +488,10 @@ tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::bubbleG() const
tmp<volScalarField> bubbleG
(
Cp_
- *sqr(liquid)*liquid.rho()
+ *liquid*liquid.rho()
*(
pow3(magUr)
- + pow(fluid.drag(gas).K()*gas.d()/liquid.rho(), 4.0/3.0)
+ + pow(fluid.drag(gas).CdRe()*liquid.nu()/gas.d(), 4.0/3.0)
*pow(magUr, 5.0/3.0)
)
*gas
@@ -501,7 +501,7 @@ tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::bubbleG() const
// Simple model
// tmp<volScalarField> bubbleG
// (
- // Cp_*sqr(liquid)*gas*fluid.drag(gas).K()*sqr(magUr)
+ // Cp_*liquid*fluid.drag(gas).K()*sqr(magUr)
// );
return bubbleG;
--
GitLab