A single-phase, isothermal k-epsilon model. More...
#include <dumux/common/properties.hh>
#include <dumux/freeflow/properties.hh>
#include <dumux/freeflow/rans/model.hh>
#include <dumux/freeflow/rans/twoeq/indices.hh>
#include <dumux/freeflow/turbulencemodel.hh>
#include "problem.hh"
#include "fluxvariables.hh"
#include "localresidual.hh"
#include "volumevariables.hh"
#include "iofields.hh"
Go to the source code of this file.
A single-phase, isothermal k-epsilon model.
Single-phase Reynolds-Averaged Navier-Stokes flow For a detailed model decription see freeflow/rans/model.hh.
The 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 Launder and Sharma [38] https://doi.org/10.1016/0094-4548(74)90150-7.
The turbulent kinetic energy balance is:
\[ \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 \varepsilon = 0 \]
.
The dissipation balance is:
\[ \frac{\partial \left( \varrho \varepsilon \right)}{\partial t} + \nabla \cdot \left( \textbf{v} \varrho \varepsilon \right) - \nabla \cdot \left( \left( \mu + \frac{\mu_\text{t}}{\sigma_{\varepsilon}} \right) \nabla \varepsilon \right) - C_{1\varepsilon} \frac{\varepsilon}{k} 2 \mu_\text{t} \textbf{S} \cdot \textbf{S} + C_{2\varepsilon} \varrho \frac{\varepsilon^2}{k} = 0 \]
.
The dynamic eddy viscosity \( \mu_\text{t} \) is:
\[ \mu_\text{t} = \varrho C_\mu \frac{k^2}{\tilde{\varepsilon}} \]
.
Finally, the model is closed with the following constants:
\[ \sigma_\text{k} = 1.00 \]
\[ \sigma_\varepsilon =1.30 \]
\[ C_{1\varepsilon} = 1.44 \]
\[ C_{2\varepsilon} = 1.92 \]
\[ C_\mu = 0.09 \]
Namespaces | |
namespace | Dumux |
namespace | Dumux::Properties |
namespace | Dumux::Properties::TTag |
Type tag for numeric models. | |