Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel > Class Template Reference

Implementation of the k-omega-SST turbulence model for incompressible and compressible flows. More...

Inheritance diagram for kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >:
Inheritance graph
[legend]
Collaboration diagram for kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >:
Collaboration graph
[legend]

Public Types

typedef BasicMomentumTransportModel::alphaField alphaField
 
typedef BasicMomentumTransportModel::rhoField rhoField
 

Public Member Functions

 kOmegaSST (const word &type, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const viscosity &viscosity)
 Construct from components. More...
 
 kOmegaSST (const kOmegaSST &)=delete
 Disallow default bitwise copy construction. More...
 
virtual ~kOmegaSST ()
 Destructor. More...
 
virtual bool read ()
 Re-read model coefficients if they have changed. More...
 
tmp< volScalarField > DkEff (const volScalarField &F1) const
 Return the effective diffusivity for k. More...
 
tmp< volScalarField > DomegaEff (const volScalarField &F1) const
 Return the effective diffusivity for omega. 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 tmp< volScalarField > omega () const
 Return the turbulence kinetic energy dissipation rate. More...
 
virtual void correct ()
 Solve the turbulence equations and correct the turbulence viscosity. More...
 
void operator= (const kOmegaSST &)=delete
 Disallow default bitwise assignment. More...
 

Protected Member Functions

virtual tmp< volScalarField > F1 (const volScalarField &CDkOmega) const
 
virtual tmp< volScalarField > F2 () const
 
virtual tmp< volScalarField > F3 () const
 
virtual tmp< volScalarField > F23 () const
 
tmp< volScalarField > blend (const volScalarField &F1, const dimensionedScalar &psi1, const dimensionedScalar &psi2) const
 
tmp< volScalarField::Internal > blend (const volScalarField::Internal &F1, const dimensionedScalar &psi1, const dimensionedScalar &psi2) const
 
tmp< volScalarField > alphaK (const volScalarField &F1) const
 
tmp< volScalarField > alphaOmega (const volScalarField &F1) const
 
tmp< volScalarField::Internal > beta (const volScalarField::Internal &F1) const
 
tmp< volScalarField::Internal > gamma (const volScalarField::Internal &F1) const
 
virtual void correctNut (const volScalarField &S2, const volScalarField &F2)
 
virtual void correctNut ()
 
virtual tmp< volScalarField::Internal > Pk (const volScalarField::Internal &G) const
 Return k production rate. More...
 
virtual tmp< volScalarField::Internal > epsilonByk (const volScalarField::Internal &F1, const volScalarField::Internal &F2) const
 Return epsilon/k which for standard RAS is betaStar*omega. More...
 
virtual tmp< fvScalarMatrix > kSource () const
 
virtual tmp< fvScalarMatrix > omegaSource () const
 
virtual tmp< fvScalarMatrix > Qsas (const volScalarField::Internal &S2, const volScalarField::Internal &gamma, const volScalarField::Internal &beta) const
 

Protected Attributes

dimensionedScalar alphaK1_
 
dimensionedScalar alphaK2_
 
dimensionedScalar alphaOmega1_
 
dimensionedScalar alphaOmega2_
 
dimensionedScalar gamma1_
 
dimensionedScalar gamma2_
 
dimensionedScalar beta1_
 
dimensionedScalar beta2_
 
dimensionedScalar betaStar_
 
dimensionedScalar a1_
 
dimensionedScalar b1_
 
dimensionedScalar c1_
 
Switch F3_
 
const volScalarField & y_
 Wall distance. More...
 
volScalarField k_
 
volScalarField omega_
 

Detailed Description

template<class MomentumTransportModel, class BasicMomentumTransportModel>
class Foam::kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >

Implementation of the k-omega-SST turbulence model for incompressible and compressible flows.

Turbulence model described in: 

    Menter, F. R. & Esch, T. (2001).
    Elements of Industrial Heat Transfer Prediction.
    16th Brazilian Congress of Mechanical Engineering (COBEM).

