63 #ifndef icoPolynomial_H 64 #define icoPolynomial_H 76 template<
class Specie,
int PolySize>
79 template<
class Specie,
int PolySize>
80 icoPolynomial<Specie, PolySize>
operator+
86 template<
class Specie,
int PolySize>
93 template<
class Specie,
int PolySize>
96 const icoPolynomial<Specie, PolySize>&,
97 const icoPolynomial<Specie, PolySize>&
100 template<
class Specie,
int PolySize>
104 const icoPolynomial<Specie, PolySize>&
112 template<
class Specie,
int PolySize=8>
152 return "icoPolynomial<" +
word(Specie::typeName_()) +
'>';
165 inline scalar
rho(scalar
p, scalar
T)
const;
168 inline scalar
H(
const scalar
p,
const scalar
T)
const;
171 inline scalar
Cp(scalar
p, scalar
T)
const;
174 inline scalar
E(
const scalar
p,
const scalar
T)
const;
177 inline scalar
Cv(scalar
p, scalar
T)
const;
180 inline scalar
Sp(
const scalar
p,
const scalar
T)
const;
183 inline scalar
Sv(
const scalar
p,
const scalar
T)
const;
186 inline scalar
psi(scalar
p, scalar
T)
const;
189 inline scalar
Z(scalar
p, scalar
T)
const;
192 inline scalar
CpMCv(scalar
p, scalar
T)
const;
198 void write(Ostream& os)
const;
215 friend icoPolynomial
operator* <Specie, PolySize>
221 friend icoPolynomial
operator== <Specie, PolySize>
223 const icoPolynomial&,
230 friend Ostream& operator<< <Specie, PolySize>
scalar H(const scalar p, const scalar T) const
Return enthalpy contribution [J/kg].
scalar psi(scalar p, scalar T) const
Return compressibility [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 Sp(const scalar p, const scalar T) const
Return entropy contribution to the integral of Cp/T [J/kg/K].
scalar E(const scalar p, const scalar T) const
Return internal energy contribution [J/kg].
Incompressible, polynomial form of equation of state, using a polynomial function for density...
static const bool incompressible
Is the equation of state is incompressible i.e. rho != f(p)
scalar Z(scalar p, scalar T) const
Return compression factor [].
scalar Cp(scalar p, scalar T) const
Return Cp contribution [J/(kg K].
static const bool isochoric
Is the equation of state is isochoric i.e. rho = const.
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.
void write(Ostream &os) const
Write to Ostream.
A class for handling words, derived from string.
scalar Cv(scalar p, scalar T) const
Return Cv contribution [J/(kg K].
An Ostream is an abstract base class for all output systems (streams, files, token lists...
void operator+=(const icoPolynomial &)
scalar Sv(const scalar p, const scalar T) const
Return entropy contribution to the integral of Cv/T [J/kg/K].
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
word name(const complex &)
Return a string representation of a complex.
Polynomial templated on size (order):
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].
static word typeName()
Return the instantiated type name.