Base class implementation of the k-omega-SST turbulence model for incompressible and compressible flows. More...

Inheritance diagram for kOmegaSSTBase< BasicEddyViscosityModel >:

Collaboration diagram for kOmegaSSTBase< BasicEddyViscosityModel >:

Public Types
typedef BasicEddyViscosityModel::alphaField	alphaField

typedef BasicEddyViscosityModel::rhoField	rhoField

typedef BasicEddyViscosityModel::transportModel	transportModel

Public Member Functions
	kOmegaSSTBase (const kOmegaSSTBase &)=delete
	No copy construct. More...

void	operator= (const kOmegaSSTBase &)=delete
	No copy assignment. More...

	kOmegaSSTBase (const word &type, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName=turbulenceModel::propertiesName)
	Construct from components. More...

virtual	~kOmegaSSTBase ()=default
	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 >	omega () const
	Return the turbulence kinetic energy dissipation rate. More...

virtual void	correct ()
	Solve the turbulence equations and correct the turbulence viscosity. More...

Protected Member Functions
void	setDecayControl (const dictionary &dict)
	Set decay control with kInf and omegaInf. More...

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
	Return the blended field. More...

tmp< volScalarField::Internal >	blend (const volScalarField::Internal &F1, const dimensionedScalar &psi1, const dimensionedScalar &psi2) const
	Return the internal blended field. More...

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)

virtual void	correctNut ()

virtual tmp< volScalarField >	S2 (const volTensorField &gradU) const
	Return square of strain rate. More...

virtual tmp< volScalarField::Internal >	Pk (const volScalarField::Internal &G) const
	Return k production rate. More...

virtual tmp< volScalarField::Internal >	epsilonByk (const volScalarField &F1, const volTensorField &gradU) const
	Return epsilon/k which for standard RAS is betaStar*omega. More...

virtual tmp< volScalarField::Internal >	GbyNu0 (const volTensorField &gradU, const volScalarField &S2) const
	Return (G/nu)_0. More...

virtual tmp< volScalarField::Internal >	GbyNu (const volScalarField::Internal &GbyNu0, const volScalarField::Internal &F2, const volScalarField::Internal &S2) const
	Return G/nu. 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_
	Flag to include the F3 term. More...

const volScalarField &	y_
	Wall distance. More...

volScalarField	k_
	Turbulent kinetic energy field [m^2/s^2]. More...

volScalarField	omega_
	Specific dissipation rate field [1/s]. More...

Switch	decayControl_
	Flag to include the decay control. More...

dimensionedScalar	kInf_

dimensionedScalar	omegaInf_

Detailed Description

template<class BasicEddyViscosityModel>
class Foam::kOmegaSSTBase< BasicEddyViscosityModel >

Base class 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.

and the optional decay control from:

    Spalart, P. R. and Rumsey, C. L. (2007).
    Effective Inflow Conditions for Turbulence Models in Aerodynamic
    Calculations
    AIAA Journal, 45(10), 2544 - 2553.

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 coefficients 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

    kOmegaSSTBaseCoeffs
    {
        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;

        // Optional decay control
        decayControl    yes;
        kInf            \<far-field k value\>;
        omegaInf        \<far-field omega value\>;
    }

Source files

Definition at line 125 of file kOmegaSSTBase.H.

Member Typedef Documentation

◆ alphaField

typedef BasicEddyViscosityModel::alphaField alphaField

Definition at line 324 of file kOmegaSSTBase.H.

◆ rhoField

typedef BasicEddyViscosityModel::rhoField rhoField

Definition at line 325 of file kOmegaSSTBase.H.

◆ transportModel

typedef BasicEddyViscosityModel::transportModel transportModel

Definition at line 326 of file kOmegaSSTBase.H.

Constructor & Destructor Documentation

◆ kOmegaSSTBase() [1/2]

kOmegaSSTBase ( const kOmegaSSTBase< BasicEddyViscosityModel > & )

delete

No copy construct.

◆ kOmegaSSTBase() [2/2]

kOmegaSSTBase	(	const word &	type,
		const alphaField &	alpha,
		const rhoField &	rho,
		const volVectorField &	U,
		const surfaceScalarField &	alphaRhoPhi,
		const surfaceScalarField &	phi,
		const transportModel &	transport,
		const word &	propertiesName = `turbulenceModel::propertiesName`
	)

