twoPhaseEulerFoam: Correction to GidaspowErgunWenYu drag model

Resolves bug-report http://www.openfoam.com/mantisbt/view.php?id=205
OpenCFD · May 26, 2011 · 40e7c06 · 40e7c06
1 parent 9b4c5b5
commit 40e7c06
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/...se/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C b/...se/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C
@@ -80,11 +80,11 @@ Foam::tmp<Foam::volScalarField> Foam::GidaspowErgunWenYu::K
     // Wen and Yu (1966)
     return
         pos(beta - 0.8)
-       *(0.75*alpha_*beta*Cds*phaseb_.rho()*Ur*bp/phasea_.d())
+       *(0.75*Cds*phaseb_.rho()*Ur*bp/phasea_.d())
       + neg(beta - 0.8)
        *(
-           150.0*sqr(alpha_)*phaseb_.nu()*phaseb_.rho()/(beta*sqr(phasea_.d()))
-         + 1.75*alpha_*phaseb_.rho()*Ur/phasea_.d()
+           150.0*alpha_*phaseb_.nu()*phaseb_.rho()/(sqr(beta*phasea_.d()))
+         + 1.75*phaseb_.rho()*Ur/(beta*phasea_.d())
         );
 }