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>
229 template<
class BasicSol
idThermo,
class MixtureType>
242 mesh.time().timeName(),
253 vectorField& KappaCells = Kappa.primitiveFieldRef();
265 ).Kappa(pCells[celli], TCells[celli]);
268 volVectorField::Boundary& KappaBf = Kappa.boundaryFieldRef();
273 const scalarField& pT = this->T_.boundaryField()[patchi];
274 const scalarField& pp = this->p_.boundaryField()[patchi];
279 this->patchFaceVolMixture
285 ).Kappa(pp[facei], pT[facei]);
293 template<
class BasicSol
idThermo,
class MixtureType>
300 const scalarField& pp = this->p_.boundaryField()[patchi];
301 const scalarField& Tp = this->T_.boundaryField()[patchi];
309 this->patchFaceVolMixture
315 ).Kappa(pp[facei], Tp[facei]);