Low-Re k-epsilon model. More...
A single-phase, isothermal low-Reynolds k-epsilon model.
Single-phase Reynolds-Averaged Navier-Stokes flow.
A single-phase, isothermal Reynolds-Averaged Navier-Stokes model.
This model implements a single-phase, isothermal Reynolds-Averaged Navier-Stokes model, solving the momentum balance equation
\[\frac{\partial (\varrho \textbf{v})}{\partial t} + \nabla \cdot (\varrho \textbf{v} \textbf{v}^{\text{T}}) = \nabla \cdot (\mu_\textrm{eff} (\nabla \textbf{v} + \nabla \textbf{v}^{\text{T}})) - \nabla p + \varrho \textbf{g} - \textbf{f} \]
The effective viscosity is composed of the fluid and the eddy viscosity:
\[ \mu_\textrm{eff} = \mu + \mu_\textrm{t} \]
.
The low-Reynolds k-epsilon models calculate the eddy viscosity with two additional PDEs, one for the turbulent kinetic energy (k) and for the dissipation ( \( \varepsilon \)). The model uses the one proposed by Chien [Chien1982a]. A good overview and additional models are given in Patel et al. [Patel1985a].
The turbulent kinetic energy balance is identical with the one from the k-epsilon model, but the dissipation includes a dampening function ( \( D_\varepsilon \)): \( \varepsilon = \tilde{\varepsilon} + D_\varepsilon \):
\[ \frac{\partial \left( \varrho k \right)}{\partial t} + \nabla \cdot \left( \textbf{v} \varhho k \right) - \nabla \cdot \left( \left( \mu + \frac{\mu_\text{t}}{\sigma_\text{k}} \right) \nabla k \right) - 2 \mu_\text{t} \textbf{S} \cdot \textbf{S} + \varrho \tilde{\varepsilon} + D_\varepsilon \varrho = 0 \]
.
The dissipation balance is changed by introducing additional functions ( \( E_\text{k}\), \( f_1 \), and \( f_2 \)) to account for a dampening towards the wall:
\[ \frac{\partial \left( \varrho \tilde{\varepsilon} \right)}{\partial t} + \nabla \cdot \left( \textbf{v} \varrho \tilde{\varepsilon} \right) - \nabla \cdot \left( \left( \mu + \frac{\mu_\text{t}}{\sigma_{\varepsilon}} \right) \nabla \tilde{\varepsilon} \right) - C_{1\tilde{\varepsilon}} f_1 \frac{\tilde{\varepsilon}}{k} 2 \mu_\text{t} \textbf{S} \cdot \textbf{S} + C_{2\tilde{\varepsilon}} \varrho f_2 \frac{\tilde{\varepsilon}^2}{k} - E_\text{k} \varrho = 0 \]
.
The kinematic eddy viscosity \( \nu_\text{t} \) is dampened by \( f_\mu \):
\[\mu_\text{t} = \varrho C_\mu f_\mu \frac{k^2}{\tilde{\varepsilon}} \]
.
The auxiliary and dampening functions are defined as:
\[ D_\varepsilon = 2 \nu \frac{k}{y^2} \]
\[ E_\text{k} = -2 \nu \frac{\tilde{\varepsilon}}{y^2} \exp \left( -0.5 y^+ \right) \]
\[ f_1 = 1 \]
\[ f_2 = 1 - 0.22 \exp \left( - \left( \frac{\mathit{Re}_\text{t}}{6} \right)^2 \right) \]
\[ f_\mu = 1 - \exp \left( -0.0115 y^+ \right) \]
\[ \mathit{Re}_\text{t} = \frac{k^2}{\nu \tilde{\varepsilon}} \]
.
Finally, the model is closed with the following constants:
\[ \sigma_\text{k} = 1.00 \]
\[ \sigma_\varepsilon =1.30 \]
\[ C_{1\tilde{\varepsilon}} = 1.35 \]
\[ C_{2\tilde{\varepsilon}} = 1.80 \]
\[ C_\mu = 0.09 \]
| Files | |
| file | fluxvariables.hh | 
| The flux variables class for the low-Reynolds k-epsilon model. This is a convenience alias for that actual, discretization-specific flux variables. | |
| file | iofields.hh | 
| Adds I/O fields for the low-Re k-epsilon turbulence model. | |
| file | localresidual.hh | 
| The local residual class for the low-Reynolds k-epsilon model. This is a convenience alias for the actual, discretization-specific local residual. | |
| file | model.hh | 
| A single-phase, isothermal low-Reynolds k-epsilon model. | |
| file | problem.hh | 
| Low-Re k-epsilon turbulence problem base class. | |
| file | fluxvariables.hh | 
| The flux variables class for the low-Reynolds k-epsilon model using the staggered grid discretization. | |
| file | localresidual.hh | 
| Element-wise calculation of the residual for low-Reynolds k-epsilon models using the staggered discretization. | |
| file | volumevariables.hh | 
| Volume variables for the isothermal single-phase low-Re k-epsilons model. | |
| Classes | |
| struct | Dumux::LowReKEpsilonIOFields | 
| Adds I/O fields for the low-Re k-epsilon turbulence model.  More... | |
| struct | Dumux::Properties::LowReKEpsilonModelTraits< dimension > | 
| Traits for the low-Reynolds k-epsilon model.  More... | |
| class | Dumux::RANSProblemImpl< TypeTag, TurbulenceModel::lowrekepsilon > | 
| Low-Re k-epsilon turbulence problem base class.  More... | |
| class | Dumux::LowReKEpsilonFluxVariablesImpl< TypeTag, BaseFluxVariables, DiscretizationMethod > | 
| The flux variables class for the low-Reynolds k-epsilon model using the staggered grid discretization.  More... | |
| class | Dumux::LowReKEpsilonResidualImpl< TypeTag, BaseLocalResidual, DiscretizationMethod > | 
| Element-wise calculation of the residual for low-Reynolds k-epsilon models using the staggered discretization.  More... | |
| class | Dumux::LowReKEpsilonVolumeVariables< Traits, NSVolumeVariables > | 
| Volume variables for the isothermal single-phase low-Re k-epsilons model.  More... | |
| Typedefs | |
| template<class TypeTag, class BaseFluxVariables> | |
| using | Dumux::LowReKEpsilonFluxVariables = LowReKEpsilonFluxVariablesImpl<TypeTag, BaseFluxVariables, typename GetPropType<TypeTag, Properties::GridGeometry>::DiscretizationMethod> | 
| The flux variables class for the low-Reynolds k-epsilon model. This is a convenience alias for that actual, discretization-specific flux variables. | |
| template<class TypeTag, class BaseLocalResidual> | |
| using | Dumux::LowReKEpsilonResidual = LowReKEpsilonResidualImpl<TypeTag, BaseLocalResidual, typename GetPropType<TypeTag, Properties::GridGeometry>::DiscretizationMethod> | 
| The local residual class for the low-Reynolds k-epsilon model. This is a convenience alias for the actual, discretization-specific local residual. | |
| using Dumux::LowReKEpsilonFluxVariables = LowReKEpsilonFluxVariablesImpl<TypeTag, BaseFluxVariables, typename GetPropType<TypeTag, Properties::GridGeometry>::DiscretizationMethod> | 
| using Dumux::LowReKEpsilonResidual = LowReKEpsilonResidualImpl<TypeTag, BaseLocalResidual, typename GetPropType<TypeTag, Properties::GridGeometry>::DiscretizationMethod> |