with updated coefficients from 

    Menter, F. R., Kuntz, M., and Langtry, R. (2003).
    Ten Years of Industrial Experience with the SST Turbulence Model.
    Turbulence, Heat and Mass Transfer 4, ed: K. Hanjalic, Y. Nagano,
    & M. Tummers, Begell House, Inc., 625 - 632.

but with the consistent production terms from the 2001 paper as form in the 2003 paper is a typo, see 

    http://turbmodels.larc.nasa.gov/sst.html

and the addition of the optional F3 term for rough walls from 

    Hellsten, A. (1998).
    "Some Improvements in Menter’s k-omega-SST turbulence model"
    29th AIAA Fluid Dynamics Conference, AIAA-98-2554.

Note that this implementation is written in terms of alpha diffusion coefficients rather than the more traditional sigma (alpha = 1/sigma) so that the blending can be applied to all coefficuients in a consistent manner. The paper suggests that sigma is blended but this would not be consistent with the blending of the k-epsilon and k-omega models.

Also note that the error in the last term of equation (2) relating to sigma has been corrected.

Wall-functions are applied in this implementation by using equations (14) to specify the near-wall omega as appropriate.

The blending functions (15) and (16) are not currently used because of the uncertainty in their origin, range of applicability and that if y+ becomes sufficiently small blending u_tau in this manner clearly becomes nonsense.

The default model coefficients are 

    kOmegaSSTCoeffs
    {
        alphaK1     0.85;
        alphaK2     1.0;
        alphaOmega1 0.5;
        alphaOmega2 0.856;
        beta1       0.075;
        beta2       0.0828;
        betaStar    0.09;
        gamma1      5/9;
        gamma2      0.44;
        a1          0.31;
        b1          1.0;
        c1          10.0;
        F3          no;
    }
Source files

Definition at line 112 of file kOmegaSSTBase.H.

Member Typedef Documentation

◆ alphaField

typedef BasicMomentumTransportModel::alphaField alphaField

Definition at line 242 of file kOmegaSSTBase.H.

◆ rhoField

typedef BasicMomentumTransportModel::rhoField rhoField

Definition at line 243 of file kOmegaSSTBase.H.

Constructor & Destructor Documentation

◆ kOmegaSST() [1/2]

kOmegaSST ( const word &  type,
const alphaField &  alpha,
const rhoField &  rho,
const volVectorField &  U,
const surfaceScalarField &  alphaRhoPhi,
const surfaceScalarField &  phi,
const viscosity &  viscosity 
)

Construct from components.

Definition at line 217 of file kOmegaSSTBase.C.

◆ kOmegaSST() [2/2]

kOmegaSST ( const kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel > &  )
delete

Disallow default bitwise copy construction.

◆ ~kOmegaSST()

virtual ~kOmegaSST ( )
inlinevirtual

Destructor.

Reimplemented in kOmegaSST< BasicMomentumTransportModel >.

Definition at line 265 of file kOmegaSSTBase.H.

Member Function Documentation

◆ F1()

virtual tmp<volScalarField> F1 ( const volScalarField &  CDkOmega) const
protectedvirtual

