33 template<
class BasicPsiThermo,
class MixtureType>
48 if (doOldTimes && (
p.nOldTimes() ||
T.nOldTimes()))
74 const typename MixtureType::thermoType& mixture_ =
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();
116 const typename MixtureType::thermoType& mixture_ =
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]);
131 const typename MixtureType::thermoType& mixture_ =
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]);
150 template<
class BasicPsiThermo,
class MixtureType>
154 const word& phaseName
157 heThermo<BasicPsiThermo, MixtureType>(
mesh, phaseName)
173 template<
class BasicPsiThermo,
class MixtureType>
177 const word& phaseName,
199 template<
class BasicPsiThermo,
class MixtureType>
206 template<
class BasicPsiThermo,
class MixtureType>
virtual ~heRhoThermo()
Destructor.
volScalarField he_
Energy field.
virtual void correct()
Update properties.
const word dictName("faMeshDefinition")
Ostream & endl(Ostream &os)
Add newline and flush stream.
#define forAll(list, i)
Loop across all elements in list.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
fvPatchField< scalar > fvPatchScalarField
A class for handling words, derived from Foam::string.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
#define DebugInFunction
Report an information message using Foam::Info.
const dimensionedScalar mu
Atomic mass unit.
Enthalpy/Internal energy for a mixture.
Mesh data needed to do the Finite Volume discretisation.
Energy for a mixture based on density.
const volScalarField & psi
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
uindirectPrimitivePatch pp(UIndirectList< face >(mesh.faces(), faceLabels), mesh.points())