• # kEpsilon模型 (7.0)

\frac{\p\alpha \rho \varepsilon}{\p t} +\nabla\cdot \left( \alpha \rho \bfU \varepsilon \right)- \nabla\cdot\left(\alpha \rho D_\varepsilon\nabla \varepsilon \right) = C_1 \alpha \rho G \frac{\varepsilon}{k} - \left( \left( \frac{2}{3} C_1 - C_{3,RDT} \right) \alpha \rho \nabla\cdot \bfU \varepsilon \right) - \left( C_2 \alpha \rho \frac{\varepsilon}{k} \varepsilon \right) + S_\varepsilon + S_{\text{fvOptions}}

\frac{\p \alpha \rho k}{\p t} + \nabla\cdot \left( \alpha \rho \bfU k \right) - \nabla\cdot \left( \alpha \rho D_k \nabla k \right) = \alpha \rho G - \left( \frac{2}{3} \alpha \rho \nabla\cdot\bfU k \right) - \left( \alpha \rho \frac{\epsilon}{k} k \right) + S_k + S_{\text{fvOptions}}

G=\nu_t\left(\nabla\bfU+\nabla\bfU^{T}-\frac{1}{3}(\nabla\cdot\bfU)\bfI\right) : \nabla\bfU

\nu_t = C_{\mu} \frac{k^2}{\epsilon}

C_\mu=0.09,C_1=1.44,C_2=1.92,C_{3,RDT}=0,\sigma_k=1,\sigma_\varepsilon=1.3,D_k=\frac{\nu_t}{\sigma_k}+\nu,D_\varepsilon=\frac{\nu_t}{\sigma_\varepsilon}+\nu

Launder, B. E., & Spalding, D. B. (1972). Lectures in mathematical models of turbulence. Launder, B. E., & Spalding, D. B. (1974). The numerical computation of turbulent flows. Computer methods in applied mechanics and engineering, 3(2), 269-289. El Tahry, S. H. (1983). k-epsilon equation for compressible reciprocating engine flows. Journal of Energy, 7(4), 345-353.

# kOmega模型 (7.0)

\frac{\p \alpha \rho k}{\p t} + \nabla\cdot \left( \alpha \rho \bfU k \right) - \nabla\cdot \left( \alpha \rho D_k \nabla k \right) = \alpha \rho G - \frac{2}{3} \alpha \rho \nabla\cdot\bfU k - C_\mu\alpha \rho\omega k + S_k + S_{\text{fvOptions}}

\frac{\p\alpha \rho \omega}{\p t} +\nabla\cdot \left( \alpha \rho \bfU \omega \right)- \nabla\cdot\left(\alpha \rho D_\omega\nabla \omega \right) = \gamma \alpha \rho G \frac{\omega}{k} - \frac{2}{3} \gamma \alpha \rho \nabla\cdot \bfU \omega - \beta \alpha \rho \omega^2+ S_\omega + S_{\text{fvOptions}}

G=\nu_t\left(\nabla\bfU+\nabla\bfU^{T}-\frac{1}{3}(\nabla\cdot\bfU)\bfI\right) : \nabla\bfU

\nu_t = \frac{k}{\omega}

C_\mu=0.09,\beta=0.072,\alpha_k=0.5,\alpha_\omega=0.5,\gamma=0.52,D_k=\alpha_k \nu_t+\nu,D_\varepsilon=\alpha_\omega \nu_t+\nu

符号标识：

$\alpha$:相分数
$\rho$:表示密度
$\bfU$:速度
$k$:湍流动能
$\varepsilon$:湍流动能耗散率
$S_k$:湍流动能源项
$S_\varepsilon$:湍流动能耗散率源项

• 东岳前辈，东岳流体“CFD中的LES大涡模拟”这项文档中 http://www.dyfluid.com/docs/LES.html ，亚格子偏应力是否应该是: T-1/3(tr T)I？

• 多谢，已更新。
这个文章有点过时了，我要找时间重写一下

CFD中文网 2016 - 2020 | 京ICP备15017992号-2