Reimplemented in kOmegaSSTLM< BasicMomentumTransportModel >.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaK(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaOmega(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::beta(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::blend(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::DkEff(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::DomegaEff(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::gamma().

Here is the caller graph for this function:

◆ F2()

virtual tmp<volScalarField> F2 ( ) const
protectedvirtual

◆ F3()

virtual tmp<volScalarField> F3 ( ) const
protectedvirtual

◆ F23()

virtual tmp<volScalarField> F23 ( ) const
protectedvirtual

◆ blend() [1/2]

tmp<volScalarField> blend ( const volScalarField &  F1,
const dimensionedScalar &  psi1,
const dimensionedScalar &  psi2 
) const
inlineprotected

Definition at line 161 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1().

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaK(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaOmega(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::beta(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::gamma().

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

◆ blend() [2/2]

tmp<volScalarField::Internal> blend ( const volScalarField::Internal &  F1,
const dimensionedScalar &  psi1,
const dimensionedScalar &  psi2 
) const
inlineprotected

Definition at line 171 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1().

Here is the call graph for this function:

◆ alphaK()

tmp<volScalarField> alphaK ( const volScalarField &  F1) const
inlineprotected

Definition at line 181 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaK1_, kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaK2_, kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::blend(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1().

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::DkEff().

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

◆ alphaOmega()

tmp<volScalarField> alphaOmega ( const volScalarField &  F1) const
inlineprotected

Definition at line 186 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaOmega1_, kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaOmega2_, kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::blend(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1().

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::DomegaEff().

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

◆ beta()

tmp<volScalarField::Internal> beta ( const volScalarField::Internal &  F1) const
inlineprotected

Definition at line 191 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::beta1_, kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::beta2_, kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::blend(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1().

Here is the call graph for this function:

◆ gamma()

tmp<volScalarField::Internal> gamma ( const volScalarField::Internal &  F1) const
inlineprotected

Definition at line 199 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::blend(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::gamma1_, and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::gamma2_.

Here is the call graph for this function:

◆ correctNut() [1/2]

void correctNut ( const volScalarField &  S2,
const volScalarField &  F2 
)
protectedvirtual

Reimplemented in kOmegaSSTSato< BasicMomentumTransportModel >.

Definition at line 118 of file kOmegaSSTBase.C.

References F2, Foam::max(), dictionary::New(), and Foam::sqrt().

Here is the call graph for this function:

◆ correctNut() [2/2]

void correctNut
protectedvirtual

Definition at line 133 of file kOmegaSSTBase.C.

References Foam::fvc::grad(), Foam::magSqr(), and Foam::symm().

Here is the call graph for this function:

◆ Pk()

tmp< volScalarField::Internal > Pk ( const volScalarField::Internal &  G) const
protectedvirtual

Return k production rate.

Reimplemented in kOmegaSSTLM< BasicMomentumTransportModel >.

Definition at line 142 of file kOmegaSSTBase.C.

References Foam::constant::universal::G, and Foam::min().

Here is the call graph for this function:

◆ epsilonByk()

tmp< volScalarField::Internal > epsilonByk ( const volScalarField::Internal &  F1,
const volScalarField::Internal &  F2 
) const
protectedvirtual

Return epsilon/k which for standard RAS is betaStar*omega.

Reimplemented in kOmegaSSTLM< BasicMomentumTransportModel >, and kOmegaSSTDES< BasicMomentumTransportModel >.

Definition at line 153 of file kOmegaSSTBase.C.

◆ kSource()

tmp< fvScalarMatrix > kSource
protectedvirtual

Definition at line 165 of file kOmegaSSTBase.C.

References Foam::dimTime, and Foam::dimVolume.

◆ omegaSource()

tmp< fvScalarMatrix > omegaSource
protectedvirtual

Definition at line 180 of file kOmegaSSTBase.C.

References Foam::dimTime, and Foam::dimVolume.

◆ Qsas()

tmp< fvScalarMatrix > Qsas ( const volScalarField::Internal &  S2,
const volScalarField::Internal &  gamma,
const volScalarField::Internal &  beta 
) const
protectedvirtual

Reimplemented in kOmegaSSTSAS< BasicMomentumTransportModel >.

Definition at line 196 of file kOmegaSSTBase.C.

◆ read()

bool read
virtual

Re-read model coefficients if they have changed.

Reimplemented in kOmegaSSTSato< BasicMomentumTransportModel >, kOmegaSSTSAS< BasicMomentumTransportModel >, kOmegaSSTLM< BasicMomentumTransportModel >, and kOmegaSSTDES< BasicMomentumTransportModel >.

Definition at line 392 of file kOmegaSSTBase.C.

References Foam::blockMeshTools::read().

Here is the call graph for this function:

◆ DkEff()

tmp<volScalarField> DkEff ( const volScalarField &  F1) const
inline

Return the effective diffusivity for k.

Definition at line 275 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaK(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1(), and GeometricField< Type, PatchField, GeoMesh >::New().

Here is the call graph for this function:

◆ DomegaEff()

tmp<volScalarField> DomegaEff ( const volScalarField &  F1) const
inline

Return the effective diffusivity for omega.

Definition at line 285 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaOmega(), kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::F1(), and GeometricField< Type, PatchField, GeoMesh >::New().

Here is the call graph for this function:

◆ k()

virtual tmp<volScalarField> k ( ) const
inlinevirtual

Return the turbulence kinetic energy.

Definition at line 295 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::k_.

◆ epsilon()

virtual tmp<volScalarField> epsilon ( ) const
inlinevirtual

Return the turbulence kinetic energy dissipation rate.

Definition at line 301 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::betaStar_, kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::k_, GeometricField< Type, PatchField, GeoMesh >::New(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::omega_.

Here is the call graph for this function:

◆ omega()

virtual tmp<volScalarField> omega ( ) const
inlinevirtual

Return the turbulence kinetic energy dissipation rate.

Definition at line 311 of file kOmegaSSTBase.H.

References kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::omega_.

◆ correct()

void correct
virtual

Solve the turbulence equations and correct the turbulence viscosity.

Reimplemented in kOmegaSSTSato< BasicMomentumTransportModel >, and kOmegaSSTLM< BasicMomentumTransportModel >.

Definition at line 420 of file kOmegaSSTBase.C.

References Foam::fvc::absolute(), alpha(), Foam::bound(), tmp< T >::clear(), fvConstraints::constrain(), Foam::MULES::correct(), Foam::fvm::ddt(), Foam::dev(), Foam::fvc::div(), Foam::fvm::div(), F1, fvConstraints(), fvModels(), Foam::constant::universal::G, Foam::fvc::grad(), Foam::fvm::laplacian(), Foam::magSqr(), Foam::max(), Foam::min(), dictionary::New(), nut, tmp< T >::ref(), rho, Foam::solve(), Foam::fvm::Sp(), Foam::sqrt(), Foam::fvm::SuSp(), Foam::symm(), Foam::twoSymm(), and U.

Here is the call graph for this function:

◆ operator=()

void operator= ( const kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel > &  )
delete

Disallow default bitwise assignment.

Member Data Documentation

◆ alphaK1_

dimensionedScalar alphaK1_
protected

Definition at line 122 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaK().

◆ alphaK2_

dimensionedScalar alphaK2_
protected

Definition at line 123 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaK().

◆ alphaOmega1_

dimensionedScalar alphaOmega1_
protected

Definition at line 125 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaOmega().

◆ alphaOmega2_

dimensionedScalar alphaOmega2_
protected

Definition at line 126 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::alphaOmega().

◆ gamma1_

dimensionedScalar gamma1_
protected

Definition at line 128 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::gamma().

◆ gamma2_

dimensionedScalar gamma2_
protected

Definition at line 129 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::gamma().

◆ beta1_

dimensionedScalar beta1_
protected

Definition at line 131 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::beta().

◆ beta2_

dimensionedScalar beta2_
protected

Definition at line 132 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::beta().

◆ betaStar_

dimensionedScalar betaStar_
protected

Definition at line 134 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::epsilon().

◆ a1_

dimensionedScalar a1_
protected

Definition at line 136 of file kOmegaSSTBase.H.

◆ b1_

dimensionedScalar b1_
protected

Definition at line 137 of file kOmegaSSTBase.H.

◆ c1_

dimensionedScalar c1_
protected

Definition at line 138 of file kOmegaSSTBase.H.

◆ F3_

Switch F3_
protected

Definition at line 140 of file kOmegaSSTBase.H.

◆ y_

const volScalarField& y_
protected

Wall distance.

Note: different to wall distance in parent RASModel which is for near-wall cells only

Definition at line 148 of file kOmegaSSTBase.H.

◆ k_

volScalarField k_
protected

Definition at line 150 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::epsilon(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::k().

◆ omega_

volScalarField omega_
protected

Definition at line 151 of file kOmegaSSTBase.H.

Referenced by kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::epsilon(), and kOmegaSST< MomentumTransportModel, BasicMomentumTransportModel >::omega().

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