35 namespace diameterModels
37 namespace binaryBreakupModels
55 const populationBalanceModel& popBal,
56 const dictionary& dict
59 binaryBreakupModel(popBal, dict),
60 gammaUpperReg2by11_(),
61 gammaUpperReg5by11_(),
62 gammaUpperReg8by11_(),
69 kolmogorovLengthScale_
73 "kolmogorovLengthScale",
80 "kolmogorovLengthScale",
86 List<Tuple2<scalar, scalar>> gammaUpperReg2by11Table;
87 List<Tuple2<scalar, scalar>> gammaUpperReg5by11Table;
88 List<Tuple2<scalar, scalar>> gammaUpperReg8by11Table;
90 gammaUpperReg2by11Table.append(Tuple2<scalar, scalar>(0.0, 1.0));
91 gammaUpperReg5by11Table.append(Tuple2<scalar, scalar>(0.0, 1.0));
92 gammaUpperReg8by11Table.append(Tuple2<scalar, scalar>(0.0, 1.0));
94 for (scalar z = 1
e-2; z <= 10.0; z = z + 1
e-2)
96 Tuple2<scalar, scalar> gamma2by11
102 Tuple2<scalar, scalar> gamma5by11
108 Tuple2<scalar, scalar> gamma8by11
114 gammaUpperReg2by11Table.append(gamma2by11);
115 gammaUpperReg5by11Table.append(gamma5by11);
116 gammaUpperReg8by11Table.append(gamma8by11);
119 gammaUpperReg2by11_ =
120 new Function1s::Table<scalar>
124 linearInterpolationWeights::typeName,
125 gammaUpperReg2by11Table
128 gammaUpperReg5by11_ =
129 new Function1s::Table<scalar>
133 linearInterpolationWeights::typeName,
134 gammaUpperReg5by11Table
137 gammaUpperReg8by11_ =
138 new Function1s::Table<scalar>
142 linearInterpolationWeights::typeName,
143 gammaUpperReg8by11Table
152 kolmogorovLengthScale_ =
178 pow(fi.x()/fj.x(), 2.0/3.0) +
pow((1 - fi.x()/fj.x()), 2.0/3.0) - 1
185 beta_*continuousPhase.rho()*
pow(fj.dSph(), 5.0/3.0)
190 const volScalarField xiMin(minEddyRatio_*kolmogorovLengthScale_/fj.dSph());
199 2*
pow(
b[celli], 3.0/11.0)*tgamma(5.0/11.0)
201 gammaUpperReg5by11_->value(
b[celli])
202 - gammaUpperReg5by11_->value(tMin[celli])
#define forAll(list, i)
Loop across all elements in list.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
LuoSvendsen(const populationBalanceModel &popBal, const dictionary &dict)
virtual void precompute()
Precompute diameter independent expressions.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const phaseCompressible::momentumTransportModel & continuousTurbulence() const
Return reference to turbulence model of the continuous phase.
const dimensionedScalar b
Wien displacement law constant: default SI units: [m K].
dimensionedScalar pow025(const dimensionedScalar &ds)
const dimensionSet dimless
Macros for easy insertion into run-time selection tables.
const dimensionSet dimLength
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual void addToBinaryBreakupRate(volScalarField &binaryBreakupRate, const label i, const label j)
Add to binary breakupRate.
const phaseModel & continuousPhase() const
Return continuous phase.
const rhoThermo & thermo() const
Return const-access to phase rhoThermo.
dimensionedScalar cbrt(const dimensionedScalar &ds)
scalar incGammaRatio_Q(const scalar a, const scalar x)
Normalised upper incomplete gamma function.
const populationBalanceModel & popBal_
Reference to the populationBalanceModel.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)
const UPtrList< sizeGroup > & sizeGroups() const
Return the sizeGroups belonging to this populationBalance.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const volScalarField & alphas() const
Return total void of phases belonging to this populationBalance.
virtual tmp< volScalarField > nu() const
Kinematic viscosity of mixture [m^2/s].
const tmp< volScalarField > sigmaWithContinuousPhase(const phaseModel &dispersedPhase) const
Return the surface tension coefficient between a given dispersed.
const doubleScalar e
Elementary charge.
virtual tmp< volScalarField > epsilon() const =0
Return the turbulence kinetic energy dissipation rate.