26 #include "randomCoalescence.H" 33 namespace diameterModels
50 const dictionary&
dict 56 alphaMax_(
"alphaMax",
dimless, dict)
65 tmp<volScalarField> tR
77 scalar Crc = Crc_.value();
78 scalar C = C_.value();
79 scalar alphaMax = alphaMax_.value();
83 scalar cbrtAlphaMax =
cbrt(alphaMax);
87 if (alpha[celli] < alphaMax - small)
89 scalar cbrtAlphaMaxMAlpha = cbrtAlphaMax -
cbrt(alpha[celli]);
92 (-12)*
phi()*kappai[celli]*alpha[celli]
95 *(1 -
exp(-C*
cbrt(alpha[celli]*alphaMax)/cbrtAlphaMaxMAlpha))
96 /(cbrtAlphaMax*cbrtAlphaMaxMAlpha);
tmp< volScalarField > Ut() const
Return the bubble turbulent velocity.
#define forAll(list, i)
Loop across all elements in list.
virtual tmp< volScalarField > R() const
const volScalarField & kappai() const
Return the interfacial curvature.
Macros for easy insertion into run-time selection tables.
static tmp< GeometricField< scalar, fvPatchField, volMesh > > New(const word &name, const Mesh &, const dimensionSet &, const word &patchFieldType=fvPatchField< scalar >::calculatedType())
Return a temporary field constructed from name, mesh, dimensionSet.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const volVectorField & U() const
const phaseModel & phase() const
Return the phase.
dimensionedScalar exp(const dimensionedScalar &ds)
dimensionedScalar cbrt(const dimensionedScalar &ds)
const phaseModel & phase() const
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
const Mesh & mesh() const
Return mesh.
defineTypeNameAndDebug(combustionModel, 0)
const IATE & iate_
Reference to the IATE this source applies to.
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 dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
A class for managing temporary objects.