The k-epsilon-phit-f turbulence closure model for incompressible and compressible flows. More...

Inheritance diagram for kEpsilonPhitF< BasicTurbulenceModel >:

Collaboration diagram for kEpsilonPhitF< BasicTurbulenceModel >:

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 ("kEpsilonPhitF")
	Runtime type information. More...

	kEpsilonPhitF (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	~kEpsilonPhitF ()=default
	Destructor. More...

virtual bool	read ()
	Re-read model coefficients if they have changed. More...

tmp< volScalarField >	DkEff () const
	Return the effective diffusivity for k (LUU:Eq. 3) More...

tmp< volScalarField >	DepsilonEff () const
	Return the effective diffusivity for epsilon (LUU:Eq. 4) More...

tmp< volScalarField >	DphitEff () const
	Return the effective diffusivity for phit (LUU:Eq. 17) More...

virtual tmp< volScalarField >	k () const
	Return the turbulent kinetic energy field. More...

virtual tmp< volScalarField >	epsilon () const
	Return the turbulent kinetic energy dissipation rate field. More...

virtual tmp< volScalarField >	phit () const
	Return the normalised wall-normal fluctuating velocity scale field. More...

virtual tmp< volScalarField >	f () const
	Return the elliptic relaxation factor field. More...

virtual void	correct ()
	Solve the transport equations and correct the turbulent 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< 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< volScalarField >	k () const=0
	Return the turbulence kinetic energy. More...

virtual tmp< volSymmTensorField >	R () 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< 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...

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

Protected Member Functions
virtual void	correctNut ()
	Update nut with the latest available k, phit, and T. More...

tmp< volScalarField >	Ts () const
	Return the turbulent time scale, T. More...

tmp< volScalarField >	Ls () const
	Return the turbulent length scale, L. More...

virtual void	correctNut ()=0

Protected Attributes
Switch	includeNu_

dimensionedScalar	Cmu_

dimensionedScalar	Ceps1a_

dimensionedScalar	Ceps1b_

dimensionedScalar	Ceps1c_

dimensionedScalar	Ceps2_

dimensionedScalar	Cf1_

dimensionedScalar	Cf2_

dimensionedScalar	CL_

dimensionedScalar	Ceta_

dimensionedScalar	CT_

dimensionedScalar	sigmaK_

dimensionedScalar	sigmaEps_

dimensionedScalar	sigmaPhit_

volScalarField	k_
	Turbulent kinetic energy [m2/s2]. More...

volScalarField	epsilon_
	Turbulent kinetic energy dissipation rate [m2/s3]. More...

volScalarField	phit_
	Normalised wall-normal fluctuating velocity scale [-]. More...

volScalarField	f_
	Elliptic relaxation factor [1/s]. More...

volScalarField	T_
	Turbulent time scale [s]. More...

dimensionedScalar	phitMin_

dimensionedScalar	fMin_

dimensionedScalar	TMin_

dimensionedScalar	L2Min_

Protected Attributes inherited from eddyViscosity< RASModel< BasicTurbulenceModel > >
volScalarField	nut_

Detailed Description

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

The k-epsilon-phit-f turbulence closure model for incompressible and compressible flows.

The model is a three-transport-equation linear-eddy-viscosity turbulence closure model alongside an elliptic relaxation equation.

Input fields

k	:	Turbulent kinetic energy [m2/s2]
epsilon	:	Turbulent kinetic energy dissipation rate [m2/s3]
phit	:	Normalised wall-normal fluctuating velocity scale [-]
f	:	Elliptic relaxation factor [1/s]

Reference:

        Standard model (Tag:LUU):
            Laurence, D. R., Uribe, J. C., & Utyuzhnikov, S. V. (2005).
            A robust formulation of the v2-f model.
            Flow, Turbulence and Combustion, 73(3-4), 169–185.
            DOI:10.1007/s10494-005-1974-8

The default model coefficients are (LUU:Eqs. 19-20):

        kEpsilonPhitFCoeffs
        {
            includeNu   true;    // include nu in (LUU: Eq. 17), see Notes
            Cmu         0.22;    // Turbulent viscosity constant
            Ceps1a      1.4;     // Model constant for epsilon
            Ceps1b      1.0;     // Model constant for epsilon
            Ceps1c      0.05;    // Model constant for epsilon
            Ceps2       1.9;     // Model constant for epsilon
            Cf1         1.4;     // Model constant for f
            Cf2         0.3;     // Model constant for f
            CL          0.25;    // Model constant for L
            Ceta        110.0;   // Model constant for L
            CT          6.0;     // Model constant for T
            sigmaK      1.0;     // Turbulent Prandtl number for k
            sigmaEps    1.3;     // Turbulent Prandtl number for epsilon
            sigmaPhit   1.0;     // Turbulent Prandtl number for phit = sigmaK
        }

Note: The name of the original variable replacing 'v2' is 'phi' (LUU:Eq. 14). However, the name 'phi' preexisted in OpenFOAM; therefore, this name was replaced by 'phit' herein.

Including nu in DphitEff even though it is not present in (LUU:Eq. 17) provided higher level of resemblance to benchmarks for the tests considered, particularly for the peak skin friction (yet, pressure-related predictions were unaffected). Users can switch off nu in DphitEff by using includeNu entry in kEpsilonPhitFCoeffs as shown above in order to follow the reference paper thereat. includeNu is left true by default. See GitLab issue #1560.

Source files

See also: kEpsilon.C

Definition at line 131 of file kEpsilonPhitF.H.

Member Typedef Documentation

◆ alphaField

typedef BasicTurbulenceModel::alphaField alphaField

Definition at line 207 of file kEpsilonPhitF.H.

◆ rhoField

typedef BasicTurbulenceModel::rhoField rhoField

Definition at line 208 of file kEpsilonPhitF.H.

◆ transportModel

typedef BasicTurbulenceModel::transportModel transportModel

Definition at line 209 of file kEpsilonPhitF.H.

Constructor & Destructor Documentation

◆ kEpsilonPhitF()

kEpsilonPhitF	(	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`
	)