33template<
class BasicPsiThermo,
class MixtureType>
36 const volScalarField&
p,
42 volScalarField&
alpha,
48 if (doOldTimes && (
p.nOldTimes() ||
T.nOldTimes()))
75 this->cellMixture(celli);
79 TCells[celli] = mixture_.THE
87 psiCells[celli] = mixture_.psi(pCells[celli], TCells[celli]);
88 rhoCells[celli] = mixture_.rho(pCells[celli], TCells[celli]);
90 muCells[celli] = mixture_.mu(pCells[celli], TCells[celli]);
91 alphaCells[celli] = mixture_.alphah(pCells[celli], TCells[celli]);
94 const volScalarField::Boundary& pBf =
p.boundaryField();
95 volScalarField::Boundary& TBf =
T.boundaryFieldRef();
96 volScalarField::Boundary& psiBf =
psi.boundaryFieldRef();
97 volScalarField::Boundary& rhoBf =
rho.boundaryFieldRef();
98 volScalarField::Boundary& heBf =
he.boundaryFieldRef();
99 volScalarField::Boundary& muBf =
mu.boundaryFieldRef();
100 volScalarField::Boundary& alphaBf =
alpha.boundaryFieldRef();
117 this->patchFaceMixture(patchi, facei);
119 phe[facei] = mixture_.HE(pp[facei], pT[facei]);
121 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
122 prho[facei] = mixture_.rho(pp[facei], pT[facei]);
123 pmu[facei] = mixture_.mu(pp[facei], pT[facei]);
124 palpha[facei] = mixture_.alphah(pp[facei], pT[facei]);
132 this->patchFaceMixture(patchi, facei);
136 pT[facei] = mixture_.THE(phe[facei], pp[facei], pT[facei]);
139 ppsi[facei] = mixture_.psi(pp[facei], pT[facei]);
140 prho[facei] = mixture_.rho(pp[facei], pT[facei]);
141 pmu[facei] = mixture_.mu(pp[facei], pT[facei]);
142 palpha[facei] = mixture_.alphah(pp[facei], pT[facei]);
150template<
class BasicPsiThermo,
class MixtureType>
154 const word& phaseName
173template<
class BasicPsiThermo,
class MixtureType>
177 const word& phaseName,
199template<
class BasicPsiThermo,
class MixtureType>
206template<
class BasicPsiThermo,
class MixtureType>
Mesh data needed to do the Finite Volume discretisation.
Energy for a mixture based on density.
virtual ~heRhoThermo()
Destructor.
virtual void correct()
Update properties.
Enthalpy/Internal energy for a mixture.
volScalarField he_
Energy field.
ThermoType thermoType
The type of thermodynamics this mixture is instantiated for.
A class for handling words, derived from Foam::string.
const volScalarField & psi
const volScalarField & mu
const word dictName("faMeshDefinition")
#define DebugInFunction
Report an information message using Foam::Info.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
Ostream & endl(Ostream &os)
Add newline and flush stream.
fvPatchField< scalar > fvPatchScalarField
#define forAll(list, i)
Loop across all elements in list.