30 template<
class Specie,
int PolySize>
44 template<
class Specie,
int PolySize>
52 rhoCoeffs_(ip.rhoCoeffs_)
56 template<
class Specie,
int PolySize>
67 template<
class Specie,
int PolySize>
80 template<
class Specie,
int PolySize>
87 return rhoCoeffs_.value(T);
91 template<
class Specie,
int PolySize>
98 return p/this->
rho(p, T);
102 template<
class Specie,
int PolySize>
113 template<
class Specie,
int PolySize>
124 template<
class Specie,
int PolySize>
135 template<
class Specie,
int PolySize>
146 template<
class Specie,
int PolySize>
157 template<
class Specie,
int PolySize>
168 template<
class Specie,
int PolySize>
181 template<
class Specie,
int PolySize>
182 inline void Foam::icoPolynomial<Specie, PolySize>::operator+=
187 scalar Y1 = this->
Y();
188 Specie::operator+=(ip);
190 if (
mag(this->
Y()) > small)
193 const scalar Y2 = ip.Y()/this->
Y();
195 rhoCoeffs_ = Y1*rhoCoeffs_ + Y2*ip.rhoCoeffs_;
200 template<
class Specie,
int PolySize>
203 Specie::operator*=(s);
209 template<
class Specie,
int PolySize>
218 static_cast<const Specie&>(ip1)
219 + static_cast<const Specie&>(ip2)
222 if (
mag(sp.Y()) < small)
232 const scalar Y1 = ip1.Y()/sp.Y();
233 const scalar Y2 = ip2.Y()/sp.Y();
238 Y1*ip1.rhoCoeffs_ + Y2*ip2.rhoCoeffs_
244 template<
class Specie,
int PolySize>
253 s*
static_cast<const Specie&
>(ip),
259 template<
class Specie,
int PolySize>
268 static_cast<const Specie&>(ip1)
269 == static_cast<const Specie&>(ip2)
272 const scalar Y1 = ip1.Y()/sp.Y();
273 const scalar Y2 = ip2.Y()/sp.Y();
278 Y2*ip2.rhoCoeffs_ - Y1*ip1.rhoCoeffs_
scalar H(const scalar p, const scalar T) const
Return enthalpy departure [J/kg].
scalar psi(scalar p, scalar T) const
Return compressibility rho/p [s^2/m^2].
autoPtr< icoPolynomial > clone() const
Construct and return a clone.
A list of keyword definitions, which are a keyword followed by any number of values (e...
scalar E(const scalar p, const scalar T) const
Return internal energy departure [J/kg].
Incompressible, polynomial form of equation of state, using a polynomial function for density...
scalar Z(scalar p, scalar T) const
Return compression factor [].
scalar Cp(scalar p, scalar T) const
Return Cp departure [J/(kg K].
static autoPtr< icoPolynomial > New(const dictionary &dict)
void operator*=(const scalar)
scalar CpMCv(scalar p, scalar T) const
Return (Cp - Cv) [J/(kg K].
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))
icoPolynomial(const Specie &sp, const Polynomial< PolySize > &rhoPoly)
Construct from components.
A class for handling words, derived from string.
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].
Polynomial templated on size (order):
PtrList< volScalarField > & Y
dimensioned< scalar > mag(const dimensioned< Type > &)
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
scalar rho(scalar p, scalar T) const
Return density [kg/m^3].