31 template<
class EquationOfState>
40 <<
"attempt to evaluate ePowerThermo<EquationOfState>" 41 " for negative temperature " << T
47 template<
class EquationOfState>
54 EquationOfState(name, jt),
64 template<
class EquationOfState>
67 const EquationOfState& st,
82 template<
class EquationOfState>
93 template<
class EquationOfState>
106 template<
class EquationOfState>
116 template<
class EquationOfState>
127 template<
class EquationOfState>
135 c0_*(
pow(T, n0_ + 1) -
pow(
Tstd, n0_ + 1))/(
pow(Tref_, n0_)*(n0_ + 1))
136 + EquationOfState::E(p, T);
140 template<
class EquationOfState>
147 return Es(p, T) + Hf();
151 template<
class EquationOfState>
158 template<
class EquationOfState>
167 + EquationOfState::Sv(p, T);
171 template<
class EquationOfState>
178 c0_*(
pow(T, n0_ + 1) -
pow(
Tstd, n0_ + 1))/(
pow(Tref_, n0_)*(n0_ + 1))
184 template<
class EquationOfState>
198 template<
class EquationOfState>
199 inline void Foam::ePowerThermo<EquationOfState>::operator+=
204 scalar Y1 = this->
Y();
206 EquationOfState::operator+=(ct);
208 if (
mag(this->
Y()) > small)
211 const scalar Y2 = ct.Y()/this->
Y();
213 Hf_ = Y1*Hf_ + Y2*ct.Hf_;
214 c0_ = Y1*c0_ + Y2*ct.c0_;
215 n0_ = Y1*n0_ + Y2*ct.n0_;
216 Tref_ = Y1*Tref_ + Y2*ct.Tref_;
223 template<
class EquationOfState>
232 static_cast<const EquationOfState&>(ct1)
233 + static_cast<const EquationOfState&>(ct2)
236 if (
mag(eofs.Y()) < small)
252 ct1.Y()/eofs.Y()*ct1.c0_
253 + ct2.Y()/eofs.Y()*ct2.c0_,
254 ct1.Y()/eofs.Y()*ct1.n0_
255 + ct2.Y()/eofs.Y()*ct2.n0_,
256 ct1.Y()/eofs.Y()*ct1.Tref_
257 + ct2.Y()/eofs.Y()*ct2.Tref_,
258 ct1.Y()/eofs.Y()*ct1.Hf_
259 + ct2.Y()/eofs.Y()*ct2.Hf_
265 template<
class EquationOfState>
274 s*
static_cast<const EquationOfState&
>(ct),
283 template<
class EquationOfState>
292 static_cast<const EquationOfState&>(ct1)
293 == static_cast<const EquationOfState&>(ct2)
299 ct2.Y()/eofs.Y()*ct2.c0_
300 - ct1.Y()/eofs.Y()*ct1.c0_,
301 ct2.Y()/eofs.Y()*ct2.n0_
302 - ct1.Y()/eofs.Y()*ct1.n0_,
303 ct2.Y()/eofs.Y()*ct2.Tref_
304 - ct1.Y()/eofs.Y()*ct1.Tref_,
305 ct2.Y()/eofs.Y()*ct2.Hf_
306 - ct1.Y()/eofs.Y()*ct1.Hf_
scalar Ea(const scalar p, const scalar T) const
Absolute internal energy [J/kg].
scalar limit(const scalar T) const
Limit the temperature to be in the range Tlow_ to Thigh_.
A list of keyword definitions, which are a keyword followed by any number of values (e...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
scalar Es(const scalar p, const scalar T) const
const dimensionedScalar & Tstd
Standard temperature.
const volScalarField & Cv
scalar S(const scalar p, const scalar T) const
Entropy [J/kg/K].
Power-function based thermodynamics package templated on EquationOfState.
static autoPtr< ePowerThermo > New(const dictionary &dict)
Selector from dictionary.
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))
autoPtr< ePowerThermo > clone() const
Construct and return a clone.
A class for handling words, derived from string.
scalar dCpdT(const scalar p, const scalar T) const
Temperature derivative of heat capacity at constant pressure.
scalar Hf() const
Enthalpy of formation [J/kg].
errorManip< error > abort(error &err)
scalar Es(const scalar p, const scalar T) const
Sensible internal energy [J/kg].
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
PtrList< volScalarField > & Y
scalar Cv(const scalar p, const scalar T) const
Heat capacity at constant volume [J/kg/K].
dimensioned< scalar > mag(const dimensioned< Type > &)
const dimensionedScalar & Pstd
Standard pressure.
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
scalar Gstd(const scalar T) const
Gibbs free energy of the mixture in the standard state [J/kg].