Go to the documentation of this file.
36 template<
class CloudType>
44 Sb_(this->coeffDict().getScalar(
"Sb")),
45 D_(this->coeffDict().getScalar(
"D")),
47 O2GlobalId_(owner.composition().carrierId(
"O2")),
48 CO2GlobalId_(owner.composition().carrierId(
"CO2")),
54 label idSolid = owner.composition().idSolid();
55 CsLocalId_ = owner.composition().localId(idSolid,
"C");
58 WO2_ = owner.thermo().carrier().W(O2GlobalId_);
59 const scalar WCO2 = owner.thermo().carrier().W(CO2GlobalId_);
62 HcCO2_ = owner.thermo().carrier().Hc(CO2GlobalId_);
67 <<
"Stoichiometry of reaction, Sb, must be greater than zero" <<
nl
71 const scalar YCloc = owner.composition().Y0(idSolid)[CsLocalId_];
72 const scalar YSolidTot = owner.composition().YMixture0()[idSolid];
73 Info<<
" C(s): particle mass fraction = " << YCloc*YSolidTot <<
endl;
77 template<
class CloudType>
86 CsLocalId_(srm.CsLocalId_),
87 O2GlobalId_(srm.O2GlobalId_),
88 CO2GlobalId_(srm.CO2GlobalId_),
97 template<
class CloudType>
122 const label idSolid = CloudType::parcelType::SLD;
123 const scalar fComb = YMixture[idSolid]*YSolid[CsLocalId_];
134 const scalar YO2 =
thermo.carrier().Y(O2GlobalId_)[celli];
140 dmC =
min(mass*fComb, dmC);
143 const scalar dmO2 = dmC/WC_*Sb_*WO2_;
146 const scalar dmCO2 = dmC + dmO2;
149 dMassSolid[CsLocalId_] += dmC;
152 dMassSRCarrier[O2GlobalId_] -= dmO2;
153 dMassSRCarrier[CO2GlobalId_] += dmCO2;
155 const scalar HsC =
thermo.solids().properties()[CsLocalId_].Hs(
T);
160 return dmC*HsC - dmCO2*HcCO2_;
Templated surface reaction model class.
COxidationDiffusionLimitedRate(const dictionary &dict, CloudType &owner)
Construct from dictionary.
Thermo package for (S)olids (L)iquids and (G)ases Takes reference to thermo package,...
virtual scalar calculate(const scalar dt, const scalar Re, const scalar nu, const label celli, const scalar d, const scalar T, const scalar Tc, const scalar pc, const scalar rhoc, const scalar mass, const scalarField &YGas, const scalarField &YLiquid, const scalarField &YSolid, const scalarField &YMixture, const scalar N, scalarField &dMassGas, scalarField &dMassLiquid, scalarField &dMassSolid, scalarField &dMassSRCarrier) const
Update surface reactions.
Different types of constants.
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
Ostream & endl(Ostream &os)
Add newline and flush stream.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
messageStream Info
Information stream (stdout output on master, null elsewhere)
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
Templated base class for dsmc cloud.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Diffusion limited rate surface reaction model for coal parcels. Limited to:
errorManipArg< error, int > exit(error &err, const int errNo=1)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
constexpr scalar pi(M_PI)
scalarField Re(const UList< complex > &cf)
Extract real component.
const Vector< label > N(dict.get< Vector< label >>("N"))