超燃冲压发动机二维热环境数值模拟

对超燃冲压发动机热环境进行了研究和计算.以二维N-S方程和一维瑞利加热规律为基础, 建立了超燃冲压发动机内流场的热环境计算模型;用MacCormark预测-校正格式编写了数值计算程序;对飞行马赫数Ma=6和Ma=8两种工况的超燃冲压发动机内部热环境进行了数值模拟, 得到了发动机内部流场的温度分布和壁面热流分布;对计算结果进行了分析.结果表明, 采用二维无反应流体计算和一维加热规律相结合的方法, 在超燃冲压发动机热环境研究中是可行且有效的.      

高超声速飞行器进气道等离子体虚拟前沿优化设计

为研究等离子体虚拟前沿在高超声速飞行器前体/发动机一体化设计中的应用,首先设计了二波系前体,然后分析了虚拟前沿对前体流场结构的影响,最后通过比较虚拟前沿结构参数及位置对总压恢复系数、流量系数等参数的改善效果,优化得出了最佳虚拟前沿。结果表明,有虚拟前沿控制前体激波时,总压恢复系数及流量系数都会增大,但压缩比不一定会增大;虚拟前沿的位置对各参数的影响很大,且存在一个最佳位置使得前体性能最好;随着虚拟前沿长度的增大,总压恢复系数和流量系数逐渐增大,压缩比逐渐减小。     

源项法模拟高超声速飞行器内外一体化流场

为了简化燃烧室流场模拟，提高内外流场数值仿真效率，为高超声速飞行器一体化构型提供设计依据和参考，发展了一种在流动控制方程组中加入源项的数值模拟方法。通过调节源项作用位置及大小，该方法可实现对超燃冲压发动机燃烧室内的添质添能流动的模拟：与直联式试车台实验数据对比表明，该方法能够较好地模拟冲压发动机内压力分布，能够满足飞行器气动／推进一体化构型方案设计及总体方案设计阶段的需求。     

吸气式高超声速飞行器机体推进一体化技术研究进展

吸气式高超声速一体化飞行器最显著的特点是子系统之间的耦合较其他类型飞行器更加强烈,这使得其设计具有挑战性。所有的子系统之间部件相互干涉,包括:气动、推进、控制、结构、装载和热防护等,特别是机体与超燃冲压发动机之间的耦合最为突出。飞行器的前体和后体下壁面既是主要的气动型面,又是超燃冲压发动机进气道外压缩型面和尾喷管的膨胀型面,在产生推力的同时也产生升力和俯仰力矩。机体与发动机的强耦合作用对飞行器的推力、升力、阻力、俯仰力矩、气动加热、机身冷却、稳定性和控制特性有直接的影响。本文介绍了国内外机体推进一体化技术的研究进展,重点介绍了中国空气动力研究与发展中心(CARDC)的相关研究工作,包括:密切曲锥曲面乘波进气道和基于双激波轴对称基准流场内转式进气道设计方法、独创的大尺度脉冲式燃烧加热风洞一体化飞行器带动力试验技术和高超声速内外流耦合数值模拟技术等。对高速飞行中激波边界层相互干扰、流动分离机理、可压缩湍流转捩及其控制、超燃冲压发动机燃烧流动机理等相关基础问题也进行了研究,强调了对高效高精度计算方法的迫切需求。     

跨声速涡轮叶栅三维流场的数值模拟

以具有试验测量数据的一个高压涡轮动叶中截面平面叶栅和一个高压涡轮导向器环形叶栅为对象,采用计算流体动力学(CFD)软件CFX-TASCflow,对两个涡轮叶栅的跨声速三维粘性流场进行了数值计算,验证了CFX-TASCflow软件在航空燃气涡轮设计计算中的可靠性,并对所研究的跨声速涡轮叶栅粘性流场进行了分析.      

超燃冲压发动机进气道不起动仿真研究

对超燃冲压发动机进气道-隔离段进行了二维稳态流场数值模拟,给出了超燃冲压发动机进气道起动到不起动整个过程的数值模拟结果。得到了不同压比下、不同攻角下进气道-隔离段内部流场等压线图,分析了进气道-隔离段壁面静压分布特性,初步给出了进气道不起动状态的判断准则。      

超燃冲压发动机磁控进气道设计影响因素分析

