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

Inheritance diagram for mixtureKEpsilon< BasicMomentumTransportModel >:

Collaboration diagram for mixtureKEpsilon< BasicMomentumTransportModel >:

[legend]

Public Types
typedef BasicMomentumTransportModel::alphaField	alphaField

typedef BasicMomentumTransportModel::rhoField	rhoField

typedef BasicMomentumTransportModel::transportModel	transportModel

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

typedef RASModel< BasicMomentumTransportModel > ::rhoField	rhoField

typedef RASModel< BasicMomentumTransportModel > ::transportModel	transportModel

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

typedef RASModel< BasicMomentumTransportModel > ::rhoField	rhoField

typedef RASModel< BasicMomentumTransportModel > ::transportModel	transportModel

Public Types inherited from RASModel< BasicMomentumTransportModel >
typedef BasicMomentumTransportModel::alphaField	alphaField

typedef BasicMomentumTransportModel::rhoField	rhoField

typedef BasicMomentumTransportModel::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 &type=typeName)
	Construct from components. More...

	mixtureKEpsilon (const mixtureKEpsilon &)=delete
	Disallow default bitwise copy construction. More...

virtual	~mixtureKEpsilon ()
	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...

void	operator= (const mixtureKEpsilon &)=delete
	Disallow default bitwise assignment. More...

Public Member Functions inherited from eddyViscosity< RASModel< BasicMomentumTransportModel > >
	eddyViscosity (const word &modelName, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport)
	Construct from components. More...

virtual	~eddyViscosity ()
	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 >	sigma () const
	Return the Reynolds stress tensor [m^2/s^2]. More...

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

Public Member Functions inherited from linearViscousStress< RASModel< BasicMomentumTransportModel > >
	linearViscousStress (const word &modelName, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport)
	Construct from components. More...

virtual	~linearViscousStress ()
	Destructor. More...

virtual tmp< volSymmTensorField >	devTau () const
	Return the effective stress tensor. More...

virtual tmp< fvVectorMatrix >	divDevTau (volVectorField &U) const
	Return the source term for the momentum equation. More...

virtual tmp< fvVectorMatrix >	divDevTau (const volScalarField &rho, volVectorField &U) const
	Return the source term for the momentum equation. More...

Public Member Functions inherited from RASModel< BasicMomentumTransportModel >
	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),(alpha, rho, U, alphaRhoPhi, phi, transport))

	RASModel (const word &type, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport)
	Construct from components. More...

	RASModel (const RASModel &)=delete
	Disallow default bitwise copy construction. More...

virtual	~RASModel ()
	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...

void	operator= (const RASModel &)=delete
	Disallow default bitwise assignment. 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< BasicMomentumTransportModel >
virtual void	printCoeffs (const word &type)
	Print model coefficients. More...

Protected Attributes
dimensionedScalar	Cmu_

dimensionedScalar	C1_

dimensionedScalar	C2_

dimensionedScalar	C3_

dimensionedScalar	Cp_

dimensionedScalar	alphap_

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< BasicMomentumTransportModel > >
volScalarField	nut_

Protected Attributes inherited from RASModel< BasicMomentumTransportModel >
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< BasicMomentumTransportModel >
static autoPtr< RASModel >	New (const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport)
	Return a reference to the selected RAS model. More...

Detailed Description

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

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;
        Cp          0.25; // Bubble generated turbulence
        alphap      0.3;  // Gas phase fraction below which
                          // bubble generated turbulence is included
        sigmak      1.0;
        sigmaEps    1.3;
    }

Source files

Definition at line 83 of file mixtureKEpsilon.H.

Member Typedef Documentation

◆ alphaField

typedef BasicMomentumTransportModel::alphaField alphaField

Definition at line 191 of file mixtureKEpsilon.H.

◆ rhoField

typedef BasicMomentumTransportModel::rhoField rhoField

Definition at line 192 of file mixtureKEpsilon.H.

◆ transportModel

typedef BasicMomentumTransportModel::transportModel transportModel

Definition at line 193 of file mixtureKEpsilon.H.

Constructor & Destructor Documentation

◆ mixtureKEpsilon() [1/2]

mixtureKEpsilon	(	const alphaField &	alpha,
		const rhoField &	rho,
		const volVectorField &	U,
		const surfaceScalarField &	alphaRhoPhi,
		const surfaceScalarField &	phi,
		const transportModel &	transport,
		const word &	type = `typeName`
	)

