Homogeneous shear flow
If the turbulence is homogeneous, the equations for the TKE and dissipation can be reduced to
The closure coefficient Cϵ2 is already known. Therefore, Cϵ1 can be estimated from experimental data. According to experimental data [3], a self similar state develops with k ⁄ ϵ ~ const and ϵ ⁄ P ~ 0.68. It follows that
(10)
(d)/(dt)⎛⎝(k)/(ϵ)⎞⎠ = 0
= (1)/(ϵ)(dk)/(dt) − (k)/(ϵ2)(dϵ)/(dt)
= (1)/(ϵ)(P − ϵ) − (k)/(ϵ2)⎛⎝Cϵ1P(ϵ)/(k) − Cϵ2(ϵ2)/(k)⎞⎠
= (P)/(ϵ)(1 − Cϵ1) + Cϵ2 − 1
or
Keeping the value Cϵ2 = 1.77 determined in the previous section, it follows that Cϵ1 = 1.52 to satisfy ϵ ⁄ P = 0.68. The reader probably realized that this process of determining the closure coefficients is questionable. The turbulence model will only be calibrated for some given flows. If the flow changes, so should the closure coefficients. Turbulence models are typically calibrated for some particular applications. For example, the Spalart-Allmaras model has been tuned specifically for simple aerospace flows (boundary layer growth over a wing, etc). Typically, some compromises are made to ensure that the model performs reasonably well in a wide range of flows. The standard k − ϵ model has Cϵ1 = 1.44 and Cϵ2 = 1.92, which suggests n ~ 1.1 instead of 1.3 for decaying isotropic turbulence, and ϵ ⁄ P ~ 0.67 instead of 0.68 for homogeneous shear flows. These coefficients could be modified to capture more accurately these two flows. However, the model accuracy would suffer for some other situations.
References
[1] L. Mydlarski and Z. Warhaft. On the onset of high-reynolds-number grid-generated wind tunnel turbulence. Journal of Fluid Mechanics, 320:331--368, 1996.
[2] W. Rodi. Turbulence Models and Their Application in Hydraulics: A State-Of-The-Art Review. Taylor & Francis, 1993.
[3] S. Tavoularis and U. Karnik. Further experiments on the evolution of turbulent stresses and scales in uniformly sheared turbulence. Journal of Fluid Mechanics, 204(1):457--478, 1989.
← Previous Page
← Previous Page