Go to the documentation of this file.
30 #include "phaseCompressibleTurbulenceModel.H"
37 namespace diameterModels
39 namespace binaryBreakupModels
62 gammaUpperReg2by11_(),
63 gammaUpperReg5by11_(),
64 gammaUpperReg8by11_(),
65 C4_(dimensionedScalar::lookupOrDefault(
"C4",
dict,
dimless, 0.923)),
66 beta_(dimensionedScalar::lookupOrDefault(
"beta",
dict,
dimless, 2.05)),
69 dimensionedScalar::lookupOrDefault(
"minEddyRatio",
dict,
dimless, 11.4)
71 kolmogorovLengthScale_
75 "kolmogorovLengthScale",
76 popBal_.time().timeName(),
82 "kolmogorovLengthScale",
96 for (scalar z = 1
e-2; z <= 10.0; z = z + 1
e-2)
116 gammaUpperReg2by11Table.
append(gamma2by11);
117 gammaUpperReg5by11Table.
append(gamma5by11);
118 gammaUpperReg8by11Table.
append(gamma8by11);
121 gammaUpperReg2by11_.reset
125 gammaUpperReg2by11Table,
126 bounds::repeatableBounding::CLAMP,
131 gammaUpperReg5by11_.reset
135 gammaUpperReg5by11Table,
136 bounds::repeatableBounding::CLAMP,
141 gammaUpperReg8by11_.reset
145 gammaUpperReg8by11Table,
146 bounds::repeatableBounding::CLAMP,
157 kolmogorovLengthScale_ =
177 const phaseModel& continuousPhase = popBal_.continuousPhase();
178 const sizeGroup& fi = popBal_.sizeGroups()[i];
179 const sizeGroup& fj = popBal_.sizeGroups()[j];
183 pow(fi.
x()/fj.
x(), 2.0/3.0) +
pow((1 - fi.
x()/fj.
x()), 2.0/3.0) - 1
188 12.0*cf*popBal_.sigmaWithContinuousPhase(fi.
phase())
190 beta_*continuousPhase.
rho()*
pow(fj.
d(), 5.0/3.0)
191 *
pow(popBal_.continuousTurbulence().epsilon(), 2.0/3.0)
195 const volScalarField xiMin(minEddyRatio_*kolmogorovLengthScale_/fj.
d());
204 2.0*
pow(
b[celli], 3.0/11.0)*tgamma(5.0/11.0)
206 gammaUpperReg5by11_()(
b[celli])
207 - gammaUpperReg5by11_()(tMin[celli])
212 C4_*(1 - popBal_.alphas())/fj.
x()
215 popBal_.continuousTurbulence().epsilon()
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
Defines the attributes of an object for which implicit objectRegistry management is supported,...
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
Dimensionless.
const phaseModel & phase() const
Return const-reference to the phase.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
static constexpr const zero Zero
Global zero.
virtual tmp< volScalarField > rho() const =0
Return the density field.
virtual void addToBinaryBreakupRate(volScalarField &binaryBreakupRate, const label i, const label j)
Add to binary breakupRate.
addToRunTimeSelectionTable(binaryBreakupModel, LehrMilliesMewes, dictionary)
void append(const T &val)
Append an element at the end of the list.
scalar incGammaRatio_Q(const scalar a, const scalar x)
Normalized upper incomplete gamma function.
#define forAll(list, i)
Loop across all elements in list.
dimensionedScalar pow025(const dimensionedScalar &ds)
virtual tmp< volScalarField > nu() const =0
Return the laminar kinematic viscosity.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
This class represents a single sizeGroup belonging to a velocityGroup. The main property of a sizeGro...
LuoSvendsen(const populationBalanceModel &popBal, const dictionary &dict)
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Class that solves the univariate population balance equation by means of a class method (also called ...
Macros for easy insertion into run-time selection tables.
const phaseModel & continuousPhase() const
Return continuous phase.
defineTypeNameAndDebug(LehrMilliesMewes, 0)
const dimensionedScalar & d() const
Return representative diameter of the sizeGroup.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const populationBalanceModel & popBal_
Reference to the populationBalanceModel.
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
Base class for binary breakup models which give the breakup rate between a sizeGroup pair directly,...
const phaseCompressibleTurbulenceModel & continuousTurbulence() const
Return reference to turbulence model of the continuous phase.
const dimensionedScalar e
Elementary charge.
dimensionedScalar cbrt(const dimensionedScalar &ds)
const dimensionedScalar & x() const
Return representative volume of the sizeGroup.
virtual void correct()
Correct diameter independent expressions.