28 #include "phaseCompressibleTurbulenceModel.H" 34 namespace diameterModels
36 namespace binaryBreakupModels
54 const populationBalanceModel& popBal,
55 const dictionary& dict
58 binaryBreakupModel(popBal, dict),
59 gammaUpperReg2by11_(),
60 gammaUpperReg5by11_(),
61 gammaUpperReg8by11_(),
68 kolmogorovLengthScale_
72 "kolmogorovLengthScale",
79 "kolmogorovLengthScale",
85 List<Tuple2<scalar, scalar>> gammaUpperReg2by11Table;
86 List<Tuple2<scalar, scalar>> gammaUpperReg5by11Table;
87 List<Tuple2<scalar, scalar>> gammaUpperReg8by11Table;
89 gammaUpperReg2by11Table.append(Tuple2<scalar, scalar>(0.0, 1.0));
90 gammaUpperReg5by11Table.append(Tuple2<scalar, scalar>(0.0, 1.0));
91 gammaUpperReg8by11Table.append(Tuple2<scalar, scalar>(0.0, 1.0));
93 for (scalar z = 1
e-2; z <= 10.0; z = z + 1
e-2)
95 Tuple2<scalar, scalar> gamma2by11
101 Tuple2<scalar, scalar> gamma5by11
107 Tuple2<scalar, scalar> gamma8by11
113 gammaUpperReg2by11Table.append(gamma2by11);
114 gammaUpperReg5by11Table.append(gamma5by11);
115 gammaUpperReg8by11Table.append(gamma8by11);
118 gammaUpperReg2by11_ =
119 new interpolationTable<scalar>
121 gammaUpperReg2by11Table,
126 gammaUpperReg5by11_ =
127 new interpolationTable<scalar>
129 gammaUpperReg5by11Table,
130 interpolationTable<scalar>::CLAMP,
134 gammaUpperReg8by11_ =
135 new interpolationTable<scalar>
137 gammaUpperReg8by11Table,
138 interpolationTable<scalar>::CLAMP,
148 kolmogorovLengthScale_ =
174 pow(fi.x()/fj.x(), 2.0/3.0) +
pow((1 - fi.x()/fj.x()), 2.0/3.0) - 1
181 beta_*continuousPhase.rho()*
pow(fj.d(), 5.0/3.0)
186 const volScalarField xiMin(minEddyRatio_*kolmogorovLengthScale_/fj.d());
195 2.0*
pow(
b[celli], 3.0/11.0)*tgamma(5.0/11.0)
197 gammaUpperReg5by11_()(
b[celli])
198 - gammaUpperReg5by11_()(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)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
virtual void correct()
Correct diameter independent expressions.
dimensionedScalar pow025(const dimensionedScalar &ds)
Macros for easy insertion into run-time selection tables.
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 dimensionedScalar b
Wien displacement law constant: default SI units: [m K].
dimensionedScalar cbrt(const dimensionedScalar &ds)
Clamp value to the start/end value.
scalar incGammaRatio_Q(const scalar a, const scalar x)
Normalized upper incomplete gamma function.
const dimensionedScalar & nu() const
Return the laminar viscosity.
const populationBalanceModel & popBal_
Reference to the populationBalanceModel.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
const phaseCompressibleTurbulenceModel & continuousTurbulence() const
Return reference to turbulence model of the continuous phase.
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)
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual tmp< volScalarField > epsilon() const =0
Return the turbulence kinetic energy dissipation rate.
const volScalarField & alphas() const
Return total void of phases belonging to this populationBalance.
const tmp< volScalarField > sigmaWithContinuousPhase(const phaseModel &dispersedPhase) const
Return the surface tension coefficient between a given dispersed.
const doubleScalar e
Elementary charge.