30 template<
class BasicPsiThermo,
class MixtureType>
37 scalarField& CpCells = this->Cp_.primitiveFieldRef();
38 scalarField& CvCells = this->Cv_.primitiveFieldRef();
39 scalarField& psiCells = this->psi_.primitiveFieldRef();
40 scalarField& muCells = this->mu_.primitiveFieldRef();
41 scalarField& alphaCells = this->alpha_.primitiveFieldRef();
45 const typename MixtureType::thermoMixtureType& thermoMixture =
46 this->cellThermoMixture(celli);
48 const typename MixtureType::transportMixtureType& transportMixture =
49 this->cellTransportMixture(celli, thermoMixture);
51 TCells[celli] = thermoMixture.THE
58 CpCells[celli] = thermoMixture.Cp(pCells[celli], TCells[celli]);
59 CvCells[celli] = thermoMixture.Cv(pCells[celli], TCells[celli]);
60 psiCells[celli] = thermoMixture.psi(pCells[celli], TCells[celli]);
62 muCells[celli] = transportMixture.mu(pCells[celli], TCells[celli]);
64 transportMixture.kappa(pCells[celli], TCells[celli])
65 /thermoMixture.Cp(pCells[celli], TCells[celli]);
68 volScalarField::Boundary& pBf =
69 this->p_.boundaryFieldRef();
71 volScalarField::Boundary& TBf =
72 this->T_.boundaryFieldRef();
74 volScalarField::Boundary& CpBf =
75 this->Cp_.boundaryFieldRef();
77 volScalarField::Boundary& CvBf =
78 this->Cv_.boundaryFieldRef();
80 volScalarField::Boundary& psiBf =
81 this->psi_.boundaryFieldRef();
83 volScalarField::Boundary& heBf =
84 this->
he().boundaryFieldRef();
86 volScalarField::Boundary& muBf =
87 this->mu_.boundaryFieldRef();
89 volScalarField::Boundary& alphaBf =
90 this->alpha_.boundaryFieldRef();
107 const typename MixtureType::thermoMixtureType&
108 thermoMixture = this->patchFaceThermoMixture(
patchi, facei);
110 const typename MixtureType::transportMixtureType&
112 this->patchFaceTransportMixture
113 (
patchi, facei, thermoMixture);
115 phe[facei] = thermoMixture.HE(pp[facei], pT[facei]);
117 pCp[facei] = thermoMixture.Cp(pp[facei], pT[facei]);
118 pCv[facei] = thermoMixture.Cv(pp[facei], pT[facei]);
119 ppsi[facei] = thermoMixture.psi(pp[facei], pT[facei]);
121 pmu[facei] = transportMixture.mu(pp[facei], pT[facei]);
123 transportMixture.kappa(pp[facei], pT[facei])
124 /thermoMixture.Cp(pp[facei], pT[facei]);
131 const typename MixtureType::thermoMixtureType& thermoMixture =
132 this->patchFaceThermoMixture(
patchi, facei);
134 const typename MixtureType::transportMixtureType&
136 this->patchFaceTransportMixture
137 (
patchi, facei, thermoMixture);
139 pT[facei] = thermoMixture.THE(phe[facei], pp[facei], pT[facei]);
141 pCp[facei] = thermoMixture.Cp(pp[facei], pT[facei]);
142 pCv[facei] = thermoMixture.Cv(pp[facei], pT[facei]);
143 ppsi[facei] = thermoMixture.psi(pp[facei], pT[facei]);
145 pmu[facei] = transportMixture.mu(pp[facei], pT[facei]);
147 transportMixture.kappa(pp[facei], pT[facei])
148 /thermoMixture.Cp(pp[facei], pT[facei]);
157 template<
class BasicPsiThermo,
class MixtureType>
161 const word& phaseName
169 this->psi_.oldTime();
175 template<
class BasicPsiThermo,
class MixtureType>
182 template<
class BasicPsiThermo,
class MixtureType>
191 this->psi_.oldTime();
#define forAll(list, i)
Loop across all elements in list.
Ostream & endl(Ostream &os)
Add newline and flush stream.
fvPatchField< scalar > fvPatchScalarField
A class for handling words, derived from string.
volScalarField scalarField(fieldObject, mesh)
hePsiThermo(const fvMesh &, const word &phaseName)
Construct from mesh and phase name.
virtual ~hePsiThermo()
Destructor.
Mesh data needed to do the Finite Volume discretisation.
Energy for a mixture based on compressibility.
virtual void correct()
Update properties.
#define InfoInFunction
Report an information message using Foam::Info.