Supriya Ghosh (Editor)

Realizable k ε Model

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

The Standard k−ɛ is a well-established model capable of resolving through the boundary layer. The second model is Realizable k−ɛ, an improvement over the standard k−ɛ model. It is a relatively recent development and differs from the standard k−ɛ model in two ways.The realizable k−ɛ model contains a new formulation for the turbulent viscosity and a new transport equation for the dissipation rate, ɛ, that is derived from an exact equation for the transport of the mean-square vorticity fluctuation. The term "realizable" means that the model satisfies certain mathematical constraints on the Reynolds stresses, consistent with the physics of turbulent flows. Neither the standard k-ɛ model nor the RNG k-ɛ model is realizable. It introduces a Variable Cμ instead of constant. 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.

Contents

Transport equations

∂ ∂ t ( ρ k ) + ∂ ∂ x j ( ρ k u j ) = ∂ ∂ x j [ ( μ + μ t σ k ) ∂ k ∂ x j ] + P k + P b − ρ ϵ − Y M + S k ∂ ∂ t ( ρ ϵ ) + ∂ ∂ x j ( ρ ϵ u j ) = ∂ ∂ x j [ ( μ + μ t σ ϵ ) ∂ ϵ ∂ x j ] + ρ C 1 S ϵ − ρ C 2 ϵ 2 k + ν ϵ + C 1 ϵ ϵ k C 3 ϵ P b + S ϵ

Where

C 1 = max [ 0.43 , η η + 5 ] , η = S k ϵ , S = 2 S i j S i j

In these equations, P k represents the generation of turbulence kinetic energy due to the mean velocity gradients, calculated in same manner as standard k-epsilon model. P b is the generation of turbulence kinetic energy due to buoyancy, calculated in same way as standard k-epsilon model.

Modelling turbulent viscosity

μ t = ρ C μ k 2 ϵ

where
C μ = 1 A 0 + A s k U ∗ ϵ
U ∗ ≡ S i j S i j + Ω ~ i j Ω ~ i j  ;
Ω ~ i j = Ω i j − 2 ϵ i j k ω k  ;
Ω i j = Ω i j ¯ − ϵ i j k ω k

where Ω i j ¯ is the mean rate-of-rotation tensor viewed in a rotating reference frame with the angular velocity ω k . The model constants A 0 and A s are given by:
A 0 = 4.04 , A s = 6 cos ⁡ ϕ

ϕ = 1 3 cos − 1 ⁡ ( 6 W ) , W = S i j S j k S k i S ~ 3 , S ~ = S i j S i j , S i j = 1 2 ( ∂ u j ∂ x i + ∂ u i ∂ x j )

Model constants

C 1 ϵ = 1.44 , C 2 = 1.9 , σ k = 1.0 , σ ϵ = 1.2

The realizable model is unable to satisfactorily predict the radial velocity; it is also the most computationally-expensive model.

References

Realizable k- ε Model Wikipedia
(Text) CC BY-SA