• 上传我的文档

fluent教程-第三章-湍流模拟解析.ppt

收藏 版权申诉 举报 下载
fluent教程-第三章-湍流模拟解析.ppt_第1页
第1页 / 共53页
fluent教程-第三章-湍流模拟解析.ppt_第2页
第2页 / 共53页
fluent教程-第三章-湍流模拟解析.ppt_第3页
第3页 / 共53页
fluent教程-第三章-湍流模拟解析.ppt_第4页
第4页 / 共53页
fluent教程-第三章-湍流模拟解析.ppt_第5页
第5页 / 共53页
资源描述:

《fluent教程-第三章-湍流模拟解析.ppt》由会员分享,可在线阅读,更多相关《fluent教程-第三章-湍流模拟解析.ppt(53页珍藏版)》请在人人文库网上搜索。

1、湍流模拟,湍流非定常,非周期性的三维速度脉动、 强化物质、动量和能量的输运. 瞬时速度分解为平均速度和脉动速度 Uit Ui uit 压力、温度、组分浓度值具有类似的脉动,什么是湍流,,,,,Time,,,U i t,Ui,uit,平均量输运方程,雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度,有 其中, 和 分别是平均速度和脉动速度(i1,2,3) 类似地,对于压力等其它标量,我们也有 把上面的表达式代入瞬时的连续与动量方程,并取平均 (去掉平均速度 上的横线),我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式,,,,,,平均量输。

2、运方程(续),,,,上面两个方程称为雷诺平均的Navier-Stokes(RANS)方程。,,如果要求解该方程,必须模拟该项以封闭方程。,湍流模拟的方法,直接数值模拟(DNS)只适合低雷诺数流动。 求解雷诺平均的 Navier-Stokes RANS 方程 其中 雷诺应力 时间平均的湍流脉动用实验常数和流场平均速度信息来模拟. 大涡模拟(Large Eddy Simulation,LES)数值求解大涡,小涡用模型.,定常, 不可压缩流动 有/无 体积力,湍流输运模型,涡旋粘性系数模型 涡旋粘性系数模型是根据湍流应力和分子运动引起的粘性应力相似而提出来的 雷诺应力模型 雷诺应力模型则是从基本方程。

3、出发,直接推出雷诺应力的方程,但在雷诺应力的方程中包含有更高阶的相关项,对这些更高阶的相关项再建立相应的模型。,湍流应力,湍流应力 在三维空间中,下标和分别可取为1,2和3,所以湍流应力有9个分量组成,是一个二阶张量,若用矩阵形式表示可写为 主对角线上的三个分量,和称为湍流正应力,其余的六个分量称为湍流切应力。三个湍流正应力之和是湍流脉动动能的两倍,,,,湍流应力,若以主对角线上的三个分量作为对称轴,则对称的两个切应力分量是相等的,很显然,这是一个对称的二阶张量。很容易可以证明,在各向同性湍流中,湍流正应力的三个分量相等,即,,判断湍流,外流,内流,自然对流,沿表面,绕流,其中,其中,其它因素。

4、,如自由流湍流度, 表面条件, 扰动 可能导致流动从层流向湍流转捩,L x, D, Dh, etc.,,,,,湍流特点,额外应变率 流向曲率 测向分离 加速或减速 有旋 回流 或分离 二次流 3D振荡流动 Transpiration 吹风/吸气 自由湍流 剪切层相互作用,需要作出选择,湍流模型 和 近壁处理,物理流体,精度要求,计算资源,计算时间要求,计算网格,,,,,,,,,,Zero-Equation Models One-Equation Models Spalart-Allmaras Two-Equation Models Standard k-e RNG k-e Realizabl。

5、e k-e Reynolds-Stress Model Large-Eddy Simulation Direct Numerical Simulation,湍流模型,,,Include More Physics,每次迭代 增加 计算时间,,,,,FLUENT 5 的湍流模型,,基于雷诺平均 (RANS)的模型,,,Fluent 6,涡旋粘性系数模型,涡旋粘性系数模型是Boussinesq于1877年最早提出来的。他假定在近似平行的剪切流中,湍流应力张量中的切应力分量和平均速度在横向方向的梯度成正比,其比例系数称为涡旋粘性系数,以后把这个假定再推广到三维的流动。即湍流应力张量和平均流场应变率之间。

