31#include "phaseCompressibleTurbulenceModel.H"
39namespace diameterModels
41namespace binaryBreakupModels
64 gammaUpperReg2by11_(),
65 gammaUpperReg5by11_(),
66 gammaUpperReg8by11_(),
73 kolmogorovLengthScale_
77 "kolmogorovLengthScale",
84 "kolmogorovLengthScale",
98 for (scalar z = 1
e-2; z <= 10.0; z = z + 1
e-2)
118 gammaUpperReg2by11Table.
append(gamma2by11);
119 gammaUpperReg5by11Table.
append(gamma5by11);
120 gammaUpperReg8by11Table.
append(gamma8by11);
123 gammaUpperReg2by11_.reset
127 gammaUpperReg2by11Table,
133 gammaUpperReg5by11_.reset
137 gammaUpperReg5by11Table,
143 gammaUpperReg8by11_.reset
147 gammaUpperReg8by11Table,
159 kolmogorovLengthScale_ =
164 popBal_.continuousPhase().nu()
166 /popBal_.continuousTurbulence().epsilon()
179 const phaseModel& continuousPhase = popBal_.continuousPhase();
180 const sizeGroup& fi = popBal_.sizeGroups()[i];
181 const sizeGroup& fj = popBal_.sizeGroups()[j];
185 pow(fi.
x()/fj.
x(), 2.0/3.0) +
pow((1 - fi.
x()/fj.
x()), 2.0/3.0) - 1
190 12.0*cf*popBal_.sigmaWithContinuousPhase(fi.
phase())
192 beta_*continuousPhase.
rho()*
pow(fj.
d(), 5.0/3.0)
193 *
pow(popBal_.continuousTurbulence().epsilon(), 2.0/3.0)
197 const volScalarField xiMin(minEddyRatio_*kolmogorovLengthScale_/fj.
d());
206 2.0*
pow(
b[celli], 3.0/11.0)*tgamma(5.0/11.0)
208 gammaUpperReg5by11_()(
b[celli])
209 - gammaUpperReg5by11_()(tMin[celli])
214 C4_*(1 - popBal_.alphas())/fj.
x()
217 popBal_.continuousTurbulence().epsilon()
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
void append(const T &val)
Append an element at the end of the list.
A 2-tuple for storing two objects of dissimilar types. The container is similar in purpose to std::pa...
Base class for binary breakup models which give the breakup rate between a sizeGroup pair directly,...
Model of Luo and Svendsen (1996). The breakup rate is calculated by.
virtual void addToBinaryBreakupRate(volScalarField &binaryBreakupRate, const label i, const label j)
Add to binary breakupRate.
virtual void correct()
Correct diameter independent expressions.
Class that solves the univariate population balance equation by means of a class method (also called ...
This class represents a single sizeGroup belonging to a velocityGroup. The main property of a sizeGro...
const dimensionedScalar & d() const
Return representative diameter of the sizeGroup.
const dimensionedScalar & x() const
Return representative volume of the sizeGroup.
const phaseModel & phase() const
Return const-reference to the phase.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
An interpolation/look-up table of scalar vs <Type> values. The reference scalar values must be monoto...
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
const dimensionedScalar & rho() const
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
scalar incGammaRatio_Q(const scalar a, const scalar x)
Regularised upper incomplete gamma function.
@ CLAMP
Clamp value to the start/end value.
const dimensionSet dimless
Dimensionless.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar cbrt(const dimensionedScalar &ds)
dimensionedScalar pow025(const dimensionedScalar &ds)
#define forAll(list, i)
Loop across all elements in list.