31#include "phaseCompressibleTurbulenceModel.H"
37namespace diameterModels
39namespace breakupModels
55Foam::diameterModels::breakupModels::LaakkonenAlopaeusAittamaa::
56LaakkonenAlopaeusAittamaa
87 const phaseModel& continuousPhase = popBal_.continuousPhase();
88 const sizeGroup& fi = popBal_.sizeGroups()[i];
91 C1_*
cbrt(popBal_.continuousTurbulence().epsilon())
96 C2_*popBal_.sigmaWithContinuousPhase(fi.
phase())
98 continuousPhase.
rho()*
pow(fi.
d(), 5.0/3.0)
99 *
pow(popBal_.continuousTurbulence().epsilon(), 2.0/3.0)
101 + C3_*continuousPhase.
mu()
104 *
cbrt(popBal_.continuousTurbulence().epsilon())
105 *
pow(fi.
d(), 4.0/3.0)
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
Base class for breakup models which give a total breakup rate and a separate daughter size distributi...
Model of Laakkonen et al. (2006). The total breakup rate is calculated by.
virtual void setBreakupRate(volScalarField &breakupRate, const label i)
Set total breakupRate.
Class that solves the univariate population balance equation by means of a class method (also called ...
This class represents a single sizeGroup belonging to a velocityGroup. The main property of a sizeGro...
const dimensionedScalar & d() const
Return representative diameter of the sizeGroup.
const phaseModel & phase() const
Return const-reference to the phase.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Dimension set for the base types, which can be used to implement rigorous dimension checking for alge...
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
const dimensionedScalar & rho() const
virtual tmp< volScalarField > mu() const =0
Return the laminar dynamic viscosity.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
dimensionedScalar erfc(const dimensionedScalar &ds)
const dimensionSet dimless
Dimensionless.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensionedScalar cbrt(const dimensionedScalar &ds)