Physics:K-epsilon turbulence model

Short description: Model used in computational fluid dynamics

K-epsilon (k-ε) turbulence model is one of the most common models used in computational fluid dynamics (CFD) to simulate mean flow characteristics for turbulent flow conditions. It is a two equation model that gives a general description of turbulence by means of two transport equations (partial differential equations, PDEs). The original impetus for the K-epsilon model was to improve the mixing-length model, as well as to find an alternative to algebraically prescribing turbulent length scales in moderate to high complexity flows.

• The first transported variable is the turbulent kinetic energy (k).
• The second transported variable is the rate of dissipation of turbulent kinetic energy (ε).

Principle

Unlike earlier turbulence models, k-ε model focuses on the mechanisms that affect the turbulent kinetic energy. The mixing length model lacks this kind of generality. The underlying assumption of this model is that the turbulent viscosity is isotropic, in other words, the ratio between Reynolds stress and mean rate of deformations is the same in all directions.

Standard k-ε turbulence model

The exact k-ε equations contain many unknown and unmeasurable terms. For a much more practical approach, the standard k-ε turbulence model (Launder and Spalding, 1974) is used which is based on our best understanding of the relevant processes, thus minimizing unknowns and presenting a set of equations which can be applied to a large number of turbulent applications.

For turbulent kinetic energy k

[math]\displaystyle{ \frac{\partial (\rho k)}{\partial t}+ \frac {\partial (\rho k u_i)}{\partial x_i}=\frac {\partial}{\partial x_j}\left[\frac {\mu_t}{\sigma_k}\frac {\partial k}{\partial x_j}\right]+2{\mu_t}{E_{ij}}{E_{ij}}-\rho \varepsilon }[/math]

For dissipation [math]\displaystyle{ \varepsilon }[/math]

[math]\displaystyle{ \frac{\partial (\rho \varepsilon)}{\partial t}+ \frac{\partial (\rho \varepsilon u_i)}{\partial x_i} = \frac {\partial}{\partial x_j}\left[\frac {\mu_t}{\sigma_\varepsilon}\frac {\partial \varepsilon}{\partial x_j}\right] + C_{1 \varepsilon} \frac{\varepsilon}{k} 2{\mu_t}{E_{ij}}{E_{ij}}- C_{2 \varepsilon } \rho \frac{\varepsilon ^2}{k} }[/math]

Rate of change of k or ε in time + Transport of k or ε by advection = Transport of k or ε by diffusion + Rate of production of k or ε - Rate of destruction of k or ε

where

[math]\displaystyle{ u_i }[/math] represents velocity component in corresponding direction

[math]\displaystyle{ E_{ij} }[/math] represents component of rate of deformation

[math]\displaystyle{ \mu_t }[/math] represents eddy viscosity

[math]\displaystyle{ \mu _t = \rho C _{\mu} \frac{k^2}{\varepsilon} }[/math]

The equations also consist of some adjustable constants [math]\displaystyle{ \sigma_k }[/math], [math]\displaystyle{ \sigma_{\varepsilon} }[/math] , [math]\displaystyle{ C_{1\varepsilon} }[/math] and [math]\displaystyle{ C_{2\varepsilon} }[/math]. The values of these constants have been arrived at by numerous iterations of data fitting for a wide range of turbulent flows. These are as follows:

[math]\displaystyle{ C_{\mu} = 0.09 }[/math]            [math]\displaystyle{ \sigma_k = 1.00 }[/math]            [math]\displaystyle{ \sigma_{\varepsilon} = 1.30 }[/math]            [math]\displaystyle{ C_{1\varepsilon} = 1.44 }[/math]            [math]\displaystyle{ C_{2\varepsilon} = 1.92 }[/math]

Applications

The k-ε model has been tailored specifically for planar shear layers and recirculating flows. This model is the most widely used and validated turbulence model with applications ranging from industrial to environmental flows, which explains its popularity. It is usually useful for free-shear layer flows with relatively small pressure gradients as well as in confined flows where the Reynolds shear stresses are most important. It can also be stated as the simplest turbulence model for which only initial and/or boundary conditions needs to be supplied.

However it is more expensive in terms of memory than the mixing length model as it requires two extra PDEs. This model would be an inappropriate choice for problems such as inlets and compressors as accuracy has been shown experimentally to be reduced for flows containing large adverse pressure gradients. The k-ε model also performs poorly in a variety of important cases such as unconfined flows, curved boundary layers, rotating flows and flows in non-circular ducts.

Other models

Realizable k-ε Model: An immediate benefit of the realizable k-ɛ model is that it provides improved predictions for the spreading rate of both planar and round jets. It also exhibits superior performance for flows involving rotation, boundary layers under strong adverse pressure gradients, separation, and recirculation. In virtually every measure of comparison, Realizable k-ɛ demonstrates a superior ability to capture the mean flow of the complex structures.

k-ω Model: used when there are wall effects present within the case.

Reynolds stress equation model: In case of complex turbulent flows, Reynolds stress models are able to provide better predictions. Such flows include turbulent flows with high degrees of anisotropy, significant streamline curvature, flow separation, zones of recirculation and influence of mean rotation effects.

Notes

• 'An Introduction to Computational Fluid Dynamics: The Finite Volume Method (2nd Edition)', H. Versteeg, W. Malalasekera; Pearson Education Limited; 2007; ISBN:0131274988
• 'Turbulence Modeling for CFD' 2nd Ed., Wilcox C. D.; DCW Industries; 1998; ISBN:0963605100
• 'An introduction to turbulence and its measurement', Bradshaw, P.; Pergamon Press; 1971; ISBN:0080166210