《fluent教程湍流模拟解析.pptx》由会员分享，可在线阅读，更多相关《fluent教程湍流模拟解析.pptx（53页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、湍流：非定常，非周期性的三维速度脉动、强化物质、动量和能量的输运.瞬时速度分解为平均速度和脉动速度:Ui(t)Ui+ui(t)压力、温度、组分浓度值具有类似的脉动什么是湍流什么是湍流?TimeU i(t)Uiui(t)第1页/共53页平均量输运方程平均量输运方程雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有：其中，和 分别是平均速度和脉动速度（i=1,2,3）类似地，对于压力等其它标量，我们也有：把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度 上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式 第2页/共53
2、页平均量输运方程（续）平均量输运方程（续）上面两个方程称为雷诺平均的Navier-Stokes（RANS）方程。如果要求解该方程，必须模拟该项以封闭方程。第3页/共53页湍流模拟的方法湍流模拟的方法直接数值模拟（DNS）只适合低雷诺数流动。求解雷诺平均的 Navier-Stokes(RANS)方程:其中 (雷诺应力)时间平均的湍流脉动用实验常数和流场平均速度信息来模拟.大涡模拟（Large Eddy Simulation，LES）数值求解大涡，小涡用模型.(定常,不可压缩流动 有/无 体积力)第4页/共53页湍流输运模型湍流输运模型涡旋粘性系数模型 涡旋粘性系数模型是根据湍流应力和分子运动引起
3、的粘性应力相似湍流应力和分子运动引起的粘性应力相似而提出来的 雷诺应力模型 雷诺应力模型则是从基本方程出发，直接推出雷诺应力的方程，但在雷诺应力的方程中包含有更高阶的相关项，对这些更高阶的相关项再建立相应的模型。第5页/共53页湍流应力湍流应力湍流应力湍流应力 在三维空间中，下标和分别可取为1，2和3，所以湍流应力有9个分量组成，是一个二阶张量，若用矩阵形式表示可写为：主对角线上的三个分量，和称为湍流正应力，其余的六个分量称为湍流切应力。三个湍流正应力之和是湍流脉动动能的两倍 第6页/共53页湍流应力湍流应力若以主对角线上的三个分量作为对称轴，则对称的两个切应力分量是相等的，很显然，这是一个对
4、称的二阶张量。很容易可以证明，在各向同性湍流中，湍流正应力的三个分量相等，即 第7页/共53页判断湍流判断湍流外流外流内流内流自然对流自然对流沿表面绕流其中其中其它因素，如自由流湍流度,表面条件,扰动 可能导致流动从层流向湍流转捩L=x,D,Dh,etc.第8页/共53页湍流特点湍流特点额外应变率流向曲率测向分离加速或减速有旋回流(或分离)二次流3D振荡流动Transpiration(吹风/吸气)自由湍流剪切层相互作用第9页/共53页需要作出选择需要作出选择湍流模型湍流模型和和近壁处理近壁处理 物理流体精度要求计算资源计算时间要求计算网格计算网格第10页/共53页Zero-Equation M
5、odelsOne-Equation Models Spalart-AllmarasTwo-Equation Models Standard k-RNG k-Realizable k-Reynolds-Stress ModelLarge-Eddy SimulationDirect Numerical Simulation湍流模型湍流模型IncludeMorePhysics每次迭代每次迭代增加增加计算时间计算时间FLUENT 5的湍流模型基于雷诺平均（RANS）的模型Fluent 6第11页/共53页涡旋粘性系数模型涡旋粘性系数模型 涡旋粘性系数模型是Boussinesq于1877年最早提出来的。
6、他假定在近似平行的剪切流中，湍他假定在近似平行的剪切流中，湍流应力张量中的切应力分量和平均速度在横向方向的梯度成正比，其比例系数称为涡旋粘性系数流应力张量中的切应力分量和平均速度在横向方向的梯度成正比，其比例系数称为涡旋粘性系数，以后把这个假定再推广到三维的流动。即湍流应力张量和平均流场应变率之间假定再推广到三维的流动。即湍流应力张量和平均流场应变率之间有关系：仿照分子运动引起的粘性系数，我们可以假定涡旋粘性系数和湍流中含能涡旋的特征长度和速度尺度成正比，即：而根据确定速度尺度V和长度尺度L方法的不同，又可以细分为各种不同的模型。而其中最简单，V和L的确定方法又是一致的模型应该是双方程模型。第
7、12页/共53页涡旋粘性系数模型涡旋粘性系数模型在双方程模型中，假定特征速度V和湍流动能k的平方根成正比，特征长度由湍流动能k和另外一个辅助的量确定。比如在k-e模型中，辅助的量选为湍流动能的耗散率，根据量纲分析得长度尺度，模型第13页/共53页涡旋粘性系数模型涡旋粘性系数模型辅助的量是湍流频率，根据量纲分析得长度尺度 模型第14页/共53页推导我们都假定流体不可压，即Favre平均和雷诺平均方程完全相同。瞬时的动量方程减去平均的动量方程得脉动速度的方程 乘以乘以并求平均，利用湍流动能的定义并求平均，利用湍流动能的定义 第15页/共53页上式右端第二项可以重新整理成：这里湍流动能的方程，方程右
8、端各项依次为输运项，产生项和耗散项 第16页/共53页各项对求导数，乘以，并求平均，可得耗散率的方程 第17页/共53页目前采用的标准ke模型方程为 如果用湍流频率代替湍流动能耗散率，频率的模型方程为：第18页/共53页uRANS 方程需要对雷诺应力进行封闭.u对于单方程模型，从修正的粘性系数输运方程简介求解湍流粘性系数。u对于双方程模型,湍流粘性系数根据湍动能(TKE)和耗散率 TKE来确定.uTransport equations for turbulent kinetic energy and dissipation rate are solved so that turbulent v
9、iscosity can be computed for RANS equations.Reynolds Stress Terms in RANS-based ModelsTurbulent Kinetic Energy:Dissipation Rate of Turbulent Kinetic Energy:湍流粘性系数:Boussinesq Hypothesis:(isotropic stresses)第19页/共53页湍流模型湍流模型0方程模型单方程模型双方程模型雷诺应力模型大涡模拟第20页/共53页湍流粘性系数:求解 的输运方程:The additional variables are
10、 functions of the modified turbulent viscosity and velocity gradients.单方程模型单方程模型:Spalart-AllmarasGenerationDiffusionDestruction第21页/共53页单方程模型单方程模型:Spalart-Allmaras用于计算航空中考虑边界作用的高速流动问题可以考虑边界层逆压梯度叶轮机械流动壁面网格可粗可细网格细了，可以考虑低雷诺数流动，考虑边界层影响.网格粗时候，可以得到相对较好的流动结果.不能预测均匀各向同性湍流的耗散。并且，单方程模型没有考虑长度尺度的变化，这对一些流动尺度变换比较
11、大的流动问题不太适合。比如，平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显。第22页/共53页应用举例应用举例无粘流结果 单方程湍流模型结果第23页/共53页无粘流结果 单方程湍流模型结果第24页/共53页双方程模型双方程模型:标准标准 k-模型模型湍动能方程湍动能方程耗散率输运方程耗散率输运方程实验常数(equations written for steady,incompressible flow w/o body forces)ConvectionGenerationDiffusionDestructionDestructionConvectionGeneratio
12、nDiffusion第25页/共53页双方程模型双方程模型:标准标准 k-模型模型“Baseline model”(Two-equation)最广泛用于工程流动与换热计算结果比较透彻，优点、缺点明显半经验性k 方程经过严格推导得到 方程通过物理推理得到只对完全湍流适用(大雷诺数流动)对于多数的湍流流动问题有较为合理的结果一般工程流动问题换热问题第26页/共53页双方程模型双方程模型:Realizable k-与标准 k-模型的区别:湍流粘性系数 是变量，其中 (A0,As,和 U*是速度梯度的函数)雷诺正应力为正;Schwarz 不等式成立新的耗散率 输运方程:GenerationDiffus
13、ionDestructionBuoyancy第27页/共53页与标准的 k-模型具有相同的湍动能输运方程在如下流动中，比标准的双方程模型更具有优越性:平板射流、圆射流边界层有强的逆压力梯度，分离有旋（rotation）,回流（recirculation）强曲率影响双方程模型双方程模型:Realizable k-第28页/共53页双方程模型双方程模型:RNG k-湍动能方程湍动能方程耗散率输运方程耗散率输运方程ConvectionDiffusionDissipationGeneration其中用RNG理论推导得到(equations written for steady,incompressib
14、le flow w/o body forces)Additional termrelated to mean strain&turbulence quantitiesConvectionGenerationDiffusionDestruction第29页/共53页双方程模型双方程模型:RNG k-k-两个方程都从瞬时NS方程出发，采用重整化群方法(Renormalization Group Method)严格推导得到.与标准的 k-方程有类似的形式、但:在 方程中多出一项用于分析快速应变的流动问题考虑了旋流对湍流的贡献采用了分析 Prandtl数的分析公式有效粘性系数公式不同对如下流动有较好模
15、拟结果:高流向曲率与应变率流动过渡态流动壁面传热传质问题第30页/共53页雷诺应力模型雷诺应力模型Reynolds Stress Model产生项产生项压力应变项压力应变项耗散项耗散项TurbulentDiffusion(modeled)(related to)(modeled)(computed)(equations written for steady,incompressible flow w/o body forces)雷诺应力输运方程雷诺应力输运方程.Pressure/velocity fluctuationsTurbulenttransport第31页/共53页雷诺应力模型（雷诺应
16、力模型（RSM，Reynolds Stress Model）RSM 封闭了雷诺平均的 Navier-Stokes方程，求解雷诺应力输运方程。方程严格推导得到。（瞬时量方程雷诺平均方程）脉动量需要求解耗散率输运方程，用以封闭方程组无需涡粘性各向同性假设雷诺应力输运方程有需要模拟量（项）.RSM 有高精度模拟复杂湍流流动的潜质.考虑了曲率,旋流,旋转和高应变率影响气旋流动,有旋燃烧室内流动问题旋转流动通道,二次流第32页/共53页大涡模拟（大涡模拟（LES，Large Eddy Simulation）大涡:决定了动量、能量和其它标量场的特性、对平均流动场起主导作用.各向异性,可能是湍流结构历史影响
17、，也许是流动场影响，大涡受到流畅、边界调节和流动参数影响。.小涡:趋于各向同性，受流动影响较小比大涡更易于模拟.LES 直接结算(求解)大涡，小涡用模型来模拟(Subgrid-Scale Modeling).计算量较大网格点,NLES 非定常计算第33页/共53页RANS 湍流模型比较湍流模型比较第34页/共53页近壁处理近壁处理 u靠近壁面区域、多数情况下 k-和RSM 湍流模型的计算结果都不够精确.u需要近壁处理.l标准壁面函数（Standard wall functions）l非平衡壁面函数（Nonequilibrium wall functions）l双层区模型（Two-layer z
18、onal model V5）l增强的壁面处理（Enhanced wall treatment，V6）Boundary layer structure第35页/共53页标准壁面函数（标准壁面函数（Standard Wall Functions）平均速度平均速度温度温度其中whereand P is a function of the fluid and turbulent Prandtl numbers.热边界层厚度第36页/共53页非平衡壁面函数非平衡壁面函数Nonequilibrium Wall FunctionsLog-law is sensitized to pressure gradi
19、ent for better prediction of adverse pressure gradient flows and separation.Relaxed local equilibrium assumptions for TKE in wall-neighboring cells.Thermal law-of-wall unchangedwhere第37页/共53页Two-Layer Zonal Model应用于低雷诺数流动或复杂流动中的近壁流动现象模拟。用点到壁面举例y为特征长度求当地雷诺数，区分区域（层）u对于高雷诺数流动，湍流核心区域用 k-模型来模拟.u双层区模型中，在粘
20、性影响的区域只求解k 的输运方程.u 通过计算得到（与特征长度相关）.u分区是动态过程，并根据计算结果调节分区区域属性.第38页/共53页近壁处理方法的比较近壁处理方法的比较第39页/共53页计算网格划分注意事项计算网格划分注意事项Wall Function ApproachTwo-Layer Zonal Model ApproachlFirst grid point in log-law regionlAt least ten points in the BL.lBetter to use stretched quad/hex cells for economy.lFirst grid po
21、int at y+1.lAt least ten grid points within buffer&sublayers.lBetter to use stretched quad/hex cells for economy.第40页/共53页估计第一个网格点的位置估计第一个网格点的位置根据经验与实验公式，估计表面摩擦系数（skin friction coefficient）:平板-管流-计算摩擦速度:反过来计算第一网格点到壁面的距离:壁面函数 Two-layer model用post-processing 检查近壁网格的空间分辨率y1=50/uy1=/u第41页/共53页边界条件设定边界条件
22、设定给出湍流 inlets&outlets(potential backflow)参数k-模型需要 k 和 Reynolds stress 模型需要 Rij 和 根据不同输入选择，给出自己熟悉或方便给出的量Turbulence intensity 和 length scale（湍流强度与长度尺度）长度尺度与包含主要能量的大涡旋结构尺度相关.对于边界层流动:l 0.4d99流体流过栅板/孔板:l opening sizeTurbulence intensity and hydraulic diameter（湍流强度与水力学直径）适合通道与管内流动Turbulence intensity and
23、turbulent viscosity ratio（湍流强度与粘性比）外流:直接给出 k 和(用户自定义函数，可以给出分布).第42页/共53页湍流模型用户界面湍流模型用户界面Define Models Viscous.湍流模型选项，湍流模型选项，Turbulence Model options近壁处理，Near Wall Treatments无粘，层流，湍流无粘，层流，湍流Inviscid,Laminar,or Turbulent其余湍流选项，Additional Turbulence options第43页/共53页Example:Channel Flow with Conjugate H
24、eat Transferadiabatic wallcold airV=50 fpmT=0 Fconstant temperature wall T=100 Finsulation1 ft1 ft10 ftPPredict the temperature at point P in the solid insulation第44页/共53页湍流模拟湍流模拟检查是否是湍流 ReDh=5,980低雷诺数下的发展湍流与壁面边界层会导致压力梯度，因此选择：RNG k-模型，用nonequilibrium wall functions处理近壁流动.网格划分简单几何形状 quadrilateral cel
25、ls估计在水平壁面的垂直方向上梯度比较大，近壁第一个网格要落在对数区（log-law region.）流向网格渐大，用于捕捉边界层增长过程.根据求解结果，调节网格，用于提高求解温度边界层（温度梯度）第45页/共53页速度等值线温度等值线上下边界层加速了核心区域的流体流动动量动量&热边界层预测结果热边界层预测结果温度边界层清晰可见P第46页/共53页举例举例:圆柱绕流圆柱绕流wallwall1 ft2 ft2 ft airV=4 fpsCompute drag coefficient of the cylinder5 ft14.5 ft第47页/共53页检查是否湍流 ReD=24,600钝体绕流
26、尾迹中，可以看到非定常的涡旋脱落现象。侧下游的分离点计算和靠近壁面的近壁影响，所以需要考虑用,RNG k-湍流模型，2-layer zonal model处理边界.网格划分简单几何形状&BLs quadrilateral cells.边界层网格，2-layer model fine mesh near surface&first cell at y+=1.湍流模拟步骤湍流模拟步骤第48页/共53页圆柱绕流的网格圆柱绕流的网格第49页/共53页湍流模拟，涡旋脱落（湍流模拟，涡旋脱落（Vortex Shedding）有效粘性系数等值线图 eff=+tCD=0.53 Strouhal Number=0.297其中第50页/共53页不同边界结果处理结果比较不同边界结果处理结果比较标准壁面函数非平衡壁面函数双层区壁面模型第51页/共53页小结小结:湍流模拟湍流模拟成功的湍流模拟需要对实际过程作出如下的正确判断:流体物理特性计算机资源项目要求精度时间可用的湍流模型与近壁处理方法可以先用 k-得到初步结果，不满意改用RNG 或Realizable k-.对于高旋流数流动，选择 RSM模型.除非是求解低雷诺数流动，或者关心复杂流动的近壁流体特性，否则壁面处理选择壁面函数.第52页/共53页感谢您的观看！第53页/共53页