SST k-omega model

From CFD-Wiki

Revision as of 13:41, 22 March 2007 by Druflin (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Kinematic Eddy Viscosity

$\nu _T = {a_1 k \over \mbox{max}(a_1 \omega, \Omega F_2) }$

Turbulence Kinetic Energy

${{\partial k} \over {\partial t}} + U_j {{\partial k} \over {\partial x_j }} = P_k - \beta ^* k\omega + {\partial \over {\partial x_j }}\left[ {\left( {\nu + \sigma_{k1} \nu _T } \right){{\partial k} \over {\partial x_j }}} \right]$

Specific Dissipation Rate

${{\partial \omega } \over {\partial t}} + U_j {{\partial \omega } \over {\partial x_j }} = \alpha S^2 - \beta \omega ^2 + {\partial \over {\partial x_j }}\left[ {\left( {\nu + \sigma_{\omega} \nu _T } \right){{\partial \omega } \over {\partial x_j }}} \right] + 2( 1 - F_1 ) \sigma_{\omega 2} {1 \over \omega} {{\partial k } \over {\partial x_i}} {{\partial \omega } \over {\partial x_i}}$

Closure Coefficients and Auxilary Relations

$F_2=\mbox{tanh} \left[ \left[ \mbox{max} \left( { 2 \sqrt{k} \over \beta^* \omega y } , { 500 \nu \over y^2 \omega } \right) \right]^2 \right]$

$P_k=\mbox{min} \left(\tau _{ij} {{\partial U_i } \over {\partial x_j }} , 20\beta^* k \omega \right)$

$F_1=\mbox{tanh} \left\{ \left\{ \mbox{min} \left[ \mbox{max} \left( {\sqrt{k} \over \beta ^* \omega y}, {500 \nu \over y^2 \omega} \right) , {4 \sigma_{\omega 2} k \over CD_{k\omega} y^2} \right] \right\} ^4 \right\}$

$CD_{k\omega}=\mbox{max} \left( 2\rho\sigma_{\omega 2} {1 \over \omega} {{\partial k} \over {\partial x_i}} {{\partial \omega} \over {\partial x_i}}, 10 ^{-10} \right )$