28 #include "twoPhaseSystem.H" 39 template<
class BasicTurbulenceModel>
40 LaheyKEpsilon<BasicTurbulenceModel>::LaheyKEpsilon
48 const word& propertiesName,
64 gasTurbulencePtr_(NULL),
106 if (type == typeName)
108 this->printCoeffs(type);
115 template<
class BasicTurbulenceModel>
120 alphaInversion_.readIfPresent(this->coeffDict());
121 Cp_.readIfPresent(this->coeffDict());
122 C3_.readIfPresent(this->coeffDict());
123 Cmub_.readIfPresent(this->coeffDict());
134 template<
class BasicTurbulenceModel>
137 typename BasicTurbulenceModel::transportModel
141 if (!gasTurbulencePtr_)
147 refCast<const twoPhaseSystem>(liquid.fluid());
162 return *gasTurbulencePtr_;
166 template<
class BasicTurbulenceModel>
170 this->gasTurbulence();
173 this->Cmu_*
sqr(this->k_)/this->epsilon_
175 *(
mag(this->U_ - gasTurbulence.
U()));
177 this->nut_.correctBoundaryConditions();
180 BasicTurbulenceModel::correctNut();
184 template<
class BasicTurbulenceModel>
188 this->gasTurbulence();
201 +
pow(fluid.
drag(gas).
CdRe()*liquid.nu()/gas.d(), 4.0/3.0)
212 template<
class BasicTurbulenceModel>
224 max(alphaInversion_ - alpha, scalar(0))
231 template<
class BasicTurbulenceModel>
238 this->gasTurbulence();
240 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
244 + phaseTransferCoeff*gasTurbulence.
k()
245 -
fvm::Sp(phaseTransferCoeff, this->k_);
249 template<
class BasicTurbulenceModel>
256 this->gasTurbulence();
258 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
261 alpha*rho*this->C3_*this->epsilon_*bubbleG()/this->k_
262 + phaseTransferCoeff*gasTurbulence.
epsilon()
263 -
fvm::Sp(phaseTransferCoeff, this->epsilon_);
267 template<
class BasicTurbulenceModel>
const transportModel & transport() const
Access function to incompressible transport model.
BasicTurbulenceModel::rhoField rhoField
const objectRegistry & db() const
Return the local objectRegistry.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
tmp< fvMatrix< Type > > Sp(const DimensionedField< scalar, volMesh > &, const GeometricField< Type, fvPatchField, volMesh > &)
virtual bool read()
Read model coefficients if they have changed.
void correct(GeometricField< Type, fvPatchField, volMesh > &field)
Apply correction to field.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
const volVectorField & U() const
Access function to velocity field.
Generic dimensioned Type class.
tmp< volScalarField > bubbleG() const
Continuous-phase k-epsilon model including bubble-generated turbulence.
Abstract base class for turbulence models (RAS, LES and laminar).
Class which solves the volume fraction equations for two phases.
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
BasicTurbulenceModel::alphaField alphaField
Standard k-epsilon turbulence model for incompressible and compressible flows including rapid distort...
static const word propertiesName
Default name of the turbulence properties dictionary.
A class for handling words, derived from string.
const Time & time() const
Return time.
static word groupName(Name name, const word &group)
virtual tmp< fvScalarMatrix > epsilonSource() const
BasicTurbulenceModel::transportModel transportModel
tmp< volScalarField > phaseTransferCoeff() const
const phaseModel & otherPhase(const phaseModel &phase) const
Constant access the phase not given as an argument.
Templated abstract base class for multiphase compressible turbulence models.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
const alphaField & alpha() const
Access function to phase fraction.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensionedScalar deltaT() const
Return time step.
fileName::Type type(const fileName &)
Return the file type: DIRECTORY or FILE.
virtual tmp< volScalarField > epsilon() const =0
Return the turbulence kinetic energy dissipation rate.
virtual void correctNut()
dimensioned< scalar > mag(const dimensioned< Type > &)
const dragModel & drag(const phaseModel &phase) const
Return the drag model for the given phase.
A class for managing temporary objects.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
virtual tmp< volScalarField > CdRe() const =0
Drag coefficient.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
virtual tmp< fvScalarMatrix > kSource() const