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

Inheritance diagram for mixtureKEpsilon< BasicTurbulenceModel >:

Collaboration diagram for mixtureKEpsilon< BasicTurbulenceModel >:

Public Types
typedef BasicTurbulenceModel::alphaField	alphaField

typedef BasicTurbulenceModel::rhoField	rhoField

typedef BasicTurbulenceModel::transportModel	transportModel

Public Types inherited from eddyViscosity< RASModel< BasicTurbulenceModel > >
typedef RASModel< BasicTurbulenceModel > ::alphaField	alphaField

typedef RASModel< BasicTurbulenceModel > ::rhoField	rhoField

typedef RASModel< BasicTurbulenceModel > ::transportModel	transportModel

Public Types inherited from linearViscousStress< RASModel< BasicTurbulenceModel > >
typedef RASModel< BasicTurbulenceModel > ::alphaField	alphaField

typedef RASModel< BasicTurbulenceModel > ::rhoField	rhoField

typedef RASModel< BasicTurbulenceModel > ::transportModel	transportModel

Public Types inherited from RASModel< 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< volScalarField >	k () const
	Return the turbulence kinetic energy. More...

virtual tmp< volScalarField >	epsilon () 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 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< volSymmTensorField >	R () const
	Return the Reynolds stress tensor. More...

virtual void	validate ()
	Validate the turbulence fields after construction. More...

Public Member Functions inherited from linearViscousStress< RASModel< 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 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...

