mixtureKEpsilon< BasicTurbulenceModel > Class Template Reference

Mixture k-epsilon turbulence model for two-phase gas-liquid systems. More...

Inheritance diagram for mixtureKEpsilon< BasicTurbulenceModel >:
[legend]
Collaboration diagram for mixtureKEpsilon< BasicTurbulenceModel >:
[legend]

Public Types

typedef BasicTurbulenceModel::alphaField alphaField
 
typedef BasicTurbulenceModel::rhoField rhoField
 
typedef BasicTurbulenceModel::transportModel transportModel
 
- Public Types inherited from eddyViscosity< RASModel< BasicTurbulenceModel > >
typedef BasicTurbulenceModel::alphaField alphaField
 
typedef BasicTurbulenceModel::rhoField rhoField
 
typedef BasicTurbulenceModel::transportModel transportModel
 
- Public Types inherited from linearViscousStress< 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< volScalarFieldk () const
 Return the turbulence kinetic energy. More...
 
virtual tmp< volScalarFieldepsilon () 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 bool read ()=0
 Re-read model coefficients if they have changed. More...
 
virtual tmp< volScalarFieldnut () const
 Return the turbulence viscosity. More...
 
virtual tmp< scalarFieldnut (const label patchi) const
 Return the turbulence viscosity on patch. More...
 
virtual tmp< volScalarFieldk () const=0
 Return the turbulence kinetic energy. More...
 
virtual tmp< volSymmTensorFieldR () const
 Return the Reynolds stress tensor. More...
 
virtual void validate ()
 Validate the turbulence fields after construction. More...
 
virtual void correct ()=0
 Solve the turbulence equations and correct the turbulence viscosity. More...
 
- Public Member Functions inherited from linearViscousStress< 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 bool read ()=0
 Re-read model coefficients if they have changed. More...
 
virtual tmp< volSymmTensorFielddevRhoReff () const
 Return the effective stress tensor. More...
 
virtual tmp< volSymmTensorFielddevRhoReff (const volVectorField &U) const
 Return the effective stress tensor based on a given velocity field. More...
 
virtual tmp< fvVectorMatrixdivDevRhoReff (volVectorField &U) const
 Return the source term for the momentum equation. More...
 
virtual tmp< fvVectorMatrixdivDevRhoReff (const volScalarField &rho, volVectorField &U) const
 Return the source term for the momentum equation. More...
 
virtual void correct ()=0
 Solve the turbulence equations and correct the turbulence viscosity. More...
 

Protected Member Functions

wordList epsilonBoundaryTypes (const volScalarField &epsilon) const
 
void correctInletOutlet (volScalarField &vsf, const volScalarField &refVsf) const
 
void initMixtureFields ()
 
virtual void correctNut ()
 
tmp< volScalarFieldCt2 () const
 
tmp< volScalarFieldrholEff () const
 
tmp< volScalarFieldrhogEff () const
 
tmp< volScalarFieldrhom () const
 
tmp< volScalarFieldmix (const volScalarField &fc, const volScalarField &fd) const
 
tmp< volScalarFieldmixU (const volScalarField &fc, const volScalarField &fd) const
 
tmp< surfaceScalarFieldmixFlux (const surfaceScalarField &fc, const surfaceScalarField &fd) const
 
tmp< volScalarFieldbubbleG () const
 
virtual tmp< fvScalarMatrixkSource () const
 
virtual tmp< fvScalarMatrixepsilonSource () const
 
tmp< volScalarFieldDkEff (const volScalarField &nutm) const
 Return the effective diffusivity for k. More...
 
tmp< volScalarFieldDepsilonEff (const volScalarField &nutm) const
 Return the effective diffusivity for epsilon. More...
 
virtual void correctNut ()=0
 

Protected Attributes

dimensionedScalar Cmu_
 
dimensionedScalar C1_
 
dimensionedScalar C2_
 
dimensionedScalar C3_
 