Construct from components.

Definition at line 233 of file kOmegaSSTBase.C.

◆ ~kOmegaSSTBase()

virtual ~kOmegaSSTBase ( )

virtualdefault

Destructor.

Member Function Documentation

◆ setDecayControl()

void setDecayControl ( const dictionary & dict )

protected

Set decay control with kInf and omegaInf.

Definition at line 442 of file kOmegaSSTBase.C.

◆ F1()

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

protectedvirtual

Reimplemented in kOmegaSSTLM< BasicTurbulenceModel >.

Definition at line 37 of file kOmegaSSTBase.C.

Referenced by kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::alphaK(), kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::alphaOmega(), kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::blend(), kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::DkEff(), and kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::DomegaEff().

Here is the caller graph for this function:

◆ F2()

tmp< volScalarField > F2 ( ) const

protectedvirtual

Definition at line 66 of file kOmegaSSTBase.C.

◆ F3()

tmp< volScalarField > F3 ( ) const

protectedvirtual

Definition at line 83 of file kOmegaSSTBase.C.

◆ F23()

tmp< volScalarField > F23 ( ) const

protectedvirtual

Definition at line 96 of file kOmegaSSTBase.C.

◆ blend() [1/2]

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

inlineprotected

Return the blended field.

Definition at line 206 of file kOmegaSSTBase.H.

Referenced by kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::alphaK(), kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::alphaOmega(), kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::beta(), and kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::gamma().

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

Return the internal blended field.

Definition at line 219 of file kOmegaSSTBase.H.

◆ alphaK()

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

inlineprotected

Definition at line 228 of file kOmegaSSTBase.H.

Referenced by kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::DkEff().

Here is the caller graph for this function:

◆ alphaOmega()

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

inlineprotected

Definition at line 233 of file kOmegaSSTBase.H.

Referenced by kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::DomegaEff().

Here is the caller graph for this function:

◆ beta()

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

inlineprotected

Definition at line 239 of file kOmegaSSTBase.H.

◆ gamma()

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

inlineprotected

Definition at line 251 of file kOmegaSSTBase.H.

◆ correctNut() [1/2]

void correctNut ( const volScalarField & S2 )

protectedvirtual

Reimplemented in kOmegaSSTSato< BasicTurbulenceModel >, kOmegaSST< BasicTurbulenceModel >, and kOmegaSSTDES< BasicTurbulenceModel >.

Definition at line 111 of file kOmegaSSTBase.C.

◆ correctNut() [2/2]

void correctNut ( )

protectedvirtual

Reimplemented in kOmegaSST< BasicTurbulenceModel >, and kOmegaSSTDES< BasicTurbulenceModel >.

Definition at line 123 of file kOmegaSSTBase.C.

◆ S2()

Foam::tmp< Foam::volScalarField > S2 ( const volTensorField & gradU ) const

protectedvirtual

Return square of strain rate.

Reimplemented in kOmegaSSTDES< BasicTurbulenceModel >, and kOmegaSSTDDES< BasicTurbulenceModel >.

Definition at line 131 of file kOmegaSSTBase.C.

◆ Pk()

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

protectedvirtual

Return k production rate.

Definition at line 141 of file kOmegaSSTBase.C.

◆ epsilonByk()

tmp< volScalarField::Internal > epsilonByk	(	const volScalarField &	F1,
		const volTensorField &	gradU
	)		const

protectedvirtual

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

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

Definition at line 151 of file kOmegaSSTBase.C.

◆ GbyNu0()

tmp< volScalarField::Internal > GbyNu0	(	const volTensorField &	gradU,
		const volScalarField &	S2
	)		const

protectedvirtual

Return (G/nu)_0.

Reimplemented in kOmegaSSTDES< BasicTurbulenceModel >, and kOmegaSSTDDES< BasicTurbulenceModel >.

Definition at line 162 of file kOmegaSSTBase.C.

◆ GbyNu()

tmp< volScalarField::Internal > GbyNu	(	const volScalarField::Internal &	GbyNu0,
		const volScalarField::Internal &	F2,
		const volScalarField::Internal &	S2
	)		const