26 #include "randomCoalescence.H" 34 namespace diameterModels
51 const dictionary&
dict 57 alphaMax_(
"alphaMax",
dimless, dict)
74 "randomCoalescence:R",
82 const scalar Crc = Crc_.value();
83 const scalar C = C_.value();
84 const scalar alphaMax = alphaMax_.value();
87 const scalar cbrtAlphaMax =
cbrt(alphaMax);
91 if (alpha[celli] < alphaMax - SMALL)
93 const scalar cbrtAlphaMaxMAlpha = cbrtAlphaMax -
cbrt(alpha[celli]);
96 12*
phi()*kappai[celli]*alpha[celli]
99 *(1 -
exp(-C*
cbrt(alpha[celli]*alphaMax)/cbrtAlphaMaxMAlpha))
100 /(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
static word timeName(const scalar, const int precision=precision_)
Return time name of given scalar time.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
DimensionedField< scalar, volMesh > Internal
Type of the internal field from which this GeometricField is derived.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
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 Time & time() const
Return time.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
A class for managing temporary objects.
Calculate the matrix for implicit and explicit sources.
virtual tmp< volScalarField > R() const =0