Skip to main content

turbulentMixingLengthDissipationRateInlet

Description

The turbulentMixingLengthDissipationRateInlet is a boundary condition that calculates turbulent kinetic energy dissipation rate, i.e. epsilon, based on a specified mixing length. The patch values are calculated using:

ϵp=Cμ0.75k1.5L \epsilon_p = \frac{C_{\mu}^{0.75} k^{1.5}}{L}

where:

PropertyDescription
ϵp\epsilon_pPatch epsilon values [m^2/s^3]
CμC_\muEmpirical model constant retrieved from turbulence model
kkTurbulent kinetic energy [m^2/s^2]
LLMixing length scale [m]

Usage

The condition requires entries in both the boundary and field files.

Boundary file

<patchName>
{
type patch;
...
}

Field file

<patchName>
{
// Mandatory entries
type turbulentMixingLengthDissipationRateInlet;
mixingLength <scalar>;

// Optional entries
Cmu <scalar>;
k <word>;
phi <word>;

// Inherited entries
...
}

where:

PropertyDescriptionTypeRequiredDefault
typeType name: turbulentMixingLengthDissipationRateInletwordyes-
mixingLengthMixing length scale [m]scalaryes-
CmuEmpirical model constantscalarno0.09
phiName of flux fieldwordnophi
kName of turbulent kinetic energy fieldwordnok

The inherited entries are elaborated in:

  • inletOutletFvPatchFields.H
  • The boundary condition is derived from inletOutlet condition. Therefore, in the event of reverse flow, a zero-gradient condition is applied.
  • The order of precedence to input the empirical model constant Cmu is: turbulence model, boundary condition dictionary, and default value=0.09.
  • The empirical model constant Cmu is not a spatiotemporal variant field. Therefore, the use of the boundary condition may not be fully consistent with the turbulence models where Cmu is a variant field, such as realizableKE closure model in this respect. Nevertheless, workflow observations suggest that the matter poses no importance.

Further information

Tutorial:

Source code:

API:

History: Introduced in version 2.2.2