diff --git a/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.C b/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.C index 346a2dccc2..fa1eba5b90 100644 --- a/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.C +++ b/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.C @@ -315,8 +315,16 @@ void Foam::SprayCloud::motion(TrackData& td) scalarField Xq(liqMix.X(q.Y())); q.rho() = liqMix.rho(q.pc(), q.T(), Xq); q.Cp() = liqMix.Cp(q.pc(), q.T(), Xq); - scalar rhs = 6.0*mq/(q.nParticle()*q.rho()*constant::mathematical::pi); - q.d() = pow(rhs, 1.0/3.0); + q.d() = + cbrt + ( + 6.0*mq + /( + q.nParticle() + *q.rho() + *constant::mathematical::pi + ) + ); } } }