Commit 556ae647 authored by andy's avatar andy
Browse files

STYLE: Minor code changes

parent f7aa88ba
......@@ -99,7 +99,7 @@ void Foam::SprayParcel<ParcelType>::calc
scalar rho1 = td.cloud().composition().liquids().rho(this->pc_, T1, X1);
this->rho() = rho1;
scalar d1 = this->d()*pow(rho0/rho1, 1.0/3.0);
scalar d1 = this->d()*cbrt(rho0/rho1);
this->d() = d1;
if (liquidCore() > 0.5)
......
......@@ -189,7 +189,7 @@ bool Foam::ReitzKHRT<CloudType>::update
d = (fraction*dc + d)/(1.0 + fraction);
//scalar ms0 = rho*pow3(dc)*mathematicalConstant::pi/6.0;
scalar ms0 = mass0*(1.0 - pow(d/d0,3.0));
scalar ms0 = mass0*(1.0 - pow3(d/d0));
ms += ms0;
if (ms/averageParcelMass > msLimit_)
......@@ -200,10 +200,10 @@ bool Foam::ReitzKHRT<CloudType>::update
// drops after breakup, see Eq. 18 in
// Patterson & Reitz, SAE 980131
bool br3 = true;
scalar ae3 = 1.;
scalar ae3 = 1.0;
scalar be3 = -dc;
scalar ce3 = 0.;
scalar de3 = d*d*(dc-d);
scalar ce3 = 0.0;
scalar de3 = d*d*(dc - d);
scalar qe3 =
pow3(be3/(3.0*ae3)) - be3*ce3/(6.0*ae3*ae3) + de3/(2.0*ae3);
scalar pe3 = (3.0*ae3*ce3 - be3*be3)/(9.0*ae3*ae3);
......@@ -214,8 +214,8 @@ bool Foam::ReitzKHRT<CloudType>::update
if (br3)
{
D3 = sqrt(D3);
scalar ue3 = cbrt(-qe3+D3);
scalar ve3 = cbrt(-qe3-D3);
scalar ue3 = cbrt(-qe3 + D3);
scalar ve3 = cbrt(-qe3 - D3);
scalar dParenDrops = ue3 + ve3 - be3/3.;
scalar mc = nParticle*(pow3(d) - pow3(dParenDrops));
scalar nChildDrops = mc/pow3(dc);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment