湍流数值模拟.pptVIP

参考书 岑可法 樊建人《燃烧流体力学》 水利电力出版社 周力行 《湍流两相流动和燃烧的理论与数值模拟》科学出版社 范维澄 《流动及燃烧的模型与计算》中国科学技术大学出版社 湍流定性认识 固壁湍流（Wall turbulence）/自由湍流（Free turbulence）； 各向同性湍流（Isotropic turbulence）/异向湍流（heterogeneous turbulence）； 湍流精细结构认识的可行性 Kolmogorov尺度 一般为Re-3/4 例： 如：Re=105，反应特征尺度为0.2m； 则：lk=Re-3/4 =0.2（mm） 一个小涡团每个坐标方向上分5个网格（很粗糙！），涡团上为125个网格， 那么：为了研究该体系湍流精细结构，所需网格数为： 103× 103 × 103 ×125=1011 显然计算量非常巨大！ 混合长度的定义 对于自由剪切流 充分发展的湍流管流 大涡模拟 设想：在各种不同类型的流动中，大涡运动是不同的，小涡运动是类似的； 对大涡直接求解，小涡湍流模型求解 大涡定义 过滤速度/过滤函数 过滤速度： 其中：G 为过滤函数；△ 过滤宽度，通常为网格的2倍。 大涡模型的应用分析 ? LES在复杂流动的模拟中可以得到很多雷诺时均方法无法获得的紊流运动的细微结构和流动图像。数值实验证明雷诺时均方法在模拟复杂流动现象如涡脱落、浮力影响、流线弯曲、旋转和压缩运动时会遇到难以克服的困难; ? 与雷诺平均模型相比，大涡数值模拟的亚格子模型具有较大的普适性。湍流大涡数值模拟方法中需要封闭的量是亚格子应力，它和大尺度脉动的相关微弱。亚格子应力是不可解小尺度脉动和可解尺度之间的动量交换，它和强烈依赖于流动边界的大尺度脉动相关性很小，因此合理的亚格子模型将有较大的普适性; 湍流大涡数值模拟可以获得流动的动态特性，而雷诺平均模型只能提供定常的气动力特性。湍流大涡数值模拟的解包含大于过滤尺度的所有脉动，由此可以获得速度谱以及气动力谱等,这些动态气动力特性对于近代航天器设计是十分重要的。 模型计算结果比较 模型计算结果比较 模型计算结果比较 需注意的地方 需注意的地方 需注意的地方 重要提示 * * 计算网格 采用的是由Harlow和Welsh提出的交错网格，通常也称为Marker-And-Cell网格，或简称为MAC网格。 所谓交错网格，是指把速度u、v、w和压力p（包括其它所有标量场及物性参数）分别存储在四套不同的网格系统上，其中u、v和w的控制容积沿各自的方向分别与p的控制容积有半个步长的错位。 交错网格具有以下两个主要优点： 1．可以避免得到不合理的速度场和压力场，有利 于提高解的稳定性和精度。 2．对位于主控制容积节点的变量求解时，控制容 积表面上的几个主要流场变量无需进行插值计 算，从而可以提高流场的求解精度。 计算网格 根据Moin等的研究结果，当DNS中采用的网格尺度小于或者跟Kolmogorov微尺度处于同一个量级时，即可得到足够精确的结果。 对于流动雷诺数Re=5000，由各向同性假设可以估算得到Kolmogorov微尺度 ＝0.043d。故在模拟中，沿x方向安排均匀的交错网格为1.45，沿z方向安排均匀的交错网格为0.725。为了能捕捉到射流强剪切核心区域内的小涡结构，沿y方向在－4.5d y 4.5d的区域内安排均匀的交错网格为0.93，在其它区域安排由中心向上下两侧拉伸的网格。总的计算网格节点数沿x、y、z三个方向布置为400×300×128，共计为15.36×106个网格 。 Kolmogorow长度尺度 叠加在大尺度气流上有许许多多很小尺寸的涡，这些涡是大涡不断破裂所形成的。这些小涡使当地气流的流型迅速发生变化，湍流能量的耗散即发生在这些最小的结构中。 而分子粘性起着把小尺寸涡的动能耗散为热能的作用，这一最小的湍流流动结构尺寸可用Kolmogorow长度尺度 表示， 它决定于单位质量能量的耗散速率 和流体运动速度 。 为此，定义Kolmogorow长度尺度 为 Kolmogorow时间尺度 Kolmogorow时间尺度表示最小湍流结构的动量扩散时间，它的定义为 有限容积方法 分步投影方法 时间积分 直接数值模拟的时间步长Δt应该小于最小涡的时间尺度，积分时间又要大于最大涡的时间尺度。同时为了保证计算