32 template<
class ChemistryModel>
36 coeffsDict_(this->subDict(
"odeCoeffs")),
44 template<
class ChemistryModel>
51 template<
class ChemistryModel>
63 if (odeSolver_->resize())
65 odeSolver_->resizeField(cTp_);
78 odeSolver_->solve(0, deltaT, cTp_, subDeltaT);
82 c[i] =
max(0.0, cTp_[i]);
An abstract base class for solving chemistry.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
virtual void solve(scalarField &c, scalar &T, scalar &p, scalar &deltaT, scalar &subDeltaT) const
Update the concentrations and return the chemical time.
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
ode(typename ChemistryModel::reactionThermo &thermo)
Construct from thermo.
Abstract base-class for ODE system solvers.
const dimensionedScalar c
Speed of light in a vacuum.
virtual ~ode()
Destructor.