dimensionedScalar Cp_
 
dimensionedScalar sigmak_
 
dimensionedScalar sigmaEps_
 
volScalarField k_
 
volScalarField epsilon_
 
autoPtr< volScalarFieldCt2_
 
autoPtr< volScalarFieldrhom_
 
autoPtr< volScalarFieldkm_
 
autoPtr< volScalarFieldepsilonm_
 
- Protected Attributes inherited from eddyViscosity< RASModel< BasicTurbulenceModel > >
volScalarField nut_
 

Detailed Description

template<class BasicTurbulenceModel>
class Foam::RASModels::mixtureKEpsilon< BasicTurbulenceModel >

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;
    }
Source files

Definition at line 86 of file mixtureKEpsilon.H.

Member Typedef Documentation

◆ alphaField

typedef BasicTurbulenceModel::alphaField alphaField

Definition at line 204 of file mixtureKEpsilon.H.

◆ rhoField

typedef BasicTurbulenceModel::rhoField rhoField

Definition at line 205 of file mixtureKEpsilon.H.

◆ transportModel

typedef BasicTurbulenceModel::transportModel transportModel

Definition at line 206 of file mixtureKEpsilon.H.

Constructor & Destructor Documentation

◆ mixtureKEpsilon()

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 49 of file mixtureKEpsilon.C.

References Foam::bound(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon_, mixtureKEpsilon< BasicTurbulenceModel >::k_, and Foam::type().

Here is the call graph for this function:

◆ ~mixtureKEpsilon()

virtual ~mixtureKEpsilon ( )
virtualdefault

Destructor.

Member Function Documentation

◆ epsilonBoundaryTypes()

wordList epsilonBoundaryTypes ( const volScalarField epsilon) const
protected

Definition at line 175 of file mixtureKEpsilon.C.

References epsilon, forAll, and GeometricBoundaryField< Type, PatchField, GeoMesh >::types().

Here is the call graph for this function:

◆ correctInletOutlet()

void correctInletOutlet ( volScalarField vsf,
const volScalarField refVsf 
) const
protected

Definition at line 197 of file mixtureKEpsilon.C.

References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), and forAll.

Here is the call graph for this function:

◆ initMixtureFields()

void initMixtureFields
protected

Definition at line 225 of file mixtureKEpsilon.C.

