不同亚格子模式在后台阶湍流流动大涡模拟中的应用

本文用大涡模拟方法研究了湍流后台阶流场中的大涡演变过程,并在此基础上研究了目前大涡模拟中比较常用的六种亚格子模式。在相同的流动几何参数以及计算条件下,给出了不同亚格子模式下湍流流动瞬时压力场以及流场中瞬时粘性大小的分布,尤其给出了流场瞬时演变的大尺度涡结构。将六种亚格子模式从回流区长度、计算时间、计算结果的准确性和流场稳定性等不同角度进行了比较。为选取合适的亚格子模式深入研究湍流耗散机理奠定了基础。     

基于有限反应速率的扩散燃烧大涡模拟

本文发展了一种基于有限反应速率的湍流燃烧大涡模拟方法。基于此方法,对常压下,甲烷/空气平板射流扩散燃烧进行了大涡模拟.甲烷/氧气反应采用包含七种组分的四步反应简化反应机理,考虑了基于Arrienius定律的有限化学反应速率,瞬态结果和时均结果与DNS结果都符合很好。为了解决考虑有限反应速率后计算量急剧增加的问题,本文引入等效化学反应速率模型,在压力泊松方程的求解中采用了求解效率较高的快速傅里叶变换方法,并对整个程序进行基于区域分解的MPI并行编程处理。     

楼梯井内火灾过程的大涡模拟

针对建筑物楼梯井内的火灾过程,采用大涡模拟方法,探讨了不同释热速率对楼梯井内火灾过程中空气卷吸量、楼梯井内平均氧气浓度、平均温度以及楼梯井内外最大压差的影响,分析了火灾烟气流动的动力学特征。研究结果表明,随着释热速率的增加,楼梯井空气卷吸量增加,但由于氧气的消耗量也在增加,导致楼梯井内平均氧气浓度降低。此外,释热速率越大,则楼梯井内温度越高,导致楼梯井内外最大压差增加,增强楼梯井的抽风拔为效应,加快了火灾烟气的蔓延。     

空间发展平板混和层流动的大涡模拟-3D与2D模拟的比较

本文利用大涡模拟方法研究了空间发展平板混合层流动,分别在三维模型和二维模型下,比较了大尺度结构涡的演化过程,以及流向速度、流向脉动速度、横向脉动速度和Reynolds应力的统计时均结果,并与实验结果对照,指出3D模拟无论在物理真实性还是预报结果准确性方面都要明显优于2D模拟。     

入口扰动频率对矩形射流影响的大涡模拟

为研究入口扰动对矩形射流的影响,对不同扰动频率下的三维矩形射流进行了大涡模拟,用分布投影法求解动量方程,对压力泊松方程采用FACR直接解法.计算结果预报的速度在横向的单峰分布和速度半宽沿流向的近似线性增长与实验结果吻合良好.计算结果表明扰动频率为射流的本征频率fe=0.22时射流的卷吸作用最强,而高频扰动则抑制卷吸作用;激励频率接近本征频率时射流提前出现速度在展向的马鞍型双峰分布,更早趋于圆射流形态.进一步的流场二维和三维瞬态涡结构图谱表明:上述现象与f=0.1和f=0.22时射流近场更频繁卷起、更丰富的大涡结构存在密切关系,正是这些大涡运动不断地在横向卷吸入外界流体,实现射流在横向的迅速扩展.     

后台阶流动再附着过程的大涡模拟研究

应用自主开发的大涡模拟程序数值模拟研究了后台阶流动中再附着过程的演变。在流动几何参数不变情况下，给出了再附长度随雷诺数的变化规律，并与实验进行了比较，二者相符得比较好。在此基础上，给出了三种典型雷诺数下，后台阶流动的回流区特征。在湍流情况下，研究了突扩比对再附长度的影响，与实验结果吻合的比较好。详细讨论了湍流情况下大涡拟序结构的瞬时再附着过程。 这些研究结果对具有再附着现象的流动结构的工程应用具有指导意义。     

后台阶分离流动中大涡结构演变的数值模拟

本文对后台阶分离流动中涡结构的演变进行了大涡模拟，研究了流场结构的变化规律。详细讨论了随着雷诺数的增加流场结构的典型特征的变化规律，指出流场中的涡结构随着雷诺数的增大变得十分复杂和丰富，回流区的数目、大小及其出现的位置也显著地不同。这些结果与已有的一些实验值和流场显示结果是吻合的。在此基础上，进一步研究了高雷诺数时流场中大尺度涡结构的瞬时发展和演化过程，展示了其中大涡的产生、追随、吸引、合并和破碎等过程。对于高雷诺数情况，对大涡模拟得到的数值结果进行了统计，得到的时均速度分布以及台阶后方的回流区长度与现有的其他实验结果符合得很好。本研究是针对后台阶分离流动深入开展湍流控制以及两相流动研究的基础。     

三维煤粉燃烧全双流体模型的数值模拟

1引言目前煤粉燃烧综合模型最为广泛的是以轨道模型为基础，美国杨伯翰大学的先进燃烧工程中心（CERC）自1980年起研制二维煤粉燃烧程序PCGC－2，从1990年起研制三维煤粉燃烧的PCGC－3程序[1]。与轨道模型发展的同时，美国Rabcock＆Wilcox公司Fiveland[2]等人研制了FURMO程序，用无滑移模型对560MW侧墙喷燃煤粉炉进行了三维全模拟。其特点是首次用全欧拉的处理方法计算三维煤粉燃烧过程，其不足之处是不考虑气粒两相间的速度滑移和温度滑移．总的看来，用轨道模型模拟煤粉燃烧，易于考虑颗粒反应经历，也可给出两相之间的速度及温…     

混合层流动的格子涡方法数值模拟

采用格子涡方法对二维平板混合层流动进行了数值模拟,重点考察了格子涡方法定量模拟混合层流动流场下游不同流向位置各湍流统计平均量的能力。模拟结果表明:除横向脉动速度均方根比实验值大外,混合层流动在下游不同流向位置的流向平均速度,流向脉动速度均方根和雷诺应力的模拟结果都与实验结果定量符合得很好。此外,模拟结果还很好地反应出混合层流动的自相似性质。从而表明,格子涡方法具有对空间发展混合层流动进行较精确的定量模拟的能力。     

混合分数PDF及其简化模型的大涡模拟研究

通过大涡模拟得到三维等密度平面射流中的混合分数的动态演化过程并统计其概率密度(PDF)的分布,对非预混燃烧雷诺时均模拟中常用的简化模型,包括β函数和截尾高斯函数,进行了定性和定量分析。混合分数的PDF分布在流场的不同区域上定性存在较大的区别,可以分为四类:追随型单峰分布、非追随型单峰分布、平板型分布和双峰分布。单峰追随型主要分布在湍流充分发展区,对不同的简化模型都有较好的适应性;单峰非追随型PDF通过β函数模拟可以得到很好的结果;平板型分布PDF,定性上难以模拟但定量模拟效果好;双峰分布主要分布于湍流外问歇性较大的区域,现有简化PDF模型难以得到满意的结果。