k-omega-SST DES turbulence model for incompressible and compressible flows. More...

#include <kOmegaSSTDES.H>

Inheritance diagram for kOmegaSSTDES< BasicTurbulenceModel >:

[legend]

Collaboration diagram for kOmegaSSTDES< BasicTurbulenceModel >:

[legend]

Public Types
typedef BasicTurbulenceModel::alphaField	alphaField
typedef BasicTurbulenceModel::rhoField	rhoField
typedef BasicTurbulenceModel::transportModel	transportModel
Public Types inherited from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >
typedef DESModel< BasicTurbulenceModel >::alphaField	alphaField
typedef DESModel< BasicTurbulenceModel >::rhoField	rhoField
typedef DESModel< BasicTurbulenceModel >::transportModel	transportModel

Public Member Functions
	TypeName ("kOmegaSSTDES")
	Runtime type information.
	kOmegaSSTDES (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.
virtual	~kOmegaSSTDES ()=default
	Destructor.
virtual bool	read ()
	Re-read model coefficients if they have changed.
virtual tmp< volScalarField >	lengthScaleRAS () const
	RAS length scale.
virtual tmp< volScalarField >	lengthScaleLES (const volScalarField &CDES) const
	LES length scale.
virtual tmp< volScalarField >	LESRegion () const
	Return the LES field indicator.
Public Member Functions inherited from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >
	kOmegaSSTBase (const kOmegaSSTBase &)=delete
	No copy construct.
void	operator= (const kOmegaSSTBase &)=delete
	No copy assignment.
virtual	~kOmegaSSTBase ()=default
	Destructor.
tmp< volScalarField >	DkEff (const volScalarField &F1) const
	Return the effective diffusivity for k.
tmp< volScalarField >	DomegaEff (const volScalarField &F1) const
	Return the effective diffusivity for omega.
virtual tmp< volScalarField >	k () const
	Return the turbulence kinetic energy.
virtual tmp< volScalarField >	omega () const
	Return the turbulence kinetic energy dissipation rate.
virtual void	correct ()
	Solve the turbulence equations and correct the turbulence viscosity.

Protected Member Functions
virtual tmp< volScalarField >	CDES (const volScalarField &F1) const
	Blending for CDES parameter.
virtual void	correctNut (const volScalarField &S2)
virtual void	correctNut ()
tmp< volScalarField >	r (const volScalarField &nur, const volScalarField &magGradU) const
virtual tmp< volScalarField >	S2 (const volTensorField &gradU) const
	Return square of strain rate.
virtual tmp< volScalarField >	dTilda (const volScalarField &magGradU, const volScalarField &CDES) const
	Return length scale.
virtual tmp< volScalarField::Internal >	epsilonByk (const volScalarField &F1, const volTensorField &gradU) const
	Return epsilon/k.
virtual tmp< volScalarField::Internal >	GbyNu0 (const volTensorField &gradU, const volScalarField &S2) const
	Return (G/nu)_0.
virtual tmp< volScalarField::Internal >	GbyNu (const volScalarField::Internal &GbyNu0, const volScalarField::Internal &F2, const volScalarField::Internal &S2) const
	Return G/nu.
Protected Member Functions inherited from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >
void	setDecayControl (const dictionary &dict)
	Set decay control with kInf and omegaInf.
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.
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 tmp< volScalarField::Internal >	Pk (const volScalarField::Internal &G) const
	Return k production rate.
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
Switch	useSigma_
	Switch to activate grey-area enhanced sigma-(D)DES.
dimensionedScalar	kappa_
dimensionedScalar	CDESkom_
	DES coefficients.
dimensionedScalar	CDESkeps_
Protected Attributes inherited from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >
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.
const volScalarField &	y_
	Wall distance.
volScalarField	k_
	Turbulent kinetic energy field [m^2/s^2].
volScalarField	omega_
	Specific dissipation rate field [1/s].
Switch	decayControl_
	Flag to include the decay control.
dimensionedScalar	kInf_
dimensionedScalar	omegaInf_

Detailed Description

template<class BasicTurbulenceModel>
class Foam::LESModels::kOmegaSSTDES< BasicTurbulenceModel >

k-omega-SST DES turbulence model for incompressible and compressible flows.

Reference:

    Strelets, M. (2001).
    Detached Eddy Simulation of Massively Separated Flows.
    39th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV.

Note: The default values of the DES constants implemented are code-specific values calibrated for OpenFOAM using decaying isotropic turbulence, and hence deviate slightly from the values suggested in the reference publication.

Source files

Definition at line 67 of file kOmegaSSTDES.H.

Member Typedef Documentation

◆ alphaField

template<class BasicTurbulenceModel>

typedef BasicTurbulenceModel::alphaField alphaField

Definition at line 171 of file kOmegaSSTDES.H.

◆ rhoField

template<class BasicTurbulenceModel>

typedef BasicTurbulenceModel::rhoField rhoField

Definition at line 172 of file kOmegaSSTDES.H.

◆ transportModel

template<class BasicTurbulenceModel>

typedef BasicTurbulenceModel::transportModel transportModel

Definition at line 173 of file kOmegaSSTDES.H.

Constructor & Destructor Documentation

◆ kOmegaSSTDES()

template<class BasicTurbulenceModel>

kOmegaSSTDES	(	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 162 of file kOmegaSSTDES.C.

References alpha, CDESkeps_, CDESkom_, Foam::endl(), dimensioned< Type >::getOrAddToDict(), Foam::GlobalIOList< Tuple2< scalar, vector > >::typeName, kappa_, kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::kOmegaSSTBase(), phi, rho, U, useSigma_, and WarningInFunction.

Here is the call graph for this function:

◆ ~kOmegaSSTDES()

template<class BasicTurbulenceModel>

virtual ~kOmegaSSTDES ( )

virtualdefault

Destructor.

Member Function Documentation

◆ CDES()

template<class BasicTurbulenceModel>

virtual tmp< volScalarField > CDES ( const volScalarField & F1 ) const

inlineprotectedvirtual

Blending for CDES parameter.

Definition at line 109 of file kOmegaSSTDES.H.

Referenced by kOmegaSSTDDES< BasicTurbulenceModel >::dTilda(), dTilda(), kOmegaSSTIDDES< BasicTurbulenceModel >::dTilda(), epsilonByk(), lengthScaleLES(), LESRegion(), kOmegaSSTDDES< BasicTurbulenceModel >::S2(), and S2().

Here is the caller graph for this function:

◆ correctNut() [1/2]

template<class BasicTurbulenceModel>

void correctNut ( const volScalarField & S2 )

protectedvirtual

Reimplemented from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >.

Definition at line 35 of file kOmegaSSTDES.C.

References correctNut(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::kOmegaSSTBase(), and S2().

Referenced by correctNut(), and correctNut().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ correctNut() [2/2]

template<class BasicTurbulenceModel>

void correctNut ( )

protectedvirtual

Definition at line 46 of file kOmegaSSTDES.C.

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

Here is the call graph for this function:

◆ r()

template<class BasicTurbulenceModel>

tmp< volScalarField > r	(	const volScalarField &	nur,
		const volScalarField &	magGradU ) const

protected

Definition at line 53 of file kOmegaSSTDES.C.

References DimensionedField< Type, GeoMesh >::dimensions(), kappa_, Foam::max(), Foam::min(), Foam::sqr(), Foam::tr(), and kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::y_.

Here is the call graph for this function:

◆ S2()

template<class BasicTurbulenceModel>

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

protectedvirtual

Return square of strain rate.

Reimplemented from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >.

Reimplemented in kOmegaSSTDDES< BasicTurbulenceModel >, kOmegaSSTDDES< EddyDiffusivity< fluidThermoCompressibleTurbulenceModel > >, and kOmegaSSTDDES< Foam::EddyDiffusivity< Foam::fluidThermoCompressibleTurbulenceModel > >.

Definition at line 71 of file kOmegaSSTDES.C.

References kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::alphaOmega2_, CDES(), dTilda(), F1, kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::F1(), Foam::fvc::grad(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::k(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::k_, kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::kOmegaSSTBase(), lengthScaleRAS(), Foam::mag(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::omega(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::omega_, Foam::pos(), tmp< T >::ref(), S2(), Foam::sqr(), and useSigma_.

Referenced by correctNut(), GbyNu(), GbyNu0(), and S2().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ dTilda()

template<class BasicTurbulenceModel>

tmp< volScalarField > dTilda	(	const volScalarField &	magGradU,
		const volScalarField &	CDES ) const

protectedvirtual

Return length scale.

Reimplemented in kOmegaSSTDDES< BasicTurbulenceModel >, kOmegaSSTDDES< EddyDiffusivity< fluidThermoCompressibleTurbulenceModel > >, kOmegaSSTDDES< Foam::EddyDiffusivity< Foam::fluidThermoCompressibleTurbulenceModel > >, kOmegaSSTIDDES< BasicTurbulenceModel >, kOmegaSSTIDDES< EddyDiffusivity< fluidThermoCompressibleTurbulenceModel > >, and kOmegaSSTIDDES< Foam::EddyDiffusivity< Foam::fluidThermoCompressibleTurbulenceModel > >.

Definition at line 108 of file kOmegaSSTDES.C.

References CDES(), lengthScaleLES(), lengthScaleRAS(), and Foam::min().

Referenced by epsilonByk(), LESRegion(), and S2().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ epsilonByk()

template<class BasicTurbulenceModel>

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

protectedvirtual

Return epsilon/k.

Reimplemented from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >.

Definition at line 119 of file kOmegaSSTDES.C.

References CDES(), dTilda(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::F1(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::k_, Foam::mag(), and Foam::sqrt().

Here is the call graph for this function:

◆ GbyNu0()

template<class BasicTurbulenceModel>

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

protectedvirtual

Return (G/nu)_0.

Reimplemented from kOmegaSSTBase< DESModel< BasicTurbulenceModel > >.

Reimplemented in kOmegaSSTDDES< BasicTurbulenceModel >, kOmegaSSTDDES< EddyDiffusivity< fluidThermoCompressibleTurbulenceModel > >, and kOmegaSSTDDES< Foam::EddyDiffusivity< Foam::fluidThermoCompressibleTurbulenceModel > >.

Definition at line 131 of file kOmegaSSTDES.C.

References GbyNu0(), kOmegaSSTBase< DESModel< BasicTurbulenceModel > >::kOmegaSSTBase(), S2(), and useSigma_.

Referenced by GbyNu(), and GbyNu0().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GbyNu()

template<class BasicTurbulenceModel>

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