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

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

typedef BasicMomentumTransportModel::rhoField	rhoField

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

typedef BasicMomentumTransportModel::rhoField	rhoField

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 viscosity &viscosity, 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 tmp< volScalarField >	omega () const
	Return the turbulence specific 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 viscosity &viscosity)
	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< BasicMomentumTransportModel >
	linearViscousStress (const word &modelName, const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const viscosity &viscosity)
	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...

Protected Member Functions
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 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_

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 183 of file mixtureKEpsilon.H.

◆ rhoField

typedef BasicMomentumTransportModel::rhoField rhoField

Definition at line 184 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 viscosity &	viscosity,
		const word &	type = `typeName`
	)

Construct from components.

Definition at line 48 of file mixtureKEpsilon.C.

References IOobject::AUTO_WRITE, Foam::bound(), mixtureKEpsilon< BasicMomentumTransportModel >::Ct2_, Foam::dimDensity, Foam::dimless, mixtureKEpsilon< BasicMomentumTransportModel >::epsilon_, mixtureKEpsilon< BasicMomentumTransportModel >::epsilonm_, phaseModel::index(), mixtureKEpsilon< BasicMomentumTransportModel >::k_, mixtureKEpsilon< BasicMomentumTransportModel >::km_, IOobject::MUST_READ, IOobject::READ_IF_PRESENT, mixtureKEpsilon< BasicMomentumTransportModel >::rhom_, and Foam::type().

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 210 of file mixtureKEpsilon.H.

Member Function Documentation

◆ initMixtureFields()

void initMixtureFields

protected

Definition at line 253 of file mixtureKEpsilon.C.

◆ correctNut()

void correctNut

protectedvirtual

Implements eddyViscosity< RASModel< BasicMomentumTransportModel > >.

Definition at line 291 of file mixtureKEpsilon.C.

References dictionary::New(), and Foam::sqr().

Here is the call graph for this function:

◆ Ct2()

tmp< volScalarField > Ct2

protected

Definition at line 334 of file mixtureKEpsilon.C.

References drag, Foam::exp(), phaseSystem::lookupInterfacialModel(), Foam::mag(), phaseSystem::phases(), phaseModel::rho(), liquid::rho(), Foam::sqr(), Foam::sqrt(), and U.

Here is the call graph for this function:

◆ rholEff()

tmp< volScalarField > rholEff

protected

Definition at line 365 of file mixtureKEpsilon.C.

References liquid::rho().

Here is the call graph for this function:

◆ rhogEff()

tmp< volScalarField > rhogEff

protected

Definition at line 374 of file mixtureKEpsilon.C.

References dispersedVirtualMassModel::Cvm(), phaseSystem::lookupInterfacialModel(), phaseSystem::phases(), phaseModel::rho(), and liquid::rho().

Here is the call graph for this function:

◆ rhom()

tmp< volScalarField > rhom

protected

Definition at line 390 of file mixtureKEpsilon.C.

◆ mix()

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

protected

Definition at line 400 of file mixtureKEpsilon.C.

◆ mixU()

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

protected

Definition at line 414 of file mixtureKEpsilon.C.

◆ mixFlux()

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

protected

Definition at line 430 of file mixtureKEpsilon.C.

References Foam::fvc::interpolate().

Here is the call graph for this function:

◆ bubbleG()

tmp< volScalarField > bubbleG

protected

Definition at line 453 of file mixtureKEpsilon.C.

References phaseModel::d(), drag, phaseSystem::lookupInterfacialModel(), Foam::mag(), phaseSystem::phases(), Foam::pos(), Foam::pow(), Foam::pow3(), liquid::rho(), and U.

Here is the call graph for this function:

◆ kSource()

tmp< fvScalarMatrix > kSource

protectedvirtual

Definition at line 494 of file mixtureKEpsilon.C.

References Foam::fvm::Su().

Here is the call graph for this function:

◆ epsilonSource()

tmp< fvScalarMatrix > epsilonSource

protectedvirtual

Definition at line 502 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 162 of file mixtureKEpsilon.H.

References GeometricField< Type, PatchField, GeoMesh >::New(), and mixtureKEpsilon< BasicMomentumTransportModel >::sigmak_.

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 172 of file mixtureKEpsilon.H.

References GeometricField< Type, PatchField, GeoMesh >::New(), and mixtureKEpsilon< BasicMomentumTransportModel >::sigmaEps_.

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 269 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< BasicMomentumTransportModel > >.

Definition at line 220 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 226 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicMomentumTransportModel >::epsilon_.

◆ omega()

virtual tmp<volScalarField> omega ( ) const

inlinevirtual

Return the turbulence specific dissipation rate.

Definition at line 232 of file mixtureKEpsilon.H.

References mixtureKEpsilon< BasicMomentumTransportModel >::Cmu_, mixtureKEpsilon< BasicMomentumTransportModel >::epsilon_, mixtureKEpsilon< BasicMomentumTransportModel >::k_, and GeometricField< Type, PatchField, GeoMesh >::New().

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

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

Here is the call graph for this function:

◆ operator=()

void operator= ( const mixtureKEpsilon< BasicMomentumTransportModel > & )

delete

Disallow default bitwise assignment.

Member Data Documentation

◆ Cmu_

dimensionedScalar Cmu_

protected

Definition at line 105 of file mixtureKEpsilon.H.

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

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

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

◆ sigmaEps_

dimensionedScalar sigmaEps_

protected

Definition at line 112 of file mixtureKEpsilon.H.

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

◆ k_

volScalarField k_

protected

Definition at line 116 of file mixtureKEpsilon.H.

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

◆ epsilon_

volScalarField epsilon_

protected

Definition at line 117 of file mixtureKEpsilon.H.

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

◆ Ct2_

autoPtr<volScalarField> Ct2_

protected

Definition at line 121 of file mixtureKEpsilon.H.

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

◆ rhom_

autoPtr<volScalarField> rhom_

protected

Definition at line 122 of file mixtureKEpsilon.H.

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

◆ km_

autoPtr<volScalarField> km_

protected

Definition at line 123 of file mixtureKEpsilon.H.

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

◆ epsilonm_

autoPtr<volScalarField> epsilonm_

protected

Definition at line 124 of file mixtureKEpsilon.H.

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

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

Detailed Description

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

Member Typedef Documentation

◆ alphaField

◆ rhoField

Constructor & Destructor Documentation

◆ mixtureKEpsilon() [1/2]

◆ mixtureKEpsilon() [2/2]

◆ ~mixtureKEpsilon()

Member Function Documentation

◆ initMixtureFields()

◆ correctNut()

◆ Ct2()

◆ rholEff()

◆ rhogEff()

◆ rhom()

◆ mix()

◆ mixU()

◆ mixFlux()

◆ bubbleG()

◆ kSource()

◆ epsilonSource()

◆ DkEff()

◆ DepsilonEff()

◆ TypeName()

◆ read()

◆ k()

◆ epsilon()

◆ omega()

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