The k-epsilon (k-ε) model for turbulence is the most common to simulate the mean flow characteristics for turbulent flow conditions. It is a two equation model which gives a general description of turbulence by means of two transport equations (PDEs), which accounts for the history effects like convection and diffusion of turbulent energy. The 2 transported variables are turbulent kinetic energy k, which determine the energy in turbulence, and turbulent dissipation ε, which determines the rate of dissipation of the turbulent kinetic energy.
The k-ε model is shown to be applicable for free-shear flows, such as the ones with relatively small pressure gradients [1], but might not be the best model for problems involving large adverse pressure gradients [2]. Hence this model might not be suitable for inlets and compressors.
The turbulent energy k is given by:
$$k=\frac { 3 }{ 2 } { \left( UI \right) }^{ 2 }$$
where \(U\) is the mean flow velocity and \(I\) is the turbulence intensity. The turbulence internsity gives the level of turbulence and can be defined as follows:
$$I \equiv \frac { u’ }{ U }$$
where \(u’\) is the root-mean-square of the turbulent velocity fluctuations and \(U\) is the mean velocity. The root-mean-square of the turbulent velocity fluctuations \(u’\) is given as:
$$u’ \equiv \sqrt { \frac { 1 }{ 3 } \left( { { u’ }_{ x } }^{ 2 } + { { u’ }_{ y } }^{ 2 } + { { u’ }_{ z } }^{ 2 } \right) } =\sqrt { \frac { 2 }{ 3 } k }$$
The mean velocity \(U\) can be calculated as follows:
$$U \equiv \sqrt { { { U }_{ x } }^{ 2 }+{ { U }_{ y } }^{ 2 }+{ { U }_{ z } }^{ 2 }}$$
By default the SimScale platform defines the default values of turbulence (like k, ε) depending on the domain of the problem. This is considered not just during the initial condition, but also for the specifications of various boundaries. If the problem requires distinct values of turbulence in a boundary condition then the Custom boundary condition can be used.
To realistically model a given problem, it is important to define the turbulence intensity at the inlets. Here are a few examples of common estimations of the incoming turbulence intensity:
The turbulent intensity at the core of a pipe for a fully developed pipe flow can be estimated as follows:
$$I=0.16 { { Re }_{ { d }_{ h } } }^{ -\frac { 1 }{ 8 } }$$
where \({ { Re }_{ { d }_{ h } } }\) is the reynolds number for a pipe of hydraulic diameter \({ { d }_{ h } }\). The turbulent dissipation rate can be calcaulted using the following formula:
$$\epsilon ={ { C }_{ \mu } }^{ \frac { 3 }{ 4 } }\frac { { k }^{ \frac { 3 }{ 2 } } }{ l }$$
where \({ { C }_{ \mu } }\) is the turbulence model constant which usually takes the value 0.009, \(k\) is the turbulent energy, \(l\) is the turbulent length scale. The turbulence length scale describes the size of large energy-containing eddies in a turbulent flow. For a fully developed pipe flow this can be given as follows:
$$l=0.07{ d }_{ h }$$
where \({ { d }_{ h } }\) is the hydraulic diameter.
References