Public Member Functions inherited from RASModel< BasicTurbulenceModel >
	TypeName ("RAS")
	Runtime type information. More...

	declareRunTimeSelectionTable (autoPtr, RASModel, dictionary,(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName),(alpha, rho, U, alphaRhoPhi, phi, transport, propertiesName))

	RASModel (const word &type, 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	~RASModel ()=default
	Destructor. More...

const dimensionedScalar &	kMin () const
	Return the lower allowable limit for k (default: SMALL) More...

const dimensionedScalar &	epsilonMin () const
	Return the lower allowable limit for epsilon (default: SMALL) More...

const dimensionedScalar &	omegaMin () const
	Return the lower allowable limit for omega (default: SMALL) More...

dimensionedScalar &	kMin ()
	Allow kMin to be changed. More...

dimensionedScalar &	epsilonMin ()
	Allow epsilonMin to be changed. More...

dimensionedScalar &	omegaMin ()
	Allow omegaMin to be changed. More...

virtual const dictionary &	coeffDict () const
	Const access to the coefficients dictionary. More...

virtual tmp< volScalarField >	nuEff () const
	Return the effective viscosity. More...

virtual tmp< scalarField >	nuEff (const label patchi) const
	Return the effective viscosity on patch. More...

virtual tmp< volScalarField >	omega () const
	Return the specific dissipation rate. More...

Protected Member Functions
wordList	epsilonBoundaryTypes (const volScalarField &epsilon) const

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

void	initMixtureFields ()

virtual void	correctNut ()

tmp< volScalarField >	Ct2 () const

tmp< volScalarField >	rholEff () const

tmp< volScalarField >	rhogEff () const

tmp< volScalarField >	rhom () const

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

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

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

tmp< volScalarField >	bubbleG () const

virtual tmp< fvScalarMatrix >	kSource () const

virtual tmp< fvScalarMatrix >	epsilonSource () const

tmp< volScalarField >	DkEff (const volScalarField &nutm) const
	Return the effective diffusivity for k. More...

tmp< volScalarField >	DepsilonEff (const volScalarField &nutm) const
	Return the effective diffusivity for epsilon. More...

Protected Member Functions inherited from RASModel< BasicTurbulenceModel >
virtual void	printCoeffs (const word &type)
	Print model coefficients. More...

	RASModel (const RASModel &)=delete
	No copy construct. More...

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

Protected Attributes
dimensionedScalar	Cmu_

dimensionedScalar	C1_

dimensionedScalar	C2_

dimensionedScalar	C3_

dimensionedScalar	Cp_

dimensionedScalar	sigmak_

dimensionedScalar	sigmaEps_

volScalarField	k_

volScalarField	epsilon_

autoPtr< volScalarField >	Ct2_

autoPtr< volScalarField >	rhom_

autoPtr< volScalarField >	km_

autoPtr< volScalarField >	epsilonm_

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

Protected Attributes inherited from RASModel< BasicTurbulenceModel >
dictionary	RASDict_
	RAS coefficients dictionary. More...

Switch	turbulence_
	Turbulence on/off flag. More...

Switch	printCoeffs_
	Flag to print the model coeffs at run-time. More...

dictionary	coeffDict_
	Model coefficients dictionary. More...

dimensionedScalar	kMin_
	Lower limit of k. More...

dimensionedScalar	epsilonMin_
	Lower limit of epsilon. More...

dimensionedScalar	omegaMin_
	Lower limit for omega. More...

Additional Inherited Members
Static Public Member Functions inherited from RASModel< BasicTurbulenceModel >
static autoPtr< RASModel >	New (const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName=turbulenceModel::propertiesName)
	Return a reference to the selected RAS model. More...

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

References Foam::bound(), 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 176 of file mixtureKEpsilon.C.

References epsilon, and forAll.

◆ correctInletOutlet()

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

protected

Definition at line 198 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_, and IOobject::READ_IF_PRESENT.

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 optionList::correct(), options::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 ( ) const

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(), Foam::sqrt(), and U.

Here is the call graph for this function:

◆ rholEff()

tmp< volScalarField > rholEff ( ) const

protected

Definition at line 405 of file mixtureKEpsilon.C.

References fluid.

◆ rhogEff()

tmp< volScalarField > rhogEff ( ) const

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 ( ) const

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

References alphal.

◆ mixU()

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

protected

Definition at line 452 of file mixtureKEpsilon.C.

References alphal.

◆ mixFlux()

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

protected

Definition at line 468 of file mixtureKEpsilon.C.

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

Here is the call graph for this function:

◆ bubbleG()

tmp< volScalarField > bubbleG ( ) const

protected

Definition at line 489 of file mixtureKEpsilon.C.

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

Here is the call graph for this function:

◆ kSource()

tmp< fvScalarMatrix > kSource ( ) const

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 ( ) const

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.

Reimplemented from RASModel< BasicTurbulenceModel >.

Definition at line 245 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicTurbulenceModel >::epsilon_.

◆ correct()

void correct ( )

virtual

Solve the turbulence equations and correct the turbulence viscosity.

Implements eddyViscosity< RASModel< BasicTurbulenceModel > >.

Definition at line 541 of file mixtureKEpsilon.C.

References Foam::fvc::absolute(), alphal, Foam::bound(), GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), tmp< T >::clear(), correct(), GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), mixtureKEpsilon< BasicTurbulenceModel >::correctNut(), Foam::fvm::ddt(), Foam::dev(), Foam::fvm::div(), Foam::fvc::div(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon_, fluid, fvOptions, Foam::fvc::grad(), mixtureKEpsilon< BasicTurbulenceModel >::k_, Foam::fvm::laplacian(), options::New(), nu, eddyViscosity< RASModel< BasicTurbulenceModel > >::nut_, phi, phig(), tmp< T >::ref(), Foam::solve(), Foam::fvm::Sp(), Foam::fvm::SuSp(), and Foam::twoSymm().

Here is the call graph for this function:

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

Definition at line 119 of file mixtureKEpsilon.H.

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

◆ k_

volScalarField k_

protected

Definition at line 123 of file mixtureKEpsilon.H.

Referenced by mixtureKEpsilon< BasicTurbulenceModel >::correct(), mixtureKEpsilon< BasicTurbulenceModel >::Ct2(), mixtureKEpsilon< BasicTurbulenceModel >::initMixtureFields(), and mixtureKEpsilon< BasicTurbulenceModel >::k().

◆ epsilon_

volScalarField epsilon_

protected

Definition at line 124 of file mixtureKEpsilon.H.

Referenced by mixtureKEpsilon< BasicTurbulenceModel >::correct(), mixtureKEpsilon< BasicTurbulenceModel >::Ct2(), mixtureKEpsilon< BasicTurbulenceModel >::epsilon(), and mixtureKEpsilon< BasicTurbulenceModel >::initMixtureFields().

◆ 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:

src/TurbulenceModels/phaseCompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.H
src/TurbulenceModels/phaseCompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.C

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

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

Member Typedef Documentation

◆ alphaField

◆ rhoField

◆ transportModel

Constructor & Destructor Documentation

◆ mixtureKEpsilon()

◆ ~mixtureKEpsilon()

Member Function Documentation

◆ epsilonBoundaryTypes()

◆ correctInletOutlet()

◆ initMixtureFields()

◆ correctNut()

◆ Ct2()

◆ rholEff()

◆ rhogEff()

◆ rhom()

◆ mix()

◆ mixU()

◆ mixFlux()

◆ bubbleG()

◆ kSource()

◆ epsilonSource()

◆ DkEff()

◆ DepsilonEff()

◆ TypeName()

◆ read()

◆ k()

◆ epsilon()

◆ correct()

Member Data Documentation

◆ Cmu_

◆ C1_

◆ C2_

◆ C3_

◆ Cp_

◆ sigmak_

◆ sigmaEps_

◆ k_

◆ epsilon_

◆ Ct2_

◆ rhom_

◆ km_

◆ epsilonm_

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