Renormalization Group (RNG) k-ε Model

The RNG k-ε turbulence model (Yakhot and Orszag, 1986) deduces the behavior of large scale eddies from that of the smaller ones by utilizing the scale similarity properties that are inherent in the energy cascade (Bradshaw, 1997).

This model employs a modified coefficient in the dissipation rate equation to account for the interaction between the turbulent dissipation and mean shear. It results in better prediction of flows containing high streamline curvature, flows over a backward facing step (Yakhot et al, 1992), and flows in an expansion duct than the standard k-ε turbulence model. However, its performance worsens when predicting flows in a contraction duct (Hanjalic, 2004).

Transport Equations

Turbulent Kinetic Energy k(1)
(ρk)t+(ρ¯ujk)xj =xj[(μ+μtσk)kxj]+Pk+Dk(ρk)t+(ρ¯¯¯¯ujk)xj =xj[(μ+μtσk)kxj]+Pk+Dk
Turbulent Dissipation Rate ε(2)
(ρε)t+(ρ¯ujε)xj =xj[(μ+μtσε)εxj]+Pε+Dε

Production Modeling

Turbulent Kinetic Energy k (3)
Pk=μtS2
Turbulent Dissipation Rate ε (4)
Pε=Cε1εkμtS2=Cε1εkPk

Dissipation Modeling

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

Modeling of Turbulent Viscosity μt

(7)
μt=Cμk2ε

Model Coefficients

Cε1 =1.44, Cε2 =1.92, Cμ =0.09, σk =1.0, σε =1.3, C'ε2=˜Cε2+Cμλ3(1λ/λ0)1+βλ3 , ˜Cε2 = 1.68, Cμ = 0.085, λ = kεS , λ0 = 4.38, β = 0.012, σk = 0.72, σε = 0.72.