Go to the documentation of this file.
30 template<
class EquationOfState>
33 const EquationOfState& st,
50 template<
class EquationOfState>
57 EquationOfState(
name, ct),
65 template<
class EquationOfState>
76 template<
class EquationOfState>
89 template<
class EquationOfState>
99 template<
class EquationOfState>
110 template<
class EquationOfState>
113 const scalar
p,
const scalar
T
116 return Cv_*(
T - Tref_) + Eref_ + EquationOfState::E(
p,
T);
120 template<
class EquationOfState>
127 template<
class EquationOfState>
130 const scalar
p,
const scalar
T
133 return Es(
p,
T) + Hc();
137 template<
class EquationOfState>
140 const scalar
p,
const scalar
T
147 template<
class EquationOfState>
150 const scalar
p,
const scalar
T
157 template<
class EquationOfState>
160 const scalar
p,
const scalar
T
169 template<
class EquationOfState>
175 scalar Y1 = this->
Y();
177 EquationOfState::operator+=(ct);
179 if (
mag(this->
Y()) > SMALL)
182 const scalar Y2 = ct.Y()/this->
Y();
184 Cv_ = Y1*Cv_ + Y2*ct.Cv_;
185 Hf_ = Y1*Hf_ + Y2*ct.Hf_;
192 template<
class EquationOfState>
201 static_cast<const EquationOfState&
>(ct1)
202 +
static_cast<const EquationOfState&
>(ct2)
205 if (
mag(eofs.Y()) < SMALL)
218 return eRefConstThermo<EquationOfState>
221 ct1.Y()/eofs.Y()*ct1.Cv_
222 + ct2.Y()/eofs.Y()*ct2.Cv_,
223 ct1.Y()/eofs.Y()*ct1.Hf_
224 + ct2.Y()/eofs.Y()*ct2.Hf_,
225 ct1.Y()/eofs.Y()*ct1.Tref_
226 + ct2.Y()/eofs.Y()*ct2.Tref_,
227 ct1.Y()/eofs.Y()*ct1.Eref_
228 + ct2.Y()/eofs.Y()*ct2.Eref_
234 template<
class EquationOfState>
238 const eRefConstThermo<EquationOfState>& ct
241 return eRefConstThermo<EquationOfState>
243 s*
static_cast<const EquationOfState&
>(ct),
252 template<
class EquationOfState>
255 const eRefConstThermo<EquationOfState>& ct1,
256 const eRefConstThermo<EquationOfState>& ct2
261 static_cast<const EquationOfState&
>(ct1)
262 ==
static_cast<const EquationOfState&
>(ct2)
265 return eRefConstThermo<EquationOfState>
268 ct2.Y()/eofs.Y()*ct2.Cv_
269 - ct1.Y()/eofs.Y()*ct1.Cv_,
270 ct2.Y()/eofs.Y()*ct2.Hf_
271 - ct1.Y()/eofs.Y()*ct1.Hf_
A class for handling words, derived from Foam::string.
scalar dCpdT(const scalar p, const scalar T) const
Temperature derivative of heat capacity at constant pressure.
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
const volScalarField & Cv
scalar Cv(const scalar p, const scalar T) const
Heat capacity at constant volume [J/(kg K)].
scalar limit(const scalar T) const
Limit the temperature to be in the range Tlow_ to Thigh_.
const dimensionedScalar Tstd
Standard temperature.
scalar dGdT(const scalar p, const scalar T) const
Derivative of Gibbs free energy w.r.t. temperature.
word name(const complex &c)
Return string representation of complex.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
scalar Es(const scalar p, const scalar T) const
Sensible internal energy [J/kg].
autoPtr< eRefConstThermo > clone() const
Construct and return a clone.
scalar Hc() const
Chemical enthalpy [J/kg].
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
dimensionedScalar log(const dimensionedScalar &ds)
PtrList< volScalarField > & Y
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
scalar S(const scalar p, const scalar T) const
Entropy [J/(kg K)].
scalar Es(const scalar p, const scalar T) const
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const volScalarField & Cp
Constant properties thermodynamics package templated into the EquationOfState.
scalar Ea(const scalar p, const scalar T) const
Absolute internal energy [J/kg].
static autoPtr< eRefConstThermo > New(const dictionary &dict)
Selector from dictionary.