34 template<
class BasicSol
idThermo,
class MixtureType>
41 scalarField& rhoCells = this->rho_.primitiveFieldRef();
44 scalarField& alphaCells = this->alpha_.primitiveFieldRef();
48 const typename MixtureType::thermoType& mixture_ =
49 this->cellMixture(celli);
51 const typename MixtureType::thermoType& volMixture_ =
52 this->cellVolMixture(pCells[celli], TCells[celli], celli);
56 TCells[celli] = mixture_.THE
64 rhoCells[celli] = volMixture_.rho(pCells[celli], TCells[celli]);
69 volMixture_.kappa(pCells[celli], TCells[celli])
71 mixture_.Cpv(pCells[celli], TCells[celli]);
74 volScalarField::Boundary& pBf = this->p_.boundaryFieldRef();
75 volScalarField::Boundary& TBf = this->T_.boundaryFieldRef();
76 volScalarField::Boundary& rhoBf = this->rho_.boundaryFieldRef();
79 volScalarField::Boundary& heBf = this->
he().boundaryFieldRef();
80 volScalarField::Boundary& alphaBf = this->alpha_.boundaryFieldRef();
82 forAll(this->T_.boundaryField(), patchi)
96 const typename MixtureType::thermoType& mixture_ =
97 this->patchFaceMixture(patchi, facei);
99 const typename MixtureType::thermoType& volMixture_ =
100 this->patchFaceVolMixture
109 phe[facei] = mixture_.HE(pp[facei], pT[facei]);
110 prho[facei] = volMixture_.rho(pp[facei], pT[facei]);
115 volMixture_.kappa(pp[facei], pT[facei])
116 / mixture_.Cpv(pp[facei], pT[facei]);
123 const typename MixtureType::thermoType& mixture_ =
124 this->patchFaceMixture(patchi, facei);
126 const typename MixtureType::thermoType& volMixture_ =
127 this->patchFaceVolMixture
137 pT[facei] = mixture_.THE(phe[facei], pp[facei] ,pT[facei]);
140 prho[facei] = volMixture_.rho(pp[facei], pT[facei]);
145 volMixture_.kappa(pp[facei], pT[facei])
146 / mixture_.Cpv(pp[facei], pT[facei]);
151 this->alpha_.correctBoundaryConditions();
156 template<
class BasicSol
idThermo,
class MixtureType>
161 const word& phaseName
172 template<
class BasicSol
idThermo,
class MixtureType>
178 const word& phaseName
189 template<
class BasicSol
idThermo,
class MixtureType>
194 const word& phaseName,
211 template<
class BasicSol
idThermo,
class MixtureType>
218 template<
class BasicSol
idThermo,
class MixtureType>
235 template<
class BasicSol
idThermo,
class MixtureType>
248 mesh.time().timeName(),
259 vectorField& KappaCells = Kappa.primitiveFieldRef();
271 ).Kappa(pCells[celli], TCells[celli]);
274 volVectorField::Boundary& KappaBf = Kappa.boundaryFieldRef();
279 const scalarField& pT = this->T_.boundaryField()[patchi];
280 const scalarField& pp = this->p_.boundaryField()[patchi];
285 this->patchFaceVolMixture
291 ).Kappa(pp[facei], pT[facei]);
299 template<
class BasicSol
idThermo,
class MixtureType>
306 const scalarField& pp = this->p_.boundaryField()[patchi];
307 const scalarField& Tp = this->T_.boundaryField()[patchi];
315 this->patchFaceVolMixture
321 ).Kappa(pp[facei], Tp[facei]);