Construct from components.

Definition at line 47 of file mixtureKEpsilon.C.

References Foam::bound(), and mixtureKEpsilon< BasicMomentumTransportModel >::epsilonBoundaryTypes().

Here is the call graph for this function:

◆ mixtureKEpsilon() [2/2]

mixtureKEpsilon ( const mixtureKEpsilon< BasicMomentumTransportModel > & )

delete

Disallow default bitwise copy construction.

◆ ~mixtureKEpsilon()

virtual ~mixtureKEpsilon ( )

inlinevirtual

Destructor.

Definition at line 219 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicMomentumTransportModel >::read().

Here is the call graph for this function:

Member Function Documentation

◆ epsilonBoundaryTypes()

wordList epsilonBoundaryTypes ( const volScalarField & epsilon ) const

protected

Definition at line 180 of file mixtureKEpsilon.C.

References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), mixtureKEpsilon< BasicMomentumTransportModel >::correctInletOutlet(), forAll, and patchi.

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::mixtureKEpsilon().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ correctInletOutlet()

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

protected

Definition at line 202 of file mixtureKEpsilon.C.

References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), forAll, and patchi.

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::epsilonBoundaryTypes().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ initMixtureFields()

void initMixtureFields ( )

protected

Definition at line 229 of file mixtureKEpsilon.C.

References IOobject::AUTO_WRITE, GeometricField< Type, PatchField, GeoMesh >::boundaryField(), mixtureKEpsilon< BasicMomentumTransportModel >::epsilon_, mixtureKEpsilon< BasicMomentumTransportModel >::k_, and IOobject::READ_IF_PRESENT.

Here is the call graph for this function:

◆ correctNut()

void correctNut ( )

protectedvirtual

Implements eddyViscosity< RASModel< BasicMomentumTransportModel > >.

Definition at line 343 of file mixtureKEpsilon.C.

References IOobject::db(), fluid, phaseModel::fluid(), IOobject::groupName(), objectRegistry::lookupObject(), phaseModel::name(), dictionary::New(), and Foam::sqr().

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::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 385 of file mixtureKEpsilon.C.

References drag, mixtureKEpsilon< BasicMomentumTransportModel >::epsilon_, Foam::exp(), fluid, phaseModel::fluid(), dragModel::K(), mixtureKEpsilon< BasicMomentumTransportModel >::k_, Foam::mag(), phaseModel::rho(), Foam::sqr(), Foam::sqrt(), and U.

Here is the call graph for this function:

◆ rholEff()

tmp< volScalarField > rholEff ( ) const

protected

Definition at line 415 of file mixtureKEpsilon.C.

References fluid, phaseModel::fluid(), phaseSystem::otherPhase(), and phaseModel::rho().

Here is the call graph for this function:

◆ rhogEff()

tmp< volScalarField > rhogEff ( ) const

protected

Definition at line 425 of file mixtureKEpsilon.C.

References virtualMassModel::Cvm(), fluid, and phaseModel::fluid().

Here is the call graph for this function:

◆ rhom()

tmp< volScalarField > rhom ( ) const

protected

Definition at line 436 of file mixtureKEpsilon.C.

References alphal, and mixtureKEpsilon< BasicMomentumTransportModel >::mix().

Here is the call graph for this function:

◆ mix()

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

protected

Definition at line 447 of file mixtureKEpsilon.C.

References alphal, and mixtureKEpsilon< BasicMomentumTransportModel >::mixU().

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::rhom().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ mixU()

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

protected

Definition at line 461 of file mixtureKEpsilon.C.

References alphal, and mixtureKEpsilon< BasicMomentumTransportModel >::mixFlux().

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::mix().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ mixFlux()

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

protected

Definition at line 477 of file mixtureKEpsilon.C.

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

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::mixU().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ bubbleG()

tmp< volScalarField > bubbleG ( ) const

protected

Definition at line 499 of file mixtureKEpsilon.C.

References dragModel::CdRe(), drag, fluid, phaseModel::fluid(), Foam::mag(), fluidThermo::nu(), Foam::pos(), Foam::pow(), Foam::pow3(), phaseModel::rho(), phaseModel::thermo(), and U.

Here is the call graph for this function:

◆ kSource()

tmp< fvScalarMatrix > kSource ( ) const

protectedvirtual

