剪切层与边界层组合流动的线性不稳定性分析

对不可压缩剪切层与边界层的组合流动完成了线性稳定性研究.组合流动的数学模型为Blasius边界层相似解与双曲正切函数的叠加,采用整体数值方法求解组合流动的稳定性方程,并验证了程序的准确性及网格无关性.研究给出组合流动的不稳定模态的辨识,即边界层模态和剪切层模态.在此基础上研究了剪切层对边界层模态不稳定性的影响以及壁面对剪切层模态的影响.由于剪切层的存在,使边界层模态中性曲线向左上方平移,临界雷诺数减小.此外,边界层模态不稳定性得到增强或抑制的影响,取决于扰动频率以及剪切层速度比的变化.组合流动中壁面边界层促使剪切层不稳定性减弱,主要表现在低频区域;而在高频区域,剪切层不稳定性几乎不受壁面边界层的影响.      

变循环发动机过渡态性能直接模拟方法

为了实现多变量可调的变循环发动机（VCE）过渡态性能模拟和控制规律设计,在稳态逆算法和过渡态性能隐式格式计算方法的基础上开发了变循环发动机过渡态直接模拟方法。在稳态逆算法模型中加入容积效应和转子动力学方程,实现了同时给定加速率、涡轮前温度和压缩部件工作点的条件下直接模拟计算过渡态过程中变几何参数和燃油流量的调节规律,验证了方法的精度和可行性,并将该方法用于变循环发动机转模态性能模拟和控制规律设计。计算结果表明,过渡态直接模拟方法的误差在0.58%以内,超声速巡航状态下由单外涵模态转换到双外涵模态的时间约为1 s,海平面静止状态下双外涵转单外涵的时间约为3 s,且推力、转速、涡轮前温度、喘振裕度等参数过渡平稳。该方法可简化转模态控制规律设计流程,并提高设计精度。     

高热流密度燃烧室壁面阵列空气射流-燃油组合冷却结构研究（爆震专刊）

针对高热流密度燃烧室壁面热防护需求提出了一种空气阵列射流冲击和燃油冷却肋板的集成冷却方式,在射流平均雷诺数(Re_j)为10000至30000、燃油进口流速(v_f)为2.33m/s至5.23m/s的范围内,采用数值模拟方法对其传热特性进行了研究,并基于壁面加热侧当量对流换热系数的概念,分析了基准肋板以及燃油冷却肋板的传热增强作用。与无肋板靶面的阵列射流冲击相比,带肋板阵列射流冲击的面积平均当量对流换热系数是前者的1.6倍,压力损失系数相对提高了约25%;采用燃油冷却肋板,加热壁面综合传热能力进一步增强,在Re_j=10000时,采用燃油冷却肋板的面积平均当量对流换热系数是基准肋板的1.5倍以上,即使在Re_j=30000时,燃油冷却肋板的传热增强比也可以达到1.2;燃油冷却肋板的出口温度相对进口温度的提升在20K~50K范围内,其提升幅度随着射流雷诺数或燃油进口流速的增大而减小。     

基于数值模拟的流场附面层边缘识别方法

附面层边缘通常取在速度达到主流速度0.99倍的位置,而复杂流场中主流流动往往并不均匀,给附面层边缘的准确识别造成了困难。为解决此问题,提出了用"参考主流"代替实际主流识别附面层边缘的方法:通过零剪切力滑移壁面边界条件下数值模拟得到不受附面层干扰的参考主流,在根据附面层定义确定附面层边缘时以该参考主流中的速度代替实际的主流速度。通过斜楔压缩和弯曲压缩两个超声速压缩流场对该识别方法进行了验证,所得到的斜楔压缩出口截面上附面层厚度与采用实际主流速度判断得到的厚度相对误差仅4.1%。根据该方法的识别结果对弯曲压缩型面设计进行附面层修正后,弯曲激波高度与无黏设计值之间的误差从修正前的2.0%降低至0.3%,压缩面末端压力的相对误差从修正前的6.6%降低至2.3%。该方法避免了指定主流速度的主观性,识别结果较为准确。     

叶根倒角模拟件设计

为了避免涡轮叶片叶根倒角低周疲劳开裂故障的发生,需借助叶根倒角特征模拟件对叶根低周疲劳强度进行考核。基于几何等效相似和载荷工况等效原则,设计了一种真实叶根倒角的特征模拟件。特征模拟件的榫头/叶身沿着周向投影宽度比例、缘板外侧与榫头外侧距离、缘板厚度、倒角半径等重要几何参数均与真实叶片一致。基于线弹性本构,采用Abaqus软件计算了特征模拟件在等效载荷工况下的应力分布。计算结果表明,特征模拟件的最大应力为187.6 MPa,出现在凸台过渡区倒角处,最危险点第一主应力方向为l₁=0.1141、m₁=0.9873、n₁=-0.1103,均与真实叶片对应部位的应力情况吻合,说明该叶根倒角特征模拟件设计合理,可用于考核真实涡轮叶片倒角的低周疲劳强度。     

