38 #ifndef chemistryModel_H 39 #define chemistryModel_H 58 template<
class CompType,
class ThermoType>
68 template<
class DeltaTType>
69 scalar solve(
const DeltaTType& deltaT);
150 inline scalar
Treact()
const;
219 const label reactionI,
225 virtual scalar solve(
const scalar deltaT);
virtual tmp< volScalarField::Internal > calculateRR(const label reactionI, const label speciei) const
Return reaction rate of the speciei in reactionI.
virtual scalar omegaI(label iReaction, const scalarField &c, const scalar T, const scalar p, scalar &pf, scalar &cf, label &lRef, scalar &pr, scalar &cr, label &rRef) const
Return the reaction rate for iReaction and the reference.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
PtrList< volScalarField > & Y_
Reference to the field of specie mass fractions.
Abstract base class for the systems of ordinary differential equations.
virtual tmp< volScalarField > tc() const
Return the chemical time scale.
scalar Treact_
Temperature below which the reaction rates are assumed 0.
virtual void jacobian(const scalar t, const scalarField &c, scalarField &dcdt, scalarSquareMatrix &dfdc) const
Calculate the Jacobian of the system.
PtrList< volScalarField::Internal > & RR()
Write access to chemical source terms.
const PtrList< ThermoType > & specieThermo() const
Thermodynamic data of the species.
virtual ~chemistryModel()
Destructor.
Simple extension of ReactionThermo to handle reaction kinetics in addition to the equilibrium thermod...
scalarField dcdt_
Temporary rate-of-change of concentration field.
virtual label nSpecie() const
The number of species.
label nReaction_
Number of reactions.
label nSpecie_
Number of species.
A class for handling words, derived from string.
scalarField c_
Temporary concentration field.
scalar Treact() const
Temperature below which the reaction rates are assumed 0.
TypeName("chemistryModel")
Runtime type information.
virtual label nEqns() const
Number of ODE's to solve.
const PtrList< ThermoType > & specieThermo_
Thermodynamic data of the species.
virtual label nReaction() const
The number of reactions.
PtrList< volScalarField::Internal > RR_
List of reaction rate per specie [kg/m3/s].
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
virtual void calculate()
Calculates the reaction rates.
Extends base chemistry model by adding a thermo package, and ODE functions. Introduces chemistry equa...
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
Mesh data needed to do the Finite Volume discretisation.
const dimensionedScalar c
Speed of light in a vacuum.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
virtual void derivatives(const scalar t, const scalarField &c, scalarField &dcdt) const
Calculate the derivatives in dydx.
virtual tmp< volScalarField > Qdot() const
Return the heat release rate [kg/m/s3].
virtual void omega(const scalarField &c, const scalar T, const scalar p, scalarField &dcdt) const
dc/dt = omega, rate of change in concentration, for each species
A class for managing temporary objects.
const PtrList< Reaction< ThermoType > > & reactions() const
The reactions.
const PtrList< Reaction< ThermoType > > & reactions_
Reactions.