36 #ifndef adiabaticPerfectFluid_H 37 #define adiabaticPerfectFluid_H 50 template<
class Specie>
54 const adiabaticPerfectFluid<Specie>&
57 template<
class Specie>
58 inline adiabaticPerfectFluid<Specie>
operator*
61 const adiabaticPerfectFluid<Specie>&
64 template<
class Specie>
65 inline adiabaticPerfectFluid<Specie>
operator==
67 const adiabaticPerfectFluid<Specie>&,
68 const adiabaticPerfectFluid<Specie>&
71 template<
class Specie>
75 const adiabaticPerfectFluid<Specie>&
83 template<
class Specie>
142 return "adiabaticPerfectFluid<" +
word(Specie::typeName_()) +
'>';
155 inline scalar
rho(scalar
p, scalar
T)
const;
158 inline scalar
H(
const scalar
p,
const scalar
T)
const;
161 inline scalar
Cp(scalar
p, scalar
T)
const;
164 inline scalar
S(
const scalar
p,
const scalar
T)
const;
167 inline scalar
psi(scalar
p, scalar
T)
const;
170 inline scalar
Z(scalar
p, scalar
T)
const;
173 inline scalar
CpMCv(scalar
p, scalar
T)
const;
179 void write(Ostream& os)
const;
193 const adiabaticPerfectFluid&
196 friend adiabaticPerfectFluid
operator* <Specie>
199 const adiabaticPerfectFluid&
202 friend adiabaticPerfectFluid
operator== <Specie>
204 const adiabaticPerfectFluid&,
205 const adiabaticPerfectFluid&
211 friend Ostream& operator<< <Specie>
214 const adiabaticPerfectFluid&
scalar Z(scalar p, scalar T) const
Return compression factor [].
static autoPtr< adiabaticPerfectFluid > New(const dictionary &dict)
A list of keyword definitions, which are a keyword followed by any number of values (e...
adiabaticPerfectFluid(const Specie &sp, const scalar p0, const scalar rho0, const scalar gamma, const scalar B)
Construct from components.
static const bool incompressible
Is the equation of state is incompressible i.e. rho != f(p)
scalar psi(scalar p, scalar T) const
Return compressibility rho/p [s^2/m^2].
static word typeName()
Return the instantiated type name.
static const bool isochoric
Is the equation of state is isochoric i.e. rho = const.
void operator+=(const adiabaticPerfectFluid &)
scalar S(const scalar p, const scalar T) const
Return entropy [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))
scalar Cp(scalar p, scalar T) const
Return Cp departure [J/(kg K].
A class for handling words, derived from string.
scalar CpMCv(scalar p, scalar T) const
Return (Cp - Cv) [J/(kg K].
scalar H(const scalar p, const scalar T) const
Return enthalpy departure [J/kg].
An Ostream is an abstract base class for all output systems (streams, files, token lists...
autoPtr< adiabaticPerfectFluid > clone() const
Construct and return a clone.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
word name(const complex &)
Return a string representation of a complex.
scalar rho(scalar p, scalar T) const
Return density [kg/m^3].
Adiabatic perfect fluid equation of state.
void operator*=(const scalar)
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
void write(Ostream &os) const
Write to Ostream.