复合固体推进剂细观结构建模及脱黏过程数值模拟

为研究复合固体推进剂损伤演化规律,基于分子动力学颗粒填充算法构建了HTPB(hydroxyl terminated polybutadiene)推进剂细观结构模型,通过在AP(ammonium perchlorate)颗粒/HTPB基体界面处引入黏接接触替代传统的黏接单元,并基于Hooke Jeeves的参数优化算法反演得到颗粒/基体界面处内聚力模型参数,利用双线性和自定义指数型损伤内聚力模型模拟了AP颗粒和HTPB基体黏接界面处损伤的萌生、发展、聚合直至宏观裂纹破坏的过程。通过数值仿真与实验结果对比发现,指数型损伤内聚力模型比双线性模型能更准确描述推进剂单轴拉伸过程中颗粒与HTPB基体界面间脱黏过程。最后对比了多阶段加载实验结果与仿真结果曲线,发现两者变化趋势基本一致,最大偏差仅为10%,验证了所建细观模型的可靠性及反演所得界面参数的准确性。      

基于正激波位置计算的进气道/发动机系统实时仿真模型

研究了超声速条件下,基于进气道内正激波位置计算的进气道/发动机系统综合仿真问题。提出了将超声速进气道分为内、外两部分独立计算的方案,外部通过斜激波计算公式获得进气道内部的边界条件,对进气道内部采用准一维CFD模型模拟背压对正激波位置的影响。相比二维CFD计算,在保证参数精度要求的前提下有效缩短了计算时间,提高了仿真实时性。进一步将进气道模型与发动机部件级模型匹配运行,构建了进气道与发动机综合实时仿真模型。进气道/发动机综合模型的仿真表明,其能够准确模拟进气道正激波位置对涡扇发动机的动稳态影响,当来流发生扰动时,通过喷管喉道面积和压缩部件导叶角的快速调节,能够抑制正激波位置的变化,并使推力下降量减少50%。      

低温液氮贮箱增压及排气流量控制方法

在对以液氮为工质的低温贮箱进行热分析的基础上推导了增压速率计算模型,并通过两种增压工况进行了验证。主要开展了直排模式、排气模式及并行模式下的气枕压力及蒸气排气流量的对比实验,通过后两者模式下的蒸气排放质量与直接排放相比,分析了压力控制带内单循环各个模式下热力排气系统（TVS）的效能,结果表明:采用低温蒸气节流的并行模式在蒸气排放质量和气枕降压速率上具有优于其他模式的明显优势,气相并行模式比直排模式质量减少97%。      

氧化层非均匀增长对热障涂层应力分布的影响

利用扫描电镜对不同氧化程度的热障涂层(TBCs)进行观察,分析发现热生长氧化层(TGO)各部位并非均匀增长。假设热生长氧化层的形状为正弦曲线,用统计方法处理扫描电镜图片可以得到热生长氧化层的基本几何参数。由此建立单波形二维轴对称有限元模型,计算并分析了非均匀增长模型的应力状况。结果表明:相同温度下,非均匀增长模型的轴向、剪切和Mises应力分量的变化趋势与均匀增长模型的对应趋势一致;但是非均匀增长模型的各个应力分量最大值均小于均匀增长模型的应力分量最大值。      

航空发动机双转子系统“临界跟随”现象的机理及影响

为了探明“临界跟随”现象导致航空发动机振动居高不下的产生机理,建立了考虑中介轴承的双转子动力学模型,推导了双转子发生“临界跟随”现象时参数之间的关系,分析了“临界跟随”状态下转子系统的动力学特性。结果表明:“临界跟随”现象会导致转子系统无法越过盘偏摆振型临界转速,造成转子振动对不平衡质量分布极其敏感,其主要影响因素为盘极转动惯量与直径转动惯量的比值(简称惯量比)以及双转子增速比;当高压盘等效惯量比等于1,或增速比与任意低压盘惯量比相等时,转子系统表现出高压激励条件下的“临界跟随”特征;当低压盘惯量比等于1,或增速比与高压盘等效惯量比的乘积等于1时,转子系统表现出低压激励条件下的“临界跟随”特征。提出了必须严格控制航空发动机叶盘的等效惯量比以及双转子转速控制律的设计建议。      

切削过程有限元仿真研究进展

得益于高性能计算机和软件技术的发展,数值仿真为深入理解切削过程、研究刀具磨损、提高加工表面质量等提供了强大的技术支持。总结了现有几种本构模型的特点,对构建切削仿真有限元模型中的几个关键问题进行了分析和总结,综述了有限元仿真技术在金属切削过程中切屑形成研究、切削温度和切削力、刀具磨损、刀具优化和残余应力预测以及加工表面显微组织演变模拟等方面的研究成果。最后,系统地分析了切削过程有限元仿真研究中存在的问题,探讨了该研究领域未来的发展趋势和需要进一步探索的热点问题。     

