氢/空气超声速燃烧研究

Ｈ２／Ａｉｒ在两种不同的燃烧室尺寸、七种燃烧喷注方式下进行了系统的超声速燃烧实验。实验空气的滞止温度在２０００Ｋ左右,滞止压力１～１．４ＭＰａ,总流量２ｋｇ／ｓ,燃烧室进口马赫数２．５,可以模拟飞行Ｍ数为７的超燃冲压发动机中的燃烧工况。新开发的一维超声速燃烧程序ＳＳＣ－１可以估算出燃烧室内的流场参数、燃烧效率和总压损失。计算结果与实验进行了比较,发现较好的一致。实验结果表明,利用垂直喷射,燃烧效率可以超过８０％,同时不引起严重的总压损失。由燃烧室壁面静压分布与燃烧效率的分析发现,燃烧室燃料注射位置应避免过于集中,宜分散按规律分布,使燃烧室静压分布尽量平直以获得高燃烧效率。     

三维横向流喷与超声速主流干扰流场的数值模拟

本文对三维隅角机翼声速喷流与超声速主流的干扰流扬进行了数值模拟。三维欧拉方程的求解采用非结构网络有限体积伽辽金法（Finite Volume Galerkin Method）。引入了总体结点积分域的概念,简化了从单元矩阵到总体矩阵的汇总过程。通量的分裂采用Osher格式,通过外差使其由一阶精度上升为二阶精度。发展了一种基于线化流量的逆风非结构网格隐式有限元格式以提高求解精度及效率。最后给出了三维隅角机翼流场的算例。     

一种利用泛函网络进行导航卫星钟差预报的方法研究

为提高导航卫星钟差预报的精度和实时性,提出一种基于多项式和泛函网络相结合的钟差预报方法。该方法首先根据卫星钟的物理特性用多项式进行拟合,并用泛函网络对多项式拟合残差进行建模,其中泛函网络的结构是通过小波降噪及相空间重构的方法确定。最后以GPS卫星为例,进行四组短期预报实验。结果表明,所提方法能够实时有效地对导航卫星钟差进行拟合和预报,且精度优于IGUP星历。     

纹影流动显示技术在超声速射流数值模拟验证中的应用

为了研究总压探针测量方法在欠膨胀超声速射流中的可行性,论文采用非结构网格有限体积法对含探针和不含探针的射流流场进行数值模拟,并通过对两个计算结果的比较,给出这种测量方法的参考性评价。为了证明数值方法及计算结果的可靠性,给出了口径6mm的喷嘴、粗1.3mm总压探针在压比4的流动条件下纹影流场显示和总压探针测量等实验结果。     

挠性壁喷管撑杆单位影响曲线相关性研究及其在2m×2m超声速风洞中的应用

对挠性壁喷管撑杆单位影响曲线的相关性进行了研究,提出了一种通过已知单位影响曲线计算同一马赫数相邻撑杆及相邻马赫数同一撑杆单位影响曲线的方法,结合实验影响法成功用于2m×2m超声速风洞二元挠性壁喷管的型面优化,有效提升了第一菱形区核心流速度场均匀性,大大提高了调试效率,降低了调试能耗,缩短了调试周期.该方法对类似喷管型面优化具有借鉴意义.     

化学动力学模型对H2/AIR超燃模拟的影响

为了研究超燃流动中不同化学反应模型的影响，提出流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超声速H2／Air 预混合气体中钝体激波诱导的振荡燃烧流场，针对飞行速度大于和接近爆轰速度两种不同流动状态，比较Evans-Schexnayder提出的7组分8反应机理（ES模型），Jachimowski提出的13组分32反应机理（J模型）及其简化为9组分19反应机理（J－1模型）对流场参数的影响，计算表明反应机理引起压力场变化较小，对温度等参数影响较大，对于飞行马赫数低于7的问题，取不包括N原子反应的Jachimowski的9组分19反应机理较合适。     

氢气引导乙烯火焰非定常燃烧过程

基于脉冲燃烧直连式试验台,开展了超燃冲压发动机氢气引导乙烯火焰的非定常燃烧过程研究。燃烧室入口条件为马赫数2、总温950 K和总压1.0 MPa。试验过程分为4个阶段：冷流、引导氢气单独燃烧、引导氢气点燃乙烯、乙烯单独燃烧。基于高频壁面压力测量和火焰荧光高速摄影,获得了代表性测点的压力时间曲线及燃烧室内火焰发展历程,提取了压力平均值、振荡幅度和频率、着火时间及反应位置等重要信息,分析了不同燃烧阶段的非定常特性。试验结果表明：在氢气单独燃烧阶段,非定常特性源于凹槽后斜坡区域氢气反应强度的变化。在氢气点燃乙烯阶段,非定常特性由氢气和乙烯火焰的“交接”引起。在乙烯单独燃烧阶段,非定常特性由燃烧和超声速流动之间的耦合引起。     

典型超声速/高超声速动导数计算方法研究

以当地流活塞理论为基础完成超声速/高超声速动导数计算公式推导,发展了一种高效的动导数计算方法.选取1个超声速、2个高超声速动导数算例标模,计算了超声速/高超声速下飞行器动导数变化规律.研究结果表明:BFM超声速流动算例的俯仰组合动导数绝对值随马赫数的增大而减小,计算重心位置后移亦会导致俯仰组合动导数绝对值减小,纵向动稳定性降低;在0° ～20°迎角范围内,高超声速流动中的尖锥和钝头旋成体俯仰组合动导数绝对值随迎角的增大而增大,纵向动稳定性增强;所提方法计算时间约为双时间动导数计算方法的1/27.