References IOobject::AUTO_WRITE, GeometricField< Type, PatchField, GeoMesh >::boundaryField(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon_, mixtureKEpsilon< BasicTurbulenceModel >::k_, IOobject::READ_IF_PRESENT, and GeometricBoundaryField< Type, PatchField, GeoMesh >::types().

Here is the call graph for this function:

◆ correctNut()

void correctNut
protectedvirtual

Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.

Definition at line 336 of file mixtureKEpsilon.C.

References Time::New(), and Foam::sqr().

Referenced by mixtureKEpsilon< BasicTurbulenceModel >::correct().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ Ct2()

tmp< volScalarField > Ct2
protected

Definition at line 378 of file mixtureKEpsilon.C.

References beta(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon_, Foam::exp(), fluid, mixtureKEpsilon< BasicTurbulenceModel >::k_, Foam::mag(), liquid::rho(), Foam::sqr(), and Foam::sqrt().

Here is the call graph for this function:

◆ rholEff()

tmp< volScalarField > rholEff
protected

Definition at line 405 of file mixtureKEpsilon.C.

References fluid.

◆ rhogEff()

tmp< volScalarField > rhogEff
protected

Definition at line 414 of file mixtureKEpsilon.C.

References virtualMassModel::Cvm(), and fluid.

Here is the call graph for this function:

◆ rhom()

tmp< volScalarField > rhom
protected

Definition at line 427 of file mixtureKEpsilon.C.

References alphal.

◆ mix()

tmp< volScalarField > mix ( const volScalarField fc,
const volScalarField fd 
) const
protected

Definition at line 437 of file mixtureKEpsilon.C.

References alphal.

◆ mixU()

tmp< volScalarField > mixU ( const volScalarField fc,
const volScalarField fd 
) const
protected

Definition at line 451 of file mixtureKEpsilon.C.

References alphal.

◆ mixFlux()

tmp< surfaceScalarField > mixFlux ( const surfaceScalarField fc,
const surfaceScalarField fd 
) const
protected

Definition at line 467 of file mixtureKEpsilon.C.

References alphal, and Foam::fvc::interpolate().

Here is the call graph for this function:

◆ bubbleG()

tmp< volScalarField > bubbleG
protected

Definition at line 489 of file mixtureKEpsilon.C.

References drag, fluid, Foam::mag(), Foam::pow(), Foam::pow3(), and liquid::rho().

Here is the call graph for this function:

◆ kSource()

tmp< fvScalarMatrix > kSource
protectedvirtual

Definition at line 527 of file mixtureKEpsilon.C.

References Foam::fvm::Su().

Here is the call graph for this function:

◆ epsilonSource()

tmp< fvScalarMatrix > epsilonSource
protectedvirtual

Definition at line 534 of file mixtureKEpsilon.C.

References Foam::fvm::Su().

Here is the call graph for this function:

◆ DkEff()

tmp< volScalarField > DkEff ( const volScalarField nutm) const
inlineprotected

Return the effective diffusivity for k.

Definition at line 177 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicTurbulenceModel >::sigmak_.

◆ DepsilonEff()

tmp< volScalarField > DepsilonEff ( const volScalarField nutm) const
inlineprotected

Return the effective diffusivity for epsilon.

Definition at line 190 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicTurbulenceModel >::sigmaEps_.

◆ TypeName()

TypeName ( "mixtureKEpsilon< BasicTurbulenceModel >"  )

Runtime type information.

◆ read()

bool read
virtual

Re-read model coefficients if they have changed.

Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.

Definition at line 316 of file mixtureKEpsilon.C.

References Foam::read().

Here is the call graph for this function:

◆ k()

virtual tmp< volScalarField > k ( ) const
inlinevirtual

Return the turbulence kinetic energy.

Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.

Definition at line 239 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicTurbulenceModel >::k_.

◆ epsilon()

virtual tmp< volScalarField > epsilon ( ) const
inlinevirtual

Return the turbulence kinetic energy dissipation rate.

Definition at line 245 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicTurbulenceModel >::epsilon_.

◆ correct()

Member Data Documentation

◆ Cmu_

dimensionedScalar Cmu_
protected

Definition at line 113 of file mixtureKEpsilon.H.

◆ C1_

dimensionedScalar C1_
protected

Definition at line 114 of file mixtureKEpsilon.H.

◆ C2_

dimensionedScalar C2_
protected

Definition at line 115 of file mixtureKEpsilon.H.

◆ C3_

dimensionedScalar C3_
protected

Definition at line 116 of file mixtureKEpsilon.H.

◆ Cp_

dimensionedScalar Cp_
protected

Definition at line 117 of file mixtureKEpsilon.H.

◆ sigmak_

dimensionedScalar sigmak_
protected

Definition at line 118 of file mixtureKEpsilon.H.

Referenced by mixtureKEpsilon< BasicTurbulenceModel >::DkEff().

◆ sigmaEps_

dimensionedScalar sigmaEps_
protected

◆ k_

◆ epsilon_

◆ Ct2_

autoPtr<volScalarField> Ct2_
protected

Definition at line 128 of file mixtureKEpsilon.H.

◆ rhom_

autoPtr<volScalarField> rhom_
protected

Definition at line 129 of file mixtureKEpsilon.H.

◆ km_

autoPtr<volScalarField> km_
protected

Definition at line 130 of file mixtureKEpsilon.H.

◆ epsilonm_

autoPtr<volScalarField> epsilonm_
protected

Definition at line 131 of file mixtureKEpsilon.H.


The documentation for this class was generated from the following files: