Realizable k-ε Model

The standard k-ε turbulence model and RNG k-ε turbulence model do not satisfy mathematical constraints on the Reynolds stresses for the consistency with physics of turbulence.

These constraints include the positivity of normal stresses ( u α ' u α ' ¯ 0 ) and Schwarz's inequality ( u α ' u α ' ¯     u β ' u β ' ¯ u α ' u β ' ¯   u α ' u β ' ¯ ). The realizable k-ε model proposed by Shih et al. (1995) employs a formulation of turbulent viscosity with a variation of C μ to satisfy these constraints. This is why this model is referred to as realizable. The second improvement made with the realizable model is to have a source term in the dissipation rate transport equation, which is derived from the transport equation of the mean square vorticity fluctuation. Compared to the standard k-ε turbulence model, it performs better when simulating planar jet flows, round jet flows and flows under adverse pressure gradients. Although the realizable model has significant performance improvement over the standard model, it has a problem in predicting impinging flow, swirling flow and secondary flows in a square duct since this model is based on the Boussinesq assumption. In addition, this model needs additional functions to simulate the near wall effects.

Transport Equations

Turbulent Kinetic Energy k(1)
( ρ k ) t + ( ρ u j ¯ k ) x j   = x j [ ( μ + μ t σ k ) k x j ] + P k + D k
Turbulent Dissipation Rate ε(2)
( ρ ε ) t + ( ρ u j ¯ ε ) x j   = x j [ ( μ + μ t σ ε ) ε x j ] + P ε + D ε + ρ C 1 S ε

Production Modeling

Turbulent Kinetic Energy k(3)
P k = μ t S 2
Turbulent Dissipation Rate ε (4)
P ε = C ε 1 ε k P k

Dissipation Modeling

Turbulent Kinetic Energy k (5)
D k = ρ ε
Turbulent Dissipation Rate ε (6)
D ε = ρ C ε 2 ε 2 k + ν ε

Modeling of Turbulent Viscosity μ t

(7)
μ t = C μ k 2 ε

Model Coefficients

C ε 1 = 1.44, C ε 2 = 1.9, σ k = 1.0, σ ε = 1.22.

where U * = S i j S i j + Ω i j ˜ Ω i j ˜ , Ω i j ˜ = Ω i j 2 ε i j k ω k , Ω i j = Ω i j ¯ 2 ε i j k ω k , A 0 = 4.04 , A s = 6   c o s ϕ , ϕ = 1 3 c o s 1 ( 6   W ) , W = S i j S j k S k j S i j S i j ,

C 1 = m a x [ 0.43 , η η + 5 ]

where η = S k ε , S = 2 S i j S i j is the strain rate magnitude.