26 #include "randomCoalescence.H" 33 namespace diameterModels
50 const dictionary&
dict 56 alphaMax_(
"alphaMax",
dimless, dict)
65 tmp<volScalarField> tR
82 scalar Crc = Crc_.value();
83 scalar C = C_.value();
84 scalar alphaMax = alphaMax_.value();
88 scalar cbrtAlphaMax =
cbrt(alphaMax);
92 if (alpha[celli] < alphaMax - SMALL)
94 scalar cbrtAlphaMaxMAlpha = cbrtAlphaMax -
cbrt(alpha[celli]);
97 (-12)*
phi()*kappai[celli]*alpha[celli]
100 *(1 -
exp(-C*
cbrt(alpha[celli]*alphaMax)/cbrtAlphaMaxMAlpha))
101 /(cbrtAlphaMax*cbrtAlphaMaxMAlpha);
#define forAll(list, i)
Loop across all elements in list.
const phaseModel & phase() const
virtual tmp< volScalarField > R() const
static word timeName(const scalar, const int precision=precision_)
Return time name of given scalar time.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensionedScalar exp(const dimensionedScalar &ds)
const Time & time() const
Return time.
dimensionedScalar cbrt(const dimensionedScalar &ds)
const volVectorField & U() const
tmp< volScalarField > Ut() const
Return the bubble turbulent velocity.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)
const IATE & iate_
Reference to the IATE this source applies to.
const volScalarField & kappai() const
Return the interfacial curvature.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
randomCoalescence(const IATE &iate, const dictionary &dict)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const Mesh & mesh() const
Return mesh.
const phaseModel & phase() const
Return the phase.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
A class for managing temporary objects.