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_.TH(hCells[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 ph[facei] = mixture_.H(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_.TH(ph[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 hCells[celli] = this->cellMixture(celli).H(TCells[celli]);
144 template<
class MixtureType>
151 template<
class MixtureType>
156 Info<<
"entering hPsiThermo<MixtureType>::correct()" <<
endl;
160 this->psi_.oldTime();
166 Info<<
"exiting hPsiThermo<MixtureType>::correct()" <<
endl;
171 template<
class MixtureType>
183 h[celli] = this->cellMixture(cells[celli]).H(T[celli]);
190 template<
class MixtureType>
202 h[facei] = this->patchFaceMixture(patchi, facei).H(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]);
330 template<
class MixtureType>
335 MixtureType::read(*
this);