支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。     

预混燃烧大涡模拟和燃烧模型的检验

用代数二阶矩亚网格(ASOM-SGS)燃烧模型对文献中测量的钝体后方丙烷-空气预混燃烧进行了大涡模拟,模拟统计的时平均速度、速度脉动均方根值和温度分布与实验数据结果吻合很好,表明所采用的ASOM-SGS亚网格燃烧模型是合理的。模拟的瞬态结果显现了钝体后方湍流流动和火焰结构。将大涡模拟数据统计得到的反应率系数-浓度关联量的分布规律,与代数二阶矩RANS(ASOM-RANS)燃烧模型的模拟值进行对比,结果发现,大涡模拟统计值和ASOM-RANS模型的模拟值很接近,从而证明了湍流燃烧代数二阶矩RANS模型的合理性。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

不同湍流预混燃烧模型在本生灯火焰中的比较

湍流预混燃烧的大涡模拟数值计算近来受到很大的关注,为了进一步发展湍流预混燃烧模型,利用正在发展的湍流预混燃烧模型等值面集合的G函数方法、层流预混火焰生成流型的化学热力学建表FGM方法、G函数/FGM耦合方法三种燃烧模型,对甲烷/空气湍流本生灯火焰进行了大涡模拟计算,并对上述燃烧模型进行了比较。大涡模拟结果表明,G函数方法能够较为准确捕捉火焰的传播特性,FGM方法能够捕捉到火焰和湍流的相互作用,但无法准确预测火焰的传播特性。而耦合的方法既可以捕捉火焰的传播特性,又能够捕捉到火焰和湍流的相互作用。     

同轴双剪切气-气喷嘴数值模拟

通过求解k-ε湍流模型的Navier-Stokes方程组,对以气氢/气氧为推进剂的同轴双剪切喷嘴燃烧室流场进行数值模拟研究.结果表明:气-气推进剂在燃烧室内形成两个燃烧面,能在大流量下实现较高的燃烧效率;同轴双剪切喷嘴的氧喷注速度和氢氧速度比是同轴双剪切喷嘴设计的关键参数,氧喷注速度越小使推进剂的燃烧位置越远离喷注面,而氢氧速度比越大使燃烧位置越靠前;与气-液同轴剪切喷嘴相比,同轴双剪切气-气喷嘴的设计可以取较大的氧喷注速度和适合的氢氧速度比.      

大涡模拟模型燃烧室燃烧性能计算

对带双级扩压器的模型燃烧室气液两相瞬态喷雾燃烧过程,在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法进行数值研究,同时采用多维经验分析法预估燃烧性能.采用 k 方程亚网格尺度模型模拟亚网格湍流黏性;亚网格EBU(eddy-break-up)燃烧模型预估化学反应速率;多维经验分析法计算燃烧性能;并在非交错网格体系下气相采用SIMPLE(semi-implicit method for pressure-linked equations)算法对控制方程进行求解,液相采用随机离散模型,两相之间的耦合采用PSIC(particle-source-in-cell)算法.通过大涡模拟瞬态及时均计算结果表明:与粒子图像测速仪(PIV)测量的瞬态速度场、出口温度分布试验数据吻合,表明在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法,数值模拟模型燃烧室两相喷雾燃烧流场,所采用的亚网格模型可以用于燃烧室气液两相喷雾燃烧流场的大涡模拟;燃烧性能计算结果与试验测量结果基本一致,说明所采用多维经验分析法可以用来数值模拟航空发动机燃烧室燃烧性能的计算,特别是污染物的预估,为设计低污染高性能航空发动机燃烧室提供有用的设计依据.      

基于超大涡模拟的燃烧室气动性能仿真研究进展

航空发动机燃烧室涉及旋流、雾化蒸发、掺混、化学反应、湍流与火焰相互作用等多尺度强耦合物理化学过程,相关的高 精度建模和数值模拟面临极大的挑战。超大涡模拟是近些年发展的兼顾计算精度、计算效率和强鲁棒性的数值模拟新方法,具备 试验室尺度和复杂工程应用场景下湍流流动与燃烧仿真能力。针对航空发动机燃烧室相关流动与燃烧基本特征,阐述了超大涡 模拟的理论方法及特点,从旋流流动、湍流燃烧、液雾雾化、碳烟生成、燃烧不稳定等典型多物理过程,以及双旋流模型燃烧室和高 温升燃烧室气动性能集成仿真等方面介绍了超大涡模拟的研究进展,对涉及的物理机制进行了分析,为超大涡模拟在航空发动机 燃烧室中规模化工程应用提供了坚实支撑。超大涡模拟在较低的计算资源消耗下具备与传统大涡模拟相当的计算精度,是一种 经济可承受的燃烧室高精度气动性能仿真新方法。     

湍流燃烧模型在燃烧室数值模拟中的对比分析

以某单头部矩形燃烧室为研究对象,采用多种湍流模型和燃烧模型进行组合计算,模拟燃烧室内部的速度场和温度场,并对计算结果进行对比分析。结果表明:湍流模型主要影响火焰筒内部主燃孔横截面上游的速度分布,Standard k-ε和Realizable k-ε模型的速度场计算结果差异相对较小;湍流动能预测受湍流模型的影响较大,并具有一定规律性;不同组合模型对燃烧室内部和出口温度分布的局部细节模拟差异较大,燃烧模型影响最大;PDF模型计算的温度值较合理,另外三种燃烧模型在单步完全反应燃烧机理下的计算值偏高,计算获得的OTDF也相差较大。     

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

旋流流动对水冲压发动机性能影响

基于某药柱式水冲压发动机,建立粉末式水冲压发动机模型.利用Fluent软件模拟研究旋流流动对水冲压发动机流动与燃烧反应性能的影响.选择k ε双方程模型作为湍流模型,有限速率/涡耗散模型作为燃烧模型,基于密度隐式耦合求解进行计算.结果表明:添加旋流叶片并增大旋流角度,燃烧室最高燃烧温度、平均燃烧温度、出口速度和比冲均增大.研究证明添加旋流叶片对水冲压发动机流动与燃烧反应性能有很大改进,且旋流角度越大,效果越明显.