Theory for Buoyancy-Induced Turbulence
When the flow is compressible or weakly compressible, an additional contribution to the production of turbulent kinetic energy may be added. This contribution can be expressed in terms of Favre-averaged fluctuations or conventional time-filtered fluctuations according to,
(3-210)
where gi is the gravitational acceleration and the last approximation is applicable for small Froude numbers. Applying a gradient-diffusion modeling approach, Equation 3-210 can be recast into,
(3-211)
where σT is the turbulent Schmidt number. The production term in the ε equation is derived from the production term in the k equation. Hence, the corresponding term in the ε equation becomes,
(3-212)
where Cε1 is the corresponding constant for every model (1.44 for the Realizable k-ε model). For the v2-f model ε/k is changed to τ-1. θ is the angle between u and g, accounting for the difference between buoyant vertical shear layers and buoyant horizontal shear layers:
(3-213)
Using the relation between k, ε, and ω, the buoyancy production term in the ω equation of the k-ω model can be derived as,
(3-214)
For the SST-model, the buoyancy production term is
(3-215)
Hence, the expressions for the various ε based and ω based turbulence models differ only in terms of constants and damping functions.
The v2-f turbulence model uses the gradient of the wall distance variable lw to compute the wall-normal direction, and relating the wall-normal turbulent fluctuations to the direction of gravity results in the following expression for buoyant production of ζ,
(3-216)
This modeling assumes that the buoyancy-induced turbulence kinetic energy is pumped into the component of k aligned with g.
When the buoyancy contribution is determined from a multiphysics node, density variations with respect to the variable in the coupled physics interface are considered. For example, in nonisothermal flow the coupled interface is a Heat Transfer in Fluids interface, and, Equation 3-211 is replaced by,
(3-217)