41 template<
class BasicMomentumTransportModel>
64 liquidTurbulencePtr_(
nullptr),
71 this->runTime_.timeName(),
91 this->printCoeffs(type);
98 template<
class BasicMomentumTransportModel>
103 alphaInversion_.readIfPresent(this->coeffDict());
114 template<
class BasicMomentumTransportModel>
131 (rhodv/(18*liquid.rho()*liquid.thermo().nu()))*
sqr(gas.d())
137 exp(
min(thetal/thetag, scalar(50))),
143 nutEff_ = omega*liquidTurbulence.
nut();
148 template<
class BasicMomentumTransportModel>
152 if (!liquidTurbulencePtr_)
160 liquidTurbulencePtr_ =
165 momentumTransportModel::typeName,
171 return *liquidTurbulencePtr_;
175 template<
class BasicMomentumTransportModel>
185 (this->alpha_ - scalar(0.5))/(alphaInversion_ - 0.5),
196 + (1.0 - blend)*rhoEff()*nutEff_/this->transport().rho()
202 template<
class BasicMomentumTransportModel>
216 gas.rho() + (virtualMass.
Cvm() + 3.0/20.0)*liquid.rho()
221 template<
class BasicMomentumTransportModel>
233 max(alphaInversion_ - alpha, scalar(0))
237 liquidTurbulence.
epsilon()/liquidTurbulence.
k(),
244 template<
class BasicMomentumTransportModel>
249 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
252 phaseTransferCoeff*liquidTurbulence.
k()
253 -
fvm::Sp(phaseTransferCoeff, this->k_);
257 template<
class BasicMomentumTransportModel>
262 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
265 phaseTransferCoeff*liquidTurbulence.
epsilon()
266 -
fvm::Sp(phaseTransferCoeff, this->epsilon_);
270 template<
class BasicMomentumTransportModel>
280 tk().boundaryField().types()
Generic thermophysical properties class for a liquid in which the functions and coefficients for each...
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
virtual tmp< fvScalarMatrix > epsilonSource() const
void correct(GeometricField< Type, fvPatchField, volMesh > &field)
Apply correction to field.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const modelType & lookupSubModel(const phasePair &key) const
Return a sub model between a phase pair.
static tmp< GeometricField< scalar, fvPatchField, volMesh > > New(const word &name, const Internal &, const PtrList< fvPatchField< scalar >> &)
Return a temporary field constructed from name,.
virtual bool read()
Re-read model coefficients if they have changed.
label k
Boltzmann constant.
Generic dimensioned Type class.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
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...
Class to represent a system of phases and model interfacial transfers between them.
A class for handling words, derived from string.
static word groupName(Name name, const word &group)
const phaseModel & otherPhase(const phaseModel &phase) const
Return the phase not given as an argument in a two-phase system.
BasicMomentumTransportModel::alphaField alphaField
virtual tmp< volScalarField > nuEff() const
Return the effective viscosity.
virtual tmp< fvScalarMatrix > kSource() const
virtual tmp< volScalarField > nut() const =0
Return the turbulence viscosity.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
virtual tmp< volScalarField > Cvm() const =0
Return the virtual mass coefficient.
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
Abstract base class for turbulence models (RAS, LES and laminar).
continuousGasKEpsilon(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &type=typeName)
Construct from components.
const momentumTransportModel & liquidTurbulence() const
Return the turbulence model for the liquid phase.
tmp< volScalarField > phaseTransferCoeff() const
const Time & time() const
Return time.
BasicMomentumTransportModel::rhoField rhoField
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
virtual tmp< volSymmTensorField > sigma() const
Return the stress tensor [m^2/s^2].
BasicMomentumTransportModel::transportModel transportModel
virtual void correctNut()
virtual tmp< volScalarField > rhoEff() const
Return the effective density for the stress.
virtual void correctNut()
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
A class for managing temporary objects.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
const objectRegistry & db() const
Return the local objectRegistry.
virtual tmp< volScalarField > epsilon() const =0
Return the turbulence kinetic energy dissipation rate.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
dimensionedScalar deltaT() const
Return time step.