32 template<
class MixtureType>
36 const scalarField& pCells = this->p_.internalField();
41 scalarField& alphaCells = this->alpha_.internalField();
45 const typename MixtureType::thermoType& mixture_ =
46 this->cellMixture(celli);
48 TCells[celli] = mixture_.TE(eCells[celli], TCells[celli]);
49 psiCells[celli] = mixture_.psi(pCells[celli], TCells[celli]);
51 muCells[celli] = mixture_.mu(TCells[celli]);
52 alphaCells[celli] = mixture_.alpha(TCells[celli]);
70 const typename MixtureType::thermoType& mixture_ =
71 this->patchFaceMixture(
patchi, facei);
73 pe[facei] = mixture_.E(pT[facei]);
75 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
76 pmu[facei] = mixture_.mu(pT[facei]);
77 palpha[facei] = mixture_.alpha(pT[facei]);
84 const typename MixtureType::thermoType& mixture_ =
85 this->patchFaceMixture(
patchi, facei);
87 pT[facei] = mixture_.TE(pe[facei], pT[facei]);
89 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
90 pmu[facei] = mixture_.mu(pT[facei]);
91 palpha[facei] = mixture_.alpha(pT[facei]);
100 template<
class MixtureType>
104 MixtureType(*this, mesh),
118 this->eBoundaryTypes()
126 eCells[celli] = this->cellMixture(celli).E(TCells[celli]);
146 template<
class MixtureType>
153 template<
class MixtureType>
158 Info<<
"entering ePsiThermo<MixtureType>::correct()" <<
endl;
162 this->psi_.oldTime();
168 Info<<
"exiting ePsiThermo<MixtureType>::correct()" <<
endl;
173 template<
class MixtureType>
185 h[celli] = this->cellMixture(cells[celli]).E(T[celli]);
192 template<
class MixtureType>
204 h[facei] = this->patchFaceMixture(patchi, facei).E(T[facei]);
211 template<
class MixtureType>
223 cp[facei] = this->patchFaceMixture(patchi, facei).Cp(T[facei]);
230 template<
class MixtureType>
256 cp[celli] = this->cellMixture(celli).Cp(this->T_[celli]);
266 pCp[facei] = this->patchFaceMixture(patchi, facei).Cp(pT[facei]);
274 template<
class MixtureType>
286 cv[facei] = this->patchFaceMixture(patchi, facei).Cv(T[facei]);
293 template<
class MixtureType>
319 cv[celli] = this->cellMixture(celli).Cv(this->T_[celli]);
332 template<
class MixtureType>
337 MixtureType::read(*
this);