28 #include "twoPhaseSystem.H" 39 template<
class BasicTurbulenceModel>
40 NicenoKEqn<BasicTurbulenceModel>::NicenoKEqn
48 const word& propertiesName,
64 gasTurbulencePtr_(NULL),
98 this->printCoeffs(type);
105 template<
class BasicTurbulenceModel>
110 alphaInversion_.readIfPresent(this->coeffDict());
111 Cp_.readIfPresent(this->coeffDict());
112 Cmub_.readIfPresent(this->coeffDict());
123 template<
class BasicTurbulenceModel>
126 typename BasicTurbulenceModel::transportModel
130 if (!gasTurbulencePtr_)
136 refCast<const twoPhaseSystem>(liquid.fluid());
151 return *gasTurbulencePtr_;
155 template<
class BasicTurbulenceModel>
159 this->gasTurbulence();
164 *(
mag(this->U_ - gasTurbulence.
U()));
166 this->nut_.correctBoundaryConditions();
169 BasicTurbulenceModel::correctNut();
173 template<
class BasicTurbulenceModel>
177 this->gasTurbulence();
181 refCast<const twoPhaseSystem>(liquid.fluid());
188 Cp_*
sqr(magUr)*fluid.
drag(gas).
K()/liquid.rho()
195 template<
class BasicTurbulenceModel>
207 max(alphaInversion_ - alpha, scalar(0))
211 this->Ce_*
sqrt(gasTurbulence.
k())/this->
delta(),
218 template<
class BasicTurbulenceModel>
225 this->gasTurbulence();
227 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
231 + phaseTransferCoeff*gasTurbulence.
k()
232 -
fvm::Sp(phaseTransferCoeff, this->k_);
tmp< volScalarField > bubbleG() const
const transportModel & transport() const
Access function to incompressible transport model.
BasicTurbulenceModel::alphaField alphaField
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)
const volVectorField & U() const
Access function to velocity field.
dimensionedScalar sqrt(const dimensionedScalar &ds)
Generic dimensioned Type class.
BasicTurbulenceModel::rhoField rhoField
Abstract base class for turbulence models (RAS, LES and laminar).
One equation eddy-viscosity model.
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.
static const word propertiesName
Default name of the turbulence properties dictionary.
virtual bool read()
Read model coefficients if they have changed.
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 > kSource() 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 deltaT() const
Return time step.
One-equation SGS model for the continuous phase in a two-phase system including bubble-generated turb...
fileName::Type type(const fileName &)
Return the file type: DIRECTORY or FILE.
tmp< volScalarField > phaseTransferCoeff() const
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.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
BasicTurbulenceModel::transportModel transportModel
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
virtual tmp< volScalarField > K(const volScalarField &Ur) const =0
The dragfunction K used in the momentum eq.
virtual void correctNut()