30 template<
class MixtureType>
34 const scalarField& pCells = this->p_.internalField();
39 scalarField& alphaCells = this->alpha_.internalField();
43 const typename MixtureType::thermoType& mixture_ =
44 this->cellMixture(celli);
46 TCells[celli] = mixture_.THs(hsCells[celli], TCells[celli]);
47 psiCells[celli] = mixture_.psi(pCells[celli], TCells[celli]);
49 muCells[celli] = mixture_.mu(TCells[celli]);
50 alphaCells[celli] = mixture_.alpha(TCells[celli]);
68 const typename MixtureType::thermoType& mixture_ =
69 this->patchFaceMixture(
patchi, facei);
71 phs[facei] = mixture_.Hs(pT[facei]);
73 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
74 pmu[facei] = mixture_.mu(pT[facei]);
75 palpha[facei] = mixture_.alpha(pT[facei]);
82 const typename MixtureType::thermoType& mixture_ =
83 this->patchFaceMixture(
patchi, facei);
85 pT[facei] = mixture_.THs(phs[facei], pT[facei]);
87 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
88 pmu[facei] = mixture_.mu(pT[facei]);
89 palpha[facei] = mixture_.alpha(pT[facei]);
98 template<
class MixtureType>
102 MixtureType(*this, mesh),
116 this->hBoundaryTypes()
124 hsCells[celli] = this->cellMixture(celli).Hs(TCells[celli]);
144 template<
class MixtureType>
151 template<
class MixtureType>
156 Info<<
"entering hsPsiThermo<MixtureType>::correct()" <<
endl;
160 this->psi_.oldTime();
166 Info<<
"exiting hsPsiThermo<MixtureType>::correct()" <<
endl;
171 template<
class MixtureType>
183 hs[celli] = this->cellMixture(cells[celli]).Hs(T[celli]);
190 template<
class MixtureType>
202 hs[facei] = this->patchFaceMixture(patchi, facei).Hs(T[facei]);
209 template<
class MixtureType>
221 cp[facei] = this->patchFaceMixture(patchi, facei).Cp(T[facei]);
228 template<
class MixtureType>
254 cp[celli] = this->cellMixture(celli).Cp(this->T_[celli]);
264 pCp[facei] = this->patchFaceMixture(patchi, facei).Cp(pT[facei]);
272 template<
class MixtureType>
284 cv[facei] = this->patchFaceMixture(patchi, facei).Cv(T[facei]);
291 template<
class MixtureType>
317 cv[celli] = this->cellMixture(celli).Cv(this->T_[celli]);
322 cv.boundaryField()[
patchi] =
330 template<
class MixtureType>
335 MixtureType::read(*
this);