26 #include "reactionDriven.H" 37 namespace diameterModels
39 namespace nucleationModels
60 const populationBalanceModel& popBal,
61 const dictionary& dict
64 nucleationModel(popBal, dict),
65 dNuc_(
"nucleationDiameter",
dimLength, dict),
70 popBal.
mesh().lookupObject<phaseModel>
75 dict.
lookup(
"velocityGroup")
82 popBal_.
mesh().lookupObject<phaseModel>
84 IOobject::groupName(
"alpha", dict.
lookup(
"reactingPhase"))
89 popBal_.
fluid().phasePairs()
91 phasePair(velGroup_.phase(), reactingPhase_)
94 dmdtfName_(dict.
lookup(
"dmdtf")),
95 specieName_(dict.
lookup(
"specie"))
99 dNuc_.value() < velGroup_.sizeGroups().first().dSph().value()
100 || dNuc_.value() > velGroup_.sizeGroups().last().dSph().value()
104 <<
"Nucleation diameter " << dNuc_.value() <<
"m outside of range [" 105 << velGroup_.sizeGroups().first().dSph().value() <<
", " 106 << velGroup_.sizeGroups().last().dSph().value() <<
"]." <<
nl 121 const sizeGroup& fi = popBal_.sizeGroups()[i];
138 const scalar dmidtfSign =
139 velGroup_.phase().name() ==
pair_.
first() ? +1 : -1;
142 popBal_.eta(i,
pi/6.0*
pow3(dNuc_))*dmidtfSign*dmidtf/rho/fi.x();
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
errorManipArg< error, int > exit(error &err, const int errNo=1)
reactionDriven(const populationBalanceModel &popBal, const dictionary &dict)
const phasePair & pair_
Phase pair.
To & refCast(From &r)
Reference type cast template function.
virtual word name() const
Pair name.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Calculate the first temporal derivative.
stressControl lookup("compactNormalStress") >> compactNormalStress
static word groupName(Name name, const word &group)
Calculate the matrix for the first temporal derivative.
virtual void addToNucleationRate(volScalarField &nucleationRate, const label i)
Add to nucleationRate.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
const Mesh & mesh() const
Return mesh.
defineTypeNameAndDebug(combustionModel, 0)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
const Type & first() const
Return first.