28 template<
class EquationOfState>
31 const EquationOfState& st,
44 template<
class EquationOfState>
51 EquationOfState(name, ct),
57 template<
class EquationOfState>
68 template<
class EquationOfState>
81 template<
class EquationOfState>
91 template<
class EquationOfState>
98 return Cp_ + EquationOfState::Cp(p, T);
102 template<
class EquationOfState>
105 const scalar p,
const scalar T
108 return Cp_*T + Hf_ + EquationOfState::H(p, T);
112 template<
class EquationOfState>
115 const scalar p,
const scalar T
118 return Cp_*T + EquationOfState::H(p, T);
122 template<
class EquationOfState>
129 template<
class EquationOfState>
132 const scalar p,
const scalar T
135 return Cp_*
log(T/
Tstd) + EquationOfState::S(p, T);
141 template<
class EquationOfState>
142 inline void Foam::hConstThermo<EquationOfState>::operator+=
147 scalar
Y1 = this->
Y();
149 EquationOfState::operator+=(ct);
151 if (
mag(this->
Y()) > SMALL)
154 scalar
Y2 = ct.Y()/this->
Y();
156 Cp_ = Y1*Cp_ + Y2*ct.Cp_;
157 Hf_ = Y1*Hf_ + Y2*ct.Hf_;
164 template<
class EquationOfState>
173 static_cast<const EquationOfState&>(ct1)
174 + static_cast<const EquationOfState&>(ct2)
177 if (
mag(eofs.Y()) < SMALL)
191 ct1.Y()/eofs.Y()*ct1.Cp_
192 + ct2.Y()/eofs.Y()*ct2.Cp_,
193 ct1.Y()/eofs.Y()*ct1.Hf_
194 + ct2.Y()/eofs.Y()*ct2.Hf_
200 template<
class EquationOfState>
209 s*
static_cast<const EquationOfState&
>(ct),
216 template<
class EquationOfState>
225 static_cast<const EquationOfState&>(ct1)
226 == static_cast<const EquationOfState&>(ct2)
232 ct2.Y()/eofs.Y()*ct2.Cp_
233 - ct1.Y()/eofs.Y()*ct1.Cp_,
234 ct2.Y()/eofs.Y()*ct2.Hf_
235 - ct1.Y()/eofs.Y()*ct1.Hf_
PtrList< volScalarField > & Y1
scalar Hs(const scalar p, const scalar T) const
Sensible enthalpy [J/kg].
Constant properties thermodynamics package templated into the EquationOfState.
scalar limit(const scalar T) const
Limit the temperature to be in the range Tlow_ to Thigh_.
dimensionedScalar log(const dimensionedScalar &ds)
A list of keyword definitions, which are a keyword followed by any number of values (e...
autoPtr< hConstThermo > clone() const
Construct and return a clone.
static autoPtr< hConstThermo > New(const dictionary &dict)
Selector from dictionary.
PtrList< volScalarField > & Y2
scalar Ha(const scalar p, const scalar T) const
Absolute Enthalpy [J/kg].
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().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))
A class for handling words, derived from string.
scalar Cp(const scalar p, const scalar T) const
Heat capacity at constant pressure [J/(kg K)].
const dimensionedScalar Tstd
Standard temperature.
const volScalarField & cp
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
PtrList< volScalarField > & Y
dimensioned< scalar > mag(const dimensioned< Type > &)
scalar S(const scalar p, const scalar T) const
Entropy [J/(kg K)].
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
scalar Hc() const
Chemical enthalpy [J/kg].