k-epsilon turbulence model

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.{{citation |url=http://www.cfd-online.com/Wiki/K-epsilon_models| title=K-epsilon models}}

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

Principle

Unlike earlier turbulence models, k-ε model focuses on the mechanisms that affect the turbulence kinetic energy. The mixing length model lacks this kind of generality.{{cite book |url=https://books.google.com/books?id=RvBZ-UMpGzIC |title= An Introduction to Computational Fluid Dynamics: The Finite Volume Method|author=Henk Kaarle Versteeg, Weeratunge Malalasekera |publisher=Pearson Education Limited |year=2007 |isbn=9780131274983}} 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{{cite journal | last1 = Launder | first1 = B.E. | last2 = Spalding | first2 = D.B. | date = March 1974 | title = The numerical computation of turbulent flows | journal = Computer Methods in Applied Mechanics and Engineering | volume = 3 | issue = 2 | pages = 269–289 | doi = 10.1016/0045-7825(74)90029-2 | bibcode = 1974CMAME...3..269L }}) 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 turbulence kinetic energy k{{cite book | last1 = Versteeg | first1 = Henk Kaarle | last2 = Malalasekera | first2 = Weeratunge | title = An introduction to Computational Fluid Dynamics: The Finite Volume Method | date = 2007 | publisher = Pearson Education }}

: $\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$

For dissipation $\varepsilon$

: $\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}$

class="wikitable"

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

: $u_i$ represents velocity component in corresponding direction

: $E_{ij}$ represents component of rate of deformation

: $\mu_t$ represents eddy viscosity

: $\mu _t = \rho C _{\mu} \frac{k^2}{\varepsilon}$

The equations also consist of some adjustable constants $\sigma_k$ , $\sigma_{\varepsilon}$ , $C_{1\varepsilon}$ and $C_{2\varepsilon}$ . 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:

$C_{\mu} = 0.09$ {{spaces|10}} $\sigma_k = 1.00$ {{spaces|10}} $\sigma_{\varepsilon} = 1.30$ {{spaces|10}} $C_{1\varepsilon} = 1.44$ {{spaces|10}} $C_{2\varepsilon} = 1.92$

Applications

The k-ε model has been tailored specifically for planar shear layers[http://www.thermopedia.com/content/1118/?tid=104&sn=1159 usage of k-e to model shear layers] and recirculating flows.[http://samlab.epfl.ch/page-15502-en.html usage of k-e approach for modelling 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.{{citation |url=http://www.innovative-cfd.com/turbulence-model.html |title=The Turbulence Model Can Make a Big Difference in Your Results |access-date=2013-11-10 |archive-date=2019-08-07 |archive-url=https://web.archive.org/web/20190807073608/http://www.innovative-cfd.com/turbulence-model.html |url-status=dead }} 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{{citation needed|reason=Why?|date=August 2016}}. The k-ε model also performs poorly in a variety of important cases such as unconfined flows,{{citation |title=Turbulent Secondary Flows |author=P Bradshaw |journal=Annual Review of Fluid Mechanics |year=1987 |doi=10.1146/annurev.fl.19.010187.000413 |volume=19 |issue=1 |pages=53–74|bibcode=1987AnRFM..19...53B }} curved boundary layers, rotating flows and flows in non-circular ducts.{{citation |title=Secondary Flow in Semi Circular Ducts |journal=Journal of Fluids Engineering |year=2011 |doi=10.1115/1.4004991 |volume=133 |issue=10 |pages=101206–101214|url=https://repository.up.ac.za/bitstream/2263/42958/1/Larsson_Secondary_2014.pdf |last1=Larsson |first1=I. A. S. |last2=Lindmark |first2=E. M. |last3=Lundström |first3=T. S. |last4=Nathan |first4=G. J. |hdl=2263/42958 |hdl-access=free }}

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.Pope, Stephen. "Turbulent Flows". Cambridge University Press, 2000. Such flows include turbulent flows with high degrees of anisotropy, significant streamline curvature, flow separation, zones of recirculation and influence of mean rotation effects.

References

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}}

Category:Turbulence models