Mixture k-epsilon turbulence model for two-phase gas-liquid systems. More...
Public Types | |
typedef BasicTurbulenceModel::alphaField | alphaField |
typedef BasicTurbulenceModel::rhoField | rhoField |
typedef BasicTurbulenceModel::transportModel | transportModel |
Public Types inherited from eddyViscosity< RASModel< BasicTurbulenceModel > > | |
typedef RASModel< BasicTurbulenceModel > ::alphaField | alphaField |
typedef RASModel< BasicTurbulenceModel > ::rhoField | rhoField |
typedef RASModel< BasicTurbulenceModel > ::transportModel | transportModel |
Public Types inherited from linearViscousStress< RASModel< BasicTurbulenceModel > > | |
typedef RASModel< BasicTurbulenceModel > ::alphaField | alphaField |
typedef RASModel< BasicTurbulenceModel > ::rhoField | rhoField |
typedef RASModel< BasicTurbulenceModel > ::transportModel | transportModel |
Public Types inherited from RASModel< BasicTurbulenceModel > | |
typedef BasicTurbulenceModel::alphaField | alphaField |
typedef BasicTurbulenceModel::rhoField | rhoField |
typedef BasicTurbulenceModel::transportModel | transportModel |
Public Member Functions | |
TypeName ("mixtureKEpsilon") | |
Runtime type information. More... | |
mixtureKEpsilon (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. More... | |
virtual | ~mixtureKEpsilon ()=default |
Destructor. More... | |
virtual bool | read () |
Re-read model coefficients if they have changed. More... | |
virtual tmp< volScalarField > | k () const |
Return the turbulence kinetic energy. More... | |
virtual tmp< volScalarField > | epsilon () const |
Return the turbulence kinetic energy dissipation rate. More... | |
virtual void | correct () |
Solve the turbulence equations and correct the turbulence viscosity. More... | |
Public Member Functions inherited from eddyViscosity< RASModel< BasicTurbulenceModel > > | |
eddyViscosity (const word &modelName, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName) | |
Construct from components. More... | |
virtual | ~eddyViscosity ()=default |
Destructor. More... | |
virtual tmp< volScalarField > | nut () const |
Return the turbulence viscosity. More... | |
virtual tmp< scalarField > | nut (const label patchi) const |
Return the turbulence viscosity on patch. More... | |
virtual tmp< volSymmTensorField > | R () const |
Return the Reynolds stress tensor. More... | |
virtual void | validate () |
Validate the turbulence fields after construction. More... | |
Public Member Functions inherited from linearViscousStress< RASModel< BasicTurbulenceModel > > | |
linearViscousStress (const word &modelName, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName) | |
Construct from components. More... | |
virtual | ~linearViscousStress ()=default |
Destructor. More... | |
virtual tmp< volSymmTensorField > | devRhoReff () const |
Return the effective stress tensor. More... | |
virtual tmp< volSymmTensorField > | devRhoReff (const volVectorField &U) const |
Return the effective stress tensor based on a given velocity field. More... | |
virtual tmp< fvVectorMatrix > | divDevRhoReff (volVectorField &U) const |
Return the source term for the momentum equation. More... | |
virtual tmp< fvVectorMatrix > | divDevRhoReff (const volScalarField &rho, volVectorField &U) const |
Return the source term for the momentum equation. More... | |
Public Member Functions inherited from RASModel< BasicTurbulenceModel > | |
TypeName ("RAS") | |
Runtime type information. More... | |
declareRunTimeSelectionTable (autoPtr, RASModel, dictionary,(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName),(alpha, rho, U, alphaRhoPhi, phi, transport, propertiesName)) | |
RASModel (const word &type, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName) | |
Construct from components. More... | |
virtual | ~RASModel ()=default |
Destructor. More... | |
const dimensionedScalar & | kMin () const noexcept |
Return the lower allowable limit for k (default: SMALL) More... | |
const dimensionedScalar & | epsilonMin () const noexcept |
Return the lower allowable limit for epsilon (default: SMALL) More... | |
const dimensionedScalar & | omegaMin () const noexcept |
Return the lower allowable limit for omega (default: SMALL) More... | |
dimensionedScalar & | kMin () noexcept |
Allow kMin to be changed. More... | |
dimensionedScalar & | epsilonMin () noexcept |
Allow epsilonMin to be changed. More... | |
dimensionedScalar & | omegaMin () noexcept |
Allow omegaMin to be changed. More... | |
virtual const dictionary & | coeffDict () const |
Const access to the coefficients dictionary. More... | |
virtual tmp< volScalarField > | nuEff () const |
Return the effective viscosity. More... | |
virtual tmp< scalarField > | nuEff (const label patchi) const |
Return the effective viscosity on patch. More... | |
virtual tmp< volScalarField > | omega () const |
Return the specific dissipation rate. More... | |
Public Member Functions inherited from RASModelBase | |
ClassName ("RASModelBase") | |
Runtime type information. More... | |
RASModelBase () noexcept=default | |
Constructor. More... | |
virtual | ~RASModelBase ()=default |
Destructor. More... | |
Protected Member Functions | |
wordList | epsilonBoundaryTypes (const volScalarField &epsilon) const |
void | correctInletOutlet (volScalarField &vsf, const volScalarField &refVsf) const |
void | initMixtureFields () |
virtual void | correctNut () |
tmp< volScalarField > | Ct2 () const |
tmp< volScalarField > | rholEff () const |
tmp< volScalarField > | rhogEff () const |
tmp< volScalarField > | rhom () const |
tmp< volScalarField > | mix (const volScalarField &fc, const volScalarField &fd) const |
tmp< volScalarField > | mixU (const volScalarField &fc, const volScalarField &fd) const |
tmp< surfaceScalarField > | mixFlux (const surfaceScalarField &fc, const surfaceScalarField &fd) const |
tmp< volScalarField > | bubbleG () const |
virtual tmp< fvScalarMatrix > | kSource () const |
virtual tmp< fvScalarMatrix > | epsilonSource () const |
tmp< volScalarField > | DkEff (const volScalarField &nutm) const |
Return the effective diffusivity for k. More... | |
tmp< volScalarField > | DepsilonEff (const volScalarField &nutm) const |
Return the effective diffusivity for epsilon. More... | |
Protected Member Functions inherited from RASModel< BasicTurbulenceModel > | |
virtual void | printCoeffs (const word &type) |
Print model coefficients. More... | |
RASModel (const RASModel &)=delete | |
No copy construct. More... | |
void | operator= (const RASModel &)=delete |
No copy assignment. More... | |
Protected Attributes | |
dimensionedScalar | Cmu_ |
dimensionedScalar | C1_ |
dimensionedScalar | C2_ |
dimensionedScalar | C3_ |
dimensionedScalar | Cp_ |
dimensionedScalar | sigmak_ |
dimensionedScalar | sigmaEps_ |
volScalarField | k_ |
volScalarField | epsilon_ |
autoPtr< volScalarField > | Ct2_ |
autoPtr< volScalarField > | rhom_ |
autoPtr< volScalarField > | km_ |
autoPtr< volScalarField > | epsilonm_ |
Protected Attributes inherited from eddyViscosity< RASModel< BasicTurbulenceModel > > | |
volScalarField | nut_ |
Protected Attributes inherited from RASModel< BasicTurbulenceModel > | |
dictionary | RASDict_ |
RAS coefficients dictionary. More... | |
Switch | turbulence_ |
Turbulence on/off flag. More... | |
Switch | printCoeffs_ |
Flag to print the model coeffs at run-time. More... | |
dictionary | coeffDict_ |
Model coefficients dictionary. More... | |
dimensionedScalar | kMin_ |
Lower limit of k. More... | |
dimensionedScalar | epsilonMin_ |
Lower limit of epsilon. More... | |
dimensionedScalar | omegaMin_ |
Lower limit for omega. More... | |
Additional Inherited Members | |
Static Public Member Functions inherited from RASModel< BasicTurbulenceModel > | |
static autoPtr< RASModel > | New (const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName=turbulenceModel::propertiesName) |
Return a reference to the selected RAS model. More... | |
Mixture k-epsilon turbulence model for two-phase gas-liquid systems.
The basic structure of the model is based on:
Behzadi, A., Issa, R. I., & Rusche, H. (2004). Modelling of dispersed bubble and droplet flow at high phase fractions. Chemical Engineering Science, 59(4), 759-770.
but an effective density for the gas is used in the averaging and an alternative model for bubble-generated turbulence from:
Lahey Jr, R. T. (2005). The simulation of multidimensional multiphase flows. Nuclear Engineering and Design, 235(10), 1043-1060.
The default model coefficients are
mixtureKEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; C3 C2; sigmak 1.0; sigmaEps 1.3; }
Definition at line 82 of file mixtureKEpsilon.H.
typedef BasicTurbulenceModel::alphaField alphaField |
Definition at line 210 of file mixtureKEpsilon.H.
typedef BasicTurbulenceModel::rhoField rhoField |
Definition at line 211 of file mixtureKEpsilon.H.
typedef BasicTurbulenceModel::transportModel transportModel |
Definition at line 212 of file mixtureKEpsilon.H.
mixtureKEpsilon | ( | 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.
Definition at line 43 of file mixtureKEpsilon.C.
References Foam::bound(), and Foam::type().
|
virtualdefault |
Destructor.
|
protected |
Definition at line 169 of file mixtureKEpsilon.C.
References epsilon, forAll, and GeometricBoundaryField< Type, PatchField, GeoMesh >::types().
|
protected |
Definition at line 191 of file mixtureKEpsilon.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), and forAll.
|
protected |
Definition at line 218 of file mixtureKEpsilon.C.
References IOobjectOption::AUTO_WRITE, GeometricField< Type, PatchField, GeoMesh >::boundaryField(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon_, mixtureKEpsilon< BasicTurbulenceModel >::k_, IOobjectOption::READ_IF_PRESENT, and GeometricBoundaryField< Type, PatchField, GeoMesh >::types().
|
protectedvirtual |
Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.
Definition at line 329 of file mixtureKEpsilon.C.
References optionList::correct(), options::New(), and Foam::sqr().
Referenced by mixtureKEpsilon< BasicTurbulenceModel >::correct().
|
protected |
Definition at line 371 of file mixtureKEpsilon.C.
References beta(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon_, Foam::exp(), fluid, mixtureKEpsilon< BasicTurbulenceModel >::k_, Foam::mag(), liquid::rho(), Foam::sqr(), Foam::sqrt(), and U.
|
protected |
Definition at line 398 of file mixtureKEpsilon.C.
References fluid.
|
protected |
Definition at line 407 of file mixtureKEpsilon.C.
References virtualMassModel::Cvm(), and fluid.
|
protected |
Definition at line 420 of file mixtureKEpsilon.C.
References alphal.
|
protected |
Definition at line 431 of file mixtureKEpsilon.C.
References alphal.
|
protected |
Definition at line 445 of file mixtureKEpsilon.C.
References alphal.
|
protected |
Definition at line 461 of file mixtureKEpsilon.C.
References alphal, and Foam::fvc::interpolate().
|
protected |
Definition at line 482 of file mixtureKEpsilon.C.
References drag, fluid, Foam::mag(), Foam::pow(), Foam::pow3(), liquid::rho(), and U.
|
protectedvirtual |
Definition at line 520 of file mixtureKEpsilon.C.
References Foam::fvm::Su().
|
protectedvirtual |
Definition at line 527 of file mixtureKEpsilon.C.
References Foam::fvm::Su().
|
inlineprotected |
Return the effective diffusivity for k.
Definition at line 181 of file mixtureKEpsilon.H.
References mixtureKEpsilon< BasicTurbulenceModel >::sigmak_.
|
inlineprotected |
Return the effective diffusivity for epsilon.
Definition at line 196 of file mixtureKEpsilon.H.
References mixtureKEpsilon< BasicTurbulenceModel >::sigmaEps_.
TypeName | ( | "mixtureKEpsilon< BasicTurbulenceModel >" | ) |
Runtime type information.
|
virtual |
Re-read model coefficients if they have changed.
Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.
Definition at line 309 of file mixtureKEpsilon.C.
References Foam::read().
|
inlinevirtual |
Return the turbulence kinetic energy.
Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.
Definition at line 255 of file mixtureKEpsilon.H.
References mixtureKEpsilon< BasicTurbulenceModel >::k_.
|
inlinevirtual |
Return the turbulence kinetic energy dissipation rate.
Reimplemented from RASModel< BasicTurbulenceModel >.
Definition at line 263 of file mixtureKEpsilon.H.
References mixtureKEpsilon< BasicTurbulenceModel >::epsilon_.
|
virtual |
Solve the turbulence equations and correct the turbulence viscosity.
Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.
Definition at line 534 of file mixtureKEpsilon.C.
References Foam::fvc::absolute(), alphal, Foam::bound(), GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), tmp< T >::clear(), correct, GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), mixtureKEpsilon< BasicTurbulenceModel >::correctNut(), Foam::fvm::ddt(), Foam::devTwoSymm(), Foam::fvm::div(), Foam::fvc::div(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon_, GeometricBoundaryField< Type, PatchField, GeoMesh >::evaluateCoupled(), fluid, fvOptions, Foam::fvc::grad(), mixtureKEpsilon< BasicTurbulenceModel >::k_, Foam::fvm::laplacian(), options::New(), nu, eddyViscosity< RASModel< BasicTurbulenceModel > >::nut_, phi, phig(), tmp< T >::ref(), Foam::solve(), Foam::fvm::Sp(), Foam::fvm::SuSp(), and GeometricBoundaryField< Type, PatchField, GeoMesh >::updateCoeffs().
|
protected |
Definition at line 115 of file mixtureKEpsilon.H.
|
protected |
Definition at line 116 of file mixtureKEpsilon.H.
|
protected |
Definition at line 117 of file mixtureKEpsilon.H.
|
protected |
Definition at line 118 of file mixtureKEpsilon.H.
|
protected |
Definition at line 119 of file mixtureKEpsilon.H.
|
protected |
Definition at line 120 of file mixtureKEpsilon.H.
Referenced by mixtureKEpsilon< BasicTurbulenceModel >::DkEff().
|
protected |
Definition at line 121 of file mixtureKEpsilon.H.
Referenced by mixtureKEpsilon< BasicTurbulenceModel >::DepsilonEff().
|
protected |
Definition at line 125 of file mixtureKEpsilon.H.
Referenced by mixtureKEpsilon< BasicTurbulenceModel >::correct(), mixtureKEpsilon< BasicTurbulenceModel >::Ct2(), mixtureKEpsilon< BasicTurbulenceModel >::initMixtureFields(), and mixtureKEpsilon< BasicTurbulenceModel >::k().
|
protected |
|
protected |
Definition at line 130 of file mixtureKEpsilon.H.
|
protected |
Definition at line 131 of file mixtureKEpsilon.H.
|
protected |
Definition at line 132 of file mixtureKEpsilon.H.
|
protected |
Definition at line 133 of file mixtureKEpsilon.H.