34 template<
class BasicSol
idThermo,
class MixtureType>
41 scalarField& rhoCells = this->rho_.primitiveFieldRef();
42 scalarField& alphaCells = this->alpha_.primitiveFieldRef();
46 const typename MixtureType::thermoType& mixture_ =
47 this->cellMixture(celli);
49 const typename MixtureType::thermoType& volMixture_ =
50 this->cellVolMixture(pCells[celli], TCells[celli], celli);
54 TCells[celli] = mixture_.THE
62 rhoCells[celli] = volMixture_.rho(pCells[celli], TCells[celli]);
65 volMixture_.kappa(pCells[celli], TCells[celli])
67 mixture_.Cpv(pCells[celli], TCells[celli]);
70 volScalarField::Boundary& pBf = this->p_.boundaryFieldRef();
71 volScalarField::Boundary& TBf = this->T_.boundaryFieldRef();
72 volScalarField::Boundary& rhoBf = this->rho_.boundaryFieldRef();
73 volScalarField::Boundary& heBf = this->
he().boundaryFieldRef();
74 volScalarField::Boundary& alphaBf = this->alpha_.boundaryFieldRef();
76 forAll(this->T_.boundaryField(), patchi)
88 const typename MixtureType::thermoType& mixture_ =
89 this->patchFaceMixture(patchi, facei);
91 const typename MixtureType::thermoType& volMixture_ =
92 this->patchFaceVolMixture
100 phe[facei] = mixture_.HE(pp[facei], pT[facei]);
101 prho[facei] = volMixture_.rho(pp[facei], pT[facei]);
104 volMixture_.kappa(pp[facei], pT[facei])
105 / mixture_.Cpv(pp[facei], pT[facei]);
112 const typename MixtureType::thermoType& mixture_ =
113 this->patchFaceMixture(patchi, facei);
115 const typename MixtureType::thermoType& volMixture_ =
116 this->patchFaceVolMixture
126 pT[facei] = mixture_.THE(phe[facei], pp[facei] ,pT[facei]);
129 prho[facei] = volMixture_.rho(pp[facei], pT[facei]);
132 volMixture_.kappa(pp[facei], pT[facei])
133 / mixture_.Cpv(pp[facei], pT[facei]);
138 this->alpha_.correctBoundaryConditions();
143 template<
class BasicSol
idThermo,
class MixtureType>
148 const word& phaseName
159 template<
class BasicSol
idThermo,
class MixtureType>
165 const word& phaseName
176 template<
class BasicSol
idThermo,
class MixtureType>
181 const word& phaseName,
198 template<
class BasicSol
idThermo,
class MixtureType>
205 template<
class BasicSol
idThermo,
class MixtureType>
216 template<
class BasicSol
idThermo,
class MixtureType>
229 mesh.time().timeName(),
240 vectorField& KappaCells = Kappa.primitiveFieldRef();
252 ).Kappa(pCells[celli], TCells[celli]);
255 volVectorField::Boundary& KappaBf = Kappa.boundaryFieldRef();
260 const scalarField& pT = this->T_.boundaryField()[patchi];
261 const scalarField& pp = this->p_.boundaryField()[patchi];
266 this->patchFaceVolMixture
272 ).Kappa(pp[facei], pT[facei]);
280 template<
class BasicSol
idThermo,
class MixtureType>
287 const scalarField& pp = this->p_.boundaryField()[patchi];
288 const scalarField& Tp = this->T_.boundaryField()[patchi];
296 this->patchFaceVolMixture
302 ).Kappa(pp[facei], Tp[facei]);