基于代理模型的大型民机机翼气动优化设计

先进的气动优化设计思想与方法,对于提升大型民机气动与综合性能具有至关重要的意义。探讨了大型民机超临界机翼气动优化设计的基本准则和要点,并结合代理优化算法,提出了一种面向工程应用的多轮次高效全局气动优化设计方法。首先,通过一系列解析函数/翼型/机翼优化测试算例进行了验证。其次,将所提出的方法与人工修型相结合,开展了针对宽体客机超临界机翼的两轮气动优化设计,使其气动性能得到显著改善。最后,采用不同的雷诺平均Navier-Stokes（RANS）方程求解器对安装优化机翼的全机巡航构型进行了典型状态的气动性能综合评估。研究结果表明,所提出的代理优化算法具有很高的优化效率、较强的约束处理和全局优化能力;将所发展的基于代理模型的多轮次气动优化设计方法与人工修型相结合,能够获得满足设计要求的气动外形,验证了该方法在大型民机超临界机翼气动设计中的有效性和工程实用性。     

钝体稳定射流火焰的火焰面模型数值研究

基于稳态火焰面(SLFM)和交互式非稳态欧拉颗粒火焰面(EPFM)模型对Sydney大学CH4/H2钝体稳定湍流扩散火焰行了数值研究,采用修正的雷诺应力模型(RSM),同时对两种不同规模的甲烷详细化学反应动力学机理进行研究,比较了燃烧模型和反应机理对湍流火焰结构、活性自由基以及氮氧化物预测精度的影响,与实验数据对比结果表明:两种反应机理得到的温度场和主要组分分布基本相同;SLFM模型能对速度场和标量场的分布进行较为准确预测,采用EPFM模型修正后,部分区域OH预测结果更加靠近实验结果;采用EPFM模型对SLFM模型耦合GRI-Mech 211的计算结果修正后,NO量级降低近2倍,预测精度明显改善,与实验结果实现较好的符合,验证了化学反应动力学机理以及非稳态效应对氮氧化物预测精度的影响。     

气动力/热与结构多场耦合试验模型方案初步设计

以多渠道、多机制交叉耦合为热防护结构特点的新一代高超声速飞行器必须采用气动力/热与结构多场耦合计算方法进行研究。目前,国外已建立较完善的耦合分析系统并用于飞行器研制,国内的中国空气动力研究与发展中心(CARDC)也已自主研发了热环境/热响应耦合计算平台(FL-CAPTER)。为验证多场耦合计算平台所用方法的有效性和计算结果的准确性,设计并开展气动力/热与结构耦合的地面试验具有十分重要的意义。本文结合气动力/热与结构多场耦合试验设计需求,以现有材料和设备能力为依托,开展了试验风洞选取、模型尺寸估算、模型材料选择、模型气动设计与模型结构设计工作。初步研究表明,模型支撑结构附近迎风面局部高温热膨胀将有利于模型前体结构产生可观的整体变形量。本文以此设计了带压缩拐角的二级压缩面结构模型,通过短时间不锈钢模型验证试验和计算对比分析初步验证了模型设计的可行性,并以此为基础预测了高温合金模型的试验结果。为下一步开展高温合金长时间风洞试验奠定了技术基础。     

Numerical simulations of turbulent flows in aeroramp injector/gas-pilot flame scramjet

To uncover the internal flow characteristics in an ethylene-fueled aeroramp injector/gaspilot(ARI/G-P)flame scramjet,a Reynolds-averaged Navier-Stokes(RANS)solver is constructed under a hybrid polyhedral cell finite volume frame.The shear stress transport(SST)k-x model is used to predict the turbulence,while the Overmann’s compressibility corrected laminar flamelet model is adopted to simulate the turbulent combustion.Nonreactive computations for Case 1(G-P jet on),Case 2(ARI jets on),and Case 3(both ARI and G-P jets on)were conducted to analyze the mixing mechanism,while reactive Cases 4–7 at equivalent ratios of 0.380,0.278,0.199 and0.167 respectively were calculated to investigate the flame structure and combustion modes.The numerical results are compared well to those of the experiments.It is shown that the G-P jet plays significant role in both the fuel/air mixing and flame holding processes;the combustion for the four reactive cases takes place intensively in the regions downstream of the ARI/G-P unit;Cases 4 and 5are under subsonic combustion mode,whereas Cases 6 and 7 are mode transition critical and supersonic combustion cases,respectively;the mode transition equivalent ratio is approximately 0.20.