基于大涡模拟的离心式喷嘴雾化过程研究

为实现离心式喷嘴雾化过程的精确数值仿真,探究喷嘴内部流动特性与外部液膜破碎形式,采用基于大涡模拟的仿真方法,对一种典型的四进口离心式喷嘴进行研究,仿真结果揭示了喷嘴内部相界面的振荡现象与外部液膜的破碎细节,并通过耦合流体体积法（VOF）与离散相模型（DPM）,获得液滴粒径的空间分布特征。研究结果表明：在液体填充过程中,喷嘴内的气液相界面存在波动与褶皱,形状并不稳定,内部的空气芯直径呈现正弦模式的振荡变化,喷嘴出口液膜厚度沿周向分布不均,这些因素导致出口附近的液膜表面出现扰动。在不同的进口条件下,不稳定性导致液膜表面上的扰动波形式不同。进口压力为0.3MPa时,液膜破碎由开尔文-亥姆霍兹（K-H）不稳定性产生的轴向正弦波所导致,产生沿周向分布的环形液带;在0.7MPa下,液膜表面开始出现由瑞利-泰勒（R-T）不稳定性引发的周向扰动波;随着压力增加至1.1MPa,液膜的破碎则由R-T不稳定性主导,产生沿轴向分布的液带结构,随后在气动力与表面张力的作用下破碎成液滴。二次雾化破碎后,喷嘴外部截面内的粒径呈“单谷”分布,液滴平均粒径计算结果与实验的最大相对误差为5.1%,与实验数据吻合度较高。     

正十烷简化机理的构建及其在发动机燃烧数值模拟中的应用

发动机燃烧室中燃料的能量释放与燃烧特性对于发动机设计具有重要作用,为了预测发动机点火包线和贫/富油极限等关键性能,迫切需要发展航空燃料及其典型组分的高精度化学动力学模型。本文针对燃料典型组分正十烷,采用自主开发的机理生成程序ReaxGen构建了其燃烧详细机理（1499种组分、5713步反应）。为了验证机理的合理性与可靠性,在当量比Φ=0.5-2.0,压力P=1-80 atm的宽工况条件下进行了点火延迟模拟验证,结果表明本文提出的正十烷详细机理在较宽的温度、压力和当量比条件下具有较高的模拟精度。为获得适用于发动机燃烧模拟的高精度简化机理,本文基于误差传播的直接关系图方法简化了正十烷燃烧详细机理,得到包含709种组分、2793步反应的正十烷半详细机理。进一步在高温范围（1000-1500 K）,采用路径通量分析方法简化得到含77种组分、359个反应的骨架机理。获得的骨架机理能够合理描述正十烷在高温下的燃烧特性,且该骨架机理尺度规模可用于基于火焰面模型的燃烧数值模拟。基于此高精度的骨架机理模型,结合火焰面生成流形湍流燃烧模型,采用大涡模拟方法进行了航空发动机环形燃烧室单头部扇形的燃烧模拟,初步获得了非稳态流场结构,其中温度模拟结果与实验值基本符合。     

航空发动机动态姿态仿真测试台架设计

为了更真实全面地仿真航空发动机的使用状态,研制了一种用于小型航空发动机整机动态姿态仿真测试的台架。该台 架具有横滚和俯仰2个自由度,可提供位置、速率及摇摆等功能,采用伺服电机带动齿轮组件和蜗轮蜗杆进行驱动,并具备排气装 置。设计结果表明：动态姿态模拟转台速度可达30°/s,加速度可达50°/s2,运动范围为±175°,定位精度实测优于±0.12°;可以任何 位置为中心做正弦、梯形、三角波等形式的运动。该设备的研制可进一步提高地面试验的仿真精度,有利于发动机检测技术的发展。     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

低压天然气喷射的船用二冲程低速双燃料发动机爆震研究

为了解决大缸径船用预混天然气双燃料发动机的爆震问题,并拓展天然气的稀燃边界。基于三维数值模拟的方法对大缸径船机进行了仿真模拟。分析了大缸径双燃料发动机爆震的特点,并对缸内涡流强度和废气再循环（EGR,Exhaust Gas Re-circulation）率对爆震的影响进行了研究。研究结果表明：大缸径预混天然气发动机的爆震位置往往发生在气缸边缘,火焰面的传播过程是引起缸内爆震的主要因素。随着缸内涡流从无到有的增强,缸内的爆震强度随之增强;当涡流到达一定程度后,随着涡流的增强,缸内的爆震强度反而降低;缸内加入EGR可以提高天然气当量比的同时减少爆震的强度,可以拓展天然气的稀燃边界。