From 2db489b562bf87968cad66d93ec74af7348e47b5 Mon Sep 17 00:00:00 2001 From: william Date: Tue, 1 Apr 2014 09:12:03 +0100 Subject: [PATCH] BUG: twoPhaseEulerFoam: fixed definition of drag in kineticTheory turbulence model --- .../kineticTheoryModel/kineticTheoryModel.C | 11 ++--------- 1 file changed, 2 insertions(+), 9 deletions(-) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseIncompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseIncompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C index 29baaf9ced..6db46aea0a 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseIncompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseIncompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C @@ -422,19 +422,12 @@ void Foam::RASModels::kineticTheoryModel::correct() *gs0_*(1.0/da)*ThetaSqrt/sqrtPi ); - // NB, drag = K*alpha*alpha2, - // (the alpha and alpha2 has been extracted from the drag function for - // numerical reasons) + // Drag volScalarField magUr(mag(U - Uc_)); volScalarField alpha2Prim ( - max - ( - alpha*(1.0 - alpha), - residualAlpha_ - ) - *phase_.fluid().drag(phase_).K()/rho + phase_.fluid().drag(phase_).K()/rho ); // Eq. 3.25, p. 50 Js = J1 - J2