Definition at line 538 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 546 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 170 of file mixtureKEpsilon.H.

References GeometricField< scalar, fvPatchField, volMesh >::New().

Here is the call graph for this function:

◆ DepsilonEff()

tmp<volScalarField> DepsilonEff ( const volScalarField & nutm ) const

inlineprotected

Return the effective diffusivity for epsilon.

Definition at line 180 of file mixtureKEpsilon.H.

References GeometricField< scalar, fvPatchField, volMesh >::New().

Here is the call graph for this function:

◆ TypeName()

TypeName ( "mixtureKEpsilon< BasicMomentumTransportModel >" )

Runtime type information.

◆ read()

bool read ( )

virtual

Re-read model coefficients if they have changed.

Implements eddyViscosity< RASModel< BasicMomentumTransportModel > >.

Definition at line 321 of file mixtureKEpsilon.C.

References Foam::read().

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::~mixtureKEpsilon().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ k()

virtual tmp<volScalarField> k ( ) const

inlinevirtual

Return the turbulence kinetic energy.

Implements eddyViscosity< RASModel< BasicMomentumTransportModel > >.

Definition at line 229 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicMomentumTransportModel >::k_.

◆ epsilon()

virtual tmp<volScalarField> epsilon ( ) const

inlinevirtual

Return the turbulence kinetic energy dissipation rate.

Definition at line 235 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicMomentumTransportModel >::correct(), mixtureKEpsilon< BasicMomentumTransportModel >::epsilon_, and mixtureKEpsilon< BasicMomentumTransportModel >::operator=().

Here is the call graph for this function:

◆ correct()

void correct ( )

virtual

Solve the turbulence equations and correct the turbulence viscosity.

Implements eddyViscosity< RASModel< BasicMomentumTransportModel > >.

Definition at line 553 of file mixtureKEpsilon.C.

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

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::epsilon().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ operator=()

void operator= ( const mixtureKEpsilon< BasicMomentumTransportModel > & )

delete

Disallow default bitwise assignment.

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::epsilon().

Here is the caller graph for this function:

Member Data Documentation

◆ Cmu_

dimensionedScalar Cmu_

protected

Definition at line 105 of file mixtureKEpsilon.H.

◆ C1_

dimensionedScalar C1_

protected

Definition at line 106 of file mixtureKEpsilon.H.

◆ C2_

dimensionedScalar C2_

protected

Definition at line 107 of file mixtureKEpsilon.H.

◆ C3_

dimensionedScalar C3_

protected

Definition at line 108 of file mixtureKEpsilon.H.

◆ Cp_

dimensionedScalar Cp_

protected

Definition at line 109 of file mixtureKEpsilon.H.

◆ alphap_

dimensionedScalar alphap_

protected

Definition at line 110 of file mixtureKEpsilon.H.

◆ sigmak_

dimensionedScalar sigmak_

protected

Definition at line 111 of file mixtureKEpsilon.H.

◆ sigmaEps_

dimensionedScalar sigmaEps_

protected

Definition at line 112 of file mixtureKEpsilon.H.

◆ k_

volScalarField k_

protected

Definition at line 116 of file mixtureKEpsilon.H.

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

◆ epsilon_

volScalarField epsilon_

protected

Definition at line 117 of file mixtureKEpsilon.H.

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

◆ Ct2_

autoPtr<volScalarField> Ct2_

protected

Definition at line 121 of file mixtureKEpsilon.H.

◆ rhom_

autoPtr<volScalarField> rhom_

protected

Definition at line 122 of file mixtureKEpsilon.H.

◆ km_

autoPtr<volScalarField> km_

protected

Definition at line 123 of file mixtureKEpsilon.H.

◆ epsilonm_

autoPtr<volScalarField> epsilonm_

protected

Definition at line 124 of file mixtureKEpsilon.H.

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

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

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

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

Member Typedef Documentation

◆ alphaField

◆ rhoField

◆ transportModel

Constructor & Destructor Documentation

◆ mixtureKEpsilon() [1/2]

◆ mixtureKEpsilon() [2/2]

◆ ~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()

◆ operator=()

Member Data Documentation

◆ Cmu_

◆ C1_

◆ C2_

◆ C3_

◆ Cp_

◆ alphap_

◆ sigmak_

◆ sigmaEps_

◆ k_

◆ epsilon_

◆ Ct2_

◆ rhom_

◆ km_

◆ epsilonm_

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