28 template<
class EquationOfState>
31 const EquationOfState& st,
48 template<
class EquationOfState>
55 EquationOfState(name, ct),
63 template<
class EquationOfState>
74 template<
class EquationOfState>
85 template<
class EquationOfState>
98 template<
class EquationOfState>
108 template<
class EquationOfState>
119 template<
class EquationOfState>
122 const scalar p,
const scalar T
129 template<
class EquationOfState>
132 const scalar p,
const scalar T
139 template<
class EquationOfState>
146 template<
class EquationOfState>
149 const scalar p,
const scalar T
158 template<
class EquationOfState>
159 inline void Foam::hRefConstThermo<EquationOfState>::operator+=
164 scalar molr1 = this->nMoles();
166 EquationOfState::operator+=(ct);
168 molr1 /= this->nMoles();
169 scalar molr2 = ct.nMoles()/this->nMoles();
171 Cp_ = molr1*Cp_ + molr2*ct.Cp_;
172 Hf_ = molr1*Hf_ + molr2*ct.Hf_;
176 template<
class EquationOfState>
177 inline void Foam::hRefConstThermo<EquationOfState>::operator-=
182 scalar molr1 = this->nMoles();
184 EquationOfState::operator-=(ct);
186 molr1 /= this->nMoles();
187 scalar molr2 = ct.nMoles()/this->nMoles();
189 Cp_ = molr1*Cp_ - molr2*ct.Cp_;
190 Hf_ = molr1*Hf_ - molr2*ct.Hf_;
196 template<
class EquationOfState>
205 static_cast<const EquationOfState&>(ct1)
206 + static_cast<const EquationOfState&>(ct2)
212 ct1.nMoles()/eofs.nMoles()*ct1.Cp_
213 + ct2.nMoles()/eofs.nMoles()*ct2.Cp_,
214 ct1.nMoles()/eofs.nMoles()*ct1.Hf_
215 + ct2.nMoles()/eofs.nMoles()*ct2.Hf_,
216 ct1.nMoles()/eofs.nMoles()*ct1.Tref_
217 + ct2.nMoles()/eofs.nMoles()*ct2.Tref_,
218 ct1.nMoles()/eofs.nMoles()*ct1.Href_
219 + ct2.nMoles()/eofs.nMoles()*ct2.Href_
224 template<
class EquationOfState>
233 static_cast<const EquationOfState&>(ct1)
234 - static_cast<const EquationOfState&>(ct2)
240 ct1.nMoles()/eofs.nMoles()*ct1.Cp_
241 - ct2.nMoles()/eofs.nMoles()*ct2.Cp_,
242 ct1.nMoles()/eofs.nMoles()*ct1.Hf_
243 - ct2.nMoles()/eofs.nMoles()*ct2.Hf_
248 template<
class EquationOfState>
257 s*
static_cast<const EquationOfState&
>(ct),
266 template<
class EquationOfState>
dimensionedScalar log(const dimensionedScalar &ds)
scalar s(const scalar p, const scalar T) const
Entropy [J/(kmol K)].
A list of keyword definitions, which are a keyword followed by any number of values (e...
static autoPtr< hRefConstThermo > New(Istream &is)
Selector from Istream.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
scalar hs(const scalar p, const scalar T) const
Sensible enthalpy [J/kmol].
scalar hc() const
Chemical enthalpy [J/kmol].
scalar ha(const scalar p, const scalar T) const
Absolute Enthalpy [J/kmol].
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))
A class for handling words, derived from string.
const dimensionedScalar Tstd
Standard temperature.
scalar limit(const scalar T) const
Limit the temperature to be in the range Tlow_ to Thigh_.
const volScalarField & cp
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
const dimensionedScalar h
Planck constant.
autoPtr< hRefConstThermo > clone() const
Construct and return a clone.
scalar cp(const scalar p, const scalar T) const
Heat capacity at constant pressure [J/(kmol K)].
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Constant properties thermodynamics package templated into the EquationOfState.