6、有关系 仿照分子运动引起的粘性系数,我们可以假定涡旋粘性系数和湍流中含能涡旋的特征长度和速度尺度成正比,即 而根据确定速度尺度V和长度尺度L方法的不同,又可以细分为各种不同的模型。而其中最简单,V和L的确定方法又是一致的模型应该是双方程模型。,,,,涡旋粘性系数模型,在双方程模型中,假定特征速度V和湍流动能k的平方根成正比,特征长度由湍流动能k和另外一个辅助的量确定。比如在k-e模型中,辅助的量选为湍流动能的耗散率,根据量纲分析得长度尺度,,,,,模型,,涡旋粘性系数模型,辅助的量是湍流频率 , 根据量纲分析得长度尺度,,模型,,,,,,推导,我们都假定流体不可压,即,Favre平均和雷诺平均。

7、方程完全相同。瞬时的动量方程减去平均的动量方程得脉动速度的方程,,,乘以,并求平均,利用湍流动能的定义,,,上式右端第二项可以重新整理成,,这里湍流动能的方程,方程右端各项依次为输运项,产生项和耗散项,各项对,求导数,乘以,,并求平均,可得耗散率,的方程,,目前采用的标准ke模型方程为,,,如果用湍流频率,代替湍流动能耗散率,,频率,的模型方程为,RANS 方程需要对雷诺应力进行封闭. 对于单方程模型,从修正的粘性系数输运方程简介求解湍流粘性系数。 对于双方程模型, 湍流粘性系数根据湍动能 TKE和耗散率 TKE来确定. Transport equations for turbulent ki。

8、netic energy and dissipation rate are solved so that turbulent viscosity can be computed for RANS equations.,Reynolds Stress Terms in RANS-based Models,Turbulent Kinetic Energy,Dissipation Rate of Turbulent Kinetic Energy,湍流粘性系数,Boussinesq Hypothesis isotropic stresses,湍流模型,0方程模型 单方程模型 双方程模型 雷诺应力模型 。

9、大涡模拟,湍流粘性系数 求解 的输运方程 The additional variables are functions of the modified turbulent viscosity and velocity gradients.,单方程模型 Spalart-Allmaras,,,,Generation,Diffusion,Destruction,单方程模型 Spalart-Allmaras,用于计算航空中考虑边界作用的高速流动问题 可以考虑边界层逆压梯度 叶轮机械流动 壁面网格可粗可细 网格细了,可以考虑低雷诺数流动,考虑边界层影响. 网格粗时候,可以得到相对较好的流动结果.,不能预。

10、测均匀各向同性湍流的耗散。并且,单方程模型没有考虑长度尺度的变化,这对一些流动尺度变换比较大的流动问题不太适合。比如,平板射流问题,从有壁面影响流动突然变化到自由剪切流,流场尺度变化明显。,应用举例,无粘流结果 单方程湍流模型结果,无粘流结果 单方程湍流模型结果,双方程模型 标准 k- 模型,湍动能方程,耗散率输运方程,实验常数,equations written for steady, incompressible flow w/o body forces,Convection,Generation,Diffusion,Destruction,双方程模型 标准 k- 模型,“Baseline。

11、 model” Two-equation 最广泛用于工程流动与换热计算 结果比较透彻,优点、缺点明显 半经验性 k 方程经过严格推导得到 方程通过物理推理得到 只对完全湍流适用大雷诺数流动 对于多数的湍流流动问题有较为合理的结果 一般工程流动问题 换热问题,双方程模型 Realizable k-,与标准 k-模型的区别 湍流粘性系数 是变量, 其中 A0, As, 和 U* 是速度梯度的函数 雷诺正应力为正; Schwarz 不等式成立 新的耗散率 输运方程,,,Generation,Diffusion,Destruction,,,Buoyancy,与标准的 k-模型具有相同的湍动能输运方程 。

