28 #include "twoPhaseSystem.H" 29 #include "dragModel.H" 40 template<
class BasicTurbulenceModel>
49 const word& propertiesName,
65 gasTurbulencePtr_(
nullptr),
107 if (type == typeName)
109 this->printCoeffs(type);
116 template<
class BasicTurbulenceModel>
121 alphaInversion_.readIfPresent(this->coeffDict());
122 Cp_.readIfPresent(this->coeffDict());
123 C3_.readIfPresent(this->coeffDict());
124 Cmub_.readIfPresent(this->coeffDict());
135 template<
class BasicTurbulenceModel>
138 typename BasicTurbulenceModel::transportModel
142 if (!gasTurbulencePtr_)
148 refCast<const twoPhaseSystem>(liquid.fluid());
163 return *gasTurbulencePtr_;
167 template<
class BasicTurbulenceModel>
171 this->gasTurbulence();
174 this->Cmu_*
sqr(this->k_)/this->epsilon_
176 *(
mag(this->U_ - gasTurbulence.
U()));
178 this->nut_.correctBoundaryConditions();
181 BasicTurbulenceModel::correctNut();
185 template<
class BasicTurbulenceModel>
189 this->gasTurbulence();
204 +
pow(drag.
CdRe()*liquid.nu()/gas.d(), 4.0/3.0)
215 template<
class BasicTurbulenceModel>
227 max(alphaInversion_ - alpha, scalar(0))
234 template<
class BasicTurbulenceModel>
241 this->gasTurbulence();
243 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
247 + phaseTransferCoeff*gasTurbulence.
k()
248 -
fvm::Sp(phaseTransferCoeff, this->k_);
252 template<
class BasicTurbulenceModel>
259 this->gasTurbulence();
261 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
264 alpha*rho*this->C3_*this->epsilon_*bubbleG()/this->k_
265 + phaseTransferCoeff*gasTurbulence.
epsilon()
266 -
fvm::Sp(phaseTransferCoeff, this->epsilon_);
270 template<
class BasicTurbulenceModel>
const modelType & lookupSubModel(const phasePair &key) const
Access a sub model between a phase pair.
BasicTurbulenceModel::rhoField rhoField
const volVectorField & U() const
Access function to velocity field.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
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.
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.
Continuous-phase k-epsilon model including bubble-generated turbulence.
Abstract base class for turbulence models (RAS, LES and laminar).
Info<< "Reading strained laminar flame speed field Su\"<< endl;volScalarField Su(IOobject("Su", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field betav\"<< endl;volScalarField betav(IOobject("betav", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Lobs\"<< endl;volScalarField Lobs(IOobject("Lobs", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field CT\"<< endl;volSymmTensorField CT(IOobject("CT", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Nv\"<< endl;volScalarField Nv(IOobject("Nv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field nsv\"<< endl;volSymmTensorField nsv(IOobject("nsv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);IOdictionary PDRProperties(IOobject("PDRProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE));autoPtr< PDRDragModel > drag
Class which solves the volume fraction equations for two phases.
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
const transportModel & transport() const
Access function to incompressible transport model.
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.
static word groupName(Name name, const word &group)
const alphaField & alpha() const
Access function to phase fraction.
BasicTurbulenceModel::transportModel transportModel
Templated abstract base class for multiphase compressible turbulence models.
const phaseModel & otherPhase(const phaseModel &phase) const
Return the phase not given as an argument.
tmp< volScalarField > phaseTransferCoeff() const
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
virtual tmp< fvScalarMatrix > kSource() const
dimensionedScalar pow3(const dimensionedScalar &ds)
const Time & time() const
Return time.
virtual tmp< volScalarField > epsilon() const =0
Return the turbulence kinetic energy dissipation rate.
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
tmp< volScalarField > bubbleG() const
virtual void correctNut()
LaheyKEpsilon(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName=turbulenceModel::propertiesName, const word &type=typeName)
Construct from components.
dimensioned< scalar > mag(const dimensioned< Type > &)
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
virtual tmp< fvScalarMatrix > epsilonSource() const
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 objectRegistry & db() const
Return the local objectRegistry.
dimensionedScalar deltaT() const
Return time step.