大涡模拟概述

湍流数值模拟方法简介——大涡模拟专业：动力机械与工程姓名：xxxxxx学号：xxxxxxxx2主要内容湍流数值模拟方法简介大涡模拟的基本思想过滤函数亚网格模型3一、湍流数值模拟方法简介4一、湍流数值模拟方法简介)()(~)(xfxfxff)(')()(xfxfxffReynoldsaveragedmodel(RAN)LargeEddySimulation(LES)5一、湍流数值模拟方法简介LESRANS过滤脉动对象空间上进行过滤，只滤去小涡脉动保留大涡脉动，得到瞬时速度时间上取平均，抹去涡团在时间上的脉动值，得到平均速度涡团尺度分辨能力取决于网格分辨率，远远小于RANS所能求解的尺度取决于所采用的湍流模型并非网格尺寸的函数边界条件初始条件必须包含对湍流脉动在时间空间上足够的分辨率，匹配时间步长和网格尺寸只需要边界上的时均信息数值精度网格越精细时间步长越小，其数值扩散越小，数值越精确与网格密度无关，精度低于LESLES与RANS的区别61.大涡模拟的基本思路二、大涡模拟的基本思想71.大涡模拟的基本思路二、大涡模拟的基本思想大涡模拟的两大问题及解决方案一.如何实现大尺度涡和小尺度脉动的分离二.如何封闭过滤后的N-S方程通过对亚网格应力不同的简化就构成了不同的亚网格模型。过滤函数均匀过滤器非均匀过滤器盒式过滤器高斯过滤器过滤函数常系数Smagorinsky模型动态Smagorinsky模型尺度相似模型8二、大涡模拟的基本思想uxPFDtDubx2不可压缩牛顿流体x方向瞬时量N-S方程：bxFyuxPyuvtu221bxFyuxPyuvtu221过滤过滤后X方向瞬时量N-S方程：过滤过程的基本性质归一性：常量过滤前后不变可交换性：过滤运算和微分及加减运算的运算顺序可交换9三、过滤函数湍流脉动的过滤过滤的基本作用：过滤是一种数学运算。过滤的作用就是突出所关注的。过滤的标准即为过滤尺度。过滤尺度：可解尺度湍流：可直接计算的大尺度脉动亚网格尺度湍流：不直接计算的小尺度脉动过滤尺度10三、过滤函数过滤是一种数学运算，在物理空间中，过滤过程可以通过积分过程来实现。物理空间中用以分离可解尺度和不可解尺度湍流的尺度称为过滤尺度，用∆表示。10'iii过滤前瞬时变量值过滤后的变量——可解尺度变量，其湍流尺度大于∆被过滤器平均化的小尺度变量——不可解尺度变量注意：这里的不是时间上的平均，而是在空间域上的平均。11三、过滤函数湍流脉动的过滤过滤方法的种类：过滤器均匀过滤器非均匀过滤器均匀盒式过滤器高斯过滤器空间三维过滤器微分过滤器非均匀卷积型过滤器非均匀盒式过滤器2阶精度可交换盒式过滤器非均匀三维过滤器'1||/2(){0||/2xGxxx2'62'()/rrGxxAexx'1/||/2(){0||/2xGxxx12三、过滤函数湍流脉动的过滤盒式过滤器：xVxxdxGtutxu)(),(),(其中，是流场中其他流体质点的x向坐标值；为x方向过滤函数；)(xxG)22(x11均匀盒式非均匀盒式0（其他情况）其中，是过滤尺度'iiix)(xxG13•物理空间的盒式过滤器（a）物理空间图形；三、过滤函数盒式过滤器滤波方法很简单，缺点是它的傅立叶变换在某些区间里有负值，并且由于滤波函数在单元边界上的间断性，难以进行微分运算。Top-hat14•高斯过滤器将过滤函数取作高斯函数，称为高斯过滤器。高斯过滤器物理空间下的数学表达式如下：21/22266()exp()xGx在谱空间的表达式：22()exp()24kGk三、过滤函数15•高斯过滤器（a）物理空间图形；高斯过滤器在物理空间和谱空间都有很好的性能，可以任意次微分。高斯滤波器性能最好，但计算很麻烦，目前用得最多得还是盒式滤波器，因为它们简单方便，易于实现。三、过滤函数Gaussian16盒式过滤器过滤运算和求导运算可以交换，用盒式过滤器对不可压缩牛顿流体N-S方程进行过滤运算，过滤后的N-S方程如下：1()()iiijjijjuupuutxxxx1()()iiijjijjuupuutxxxx1()()()iiijijijjijjjuupuuuuuutxxxxx()iijiijjjuuuuuuuu亚网格应力ij物理意义：亚网格应力是可解尺度脉动和过滤掉的小尺度脉动（不可解尺度脉动）间的动量输运。可以看到滤波后的方程里除了有大尺度涡的未知量和，还出现了新的未知量亚网格应力，这样，方程组就不能封闭，需要构造模型使过滤后的N-S方程封闭。所构造的模型就是下面将要介绍的亚网格模型。iuipijijijuuuu四、亚网格模型17要准确理解已有的亚网格模型，和构造新的亚网格模型，有必要了解一下亚网格力的机制：1()()()iiijijijjijjjuupuuuuuutxxxxxijijijuuuu''''ijijijijjiijuuuuuuuuuuijijijLuuuu''ijijjiCuuuu''ijijRuu代数运算称为里昂纳特（Leonard）应力，它由可解尺度间的相互作用产生'iiiuuu称为交叉应力，它是可解尺度脉动和不可解尺度脉动动量交换对亚网格应力的贡献；称为亚网格雷诺应力，它是不可解尺度脉动的脉动动量交换对亚网格应力的贡献。四、亚网格模型18四、亚网格模型亚网格尺度模型的求解T常系数Smagorinsky模型亚网格尺度模型2/12T)(CijijsSS动态Smagorinsky模型),,(CCtyxss尺度相似模型（SSM）122()(2)ijijisijijjuuuuCSS19四、亚网格模型Smargorinsky模型是最早提出的亚网格应力模型，是参照雷诺平均模式的涡粘模型，以各向同性湍流为基础，认为亚网格湍流具有混合长度型涡粘系数。实际应用表明Cs应取更小的值，以减小亚网格应力的扩散影响。19123ijtijijkkS亚网格涡粘系数122(2)tsijijCSS1()2jiijjiuuSxx可解尺度的变形率张量Smagorinsky常数Cs=0.1820四、亚网格模型Smargorinsky模型的优点是，概念简单、易于实施且计算方便，只要增加一个涡粘系数的模块，就可以利用N-S方程的数值计算方法和程序；主要缺陷是耗散过大，属于唯象论模型。尤其是壁面处，该影响尤为明显，可以利用近壁阻尼系数对Smargorinsky系数Cs做修正：200.18[1exp(/)]sCyA/A26yyu近壁阻尼系数到壁面最近距离半经验常数21四、亚网格模型——尺度相似模型为了克服唯象论亚网格涡粘和涡扩散的缺陷，使亚网格模型能够适应复杂的湍流运动，提出了尺度相似模型。前面曾经指出，亚网格应力实质上是可解尺度湍流和不可解尺度湍流间的动量输运。从输运机制出发建立亚网格应力模型是正确的途径。尺度相似模型（SSM）假定可解尺度中的最小尺度脉动和不可解尺度脉动具有相似性，根据这一假定，依据Bardina尺度相似模式可以推导出亚网格应力：21()ijijijijijijLCRuuuu''iiuu根据尺度相似概念，将可解尺度脉动，再做一次过滤得到，则：'iiiuuu''iiiiuuuu将上式代入亚网格雷诺应力和交叉应力中，与Leonard应力相加，得亚网格应力公式：22相似模型的优点是能够比较准确地表达可解尺度和不可解尺度间的动量输运关系，这一点已由各向同性湍流的直接数值模拟结果证实。缺点是该模型抛弃了涡粘假设，它不是单纯耗散性的，既可以由可解尺度湍流向不可解尺度湍流输送能量，也可以由逆传。故存在严重耗散不足；此外，由于存在逆传（相当于负涡粘系数），数值计算稳定性很差。综合Smargorinsky和SSM模型各自优点，发展了混合模型，即将两个模型做线性叠加：该模型既有正确的亚网格动量输运，又有足够的亚网格耗散。22122()(2)ijijisijijjuuuuCSS四、亚网格模型——尺度相似模型23四、亚网格模型——动态模型类似于尺度相似的思想，20世纪90年代发展了动态模型。动态模型实际上是动态确定亚网格涡粘模型的系数，主要对流场做两次过滤，一次是细过滤，过滤后再做一次粗过滤。以∆1过滤的可解速度用上标“—”表示，以∆2过滤的可解速度用上标“~”表示，﹒表示在空间均匀的方向上取平均。其基本思想是：粗过滤中的小尺度脉动和细过滤的脉动相似。根据这一原则采用Germano等式确定Smargorinsky模型中的系数就叫做动态Smargorinsky模型23122(2)tsijijCSS22212[||||]ijijsijijijijijijijijMLCMMLuuuuMSSSS24采用二次滤波的方法建立亚网格应力模型小尺度G-levelF-levelGermano恒等式：F-滤波+G-滤波与FG滤波之间的关系式F-level滤波滤波尺度为，G-level滤波滤波尺度为FG-level滤波：ffkfˆjijiijuuuujijiijuuuuTˆˆ^ˆˆ(1)ijijijijijTuuuuL特点：该量无需模型，可直接计算FG滤波F滤波+G滤波^^ˆijijijuuuu四、亚网格模型——动态模型2525特点：无需模化，可“精确”算出^ˆˆijijijijTuuuuijijTFG滤波（）亚网格应力k经过G-滤波后的F-滤波（）亚网格应力Germano恒等式启发：Germano提供了亚网格模型的一个约束条件，可用来改进模型)()(22uSuSCijdij(,,)ijfuC模型系数，动态可调，需要计算ˆ(,,)ijTfukC^^ˆˆˆ(,,)(,,)ijijijijTfukCfuCuuuu仅C是未知数，可解6个方程1个未知数，通常采用最小二乘解四、亚网格模型——动态模型26动态涡粘模型F-levelFG-level22ˆˆˆˆ(()()()(()()))ijijdijijdijTCkSuSuSuSuCMijdijijijMCTLijijijijdMMLMC2()()ijdijCSuSu2ˆˆ()()()ijdijTCkSuSu涡粘系数C动态可调通过两次滤波，确定该系数FG滤波，相当于用进行滤波k可直接计算，无需模型四、亚网格模型——动态模型27•在近壁区，湍涡的长度尺度趋于零，若要分辨出所有涡团，理论上需要无限小的网格。在LES中的粗网格下是无法准确模拟近壁区域的。•因而引入壁函数，即在黏性底层之外建立的半经验边界条件。如Groetzbach-Schmann提出的公式：0nuloglog2,,wtstwuzyxu式中，表示Favre平均，和分别是垂直和平行于壁面的速度分量，是壁面切应力，log表示按照对数律得出的平均量。.nutuw四、亚网格模型——壁面条件28谢谢各位！