41 void parcel::setRelaxationTimes
46 scalar& tauHeatTransfer,
51 const scalar temperature,
52 const scalar pressure,
61 scalar mCell = rho*sDB.
mesh().
V()[
cell()];
68 scalar Tf =
T() + (temperature -
T())/3.0;
72 scalar kMixture = 0.0;
73 scalar cpMixture = 0.0;
76 for(label i=0; i<Ns; i++)
92 for(label i=0; i<Nf; i++)
99 psat[i] = fuels.
properties()[i].pv(pressure, temperature);
100 msat[i] =
min(1.0, psat[i]/pressure)*Wi/W;
103 scalar nuf = muf/
rho;
105 scalar liquidDensity = fuels.
rho(pressure,
T(),
X());
106 scalar liquidcL = fuels.
cp(pressure,
T(),
X());
107 scalar heatOfVapour = fuels.
hl(pressure,
T(),
X());
116 scalar pressureAtSurface = fuels.
pv(pressure,
T(),
X());
117 scalar rhoFuelVap = pressureAtSurface*fuels.
W(
X())/(
specie::RR*Tf);
123 for(label i=0; i<Nf; i++)
130 for(label i=0; i<Nf; i++)
134 Ys[i] = Xs[i]*Wi/Wliq;
137 scalar Reynolds =
Re(Up, nuf);
138 scalar Prandtl =
Pr(cpMixture, muf, kMixture);
179 for(label i=0; i<Nf; i++)
182 bool boiling = fuels.
properties()[i].pv(pressure, Td) >= 0.999*pressure;
183 scalar Di = fuels.
properties()[i].D(pressure, Td);
184 scalar Schmidt =
Sc(nuf, Di);
186 scalar partialPressure = Xf[i]*pressure;
189 if(partialPressure > psat[i])
191 tauEvaporation[i] = GREAT;
199 scalar dm =
max(SMALL, 0.9999*msat[i] - mfg[i]);
222 scalar tBoilingSurface = Td;
225 scalar deltaT = 10.0;
226 scalar dp0 = fuels.
properties()[i].pv(pressure, tBoilingSurface) - pressure;
227 while ((Niter < 200) && (
mag(deltaT) > 1.0
e-3))
230 scalar pBoil = fuels.
properties()[i].pv(pressure, tBoilingSurface);
231 scalar dp = pBoil - pressure;
232 if ( (dp > 0.0) && (dp0 > 0.0) )
234 tBoilingSurface -= deltaT;
238 if ( (dp < 0.0) && (dp0 < 0.0) )
240 tBoilingSurface += deltaT;
245 if ( (dp > 0.0) && (dp0 < 0.0) )
247 tBoilingSurface -= deltaT;
251 tBoilingSurface += deltaT;
258 scalar vapourSurfaceEnthalpy = 0.0;
259 scalar vapourFarEnthalpy = 0.0;
267 scalar kLiquid = fuels.
properties()[i].K(pressure, 0.5*(tBoilingSurface+
T()));
282 vapourSurfaceEnthalpy,