高雷诺数湍流非预混火焰及NO生成的大涡模拟

采用两种亚网格湍流燃烧模型,即化学建表方法结合假定概率密度模型和稀疏拉格朗日过滤密度函数方法,对高雷诺数湍流非预混火焰Flame D进行数值研究,定量比较不同亚网格模型的差异,并对火焰特征、污染生物生成特性进行分析。结果表明,两类亚网格燃烧模型预测的温度及大组分分布相近,稀疏拉格朗日过滤密度函数方法可以更好地模拟CO质量分数分布。不同的假定概率分布均可合理描述湍流与火焰的相互作用,之间的差别主要体现在NO分布,Dirac函数远高估了NO生成,而Top-hat函数则略低估了NO生成,Beta函数表现最优。Flame D的高温区及NO质量分数均主要分布在当量混合线及富燃侧附近。受高温伴流的影响,NO质量分数与温度一直保持高度正相关,峰值主要集中在标量耗散率很小的区域。不同截面上,反应物中的O_2和生成物中的H_2O均与NO高度相关。     

转静盘腔流动瞬态响应特性研究

以航空发动机转静盘腔内非稳态流动为研究对象,基于用户自定义函数编程设置进口压力按斜坡函数变化时的非稳态边界条件,利用数值模拟方法,在滞后效应、流场演化和物性参数方面,研究了不同无量纲跃升幅值下转静盘腔流动瞬态响应机理。结果表明:无量纲跃升幅值由1.05增至1.2,盘腔内部的当地响应时间延长了81.1%,相对总温的超调量增加了157.0%,其峰值时间增加了31.3%,且进口压力波动的传播方向总是由进口指向出口;瞬态过程中,转静盘腔的容积效应和旋转效应共同作用使出口流量明显滞后进口流量,且其滞留质量流量与相对总温的超调现象密切相关;冲击射流、旋转贴壁射流和静子壁面侧的涡系构成了盘腔流场演化的主体,无量纲跃升幅值越大,流场演化过程越剧烈,响应时间越长;声速和热扩散率与盘腔响应速度密切相关。     

航空发动机中冰晶结冰的研究进展

一般认为经过风扇压缩后空气温度高于冰点，压气机内不会发生结冰。然而近年来，国外研究者对若干起发动机推力损失故障的研究表明冰晶能够引起压气机内部结冰，并在试验室中再现了温度高于冰点时压气机内结冰这一现象。我国在航空发动机冰晶结冰研究方面刚刚起步，为准确掌握国内外研究现状，从数值和试验两方面对国内外已开展的压气机内结冰研究进行了分析与总结，讨论了当前冰晶结冰的主要研究成果和存在的局限，提出了开展冰晶结冰研究需重点关注的方向，为我国航空发动机冰晶摄入结冰研究和适航审定研究提供一定参考。     

高焓电弧风洞试验热化学非平衡流场数值模拟

针对高焓电弧风洞内部流动的热化学非平衡效应及气体组分和振动能量冻结效应导致的试验数据外推困难问题，基于高焓风洞喷管/试验段/试验模型一体化数值模拟的思路，通过数值求解三维热化学非平衡Navier-Stokes方程，开展了FD-15高焓电弧风洞典型运行状态下流场的数值模拟，与典型试验状态的气动热数据进行了对比验证，研究了试验数据外推飞行条件的方法及有效性问题，分析了提高驻室总压对试验数据外推的影响。研究表明:（1）风洞试验段来流离解度高，热化学非平衡效应及其冻结现象严重；（2）热流校核试验测量数据位于一体化数值模拟的完全催化热流和非催化热流之间，分布合理，验证了计算方法和程序的正确性；（3）试验模型安放位置对模型表面压力和热流存在影响，模型与喷管出口的距离越大，模型表面压力和热流越低；（4）当驻室总压较低时，通过双尺度模拟准则（模拟飞行条件总焓和双尺度参数ρ_∞L）外推热流失效，使用部分模拟准则（模拟飞行条件总焓和驻点压力）外推热流也会出现较大差异，在非催化条件下这一现象更加明显；（5）当驻室总压较高时，使用双尺度模拟准则或部分模拟准则外推飞行条件，产生的热流差异明显减小。     

改进人工蜂群算法及其轴流压气机优化应用

