38 #ifndef StandardChemistryModel_H 39 #define StandardChemistryModel_H 58 template<
class ReactionThermo,
class ThermoType>
68 template<
class DeltaTType>
69 scalar solve(
const DeltaTType& deltaT);
148 inline scalar
Treact()
const;
203 const label reactionI,
209 virtual scalar solve(
const scalar deltaT);
virtual label nSpecie() const
The number of species.
StandardChemistryModel(const ReactionThermo &thermo)
Construct from thermo.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Basic chemistry model templated on thermodynamics.
virtual tmp< volScalarField > Qdot() const
Return the heat release rate [kg/m/s^3].
Abstract base class for the systems of ordinary differential equations.
scalarField dcdt_
Temporary rate-of-change of concentration field.
virtual void jacobian(const scalar t, const scalarField &c, const label li, scalarField &dcdt, scalarSquareMatrix &J) const
Calculate the Jacobian of the system.
virtual void omega(const scalar p, const scalar T, const scalarField &c, const label li, scalarField &dcdt) const
dc/dt = omega, rate of change in concentration, for each species
label nSpecie_
Number of species.
PtrList< volScalarField::Internal > RR_
List of reaction rate per specie [kg/m^3/s].
const dimensionedScalar & c
Speed of light in a vacuum.
Extends base chemistry model by adding a thermo package, and ODE functions. Introduces chemistry equa...
const ReactionList< ThermoType > reactions_
Reactions.
void operator=(const StandardChemistryModel &)=delete
Disallow default bitwise assignment.
const PtrList< Reaction< ThermoType > > & reactions() const
The reactions.
virtual void calculate()
Calculates the reaction rates.
scalar Treact() const
Temperature below which the reaction rates are assumed 0.
scalar Treact_
Temperature below which the reaction rates are assumed 0.
const PtrList< ThermoType > & specieThermos_
Thermodynamic data of the species.
const PtrList< ThermoType > & specieThermos() const
Thermodynamic data of the species.
virtual scalar omegaI(label iReaction, const scalar p, const scalar T, const scalarField &c, const label li, scalar &pf, scalar &cf, label &lRef, scalar &pr, scalar &cr, label &rRef) const
Return the reaction rate for iReaction and the reference.
virtual ~StandardChemistryModel()
Destructor.
scalarField c_
Temporary concentration field.
label nReaction_
Number of reactions.
virtual void derivatives(const scalar t, const scalarField &c, const label li, scalarField &dcdt) const
Calculate the derivatives in dydx.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
PtrList< volScalarField::Internal > & RR()
Write access to chemical source terms.
virtual tmp< volScalarField::Internal > calculateRR(const label reactionI, const label speciei) const
Return reaction rate of the speciei in reactionI.
List of templated reactions.
virtual tmp< volScalarField > tc() const
Return the chemical time scale.
const ReactionThermo & thermo() const
Return const access to the thermo package.
virtual label nEqns() const
Number of ODE's to solve.
const PtrList< volScalarField > & Y_
Reference to the field of specie mass fractions.
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
A class for managing temporary objects.
TypeName("standard")
Runtime type information.
virtual label nReaction() const
The number of reactions.