阐述了磁控进气道的基本原理,建立了一维数学模型。对磁控进气道设计影响因素磁感应强度、负载系数和通道截面积变化率进行了特性分析,得到磁控进气道出口参数随设计参数的变化规律,为磁控进气道设计提供了依据。给出了一组磁控进气道设计参数,仿真结果表明采用磁控进气道技术可以有效地调节燃烧室进口来流状态,提高比冲,拓宽飞行马赫数范围。      

超燃冲压发动机支板热环境及热防护方案

在超燃冲压发动机工作过程中,支板前缘的热环境非常恶劣。本文研究了超燃冲压发动机支板前缘的热环境,得到了热载荷分布。提出了支板前缘金属结构再生冷却方案、耐烧蚀材料热防护方案和气体喷射热防护方案,并对这三种方案的热防护效果进行了数值模拟。数值模拟的结果表明,金属结构再生冷却方案无法对支板进行有效的热防护;耐烧蚀材料方案可以在飞行马赫数8以下起到很好的热防护效果;而当马赫数大于8时,则只有气体喷射方案可以实现有效的热防护。      

吸气式高超声速飞行器前体/进气道一体化设计与试验

针对吸气式高超声速飞行器前体/进气道一体化设计的问题,以总体指标为约束,采用数值设计手段开展了前体/进气道一体化设计,并对高超声速飞行器进行测压/测力试验,考核了前体/进气道的一些主要性能,结果表明:①设计状态下,数值计算结果表明前体/进气道性能符合总体指标要求,设计手段有效;②数值手段模拟结果和风洞试验结果吻合良好,流量系数最大误差为4％,总压恢复系数最大误差为42％,数值算法有效;③前体/进气道的附加阻力随来流马赫数的增大而减小,0°攻角下,在来流马赫数为4时,附加阻力占总阻力的172％,在总体设计时应予以考虑;④在进行吸气式高超声速飞行器通流测压/测力试验设计时,应充分考虑进气道不起动的试验预案,防止由于进气道不起动导致整个试验的失败.      

流线追踪Busemann进气道马赫数3.85实验研究

为了研究流线追踪Busemann进气道在低马赫数来流条件下的起动性能、总体性能以及抗反压能力, 对设计马赫数5.3、收缩比CR=5的6°截短后的流线追踪进气道进行了马赫数3.85风洞实验.结果表明:该进气道在马赫数3.85来流条件下能够自起动;进气道流量系数为0.83, 出口总压恢复0.78, 马赫数2.01, 增压比11.92;进气道能够承受的最大反压为34倍的来流静压.      

改进的翼型积冰数值模拟方法

基于微团概念,以雨滴运动的基本假设出发,对翼型积冰的拉格朗日方法进行了改进.新方法克服了原方法计算效率低、冰形同雨滴分布密度及网格密度相关的缺点,继承了原方法的鲁棒性,不存在欧拉方法的收敛性问题.为验证方法的正确性,计算了NACA0012翼型在0°和4°攻角下的霜冰积冰情况、表面水收集量和局部撞击系数.结果表明新方法减少了计算量,预测的冰形和实验结果一致,并分析了积冰对气动性能的影响.      

气动热数值模拟中壁面热流后处理方法

针对气动热数值模拟中壁面热流后处理方法进行了系统的数值研究,提出了一种有别于常规微分法和积分法的壁面热流后处理方法.通过3种算例进行了对比分析,结果表明:即使在流场参数相同情况下,采用不同的壁面热流后处理方法计算的热流结果也会有较大差异,且表现出不同的网格依赖性.微分法对网格的依赖性较大,单纯提高差分精度不一定能够减弱网格依赖性和提高计算精度,还与近壁处网格分布密切相关;积分法从能量平衡积分方程出发,极大地减弱了网格依赖性,具有较高的计算精度.所提出方法具有与积分法相当的网格依赖性和计算精度,尤其在驻点区域,其网格依赖性更弱于积分法,且计算简单,提高了计算效率.      

复杂积冰翼型气动性能分析

为了克服复杂外形难以划分结构化网格的问题, 采用结构/非结构混合多块网格对复杂积冰翼型的气动性能-进行了分析计算.在靠近积冰边界的内层采用了非结构网格, 在非结构网格外采用结构化网格以节省存储空间和计算时间.不同类型网格之间采用了点点对接的方式以简化通量重构的处理方式.整个流动区域采用SST湍流模型, 求解了雷诺平均N-S方程.数值计算结果表明结构/非结构混合多块网格在复杂积冰翼型气动性能计算上的可行性, 积冰对翼型的气动性能造成了严重的影响.      

