31 #include "twoPhaseSystem.H" 32 #include "dragModel.H" 43 template<
class BasicTurbulenceModel>
52 const word& propertiesName,
68 gasTurbulencePtr_(
nullptr),
110 if (
type == typeName)
112 this->printCoeffs(
type);
119 template<
class BasicTurbulenceModel>
124 alphaInversion_.readIfPresent(this->coeffDict());
125 Cp_.readIfPresent(this->coeffDict());
126 C3_.readIfPresent(this->coeffDict());
127 Cmub_.readIfPresent(this->coeffDict());
136 template<
class BasicTurbulenceModel>
137 const PhaseCompressibleTurbulenceModel
139 typename BasicTurbulenceModel::transportModel
141 LaheyKEpsilon<BasicTurbulenceModel>::gasTurbulence()
const 143 if (!gasTurbulencePtr_)
147 const transportModel& liquid = this->transport();
148 const twoPhaseSystem&
fluid =
149 refCast<const twoPhaseSystem>(liquid.fluid());
150 const transportModel& gas =
fluid.otherPhase(liquid);
154 .lookupObject<PhaseCompressibleTurbulenceModel<transportModel>>
164 return *gasTurbulencePtr_;
168 template<
class BasicTurbulenceModel>
172 this->gasTurbulence();
175 this->Cmu_*
sqr(this->k_)/this->epsilon_
177 *(
mag(this->U_ - gasTurbulence.
U()));
179 this->nut_.correctBoundaryConditions();
182 BasicTurbulenceModel::correctNut();
186 template<
class BasicTurbulenceModel>
190 this->gasTurbulence();
216 template<
class BasicTurbulenceModel>
221 const alphaField&
alpha = this->alpha_;
222 const rhoField&
rho = this->rho_;
228 max(alphaInversion_ -
alpha, scalar(0))
230 *
min(gasTurbulence.
epsilon()/gasTurbulence.
k(), 1.0/
U.time().deltaT())
235 template<
class BasicTurbulenceModel>
238 const alphaField&
alpha = this->alpha_;
239 const rhoField&
rho = this->rho_;
242 this->gasTurbulence();
244 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
248 + phaseTransferCoeff*gasTurbulence.
k()
249 -
fvm::Sp(phaseTransferCoeff, this->k_);
253 template<
class BasicTurbulenceModel>
256 const alphaField&
alpha = this->alpha_;
257 const rhoField&
rho = this->rho_;
260 this->gasTurbulence();
262 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
265 alpha*
rho*this->C3_*this->epsilon_*bubbleG()/this->k_
266 + phaseTransferCoeff*gasTurbulence.
epsilon()
267 -
fvm::Sp(phaseTransferCoeff, this->epsilon_);
271 template<
class BasicTurbulenceModel>
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.
BasicTurbulenceModel::rhoField rhoField
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const volVectorField & U() const
Access function to velocity field.
virtual bool read()
Read model coefficients if they have changed.
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 void correct()
Solve the turbulence equations and correct the turbulence viscosity.
Generic dimensioned Type class.
GeometricField< vector, fvPatchField, volMesh > volVectorField
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.
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.
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 Foam::string.
const alphaField & alpha() const
Access function to phase fraction.
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.
BasicTurbulenceModel::transportModel transportModel
Templated abstract base class for multiphase compressible turbulence models.
tmp< volScalarField > phaseTransferCoeff() const
virtual bool read()
Re-read model coefficients if they have changed.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
virtual tmp< fvScalarMatrix > kSource() const
dimensionedScalar pow3(const dimensionedScalar &ds)
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.
tmp< volScalarField > bubbleG() const
virtual void correctNut()
virtual tmp< fvScalarMatrix > epsilonSource() const
A class for managing temporary objects.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to database if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].