BrownianMotionForce: Corrected expression for s0

Patch provided by Bruno Santos
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=2036
This commit is contained in:
Henry Weller
2016-04-17 15:39:09 +01:00
parent 2e990f19d8
commit 6862e3d76e

View File

@ -187,9 +187,8 @@ Foam::forceSuSp Foam::BrownianMotionForce<CloudType>::calcCoupled
} }
else else
{ {
const scalar rhoRatio = p.rho()/p.rhoc();
const scalar s0 = const scalar s0 =
216*muc*sigma*Tc/(sqr(mathematical::pi)*pow5(dp)*(rhoRatio)*cc); 216*muc*sigma*Tc/(sqr(mathematical::pi)*pow5(dp)*sqr(p.rho())*cc);
f = eta*sqrt(mathematical::pi*s0/dt); f = eta*sqrt(mathematical::pi*s0/dt);
} }