58 coeffsDict_(dict.
subDict(typeName +
"Coeffs")),
61 Cstrip_(
readScalar(coeffsDict_.lookup(
"Cstrip"))),
93 label cellI = p.
cell();
94 scalar pressure = spray_.p()[cellI];
95 scalar temperature = spray_.T()[cellI];
96 scalar Taverage = p.
T() + (temperature - p.
T())/3.0;
98 scalar muAverage = 0.0;
100 for(label i=0; i<Ns; i++)
102 Winv += Y[i][cellI]/spray_.gasProperties()[i].W();
103 muAverage += Y[i][cellI]*spray_.gasProperties()[i].mu(Taverage);
108 scalar rhoAverage = pressure/R/Taverage;
109 scalar nuAverage = muAverage/rhoAverage;
110 scalar sigma = fuels.
sigma(pressure, p.
T(), p.
X());
116 scalar WeberNumber = p.
We(vel, rhoAverage, sigma);
117 scalar ReynoldsNumber = p.
Re(vel, nuAverage);
119 scalar sqRey =
sqrt(ReynoldsNumber);
121 if (WeberNumber > Cbag_)
123 if (WeberNumber > Cstrip_*sqRey)
126 pow(2.0*Cstrip_*sigma, 2.0)/
137 fuels.
rho(pressure, p.
T(), p.
X())
142 scalar fraction = deltaT/tauStrip;
145 p.
d() = (fraction*dStrip + p.
d())/(1.0 + fraction);
160 fuels.
rho(pressure, p.
T(), p.
X())
165 scalar fraction = deltaT/tauBag;
168 p.
d() = (fraction*dBag + p.
d())/(1.0 + fraction);