Supriya Ghosh (Editor)

K epsilon turbulence model

Updated on
Edit
Like
Comment
Share on FacebookTweet on TwitterShare on LinkedInShare on Reddit

K-epsilon (k-ε) turbulence model is the most common model used in Computational Fluid Dynamics (CFD) to simulate 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). 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.

Contents

  • The first transported variable determines the energy in the turbulence and is called turbulent kinetic energy (k).
  • The second transported variable is the turbulent dissipation (ε) which determines the rate of dissipation of the 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

    ( ρ k ) t + ( ρ k u i ) x i = x j [ μ t σ k k x j ] + 2 μ t E i j E i j ρ ϵ

    For dissipation ϵ

    ( ρ ϵ ) t + ( ρ ϵ u i ) x i = x j [ μ t σ ϵ ϵ x j ] + C 1 ϵ ϵ k 2 μ t E i j E i j C 2 ϵ ρ ϵ 2 k

    where

    u i represents velocity component in corresponding direction E i j represents component of rate of deformation μ t represents eddy viscosity μ t = ρ C μ k 2 ϵ

    The equations also consist of some adjustable constants σ k , σ ϵ , C 1 ϵ and C 2 ϵ . 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 μ = 0.09            σ k = 1.00            σ ϵ = 1.30            C 1 ϵ = 1.44            C 2 ϵ = 1.92

    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

    Realisable 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.

    References

    K-epsilon turbulence model Wikipedia