湍流的数学模型课件

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,湍流的数学模型,第六讲,流体仿真与应用,湍流的数学模型第六讲流体仿真与应用,1,湍流模型,雷诺应力模型（RSM）,雷诺应力方程是微分形式的，称为,雷诺应力方程模型,。,Reynolds Stress equation ModelRSM,若将雷诺应力方程的微分形式简化为代数方程的形式，则称为,代数应力方程模型,。,这个模型的特点是直接构建表示雷诺应力的补充方程，然后联立求解湍流时均运动控制方程组。,湍流模型雷诺应力模型（RSM）雷诺应力方程是微分形式的,2,湍流模型,Reynolds应力输运方程,第一步，建立 的输运方程。首先，将 乘以 的N-S方程与 乘以 的,N-S,方程相加，得到的 方程，再对此方程作Reynolds时均、分解，即得到 的输运方程,。,第二步，建立 的输运方程。乘以 的Reynolds时均方程与 乘以 的Reynolds时均方程相加，即得到的 输运方程。,湍流模型Reynolds应力输运方程 第一步，建立,3,湍流模型,Reynolds应力输运方程,湍流模型Reynolds应力输运方程,4,湍流模型,Reynolds应力输运方程,湍流模型Reynolds应力输运方程,5,湍流模型,Reynolds应力输运方程,湍动扩散项 的计算,推荐用下式：,Daly和Harlow所给出的广义梯度扩散模型来计算：,湍流模型Reynolds应力输运方程 湍动扩散项,6,湍流模型,Reynolds应力输运方程,浮力产生项 的计算,热膨胀系数,理想气体,不可压缩流体,湍流模型Reynolds应力输运方程 浮力产生项,7,湍流模型,Reynolds应力输运方程,压力应变项 的计算,压力应变项 的存在是雷诺应力方程与 方程的,最大区别之处,。因此，仅在湍流各分量项存在，当 时，它表示减小剪应力，使湍流趋向各向同性；当 时，它表示使湍动能在各应力分量间重新分配，对总量无影响，可见，此项并不产生脉动能量，仅起到再分配的作用。因此，在的的文献中称此项为,再分配项,。,慢的压力应变项,快的压力应变项,壁面反射项,湍流模型Reynolds应力输运方程 压力应变项,8,湍流模型,Reynolds应力输运方程,压力应变项 的计算,湍流模型Reynolds应力输运方程 压力应变项,9,湍流模型,Reynolds应力输运方程,粘性耗散项 的计算,耗散项表示分子粘性对雷诺应力产生的耗散。在建立耗散项的计算公式时，认为大尺度旋涡承担动能输运，小尺度旋涡承担粘性耗散，因此小尺度涡团可看成是各向同性的。即认为局部是各向同性的。,湍流模型Reynolds应力输运方程 粘性耗散项,10,湍流模型,Reynolds应力输运方程,FLUENT等大多数CFD软件所使用的广义的Reynolds应力输运方程,湍流模型Reynolds应力输运方程 FLUENT等大多数,11,湍流模型,Reynolds应力输运方程,若不考虑浮力的作用（即 ）及旋转的影响（），同时压力应变项不考虑壁面反射（=0）,湍流模型Reynolds应力输运方程 若不考虑浮力的作用,12,湍流模型,RSM的控制方程（补充方程）,时均连续性方程、雷诺方程、,雷诺,应力方程、方程和 方程等，共12个方程构成了封闭的三维湍流流动问题的基本控制方程组。,湍流模型RSM的控制方程（补充方程）时均连续性方程、雷,13,湍流模型,RSM的控制方程的说明,如果需要对能量或组分等进行计算，需要建立针对标量型变量（如温度、组分浓度）的脉动量的控制方程。,由于从Reynolds应力输运方程的3个正应力项可以得出脉动动能，即 。因此，曾有文献不把 作为独立变量，也不引入 方程，但大多数文献把 方程列入控制方程之一，就本书所采用的方式一样。,湍流模型RSM的控制方程的说明 如果需要对能量或组分等进,14,湍流模型,RSM适用性,RSM是针对湍流发展非常充分的湍流运动来建立的。即是针对高Re湍流模型，而当Re较低时，上述方程不再适用。常用解决壁面流动方法有：,一种是壁面函数法,；,另一种是采用低Re的RSM模型,。,低Re的RSM模型基本思想是修正高Re的RSM模型耗散函数（扩散项）及压力应变重新分配项的表达式，以使RSM模型方程可以直接应用到壁面区域。,计算实践表明，RSM虽然能考虑一些各向异性效应，但并不一定比其它模型效果好，在计算,突扩流动分离区,和,计算湍流输运各向异性,较强的流动时，RSM优于两方程模型，采用RSM意味着要多解6个雷诺应力的微分方程，计算量大，对计算机的要求高。因此，RSM不如 模型应用更广泛，但许多文献认为RSM是一种更有潜力的湍流模型。,湍流模型RSM适用性 RSM是针对湍流发展非常充分的湍流,15,湍流模型(LES),大涡模拟（LES-Large Eddy Simulation）,湍流中包含了不同时间与长度尺度的涡旋。最大长度尺度通常为平均流动的特征长度尺度。最小尺度为Komogrov尺度。,LES的基本假设是,1，动量、能量、质量及其它标量主要由大涡输运；2，流动的几何和边界条件决定了大涡的特性，而流动特性主要在大涡中体现；,3，小尺度涡旋受几何和边界条件影响较小，并且各向同性；大涡模拟过程中，直接求解大涡，不直接计算小尺度涡旋，从而使得网格要求比DNS低。,湍流模型(LES)大涡模拟（LES-Large Eddy,16,湍流模型(LES),大涡模拟（LES）控制方程,LES的控制方程是对Navier-Stokes方程在物理空间进行过滤得到的。过滤的过程是去掉比过滤宽度或者给定物理宽度小的涡旋，从而得到大涡旋的控制方程。,过滤变量（上横线）定义为,D表示流体区域；G是决定涡旋大小的过滤函数,在FLUENT中，有限控制体离散本身暗中包括了过滤运算,湍流模型(LES)大涡模拟（LES）控制方程LES的控制方,17,湍流模型(LES),大涡模拟（LES）控制方程,大涡模拟对不可压流动问题得到较多应用，但在可压缩问题中的应用还很少，因此在FLUENT中，大涡模拟只能,针对不可压流体的流动,。,过滤不可压的Navier-Stokes方程后，可以得到LES控制方程,上述方程与雷诺平均方程很相似，只不过大涡模拟中的变量是过滤过的量，而非时间平均量，并且湍流应力也不同。,亚网格应力,湍流模型(LES)大涡模拟（LES）控制方程大涡模拟对不可,18,湍流模型(LES),亚网格模型,由于LES中亚网格应力项是未知的，并且需要模拟以封闭方程。目前，采用比较多的亚网格模型为涡旋粘性模型，形式为：,亚网格湍流粘性系数,求解尺度下的应变率张量,湍流模型(LES)亚网格模型由于LES中亚网格应力项是未知,19,湍流模型(LES),亚网格模型(Smagorinsky-Lilly模型),求解亚网格湍流粘性系数,Smagorinsky-Lilly模型,是亚网格的混合长度,是Smagorinsky常数,0.23，对于有平均剪切或者过渡流动中，该系数过高估计了大尺度涡旋的阻尼作用,0.1有比较好的模拟结果,湍流模型(LES)亚网格模型(Smagorinsky-Li,20,湍流模型(LES),亚网格模型(重整化群亚网格模型),求解亚网格湍流粘性系数,V是计算控制体体积；重整化群常数 ,而常数C=100,H(x)是Heaviside函数,湍流模型(LES)亚网格模型(重整化群亚网格模型),21,湍流模型(LES),亚网格模型(重整化群亚网格模型),对于高雷诺数流动（），基于重整化群理论的亚网格模型就与Smagorinsky-Lilly模型相同，只是模型常数有区别。,在流动场的低雷诺数区域，上面的函数就小于零，从而只有分子粘性起作用。,所以，基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果。,湍流模型(LES)亚网格模型(重整化群亚网格模型)对于高,22,湍流模型总览,湍流模型总览,23,