32 const Foam::scalar Foam::liquidMixture::TrMax = 0.999;
41 components_(thermophysicalProperties.
lookup(
"liquidComponents")),
42 properties_(components_.size())
59 liquid::New(subDictPtr->
lookup(components_[i]))
67 liquid::New(thermophysicalProperties.
lookup(components_[i]))
98 scalar x1 = x[i]*properties_[i].Vc();
100 vTc += x1*properties_[i].Tc();
115 Tpc += x[i]*properties_[i].Tc();
131 Vc += x[i]*properties_[i].Vc();
132 Zc += x[i]*properties_[i].Zc();
135 return specie::RR*Zc*Tpc(x)/Vc;
147 omega += x[i]*properties_[i].omega();
168 scalar Ti =
min(TrMax*properties_[i].Tc(), Tl);
169 xs[i] = properties_[i].pv(p, Ti)*xl[i]/
p;
183 W += x[i]*properties_[i].W();
199 Y[i] *= properties_[i].W();
215 Winv += Y[i]/properties_[i].W();
216 X[i] = Y[i]/properties_[i].W();
236 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
237 scalar
rho = SMALL + properties_[i].rho(p, Ti);
238 v += x[i]*properties_[i].W()/
rho;
259 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
260 pv += x[i]*properties_[i].pv(p, Ti)*properties_[i].W();
281 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
282 hl += x[i]*properties_[i].hl(p, Ti)*properties_[i].W();
303 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
304 cp += x[i]*properties_[i].cp(p, Ti)*properties_[i].W();
326 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
327 scalar Pvs = properties_[i].pv(p, Ti);
328 scalar xs = x[i]*Pvs/
p;
337 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
338 sigma += (Xs[i]/XsSum)*properties_[i].sigma(p, Ti);
359 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
360 mu += x[i]*
log(properties_[i].
mu(p, Ti));
381 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
383 scalar Vi = properties_[i].W()/properties_[i].rho(p, Ti);
397 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
401 scalar Tj =
min(TrMax*properties_[j].Tc(), T);
406 1.0/properties_[i].K(p, Ti)
407 + 1.0/properties_[j].K(p, Tj)
409 K += phii[i]*phii[j]*Kij;
431 scalar Ti =
min(TrMax*properties_[i].Tc(), T);
432 Dinv += x[i]/properties_[i].D(p, Ti);