30 template<
class BasicPsiThermo,
class MixtureType>
37 scalarField& psiCells = this->psi_.primitiveFieldRef();
38 scalarField& rhoCells = this->rho_.primitiveFieldRef();
39 scalarField& muCells = this->mu_.primitiveFieldRef();
40 scalarField& alphaCells = this->alpha_.primitiveFieldRef();
44 const typename MixtureType::thermoType& mixture_ =
45 this->cellMixture(celli);
47 TCells[celli] = mixture_.THE
54 psiCells[celli] = mixture_.psi(pCells[celli], TCells[celli]);
55 rhoCells[celli] = mixture_.rho(pCells[celli], TCells[celli]);
57 muCells[celli] = mixture_.mu(pCells[celli], TCells[celli]);
58 alphaCells[celli] = mixture_.alphah(pCells[celli], TCells[celli]);
61 volScalarField::Boundary& pBf =
62 this->p_.boundaryFieldRef();
64 volScalarField::Boundary& TBf =
65 this->T_.boundaryFieldRef();
67 volScalarField::Boundary& psiBf =
68 this->psi_.boundaryFieldRef();
70 volScalarField::Boundary& rhoBf =
71 this->rho_.boundaryFieldRef();
73 volScalarField::Boundary& heBf =
74 this->
he().boundaryFieldRef();
76 volScalarField::Boundary& muBf =
77 this->mu_.boundaryFieldRef();
79 volScalarField::Boundary& alphaBf =
80 this->alpha_.boundaryFieldRef();
96 const typename MixtureType::thermoType& mixture_ =
97 this->patchFaceMixture(
patchi, facei);
99 phe[facei] = mixture_.HE(pp[facei], pT[facei]);
101 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
102 prho[facei] = mixture_.rho(pp[facei], pT[facei]);
103 pmu[facei] = mixture_.mu(pp[facei], pT[facei]);
104 palpha[facei] = mixture_.alphah(pp[facei], pT[facei]);
111 const typename MixtureType::thermoType& mixture_ =
112 this->patchFaceMixture(
patchi, facei);
114 pT[facei] = mixture_.THE(phe[facei], pp[facei], pT[facei]);
116 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
117 prho[facei] = mixture_.rho(pp[facei], pT[facei]);
118 pmu[facei] = mixture_.mu(pp[facei], pT[facei]);
119 palpha[facei] = mixture_.alphah(pp[facei], pT[facei]);
128 template<
class BasicPsiThermo,
class MixtureType>
132 const word& phaseName
143 template<
class BasicPsiThermo,
class MixtureType>
150 template<
class BasicPsiThermo,
class MixtureType>
#define forAll(list, i)
Loop across all elements in list.
virtual ~heRhoThermo()
Destructor.
virtual void correct()
Update properties.
Ostream & endl(Ostream &os)
Add newline and flush stream.
fvPatchField< scalar > fvPatchScalarField
A class for handling words, derived from string.
volScalarField scalarField(fieldObject, mesh)
Mesh data needed to do the Finite Volume discretisation.
Energy for a mixture based on density.
#define InfoInFunction
Report an information message using Foam::Info.