37 populationBalance_(popBal),
38 kolmogorovLengthScale_
42 "kolmogorovLengthScale",
43 populationBalance_.time().
name(),
44 populationBalance_.mesh()
46 populationBalance_.mesh(),
49 "kolmogorovLengthScale",
59 populationBalance_.time().
name(),
60 populationBalance_.mesh()
62 populationBalance_.mesh(),
75 populationBalance_.time().
name(),
76 populationBalance_.mesh()
78 populationBalance_.mesh(),
93 kolmogorovLengthScale_ =
96 pow3(populationBalance_.continuousPhase().thermo().nu())
97 /populationBalance_.continuousTurbulence().epsilon()
107 populationBalance_.continuousPhase().rho()
108 *populationBalance_.continuousTurbulence().epsilon()
109 /populationBalance_.continuousPhase().thermo().mu()
112 if (uTerminal_.empty())
114 const fvMesh& mesh = populationBalance_.mesh();
122 gAverage(populationBalance_.continuousPhase().thermo().nu()())
129 gAverage(populationBalance_.continuousPhase().rho())
136 gAverage(populationBalance_.sizeGroups()[1].phase().rho())
145 populationBalance_.sigmaWithContinuousPhase
147 populationBalance_.sizeGroups()[1].phase()
152 for(
int m = 0; m < populationBalance_.sizeGroups().size(); m++)
154 const sizeGroup&
f = populationBalance_.sizeGroups()[m];
177 Re = uTerminal*
f.dSph()/nuc;
181 +
neg(1000 -
Re)*0.44;
183 CdEllipse = 0.6666*
sqrt(Eo);
188 +
neg(CdEllipse - Cd)*Cd;
191 4.0/3.0*(rhoc - rhod)*
mag(g)*
f.dSph()
192 - rhoc*Cd*
sqr(uTerminal);
194 ReX = (uTerminal + uTerminalX)*
f.dSph()/nuc;
198 *24/ReX*(1 + 0.1*
pow(ReX, 0.75))
199 +
neg(1000 - ReX)*0.44;
202 pos0(CdEllipse - CdX)
204 +
neg(CdEllipse - CdX)*CdX;
206 dCd = (CdX - Cd)/uTerminalX;
208 dF = -(2*rhoc*uTerminal*Cd + rhoc*
sqr(uTerminal)*dCd);
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
virtual void precompute()
Precompute diameter independent expressions.
LiaoBase(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 ...
A list of keyword definitions, which are a keyword followed by any number of values (e....
const dimensionSet & dimensions() const
Return const reference to dimensions.
const Type & value() const
Return const reference to value.
Mesh data needed to do the Finite Volume discretisation.
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type and name.
Calculate the gradient of the given field.
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
const dimensionedScalar F
Faraday constant: default SI units: [C/mol].
tmp< VolField< typename outerProduct< vector, Type >::type > > grad(const SurfaceField< Type > &ssf)
const dimensionSet dimViscosity
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
dimensionedScalar pos0(const dimensionedScalar &ds)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimless
const dimensionSet dimLength
const dimensionSet dimForce
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
const dimensionSet dimDensity
dimensioned< scalar > mag(const dimensioned< Type > &)
dimensionedScalar neg(const dimensionedScalar &ds)
Type gAverage(const FieldField< Field, Type > &f)
const dimensionSet dimVelocity
const dimensionSet dimArea
word name(const complex &)
Return a string representation of a complex.
dimensionedScalar pow025(const dimensionedScalar &ds)
scalarField Re(const UList< complex > &cf)