43 standardEvaporationModel,
57 evapDict_(dict.
subDict(typeName +
"Coeffs")),
58 preReScFactor_(
readScalar(evapDict_.lookup(
"preReScFactor"))),
59 ReExponent_(
readScalar(evapDict_.lookup(
"ReExponent"))),
60 ScExponent_(
readScalar(evapDict_.lookup(
"ScExponent"))),
61 evaporationScheme_(evapDict_.lookup(
"evaporationScheme")),
64 if (evaporationScheme_ ==
"implicit")
68 else if (evaporationScheme_ ==
"explicit")
75 <<
"evaporationScheme type " << evaporationScheme_
77 <<
"Use implicit or explicit."
99 const scalar ReynoldsNumber,
100 const scalar SchmidtNumber
103 return 2.0 + preReScFactor_*
pow(ReynoldsNumber,ReExponent_)*
pow(SchmidtNumber,ScExponent_);
108 const scalar diameter,
109 const scalar liquidDensity,
110 const scalar rhoFuelVapor,
111 const scalar massDiffusionCoefficient,
112 const scalar ReynoldsNumber,
113 const scalar SchmidtNumber,
141 scalar Xratio = (Xs - Xf)/
max(SMALL, 1.0 - Xs);
145 lgExpr =
log(1.0 + Xratio);
153 6.0 * massDiffusionCoefficient
154 * Sh(ReynoldsNumber, SchmidtNumber)
155 * rhoFuelVapor * lgExpr;
157 if (denominator > SMALL)
159 time =
max(VSMALL, liquidDensity *
pow(diameter, 2.0)/denominator);
168 const scalar liquidDensity,
170 const scalar heatOfVapour,
172 const scalar Nusselt,
173 const scalar deltaTemp,
174 const scalar diameter,
192 scalar deltaT =
max(0.5, deltaTemp);
194 time = liquidDensity*cpFuel*
sqr(diameter)/
196 6.0 * kappa * Nusselt *
log(1.0 + cpFuel * deltaT/
max(SMALL, heatOfVapour))
199 time =
max(VSMALL, time);