32 template<
class Specie>
36 const coeffList& coeffs
46 template<
class Specie>
58 template<
class Specie>
66 template<
class Specie>
79 template<
class Specie>
82 return 1/(C_[0] + (C_[1] + C_[2]*
T - C_[4]*
p)*
T - C_[3]*
p);
86 template<
class Specie>
93 template<
class Specie>
100 template<
class Specie>
107 template<
class Specie>
114 template<
class Specie>
121 template<
class Specie>
128 template<
class Specie>
135 template<
class Specie>
144 template<
class Specie>
150 const scalar Y1 = this->
Y();
151 Specie::operator+=(
rp);
153 if (
mag(this->
Y()) > SMALL)
155 C_ = (Y1*C_ +
rp.Y()*
rp.C_)/this->
Y();
160 template<
class Specie>
163 Specie::operator*=(
s);
169 template<
class Specie>
172 const rPolynomial<Specie>& rp1,
173 const rPolynomial<Specie>& rp2
178 static_cast<const Specie&>(rp1)
179 + static_cast<const Specie&>(rp2)
182 if (
mag(sp.Y()) < SMALL)
184 return rPolynomial<Specie>
191 return rPolynomial<Specie>
194 (rp1.Y()*rp1.C_ + rp2.Y()*rp2.C_)/sp.Y()
199 template<
class Specie>
203 const rPolynomial<Specie>&
rp 206 return rPolynomial<Specie>
208 s*
static_cast<const Specie&
>(
rp),
214 template<
class Specie>
217 const rPolynomial<Specie>& rp1,
218 const rPolynomial<Specie>& rp2
223 static_cast<const Specie&>(rp1)
224 == static_cast<const Specie&>(rp2)
227 return rPolynomial<Specie>
230 (rp1.Y()*rp1.C_ - rp2.Y()*rp2.C_)/sp.Y()
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
A list of keyword definitions, which are a keyword followed by a number of values (eg...
scalar Z(scalar p, scalar T) const
Return compression factor [].
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Reciprocal polynomial equation of state for liquids and solids.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
scalar S(const scalar p, const scalar T) const
Return entropy [J/kg/K].
scalar Cv(scalar p, scalar T) const
Return Cv departure [J/(kg K].
autoPtr< rPolynomial > clone() const
Construct and return a clone.
scalar H(const scalar p, const scalar T) const
Return enthalpy departure [J/kg].
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
A class for handling words, derived from Foam::string.
regionProperties rp(runTime)
void operator*=(const scalar)
scalar psi(scalar p, scalar T) const
Return compressibility [s^2/m^2].
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
scalar CpMCv(scalar p, scalar T) const
Return (Cp - Cv) [J/(kg K].
PtrList< volScalarField > & Y
static autoPtr< rPolynomial > New(const dictionary &dict)
scalar rho(scalar p, scalar T) const
Return density [kg/m^3].
Pointer management similar to std::unique_ptr, with some additional methods and type checking...
scalar E(const scalar p, const scalar T) const
Return internal energy departure [J/kg].
scalar Cp(scalar p, scalar T) const
Return Cp departure [J/(kg K].
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))
rPolynomial(const Specie &sp, const coeffList &coeffs)
Construct from components.