35 namespace diameterModels
37 namespace binaryBreakupModels
60 gammaUpperReg2by11_(),
61 gammaUpperReg5by11_(),
62 gammaUpperReg8by11_(),
69 kolmogorovLengthScale_
73 "kolmogorovLengthScale",
74 popBal_.time().
name(),
80 "kolmogorovLengthScale",
94 for (scalar z = 1
e-2; z <= 10.0; z = z + 1
e-2)
114 gammaUpperReg2by11Table.
append(gamma2by11);
115 gammaUpperReg5by11Table.
append(gamma5by11);
116 gammaUpperReg8by11Table.
append(gamma8by11);
119 gammaUpperReg2by11_ =
124 linearInterpolationWeights::typeName,
126 gammaUpperReg2by11Table
129 gammaUpperReg5by11_ =
134 linearInterpolationWeights::typeName,
136 gammaUpperReg5by11Table
139 gammaUpperReg8by11_ =
144 linearInterpolationWeights::typeName,
146 gammaUpperReg8by11Table
155 kolmogorovLengthScale_ =
160 popBal_.continuousPhase().fluidThermo().nu()
162 /popBal_.continuousTurbulence().epsilon()
175 const phaseModel& continuousPhase = popBal_.continuousPhase();
176 const sizeGroup& fi = popBal_.sizeGroups()[i];
177 const sizeGroup& fj = popBal_.sizeGroups()[j];
181 pow(fi.
x()/fj.
x(), 2.0/3.0) +
pow((1 - fi.
x()/fj.
x()), 2.0/3.0) - 1
186 12*cf*popBal_.sigmaWithContinuousPhase(fi.
phase())
188 beta_*continuousPhase.
rho()*
pow(fj.
dSph(), 5.0/3.0)
189 *
pow(popBal_.continuousTurbulence().epsilon(), 2.0/3.0)
202 2*
pow(
b[celli], 3.0/11.0)*tgamma(5.0/11.0)
204 gammaUpperReg5by11_->value(
b[celli])
205 - gammaUpperReg5by11_->value(tMin[celli])
210 C4_*(1 - popBal_.alphas())/fj.
x()
213 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.
Macros for creating standard TableReader-s.
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)
word name(const bool)
Return a word representation of a bool.
const dimensionSet dimless
const dimensionSet dimLength
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar cbrt(const dimensionedScalar &ds)
scalar incGammaRatio_Q(const scalar a, const scalar x)
Normalised upper incomplete gamma function.
dimensionedScalar pow025(const dimensionedScalar &ds)