29 #include "phaseCompressibleTurbulenceModel.H" 35 namespace diameterModels
37 namespace coalescenceModels
57 const populationBalanceModel& popBal,
58 const dictionary&
dict 61 coalescenceModel(popBal, dict),
62 C1_(
"C1",
dimless, dict.lookupOrDefault<scalar>(
"C1", 0.356)),
63 h0_(
"h0",
dimLength, dict.lookupOrDefault<scalar>(
"h0", 1
e-4)),
64 hf_(
"hf",
dimLength, dict.lookupOrDefault<scalar>(
"h0", 1
e-8)),
65 turbulentCollisions_(dict.
lookup(
"turbulentCollisions")),
66 buoyantCollisions_(dict.
lookup(
"buoyantCollisions")),
67 laminarShearCollisions_(dict.
lookup(
"laminarShearCollisions"))
95 pow3(rij)*continuousPhase.rho()
103 if (turbulentCollisions_)
107 C1_*
pi*
sqr(fi.d() + fj.d())
109 *
sqrt(
pow(fi.d(), 2.0/3.0) +
pow(fj.d(), 2.0/3.0))
111 *collisionEfficiency;
114 if (buoyantCollisions_)
126 /(continuousPhase.rho()*fi.d()) + 0.505*
mag(g)*fi.d()
131 /(continuousPhase.rho()*fj.d()) + 0.505*
mag(g)*fj.d()
135 *collisionEfficiency;
138 if (laminarShearCollisions_)
141 <<
"Laminar shear collision contribution not implemented for " 142 << this->
type() <<
" coalescence model."
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
dimensionedScalar log(const dimensionedScalar &ds)
errorManipArg< error, int > exit(error &err, const int errNo=1)
UniformDimensionedField< vector > uniformDimensionedVectorField
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
stressControl lookup("compactNormalStress") >> compactNormalStress
dimensionedScalar exp(const dimensionedScalar &ds)
const phaseModel & continuousPhase() const
Return continuous phase.
PrinceBlanch(const populationBalanceModel &popBal, const dictionary &dict)
dimensionedScalar cbrt(const dimensionedScalar &ds)
const List< sizeGroup * > & sizeGroups() const
Return the sizeGroups belonging to this populationBalance.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
const phaseCompressibleTurbulenceModel & continuousTurbulence() const
Return reference to turbulence model of the continuous phase.
const populationBalanceModel & popBal_
Reference to the populationBalanceModel.
defineTypeNameAndDebug(combustionModel, 0)
fileName::Type type(const fileName &, const bool followLink=true)
Return the file type: DIRECTORY or FILE.
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.
dimensioned< scalar > mag(const dimensioned< Type > &)
const fvMesh & mesh() const
Return reference to the mesh.
const tmp< volScalarField > sigmaWithContinuousPhase(const phaseModel &dispersedPhase) const
Return the surface tension coefficient between a given dispersed.
const doubleScalar e
Elementary charge.
virtual void addToCoalescenceRate(volScalarField &coalescenceRate, const label i, const label j)
Add to coalescenceRate.