zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: twoPhaseEulerFoam: fixed definition of drag in kineticTheory turbulence model

Browse Source
This commit is contained in:
william
2014-04-01 09:12:03 +01:00
committed by Andrew Heather
parent c07e1e0486
commit 2db489b562
1 changed files with 2 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
applications/solvers/multiphase/twoPhaseEulerFoam/phaseIncompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
View File

@ -422,19 +422,12 @@ void Foam::RASModels::kineticTheoryModel::correct()
*gs0_*(1.0/da)*ThetaSqrt/sqrtPi *gs0_*(1.0/da)*ThetaSqrt/sqrtPi
); );
// NB, drag = K*alpha*alpha2, // Drag
// (the alpha and alpha2 has been extracted from the drag function for
// numerical reasons)
volScalarField magUr(mag(U - Uc_)); volScalarField magUr(mag(U - Uc_));
volScalarField alpha2Prim volScalarField alpha2Prim
( (
max phase_.fluid().drag(phase_).K()/rho
(
alpha*(1.0 - alpha),
residualAlpha_
)
*phase_.fluid().drag(phase_).K()/rho
); );
// Eq. 3.25, p. 50 Js = J1 - J2 // Eq. 3.25, p. 50 Js = J1 - J2