为了探索更优的全局优化算法并将其应用于轴流压气机气动优化上，将标准人工蜂群（ABC）算法中采蜜蜂和观察蜂探索新蜜源的方式进行改进，从而更好地利用总体蜜源的探索信息，得到改进人工蜂群（IABC）算法。经基准函数测试表明，改进算法既能提升全局寻优能力，又能加快收敛速度。采用IABC算法和经过校核的CFD数值方法搭建优化平台，对单级跨声速轴流压气机Stage35进行优化。优化变量为动叶和静叶径向6个截面的弯、掠、前缘重弯值和尾缘重弯值，以流量和压比相对变化保持在0.5%以内为约束条件，以提高绝热效率为优化目标。结果表明：在设计转速下，优化后设计点绝热效率提升0.83%，全工况范围内平均绝热效率提升2.0%，同时喘振裕度提升1.0%，验证了IABC算法在轴流压气机优化中的有效性。     

LENS风洞试验返回器模型气动热特性模拟

国家数值风洞高超声速流动模拟软件HyFLOW的研制对打破国外同类软件的技术壁垒具有重要意义。与国外DPLR软件进行了对比研究，同时系统介绍了HyFLOW软件求解器的数值方法、物理化学模型以及壁面催化特性计算模型等主要方法，采用典型算例对有限催化模型进行了数值验证，最后基于LENS风洞试验146 mm返回器标模外形开展了高超声速气动热特性数值模拟。研究结果表明，HyFLOW软件在高超声速热化学非平衡流动模拟与评估方面的气动力计算精度高，与国外同类软件DPLR相当，同时其壁面催化条件下的气动热计算精度可靠，可信度高。     

全局方向模板在高阶精度非结构有限体积方法中的推广

梯度与高阶导数重构是影响高阶精度非结构有限体积(Finite?Volume,?FV)格式计算效果的主要过程,其中,不同的模板选择方式发挥了重要作用.传统的模板选择方式往往依赖于固定的网格拓扑关系,无法有效反映流动变化特征,并且随着求解精度的提高,模板单元的数量上升明显,导致找到的模板单元包含过多冗余信息的同时,显著增大计算量,降低求解效率.基于此现状,文章将基于二阶精度FV格式发展的全局方向模板推广至高阶精度FV方法,以充分发挥模板的空间延展性优势,并减少冗余的模板单元数量.此外,文章通过基于制造解的流动与真实超声速涡流两个数值算例,测试了全局方向模板的数值表现.经检验,全局方向模板的使用可有效减少重构过程所需的模板单元数量,并且计算误差相比传统基于网格拓扑的共点、共面模板更低,计算稳定性优于局部方向模板.因此,全局方向模板选择方法在三阶精度非结构有限体积方法中具有较好的数值表现,具备进一步推广与应用的可行性.     

非恒定增压式微尺度拉瓦尔喷管流动特性

针对微尺度喷流在航天器运动状态切换时出现的非恒定增压变化,采用直接模拟蒙特卡洛(DSMC)方法对阶跃式增压和线性式增压两种模式下的微尺度拉瓦尔喷管流场进行了模拟,并对其变化过程中的流动特性进行了对比分析。结果显示：阶跃式增压会导致流动特性出现较大幅值的峰谷式波动,而线性式增压下的流动特性则呈现出线性变化的特点;黏性力对微尺度喷流的非恒定增压变化产生了重要的黏滞作用,在喉部扩张段至出口的流场中尤为明显;在设定的条件下,阶跃式增压过程中喷流产生的总冲量较线性式增压高59.5%,质量流量高74.7%,单位工质提供的冲量低约8.6%,波动性也高于线性式模型,阶跃式增压适用于系统需要较大推力改变运动状态且推进剂充足的情况,而线性式增压在系统精确微调或需要推进剂产生更高效能时具有明显的优势。     

基于欠采样的零阶优化算法

非平衡学习吸引了许多研究者的关注。一般情况下,少数类是更值得关注的,并且其误分类代价要远高于多数类。由于非平衡数据分布的非均衡性,标准的分类算法将难以适用。为了解决非平衡数据分类问题,给出了基于欠采样的零阶优化算法。首先,为了降低数据非平衡分布的影响,针对不同非平衡比的数据集给出了不同的两种采样策略。然后,采用了一种引入间隔均值项的支持向量机(Support vector machine,SVM)优化模型进行分类,并使用带有方差减小的零阶随机梯度下降算法进行求解,提高了算法的精度。在非平衡数据上进行了对比实验,实验证明提出的方法有效提高了非平衡数据的分类效果。