From 04fe88592012870be5cacb94cb5c9ce2aa84339e Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Sun, 17 Apr 2016 15:39:09 +0100
Subject: [PATCH] BrownianMotionForce: Corrected expression for s0 Patch
provided by Bruno Santos Resolves bug-report
http://www.openfoam.org/mantisbt/view.php?id=2036
---
.../ParticleForces/BrownianMotion/BrownianMotionForce.C | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/src/lagrangian/turbulence/submodels/Thermodynamic/ParticleForces/BrownianMotion/BrownianMotionForce.C b/src/lagrangian/turbulence/submodels/Thermodynamic/ParticleForces/BrownianMotion/BrownianMotionForce.C
index 9dffc11bcd..3903093628 100644
--- a/src/lagrangian/turbulence/submodels/Thermodynamic/ParticleForces/BrownianMotion/BrownianMotionForce.C
+++ b/src/lagrangian/turbulence/submodels/Thermodynamic/ParticleForces/BrownianMotion/BrownianMotionForce.C
@@ -187,9 +187,8 @@ Foam::forceSuSp Foam::BrownianMotionForce<CloudType>::calcCoupled
}
else
{
- const scalar rhoRatio = p.rho()/p.rhoc();
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);
}
--
GitLab