12、在如下流动中,比标准的双方程模型更具有优越性 平板射流、圆射流 边界层有强的逆压力梯度,分离 有旋(rotation),回流( recirculation) 强曲率影响,双方程模型 Realizable k-,双方程模型 RNG k-,湍动能方程,耗散率输运方程,其中,用RNG理论推导得到,equations written for steady, incompressible flow w/o body forces,双方程模型 RNG k-,k- 两个方程都从瞬时NS方程出发,采用重整化群方法 Renormalization Group 严格推导得到. 与标准的 k- 方程有类似的形式、。

13、但 在 方程中多出一项用于分析快速应变的流动问题 考虑了旋流对湍流的贡献 采用了分析 Prandtl数的分析公式 有效粘性系数公式不同 对如下流动有较好模拟结果 高流向曲率与应变率流动 过渡态流动 壁面传热传质问题,雷诺应力模型Reynolds Stress Model,产生项,压力应变项,耗散项,Turbulent Diffusion,modeled,related to e,modeled,computed,equations written for steady, incompressible flow w/o body forces,雷诺应力输运方程.,Pressure/velocit。

14、y fluctuations,Turbulent transport,,,雷诺应力模型(RSM,Reynolds Stress Model),RSM 封闭了雷诺平均的 Navier-Stokes方程,求解雷诺应力输运方程。 方程严格推导得到。(瞬时量方程雷诺平均方程)脉动量 需要求解耗散率输运方程,用以封闭方程组 无需涡粘性各向同性假设 雷诺应力输运方程有需要模拟量(项). RSM 有高精度模拟复杂湍流流动的潜质. 考虑了曲率, 旋流, 旋转和高应变率影响 气旋流动, 有旋燃烧室内流动问题 旋转流动通道, 二次流,大涡模拟(LES,Large Eddy Simulation),大涡 决定了动量。

15、、能量和其它标量场的特性、对平均流动场起主导作用. 各向异性, 可能是湍流结构历史影响,也许是流动场影响,大涡受到流畅、边界调节和流动参数影响。. 小涡 趋于各向同性,受流动影响较小 比大涡更易于模拟. LES 直接结算 求解 大涡 ,小涡用模型来模拟 Subgrid-Scale Modeling. 计算量较大 网格点, NLES 非定常计算,RANS 湍流模型比较,近壁处理,靠近壁面区域、多数情况下 k- 和RSM 湍流模型的计算结果都不够精确. 需要近壁处理. 标准壁面函数(Standard wall functions) 非平衡壁面函数(Nonequilibrium wall funct。

16、ions) 双层区模型(Two-layer zonal model V5) 增强的壁面处理(Enhanced wall treatment,V6),Boundary layer structure,标准壁面函数(Standard Wall Functions),平均速度,温度,其中,where,and P is a function of the fluid and turbulent Prandtl numbers.,热边界层厚度,,非平衡壁面函数Nonequilibrium Wall Functions,Log-law is sensitized to pressure gradient 。

17、for better prediction of adverse pressure gradient flows and separation. Relaxed local equilibrium assumptions for TKE in wall-neighboring cells. Thermal law-of-wall unchanged,where,Two-Layer Zonal Model,应用于低雷诺数流动或复杂流动中的近壁流动现象模拟。 用点到壁面举例y为特征长度求当地雷诺数,区分区域(层),对于高雷诺数流动,湍流核心区域用 k-e 模型来模拟. 双层区模型中,在粘性影响的区。

18、域只求解k 的输运方程. 通过计算得到(与特征长度相关). 分区是动态过程,并根据计算结果调节分区区域属性.,近壁处理方法的比较,计算网格划分注意事项,Wall Function Approach,Two-Layer Zonal Model Approach,First grid point in log-law region At least ten points in the BL. Better to use stretched quad/hex cells for economy.,First grid point at y 1. At least ten grid points wi。

