35 namespace diameterModels
37 namespace nucleationModels
56 const populationBalanceModel& popBal,
57 const dictionary&
dict 60 nucleationModel(popBal, dict),
65 popBal.
mesh().lookupObject<phaseModel>
70 dict.
lookup(
"velocityGroup")
88 const volScalarField::Boundary& alphatBf = alphat.boundaryField();
90 typedef compressible::alphatWallBoilingWallFunctionFvPatchScalarField
91 alphatWallBoilingWallFunction;
97 isA<alphatWallBoilingWallFunction>(alphatBf[patchi])
100 const alphatWallBoilingWallFunction& alphatw =
101 refCast<const alphatWallBoilingWallFunction>(alphatBf[
patchi]);
105 if (
min(dDep) < velGroup_.sizeGroups().first().dSph().value())
108 <<
"Minimum departure diameter " <<
min(dDep)
109 <<
" m outside of range [" 110 << velGroup_.sizeGroups().first().dSph().value() <<
", " 111 << velGroup_.sizeGroups().last().dSph().value() <<
"] m" 112 <<
" at patch " << alphatw.patch().name()
114 <<
" The nucleation rate in populationBalance " 116 <<
" Adjust discretization over property space to" 117 <<
" suppress this warning." 120 else if (
max(dDep) > velGroup_.sizeGroups().last().dSph().value())
123 <<
"Maximum departure diameter " <<
max(dDep)
124 <<
" m outside of range [" 125 << velGroup_.sizeGroups().first().dSph().value() <<
", " 126 << velGroup_.sizeGroups().last().dSph().value() <<
"] m" 127 <<
" at patch " << alphatw.patch().name()
129 <<
" The nucleation rate in populationBalance " 131 <<
" Adjust discretization over property space to" 132 <<
" suppress this warning." 148 const phaseModel& phase = fi.phase();
157 const volScalarField::Boundary& alphatBf = alphat.boundaryField();
159 typedef compressible::alphatWallBoilingWallFunctionFvPatchScalarField
160 alphatWallBoilingWallFunction;
166 isA<alphatWallBoilingWallFunction>(alphatBf[patchi])
169 const alphatWallBoilingWallFunction& alphatw =
170 refCast<const alphatWallBoilingWallFunction>(alphatBf[
patchi]);
175 const labelList& faceCells = alphatw.patch().faceCells();
181 if (dmdt[facei] > small)
183 const label faceCelli = faceCells[facei];
185 nucleationRate[faceCelli] +=
189 fi.x()/
pow3(fi.dSph())*
pow3(dDep[facei]*unitLength)
191 *dmdt[facei]/rho[faceCelli]/fi.x().
value();
virtual void correct()
Correct diameter independent expressions.
#define forAll(list, i)
Loop across all elements in list.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
const word & name() const
Return name.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
To & refCast(From &r)
Reference type cast template function.
Ostream & endl(Ostream &os)
Add newline and flush stream.
const word & name() const
const populationBalanceModel & popBal_
Reference to the populationBalanceModel.
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
Macros for easy insertion into run-time selection tables.
const dimensionedScalar eta(const label i, const dimensionedScalar &v) const
Return allocation coefficient.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
stressControl lookup("compactNormalStress") >> compactNormalStress
const phaseModel & continuousPhase() const
Return continuous phase.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
static word groupName(Name name, const word &group)
const Type & value() const
Return const reference to value.
List< label > labelList
A List of labels.
virtual void addToNucleationRate(volScalarField &nucleationRate, const label i)
Add to nucleationRate.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)
const UPtrList< sizeGroup > & sizeGroups() const
Return the sizeGroups belonging to this populationBalance.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
wallBoiling(const populationBalanceModel &popBal, const dictionary &dict)
const fvMesh & mesh() const
Return reference to the mesh.