28 #include "twoPhaseSystem.H" 29 #include "virtualMassModel.H" 40 template<
class BasicTurbulenceModel>
41 continuousGasKEpsilon<BasicTurbulenceModel>::continuousGasKEpsilon
49 const word& propertiesName,
65 liquidTurbulencePtr_(NULL),
72 this->runTime_.timeName(),
92 this->printCoeffs(type);
99 template<
class BasicTurbulenceModel>
104 alphaInversion_.readIfPresent(this->coeffDict());
115 template<
class BasicTurbulenceModel>
132 exp(
min(thetal/thetag, scalar(50))),
138 nutEff_ = omega*liquidTurbulence.
nut();
143 template<
class BasicTurbulenceModel>
147 if (!liquidTurbulencePtr_)
153 refCast<const twoPhaseSystem>(gas.fluid());
156 liquidTurbulencePtr_ =
167 return *liquidTurbulencePtr_;
171 template<
class BasicTurbulenceModel>
181 (this->alpha_ - scalar(0.5))/(alphaInversion_ - 0.5),
194 + (1.0 - blend)*rhoEff()*nutEff_/this->transport().rho()
201 template<
class BasicTurbulenceModel>
214 gas.rho() + (fluid.
virtualMass(gas).
Cvm() + 3.0/20.0)*liquid.rho()
220 template<
class BasicTurbulenceModel>
232 max(alphaInversion_ - alpha, scalar(0))
236 liquidTurbulence.
epsilon()/liquidTurbulence.
k(),
243 template<
class BasicTurbulenceModel>
248 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
251 phaseTransferCoeff*liquidTurbulence.
k()
252 -
fvm::Sp(phaseTransferCoeff, this->k_);
256 template<
class BasicTurbulenceModel>
261 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
264 phaseTransferCoeff*liquidTurbulence.
epsilon()
265 -
fvm::Sp(phaseTransferCoeff, this->epsilon_);
269 template<
class BasicTurbulenceModel>
282 this->runTime_.timeName(),
288 tk().boundaryField().types()
BasicTurbulenceModel::transportModel transportModel
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
const virtualMassModel & virtualMass(const phaseModel &phase) const
Return the virtual mass model for the given phase.
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 > &)
void correct(GeometricField< Type, fvPatchField, volMesh > &field)
Apply correction to field.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
virtual bool read()
Re-read model coefficients if they have changed.
label k
Boltzmann constant.
Generic dimensioned Type class.
tmp< volScalarField > phaseTransferCoeff() const
virtual tmp< fvScalarMatrix > epsilonSource() const
Abstract base class for turbulence models (RAS, LES and laminar).
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
Class which solves the volume fraction equations for two phases.
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
dimensionedScalar exp(const dimensionedScalar &ds)
static const Identity< scalar > I
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< volScalarField > rhoEff() const
Return the effective density for the stress.
virtual tmp< volScalarField > nuEff() const
Return the effective viscosity.
const phaseModel & otherPhase(const phaseModel &phase) const
Constant access the phase not given as an argument.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
virtual tmp< volScalarField > Cvm() const =0
Return the virtual mass coefficient.
BasicTurbulenceModel::alphaField alphaField
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.
const turbulenceModel & liquidTurbulence() const
Return the turbulence model for the liquid phase.
virtual tmp< fvScalarMatrix > kSource() const
virtual void correctNut()
BasicTurbulenceModel::rhoField rhoField
virtual void correctNut()
virtual tmp< volScalarField > nut() const =0
Return the turbulence viscosity.
A class for managing temporary objects.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
virtual tmp< volSymmTensorField > R() const
Return the Reynolds stress tensor.