19、thin buffer sublayers. Better to use stretched quad/hex cells for economy.,估计第一个网格点的位置,根据经验与实验公式,估计表面摩擦系数(skin friction coefficient) 平板- 管流- 计算摩擦速度 反过来计算第一网格点到壁面的距离 壁面函数 Two-layer model 用post-processing 检查近壁网格的空间分辨率,y1 50/u,y1 / u,边界条件设定,给出湍流 inlets outlets potential backflow参数 k- 模型需要 k 和 Reyno。

20、lds stress 模型需要 Rij 和 根据不同输入选择,给出自己熟悉或方便给出的量 Turbulence intensity 和 length scale (湍流强度与长度尺度) 长度尺度与包含主要能量的大涡旋结构尺度相关. 对于边界层流动 l 0.4d99 流体流过栅板/孔板l opening size Turbulence intensity and hydraulic diameter(湍流强度与水力学直径) 适合通道与管内流动 Turbulence intensity and turbulent viscosity ratio(湍流强度与粘性比) 外流 直接给出 k 和 用户自定。

21、义函数,可以给出分布.,湍流模型用户界面,Define Models Viscous...,湍流模型选项,Turbulence Model options,近壁处理,Near Wall Treatments,无粘,层流,湍流Inviscid, Laminar, or Turbulent,,,,其余湍流选项,Additional Turbulence options,,Example Channel Flow with Conjugate Heat Transfer,,,,adiabatic wall,cold air V 50 fpm T 0 F,constant temperature 。

22、wall T 100 F,insulation,1 ft,1 ft,10 ft,,,,,,,,P,Predict the temperature at point P in the solid insulation,,,湍流模拟,检查是否是湍流 ReDh 5,980 低雷诺数下的发展湍流与壁面边界层会导致压力梯度,因此选择 RNG k- 模型,用nonequilibrium wall functions处理近壁流动. 网格划分 简单几何形状 quadrilateral cells 估计在水平壁面的垂直方向上梯度比较大,近壁第一个网格要落在对数区( log-law region.) 流向网格渐。

23、大,用于捕捉边界层增长过程. 根据求解结果,调节网格,用于提高求解温度边界层(温度梯度),速度等值线,温度等值线,上下边界层加速了核心区域的流体流动,动量 热边界层预测结果,温度边界层清晰可见,,,,,P,举例 圆柱绕流,,,,wall,wall,,,,1 ft,2 ft,2 ft,,,,,air V 4 fps,Compute drag coefficient of the cylinder,,,,5 ft,14.5 ft,,,,,检查是否湍流 ReD 24,600 钝体绕流尾迹中,可以看到非定常的涡旋脱落现象。侧下游的分离点计算和靠近壁面的近壁影响,所以需要考虑用, RNG k- 湍。

24、流模型, 2-layer zonal model处理边界. 网格划分 简单几何形状 BLs quadrilateral cells. 边界层网格, 2-layer model fine mesh near surface first cell at y 1.,湍流模拟步骤,圆柱绕流的网格,湍流模拟,涡旋脱落(Vortex Shedding),有效粘性系数等值线图 eff t,CD 0.53 Strouhal Number 0.297,其中,不同边界结果处理结果比较,标准壁面函数,非平衡壁面函数,双层区壁面模型,小结 湍流模拟,成功的湍流模拟需要对实际过程作出如下的正确判断 流体物理特性 计算机资源 项目要求 精度 时间 可用的湍流模型与近壁处理方法 可以先用 k- 得到初步结果,不满意改用RNG 或Realizable k-. 对于高旋流数流动,选择 RSM模型. 除非是求解低雷诺数流动,或者关心复杂流动的近壁流体特性,否则壁面处理选择壁面函数.,。

展开阅读全文
我图网
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

网站客服QQ:2846424093   微信号:renrenwenkuwang   人人文档上传用户QQ群:460291265   

copyright@ 2020-2023  renrendoc.com 人人文库版权所有   联系电话:0512-65154990

备案号:苏ICP备12009002号-5  经营许可证:苏B2-20200052  苏公网安备:32050602011097号

           

本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知人人文库网,我们立即给予删除!