第章fluent湍流模型.docVIP

3.3.6 大涡模拟（LES） 湍流中包含了不同时间与长度尺度的涡旋。最大长度尺度通常为平均流动的特征长度尺度。最小尺度为Komogrov尺度。 LES的基本假设是1，动量、能量、质量及其它标量主要由大涡输运；2，流动的几何和边界条件决定了大涡的特性，而流动特性主要在大涡中体现；3，小尺度涡旋受几何和边界条件影响较小，并且各向同性；大涡模拟过程中，直接求解大涡，小尺度涡旋模拟，从而使得网格要求比DNS低。 3.3.6.1大涡模拟的控制方程 LES的控制方程是对Navier-Stokes方程在波数空间或者物理空间进行过滤得到的。过滤的过程是去掉比过滤宽度或者给定物理宽度小的涡旋，从而得到大涡旋的控制方程。 过滤变量（上横线）定义为： 3－55 其中，D表示流体区域；G是决定涡旋大小的过滤函数。 在FLUENT中，有限控制体离散本身暗中包括了过滤运算， ， 3－56 其中V是计算控制体体积，过滤函数为： 3－57 目前，大涡模拟对不可压流动问题得到较多应用，但在可压缩问题中的应用还很少，因此这里涉及的理论都是针对不可压流动的大涡模拟方法。在FLUENT中，大涡模拟只能针对不可压流体（当然并非说是密度是常数）的流动。 过滤不可压的Navier-Stokes方程后，可以得到LES控制方程： 3－58 3－59 其中，为亚网格应力，定义为： 3－60 很明显，上述方程与雷诺平均方程很相似，只不过大涡模拟中的变量是过滤过的量，而非时间平均量，并且湍流应力也不同。 3.3.6.2 亚网格模型 由于LES中亚网格应力项是未知的，并且需要模拟以封闭方程。目前，采用比较多的亚网格模型为涡旋粘性模型，形式为： 3－61 式中，是亚网格湍流粘性系数；是求解尺度下的应变率张量，定义为： 3－62 求解亚网格湍流粘性系数时，FLUENT提供了两种方法。第一，Smagorinsky-Lilly模型；第二，基于重整化群的亚网格模型。 最基本的亚网格模型是Smagorinsky [L145]最早提出的，Lilly [L99]把它进行了改善，这就是今天的Smagorinsky-Lilly模型。该模型的涡粘性计算方程为： 3－63 式中，是亚网格的混合长度；。是Smagorinsky常数，则亚网格混合长度可以用下式计算。 3－64 其中，k=0.42，d是到最近壁面的距离，V是计算控制体体积。 Lilly通过对均匀各向同性湍流惯性子区湍流分析，得到了＝0.23。但是研究中发现，对于有平均剪切或者过渡流动中，该系数过高估计了大尺度涡旋的阻尼作用。因此，对于比较多的流动问题，＝0.1有比较好的模拟结果，该值是FLUENT的默认设置值。 我们再来看看基于重整化群思想的亚网格模型。人们用重整化群理论推导出了亚网格涡旋粘性系数[L182]，该方法得到的是亚网格有效粘性系数，，而 3－65 式中，，H(x)是Heaviside函数， 3－66 V是计算控制体体积；重整化群常数,而常数C=100。 对于高雷诺数流动（ ），，基于重整化群理论的亚网格模型就与Smagorinsky-Lilly模型相同，只是模型常数有区别。在流动场的低雷诺数区域，上面的函数就小于零，从而只有分子粘性起作用。所以，基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果。 3.3.6.3 大涡模拟的边界条件 对于给定进口速度边界条件，速度等于各个方