湍流的数学模型简介【精心整理版】

湍流的数学模型报告人：指导老师：中南大学防灾科学与安全技术研究所2013.04Contents湍流导论1湍流的数学模型简介2湍流模型RANS3直接模拟DNS4大涡模拟LES5湍流燃烧模型简介6第1章湍流导论湍流现象描述湍流是一种高度复杂的三维非稳态、带旋转的不规则流动。湍流中流体的各个物理参数，如速度、压力、温度等都随时间与空间发生随机变化。湍流与层流自然界中的流体流动状态主要有两种形式，即层流(laminar)和湍流(trubulence)。层流是指流体在流动过程中两层之间没有相互混掺，而湍流是指流体不是处于分层流动状态。一般说来，湍流是普遍的，而层流则属于个别情况。判断流动是层流还是湍流，是看其雷诺数是否超过临界雷诺数。雷诺数的定义如下：式中：V为截面的平均速度；L为特征长度；为流体的运动粘度。当Re2000，管内流动保持稳定的层流状态。VLRe1.1、湍流的认识第1章湍流导论1.1、湍流的认识葛饰北斋的浮世绘作品《神奈川冲浪里》1.1湍流的认识湍流物理特征大尺度的涡旋小尺度的涡旋主要由流动边界条件决定，从主流获得能量，是引起低频脉动的原因。由于流体粘性的作用，不断消失，从而产生能量耗散；是引起高频脉动的原因。“随机”和“脉动”是湍流流场的重要的物理特征。1.1湍流的认识Kolmogorow尺度分布理论在描述湍流行为的理论中，Kolmogorov尺度分布理论是相当重要也是非常普适的一种。1Kolmogorow长度尺度湍流能量的耗散发生在小涡结构中，这一最小的湍流流动结构尺寸可用Kolmogorow长度尺度表示：2Kolmogorow时间尺度134Kl12KKolmogorow时间尺度表示最小湍流结构的动量扩散时间，它的定义为第1章湍流导论1.2湍流的统计平均法统计平均方法是湍流研究的开始.他将不规则的流场分解为规则的平均场和不规则的脉动场，同时也引出了封闭雷诺方程的世纪难题。湍流的随机性统计平均方法是处理湍流流动的基本手段，这是由湍流的随机性所决定的。研究湍流的统计平均方法在湍流理论中，有多种统计平均方法。例如时均法、体均法、按概率平均法（或称系综平均法）等。下面将分别予以讨论，然后在进行比较。1.2湍流的统计平均法1时均法时均法的确切定义是：001()()tTiitututdtT上式中的速度瞬时值是任一次试验结果，积分限中的下线可以任意取，即一次试验中，从任何时候开始都不能影响平均值的结果。当时间间隔T很长时，有：这时，速度时均值不再是时间的函数，这就是雷诺平均。应用时均法需满足下列要求：平均值与平均的起始时刻及时间间隔T（只要足够长）无关。而且平均值本身不再是时间的函数，因此，时均法只能用于讨论定常的湍流流动。001()lim()tTiitTututdtT1.2湍流的统计平均法2体均法湍流的随机变量不仅表现在时间上，在空间分布上也具有随机性。体均值要求与积分体积的大小及所处的坐标位置无关。因此严格说来，体均法只适用于描述对体均值而言的均匀的湍流流场。3概率平均法（系综平均法）()1()(,,,)iitVtdddV时均法和体均法只适用于两种特殊状态的湍流，前者适用于定常湍流，后者适用于均匀湍流。对于一般的不定常非均匀流，可以采用随机变量的一般平均法，即概率平均法。在相同条件下重复N次试验，再对此N次试验值取平均。若能对某种湍流找到相应的概率密度，则湍流问题就可认为已经解决。1pk1limNNkiittVNV1.2湍流的统计平均法三种平均法之间的关系及各态遍历假说tpiiiVVV时均法只适用于定常湍流，体均法只适用于均匀不定常湍流。在什么物理条件下，普遍适用的概率平均值和时均值或体均值等价？•各态遍历假说的思想：一个随机变量在重复许多次的试验中出现的所有可能状态，能够在一次试验的相当长的时间或相当大的空间范围内以相同的概率出现。•各态遍历假说的结论：对于一个满足各态遍历的系统，三种平均值相等•在各态遍历假说成立的前提下，可以用时均法研究不定常流动。脉动值随机值与平均值之差称为涨落，在湍流中称为脉动iiiVtVtVt脉动值是随机变量，平均值是统计的决定性变量，全部湍流理论就是研究脉动值和平均值之间的互相关系。第1章湍流导论1.3、湍流的基本方程湍流瞬时控制方程(包括连续方程、动量方程和能量方程)可用通用微分方程表示。一般认为，无论湍流流动多么复杂，非稳态的连续性方程和N-S方程(动量方程)仍然适用于湍流的瞬时流动。第1章湍流导论'iiiuuu1.3、湍流的基本方程（不可压）N-S方程平均值与脉动值之和为流动变量的瞬时值将非稳态N-S方程对时间作平均，即把湍流的运动看成是时间平均流动与瞬间脉动流动的叠加：'dtttttt1'iiiuuu'iiiPPP'iiiTTT1.3湍流的基本方程''iijijiijijjjuuuuuupftxxxxx''ijijuu以上为Reynolds时均方程，引入的Reynolds应力有6个未知分量，由于雷诺平均方程中未知数个数大大多于方程个数而出现了方程不封闭的问题。必须做假设引入雷诺应力的封闭模型即建立湍流模型才能求解出平均流场。Reynolds时均方程0iixu1.3湍流的基本方程ijc雷诺应力输运方程上式称为不可压缩湍流的雷诺应力输运方程，方程中各项分别用，，，，来表示。ijijPijDijEijcijijPijDijEijc雷诺应力在平均运动轨迹上的增长率。脉动压强和脉动速度变形率张量相关的平均值，称再分配项。雷诺应力与平均运动速度梯度的乘积，产生湍动能的关键，称生成项。具有扩散性质，称雷诺应力扩散项。脉动速度梯度乘积的平均值，使湍流能耗散，故称耗散项。1.3湍流的基本方程其它变量时均方程时均化的能量方程iiiiiijquTuDTccDtxxx雷诺热流二阶相关量3个未知量()()('')jjjjjuuStxxx第一章湍流导论1.4、湍流封闭问题湍流模式理论的主要任务就是研究湍流方程的封闭方法。核心问题求解雷诺应力第2章湍流的数值模拟方法简介2.1湍流数值模拟方法的分类湍流运动的数值模拟方法可以分为直接数值模拟方法和非直接数值模拟方法。所谓直接数值模拟方法是指求解瞬时湍流控制方程。非直接数值模拟方法就是不直接计算湍流的脉动特性，而是设法对湍流做某种程度的近似和简化处理。根据依赖所采用的近似和简化方法不同，非直接数值模拟方法分为大涡模拟、统计平均法和Reynolds平均法。第2章湍流的数值模拟方法简介2.2模型比较湍流模型方法(RANS方法)大涡模拟方法(LES方法)直接数值模拟(DNS方法)给出了时间平均的流动信息，易于工程应用抹去了流动的瞬态特性及细观结构，适合高雷诺数，不具普适性介于RANS与DNS之间，非常成功的应用于RANS不能满足要求的高端应用，如燃烧、混合、外部空气动力学。亚格子湍流模型有待进一步完善无需湍流模型，能精确给出湍流瞬态演变过程数值求解方法难度大，适合低雷诺数第3章湍流模型（RANS）不可压缩时均运动控制方程组之所以出现方程组不封闭（需求解的未知函数较方程数多），在于方程中出现了湍流脉动值的雷诺应力项。要使方程组封闭，必须对雷诺应力做出某些假定，即建立应力的表达式（或者引入新的湍流方程），通过这此表达式把湍流的脉动值与时均值等联系起来。基于某些假定所得出的湍流控制方程，称为湍流模型。所谓湍流模型，是依靠理论与经验的结合，引进一系列模型假设，把脉动值附加项与时均值联系起来的一些特定的关系式。3.1湍流模型的分类湍流涡粘模型雷诺应力模型1.湍流涡粘模型（Eddy-ViscosityModels，EVM）这类模型的处理方法不直接处理雷诺应力项，而是引入涡粘系数（EddyViscosity），然后把湍流应力表示成为涡粘系数的函数，整个计算关键在于确定这种湍流粘性系数。引入Boussinesq涡粘性假设，认为雷诺应力与平均速度梯度成正比，即将Reynolds应力项表示为湍流粘性系数ijijkkijjijiijkxUxUxUuu3232tt基于不同的假设，湍流模型分为湍动能：2'2'2'''212wvuii3.1湍流模型的分类t一方程模型•常系数模型•二维Prandtl混合长度理论minmaxuuCtyulmt2零方程模型一方程模型两方程模型lkCt2/1/零方程模型根据确定湍流粘性系数的微分方程数目，又可分为3.1湍流模型的分类由求解湍流特征参数的微分方程来确定湍流粘性。包括k-ε、k-ω、k-g模型等。其中，应用最普遍的是k-ε模型。两方程模型以上介绍的模型都是基于Boussinesq假设，认为湍流粘性系数各向同性，难于考虑旋转流动及流动方向表面曲率变化的影响，不适用于复杂流动。针对k-ε模型不足，许多学者对标准的模型进行了修正。重整化群k-ε模型（renormalizationgroup，RNGmodel）可实现k-ε模型（realizablek-εmodel）多尺度k-ε模型（multiscalemodelofturbulence）3.1湍流模型的分类雷诺应力方程模型（ReynoldsStressModel，RSM）由各项异性的前提出发，完全抛弃了Boussinesq表达式及的概念，直接建立以雷诺应力为因变量的微分方程，然后作适当假设使之封闭。这种模型也称为二阶封闭模型。t主要思想是设法将应力的微分方程简化为代数表达式，以减少RSM模型过分复杂的弱点，同时保留湍流各项异性的基本特点。代数应力方程模型（AlgebraicStressModel，ASM）2雷诺应力方程模型3.2湍流模型具体介绍双方程模型标准模型k可实现模型kRNG模型kReynolds应力模型（RSM）代数应力模型（ASM）零方程模型单方程模型3.2湍流模型具体介绍1零方程模型–代数涡粘模型这一假设并无物理基础，且采用各向同性的湍流动力粘度来计算湍流应力，难于考虑旋转流动和表面曲率变化的影响，但以此为基础的湍流模型目前在工程计算却应用最为广泛。所谓零方程模型就是不使用微分方程，而是基于Boussinesq1877年的假设，用代数关系式，把湍流粘度与时均值联系起来的模型，它只用湍流的时均连续方程和Reynolds方程组成方程组，把方程组中的Reynolds应力用平均速度场的局部速度梯度来表示。ijijkkijjijiijkxUxUxUuu3232tt1零方程模型零方程–Prandtal混合长度理论零方程中最著名的是Prandtl提出的混合长度模型（mixinglengthmodel）。混合长度定义：脉动微团在经历这段距离内保持有不变的脉动速度值。表示：微流微团的作用范围。混合长度模型的特点：直接用平均量梯度代数表达式来模拟Reynolds时均方程组中未知的应力或热流、物质流关联项。lm由实验或直观判断加以确定。)()(xyxlm42)1(06.0)1(08.014.0RyRyRlm对于自由剪切流充分发展的湍流管流yulmt21零方程模型普朗特混合长度模型的评价•优点：直观、简单，无须附加湍流特性的微分方程——适用于简单流动，如射流、边界层、管流、喷管流动等。另外，研究历史较长，积累了很多经验。•缺点1：在处必然是湍流粘性T为零，或剪力、热流、扩散流均为零——与实际不符。0uy混合长度模型相当于湍流能量达到局部平衡，即湍流的产生等于湍流的耗散，亦