34 namespace diameterModels
36 namespace binaryBreakupModels
59 gammaUpperReg2by11_(),
60 gammaUpperReg5by11_(),
61 gammaUpperReg8by11_(),
68 kolmogorovLengthScale_
72 "kolmogorovLengthScale",
73 popBal_.time().
name(),
79 "kolmogorovLengthScale",
93 for (scalar z = 1
e-2; z <= 10.0; z = z + 1
e-2)
113 gammaUpperReg2by11Table.
append(gamma2by11);
114 gammaUpperReg5by11Table.
append(gamma5by11);
115 gammaUpperReg8by11Table.
append(gamma8by11);
118 gammaUpperReg2by11_ =
123 linearInterpolationWeights::typeName,
125 gammaUpperReg2by11Table
128 gammaUpperReg5by11_ =
133 linearInterpolationWeights::typeName,
135 gammaUpperReg5by11Table
138 gammaUpperReg8by11_ =
143 linearInterpolationWeights::typeName,
145 gammaUpperReg8by11Table
154 kolmogorovLengthScale_ =
159 popBal_.continuousPhase().thermo().nu()
161 /popBal_.continuousTurbulence().epsilon()
174 const phaseModel& continuousPhase = popBal_.continuousPhase();
175 const sizeGroup& fi = popBal_.sizeGroups()[i];
176 const sizeGroup& fj = popBal_.sizeGroups()[j];
180 pow(fi.
x()/fj.
x(), 2.0/3.0) +
pow((1 - fi.
x()/fj.
x()), 2.0/3.0) - 1
185 12*cf*popBal_.sigmaWithContinuousPhase(fi.
phase())
187 beta_*continuousPhase.
rho()*
pow(fj.
dSph(), 5.0/3.0)
188 *
pow(popBal_.continuousTurbulence().epsilon(), 2.0/3.0)
201 2*
pow(
b[celli], 3.0/11.0)*tgamma(5.0/11.0)
203 gammaUpperReg5by11_->value(
b[celli])
204 - gammaUpperReg5by11_->value(tMin[celli])
209 C4_*(1 - popBal_.alphas())/fj.
x()
212 popBal_.continuousTurbulence().epsilon()
#define forAll(list, i)
Loop across all elements in list.
Macros for easy insertion into run-time selection tables.
Generic GeometricField class.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
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 &)
Append an element at the end of the list.
Base class to read table data for tables.
A 2-tuple for storing two objects of different types.
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Base class for binary breakup models that provide a breakup rate between a size class pair directly,...
Model of Luo and Svendsen (1996). The breakup rate is calculated by.
virtual void precompute()
Precompute diameter independent expressions.
virtual void addToBinaryBreakupRate(volScalarField &binaryBreakupRate, const label i, const label j)
Add to binary breakupRate.
LuoSvendsen(const populationBalanceModel &popBal, const dictionary &dict)
Model for tracking the evolution of a dispersed phase size distribution due to coalescence (synonymou...
Single size class fraction field representing a fixed particle volume as defined by the user through ...
const dimensionedScalar & dSph() const
Return representative spherical 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 any number of values (e....
virtual const volScalarField & rho() const =0
Return the density field.
defineTypeNameAndDebug(LehrMilliesMewes, 0)
addToRunTimeSelectionTable(binaryBreakupModel, LehrMilliesMewes, dictionary)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimless
const dimensionSet dimLength
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar cbrt(const dimensionedScalar &ds)
word name(const complex &)
Return a string representation of a complex.
scalar incGammaRatio_Q(const scalar a, const scalar x)
Normalised upper incomplete gamma function.
dimensionedScalar pow025(const dimensionedScalar &ds)