1、湍流数值模拟方法简介 大涡模拟,专业：动力机械与工程姓名：xxxxxx学号：xxxxxxxx,主要内容,湍流数值模拟方法简介 大涡模拟的基本思想 过滤函数 亚网格模型,一、湍流数值模拟方法简介,一、湍流数值模拟方法简介,f,f,Reynolds averaged model (RAN),Large Eddy Simulation (LES),一、湍流数值模拟方法简介,LES与RANS的区别,大涡模拟的基本思路,二、大涡模拟的基本思想,大涡模拟的基本思路,二、大涡模拟的基本思想,大涡模拟的两大问题及解决方案,一. 如何实现大尺度涡和小尺度脉动的分离,二. 如何封闭过滤后的N-S方程通过对亚网格应
2、力不同的简化就构成了不同的亚网格模型。,过滤函数,均匀过滤器,非均匀过滤器,盒式过滤器 高斯过滤器,过滤函数,常系数Smagorinsky模型,动态Smagorinsky模型,尺度相似模型,二、大涡模拟的基本思想,不可压缩牛顿流体x方向瞬时量N-S方程：,过滤,过滤后X方向瞬时量N-S方程：,过滤过程的基本性质,归一性：常量过滤前后不变,可交换性：过滤运算和微分及加减运算的运算顺序可交换,三、过滤函数,湍流脉动的过滤,过滤的基本作用：,过滤是一种数学运算。过滤的作用就是突出所关注的。过滤的标准即为过滤尺度。,过滤尺度：,可解尺度湍流： 可直接计算的大尺度脉动,亚网格尺度湍流： 不直接计算的小尺
3、度脉动,过滤尺度,三、过滤函数,过滤是一种数学运算，在物理空间中，过滤过程可以通过积分过程来实现。物理空间中用以分离可解尺度和不可解尺度湍流的尺度称为过滤尺度，用表示。,10,过滤 前瞬 时变 量值,过滤后的变量可解尺度变 量，其湍流尺度大于,被过滤器平均化的小尺度变量 不可解尺度变量,注意：这里的 不是时间上的平均，而是在空间域上的平均。,三、过滤函数,湍流脉动的过滤,过滤方法的种类：,过滤器,均匀过滤器,非均匀过滤器,均匀盒式过滤器,高斯过滤器,空间三维过滤器,微分过滤器,非均匀卷积型过滤器,非均匀盒式过滤器,2阶精度可交换盒式过滤器,非均匀三维过滤器,三、过滤函数,湍流脉动的过滤,盒式过
4、滤器：,其中， 是流场中其他流体质点的x向坐标值； 为x方向过滤函数；,均匀盒式,非均匀盒式,（其他情况）,其中， 是过滤尺度,物理空间的盒式过滤器,（a）物理空间图形；,三、过滤函数,盒式过滤器滤波方法很简单，缺点是它的傅立叶变换在某些区间里有负值，并且由于滤波函数在单元边界上的间断性，难以进行微分运算。,Top-hat,高斯过滤器,将过滤函数取作高斯函数，称为高斯过滤器。高斯过滤器物理空间下的数学表达式如下：,在谱空间的表达式：,三、过滤函数,高斯过滤器,（a）物理空间图形；,高斯过滤器在物理空间和谱空间都有很好的性能，可以任意次微分。高斯滤波器性能最好，但计算很麻烦，目前用得最多得还是盒
5、式滤波器，因为它们简单方便，易于实现。,三、过滤函数,Gaussian,盒式过滤器过滤运算和求导运算可以交换，用盒式过滤器对不可压缩牛顿流体N-S方程进行过滤运算，过滤后的N-S方程如下：,亚网格应力,物理意义：亚网格应力是可解尺度脉动和过滤掉的小尺度脉动（不可解尺度脉动）间的动量输运。,可以看到滤波后的方程里除了有大尺度涡的未知量 和 ，还出现了新的未知量亚网格应力 ，这样，方程组就不能封闭，需要构造模型使过滤后的N-S方程封闭。所构造的模型就是下面将要介绍的亚网格模型。,四、亚网格模型,要准确理解已有的亚网格模型，和构造新的亚网格模型，有必要了解一下亚网格力的机制：,代数运算,称为里昂纳特
6、（Leonard）应力，它由可解尺度间的相互作用产生,称为交叉应力，它是可解尺度脉动和不可解尺度脉动动量交换对亚网格应力的贡献；,称为亚网格雷诺应力，它是不可解尺度脉动的脉动动量交换对亚网格应力的贡献。,四、亚网格模型,四、亚网格模型,亚网格尺度模型,常系数Smagorinsky模型,亚网格尺度模型,动态Smagorinsky模型,尺度相似模型（SSM）,四、亚网格模型,Smargorinsky模型是最早提出的亚网格应力模型，是参照雷诺平均模式的涡粘模型， 以各向同性湍流为基础，认为亚网格湍流具有混合长度型涡粘系数 。实际应用表明Cs应取更小的值，以减小亚网格应力的扩散影响。,19,亚网格涡
7、粘系数,可解尺度的 变形率张量,Smagorinsky常数 Cs= 0.18,四、亚网格模型,Smargorinsky模型的优点是，概念简单、易于实施 且计算方便，只要增加一个涡粘系数的模块，就可以利用N-S方程的数值计算方法和程序；主要缺陷是耗散过大，属于唯象论模型。尤其是壁面处，该影响尤为明显，可以利用近壁阻尼系数对Smargorinsky 系数Cs做修正：,20,近壁阻尼 系数,到壁面 最近距 离,半经验常数,四、亚网格模型尺度相似模型,为了克服唯象论亚网格涡粘和涡扩散的缺陷，使亚网格模型能够适应复杂的湍流运动，提出了尺度相似模型。前面曾经指出，亚网格应力实质上是可解尺度湍流和不可解尺度
8、湍流间的动量输运。从输运机制出发建立亚网格应力模型是正确的途径。尺度相似模型（SSM）假定可解尺度中的最小尺度脉动和不可解尺度脉动具有相似性， 根据这一假定，依据Bardina尺度相似模式可以推导出亚网格应力：,21,根据尺度相似概念，将可解尺度脉动 ，再做一次过滤得到 ，则：,将上式代入亚网格雷诺应力和交叉应力中，与Leonard应力相加，得亚网格应力公式：,相似模型的优点是能够比较准确地表达可解尺度和不可解尺度间的动量输运关系，这一点已由各向同性湍流的直接数值模拟结果证实。缺点是该模型抛弃了涡粘假设，它不是单纯耗散性的，既可以由可解尺度湍流向不可解尺度湍流输送能量，也可以由逆传。故存在严重
9、耗散不足；此外，由于存在逆传（相当于负涡粘系数），数值计算稳定性很差。综合 Smargorinsky和 SSM模型各自优点，发展了混合模型，即将两个模型做线 性叠加：该模型既有正确的亚网格动量输运，又有足够的亚网格耗散。,22,四、亚网格模型尺度相似模型,四、亚网格模型动态模型,类似于尺度相似的思想，20世纪90年代发展了动态模型。动态模型实际上是动态确定亚网格涡粘模型的系数，主要对流场做两次过滤，一次是细过滤，过滤后再做一次粗过滤。以1过滤的可解速度用上标“”表示，以2过滤的可解速度用上标“”表示，表示在空间均匀的方向上取平均。其基本思想是：粗过滤中的小尺度脉动和细过滤的脉动相似。根据这一原
10、则采用Germano等式确定Smargorinsky模型中的系数就叫做动态Smargorinsky模型,23,采用二次滤波的方法建立亚网格应力模型,小尺度,G-level,F-level,Germano 恒等式： F-滤波+ G-滤波 与 FG滤波之间的关系式,F-level 滤波 滤波尺度为 ，G-level滤波 滤波尺度为 FG-level滤波：,特点： 该量无需模型，可直接计算,FG滤波,F滤波+ G滤波,四、亚网格模型动态模型,25,特点： 无需模化，可“精确”算出,FG滤波 （ ）亚网格应力,经过G-滤波后的 F-滤波（ ）亚网格应力,Germano恒等式,启发： Germano 提
11、供了亚网格模型的一个约束条件，可用来改进模型,模型系数，动态可调，需要计算,仅C是未知数，可解,6个方程1个未知数，通常采用最小二乘解,四、亚网格模型动态模型,动态涡粘模型,F-level,FG-level,涡粘系数C动态可调,通过两次滤波，确定该系数,FG滤波，相当于用 进行滤波,可直接计算，无需模型,四、亚网格模型动态模型,在近壁区，湍涡的长度尺度趋于零，若要分辨出所有涡团，理论上需要无限小的网格。在LES中的粗网格下是无法准确模拟近壁区域的。 因而引入壁函数，即在黏性底层之外建立的半经验边界条件。如Groetzbach-Schmann提出的公式：,式中， 表示Favre平均， 和 分别是垂直和平行于壁面的速度分量， 是壁面切应力，log表示按照对数律得出的平均量。,四、亚网格模型壁面条件,谢谢各位！,