降落伞主充气阶段的动态仿真及流场分析

建立了有顶孔的平面圆形降落伞的轴向和径向动量方程,使用数值方法计算出选定的降落伞在主充气过程中各种参数的变化,对降落伞的主充气过程进行了动态仿真.在使用动量方程计算出的的充气模型的基础上,利用降落伞在各个时间结点时的外形和下降速度等相关数据,采用计算流体力学的方法和准定常假设,对降落伞充气过程中的流场进行了模拟和分析,对降落伞充气过程的机理及流场衍变情况进行了探讨.      

夹气喷嘴瞬态喷雾的CFD仿真及试验

通过CFD仿真和试验的方法研究了夹气喷嘴的瞬态喷雾特性并验证了喷雾仿真模型的准确性.利用Fluent软件在不同环境背压和不同喷气压力的条件下进行了瞬态喷雾的CFD仿真计算,同时利用高速摄影机拍摄定容弹内部的喷雾形态,以及通过激光粒径测试仪对大气环境下的喷雾进行燃油液滴SMD的测量.结果表明:喷嘴出口处高速流出的压缩气体与燃油发生强烈的耦合,并于下游处产生强烈的涡环结构;试验中环境背压减少0.2MPa时所增加的平均贯穿速率是喷气压力增加0.2MPa时的3.62倍;喷雾平均贯穿速率的试验结果与仿真结果的平均误差为9.94%;试验得到的喷嘴外围的喷雾平均索特直径与仿真平均索特直径的平均误差为4.33%.      

高超声速二维混压式前体/进气道设计方法研究

以飞行Ma数Ma=6,H=25km为设计点,分别采用等激波角和等激波强度设计方法,并考虑变比热、激波与附面层干扰等因素的影响,分别对唇口平直和唇口带有斜楔的超燃冲压发动机二维混压式前体/进气道进行了初步设计,分析并比较了几种方案进气道的设计点和非设计点性能及二维流场。研究表明,在低飞行Ma数下,唇口带有斜楔的前体/进气道起动性能和总压恢复优于唇口平直的,在高飞行Ma数下,唇口平直的前体/进气道冲压比高、外罩阻力小,而唇口带有斜楔的前体/进气道总压恢复系数高,外罩阻力相对较大。另外,分别采用等激波角和等激波强度方法设计的前体/进气道性能相近。本文提出的方法对于二维混压式高超声速前体/进气道方案的初步筛选具有一定的适用性。      

一种带前体的高超声速矩转圆形进气道研究

对一种带前体的高超声速矩转圆形进气道开展了数值仿真及高焓脉冲风洞试验研究。结果表明：(1)来流马赫数为6时进气道出口马赫数为2.60,总压恢复系数为0.40,增压比为39.8,流量系数为0.769;来流马赫数为5时进气道出口马赫数为2.28,总压恢复系数为0.45,增压比为19.7,流量系数为0.643;(2)超声通流时,内通道上、下壁面静压大幅波动;(3)隔离段内较强的横向压力梯度使得进气道出口流场畸变较大;(4)前体横向溢流较大,对进气道性能不利;(5)数值仿真与实验结果吻合较好,验证了计算方法的可靠性。     

战斗机垂尾脉动压力数值模拟

在亚跨超计算流体力学（CFD）软件平台（TRIP）上开发了基于RANS/LES混合思路的IDDES流动模拟技术,并通过NACA0021翼型60°大迎角分离流动与串列圆柱绕流模拟对RANS/LES混合方法的精确度进行了验证,针对某战斗机外形的垂尾脉动压力开展了数值模拟研究。战斗机来流马赫数为0.1,基于全机长度的雷诺数为2×10⁶,模型迎角为20°、30°和40°。分别通过脉动压力系数、脉动压力功率谱密度、空间流动结构以及侧向力响应曲线等对战斗机的垂尾脉动压力进行了分析。脉动压力模拟结果表明：当垂尾完全沉浸在边条翼脱体涡破碎后的宽频湍流脉动气流中时,垂尾翼梢位置的脉动压力会发生明显的增大。