31 #include "twoPhaseSystem.H" 32 #include "virtualMassModel.H" 33 #include "dragModel.H" 44 template<
class BasicTurbulenceModel>
53 const word& propertiesName,
69 liquidTurbulencePtr_(
nullptr),
76 this->runTime_.timeName(),
96 this->printCoeffs(
type);
103 template<
class BasicTurbulenceModel>
108 alphaInversion_.readIfPresent(this->coeffDict());
117 template<
class BasicTurbulenceModel>
137 exp(
min(thetal/thetag, scalar(50))),
143 nutEff_ = omega*liquidTurbulence.
nut();
148 template<
class BasicTurbulenceModel>
152 if (!liquidTurbulencePtr_)
158 refCast<const twoPhaseSystem>(gas.fluid());
161 liquidTurbulencePtr_ =
172 return *liquidTurbulencePtr_;
176 template<
class BasicTurbulenceModel>
186 (this->alpha_ - scalar(0.5))/(alphaInversion_ - 0.5),
199 + (1.0 - blend)*rhoEff()*nutEff_/this->transport().rho()
206 template<
class BasicTurbulenceModel>
228 template<
class BasicTurbulenceModel>
233 const alphaField&
alpha = this->alpha_;
234 const rhoField&
rho = this->rho_;
240 max(alphaInversion_ -
alpha, scalar(0))
244 liquidTurbulence.
epsilon()/liquidTurbulence.
k(),
245 1.0/
U.time().deltaT()
251 template<
class BasicTurbulenceModel>
256 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
259 phaseTransferCoeff*liquidTurbulence.
k()
260 -
fvm::Sp(phaseTransferCoeff, this->k_);
264 template<
class BasicTurbulenceModel>
269 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
273 -
fvm::Sp(phaseTransferCoeff, this->epsilon_);
277 template<
class BasicTurbulenceModel>
290 this->runTime_.timeName(),
296 tk().boundaryField().types()
BasicTurbulenceModel::transportModel transportModel
Generic thermophysical properties class for a liquid in which the functions and coefficients for each...
void correct(GeometricField< Type, PatchField, GeoMesh > &field)
Apply correction to field.
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
virtual tmp< fvScalarMatrix > epsilonSource() const
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
virtual bool read()
Re-read model coefficients if they have changed.
label k
Boltzmann constant.
Generic dimensioned Type class.
Ignore writing from objectRegistry::writeObject()
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.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
scalar rho(scalar p, scalar T) const
Liquid density [kg/m^3].
dimensionedScalar exp(const dimensionedScalar &ds)
k-epsilon model for the gas-phase in a two-phase system supporting phase-inversion.
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 Foam::string.
Reading is optional [identical to LAZY_READ].
virtual tmp< volScalarField > nuEff() const
Return the effective viscosity.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
virtual tmp< fvScalarMatrix > kSource() const
virtual tmp< volScalarField > Cvm() const =0
Return the virtual mass coefficient.
virtual bool read()
Re-read model coefficients if they have changed.
virtual tmp< volSymmTensorField > R() const
Return the Reynolds stress tensor.
BasicTurbulenceModel::alphaField alphaField
tmp< volScalarField > phaseTransferCoeff() const
Base-class for all transport models used by the incompressible turbulence models. ...
virtual tmp< volScalarField > epsilon() const =0
Return the turbulence kinetic energy dissipation rate.
Automatically write from objectRegistry::writeObject()
virtual void correctNut()
BasicTurbulenceModel::rhoField rhoField
virtual tmp< volScalarField > rhoEff() const
Return the effective density for the stress.
const turbulenceModel & liquidTurbulence() const
Return the turbulence model for the liquid phase.
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 database if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
Defines the attributes of an object for which implicit objectRegistry management is supported...
SymmTensor< Cmpt > devTwoSymm(const SymmTensor< Cmpt > &st)
Return the deviatoric part of twice the symmetric part of a SymmTensor.