31 #include "twoPhaseSystem.H" 42 template<
class BasicTurbulenceModel>
51 const word& propertiesName,
67 gasTurbulencePtr_(
nullptr),
81 this->printCoeffs(
type);
88 template<
class BasicTurbulenceModel>
93 Cmub_.readIfPresent(this->coeffDict());
102 template<
class BasicTurbulenceModel>
103 const PhaseCompressibleTurbulenceModel
105 typename BasicTurbulenceModel::transportModel
107 kOmegaSSTSato<BasicTurbulenceModel>::gasTurbulence()
const 109 if (!gasTurbulencePtr_)
113 const transportModel& liquid = this->transport();
114 const twoPhaseSystem&
fluid =
115 refCast<const twoPhaseSystem>(liquid.fluid());
116 const transportModel& gas =
fluid.otherPhase(liquid);
120 .lookupObject<PhaseCompressibleTurbulenceModel<transportModel>>
130 return *gasTurbulencePtr_;
134 template<
class BasicTurbulenceModel>
141 this->gasTurbulence();
145 pow(this->betaStar_, 0.25)*this->y_*
sqrt(this->k_)/this->
nu()
152 this->a1_*this->omega_,
153 this->b1_*this->F23()*
sqrt(S2)
157 *(
mag(this->U_ - gasTurbulence.
U()));
159 this->nut_.correctBoundaryConditions();
162 BasicTurbulenceModel::correctNut();
166 template<
class BasicTurbulenceModel>
void correct(GeometricField< Type, PatchField, GeoMesh > &field)
Apply correction to field.
Implementation of the k-omega-SST turbulence model for dispersed bubbly flows with Sato (1981) bubble...
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const volVectorField & U() const
Access function to velocity field.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
BasicTurbulenceModel::transportModel transportModel
BasicTurbulenceModel::alphaField alphaField
Generic dimensioned Type class.
GeometricField< vector, fvPatchField, volMesh > volVectorField
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.
dimensionedScalar exp(const dimensionedScalar &ds)
static const word propertiesName
Default name of the turbulence properties dictionary.
A class for handling words, derived from Foam::string.
BasicTurbulenceModel::rhoField rhoField
const alphaField & alpha() const
Access function to phase fraction.
Implementation of the k-omega-SST turbulence model for incompressible and compressible flows...
Templated abstract base class for multiphase compressible turbulence models.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
virtual bool read()
Read model coefficients if they have changed.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
virtual bool read()
Re-read model coefficients if they have changed.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to database if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
virtual void correctNut()