三单元气氢/气氧塞式喷管发动机实验

以气氢/气氧为推进剂, 对三单元直排塞式喷管模型发动机(内喷管和燃烧室内衬采用耐烧蚀材料钨渗铜, 塞锥为碳钢)进行了热试实验研究.给出了实验结果, 进行了数据分析.在两个压比(燃烧室压强与环境反压之比)下, 获得了塞式喷管性能数据和塞锥表面5个测点处的压强分布.实验表明, 塞式喷管具有良好的高度补偿能力和较高的效率.在压比值为50附近, 效率达到92%-93.5%;在压比值为350附近, 效率达到95%-96%.数值模拟的预示结果与实验数据吻合较好.      

喷管对旋转爆震发动机性能影响的实验

为了研究喷管对旋转爆震发动机（RDE）性能的影响,设计了RDE推力测试台,RDE环形燃烧室的内径和外径分别为70mm和80mm,长度为40mm,以氢气和空气分别作为燃料和氧化剂,采用环缝-喷孔对撞式掺混方式,使用高能火花塞点火,并用压电式压力传感器测量推力,实验研究在燃烧室出口安装收敛喷管、扩张喷管以及拉瓦尔喷管对发动机工作性能的影响.结果表明,入口总质量流量小于0.126kg/s时,在发动机稳定工作的工况范围内,收敛喷管对发动机推进性能的提高最为明显.爆震波的传播速度及RDE的推进性能随着入口质量流量的增大而逐渐提高,而当量比则存在一个最佳值,使爆震波的传播速度及RDE的推力达到最大.基于混合物的比冲最高能够达到95.21s.      

主要结构参数对直排塞式喷管性能的影响

数值模拟了内喷管面积比、内喷管型式和塞锥截短对塞式喷管性能的影响,比较了不同结构参数选择下对应的塞锥尺寸。内喷管在地面处于欠膨胀和完全膨胀时塞式喷管的性能高于内喷管在地面处于过膨胀时的性能;内喷管为矩型喉部矩型出口二维喷管的塞式喷管效率最高,内喷管为圆喉方型出口和圆喉矩型出口的塞式喷管性能接近,但都比具有二维内喷管的塞式喷管效率稍低;塞锥截得越短,塞式喷管的效率越低,不同塞锥截短率塞式喷管的性能差别随着压比的增加而减小;随着塞锥长度的增加,塞锥长度对塞式喷管性能的影响逐渐降低。     

两相流对固体火箭发动机塞式喷管性能的影响

针对固体火箭发动机中气固两相流对塞式喷管性能的影响,利用颗粒轨道模型法计算了环喉型塞式喷管两相流流场,辨识颗粒在流场中的运动轨迹,分析不同粒径颗粒对流场及推力性能的影响。结合传统钟形喷管,比较不同工况下塞式喷管与钟形喷管的推力性能,分析气固两相流对塞式喷管高度补偿特性的影响。结果表明：随着颗粒粒径的增大,颗粒会出现交错式的运动,增加了流场的复杂性;单个颗粒与该截断式塞锥只会发生一次碰撞,且粒径越大,碰撞位置越靠前;单一粒径颗粒带来的推力性能损失与颗粒粒径成正比,与工作压比成反比;真实工作条件下,相比于钟形喷管,气固两相流在设计工况对塞式喷管带来约3%更多的推力效率损失,推力效率最终趋于弱于钟型喷管2.7%,但塞式喷管在低压比工况仍具有可观的性能优势。     

塞式喷管主喷管角度特性

首先建立了塞式喷管角度特性的研究模型，并运用简化的塞锥设计方法设计塞锥型面，在此基础上，通过连续改变主喷管倾角来研究主喷管倾角对塞式喷管的推力性能的影响，并确定在给定设计条件下的主喷管最佳倾角，还分析研究了底部高度，内膨胀比，总膨胀比，总压和飞行高度与主喷管最佳倾角的关系以及变化规律，针对特定的塞式喷管实验发动机，进行了初步的变角度实验，并与计算结果进行比较和分析。     

气氢/气氧塞式喷管模型发动机实验

以气氢/气氧为推进剂,对圆转方内喷管单元直排塞式喷管进行了热试实验研究。介绍了实验系统及实验发动机主要零部件的结构和设计参数,给出了实验参数测量结果、实验照片和数据分析。无再生冷却塞式喷管发动机采用耐烧蚀材料钨渗铜加工内喷管和燃烧室内衬,碳钢材料加工塞锥,成功进行了多次短时间热试实验。在三个压比下获得了塞式喷管性能数据,实验表明,塞式喷管具有良好的高度补偿能力和较高的效率。在CNPR=110附近,效率达到93%~95%;在CNPR=450附近,效率达到96%~98%;在CNPR=1000附近,效率达到93%~96%。预计在设计点的效率不低于98%。     

塞式喷管与当量钟型喷管性能的实验测定和分析

为了合理地评价塞式喷管是否具有高度补偿特性，利用空气作为工作介持，对实验塞式喷管和钟型喷管进行了冷流比较实验研究。三种不同的喷管效率计算方法对两种喷管的实验结果进行了处理，得到不同的性能曲线，分析比较两种喷管在不同高度下的效率。与当量实验钟型喷管相比，实验塞式喷管具有高度补偿特性，算法一和算法二分别比较了喷管性能的不同方面，评价喷管性能时可以综合考虑。     

结构参数对环簇式塞式喷管性能的影响

为研究环簇式塞式喷管设计中人为选择的结构参数对性能的影响,在单元间距与内喷管出口直径比、单元喷管数目、内喷管面积比分别相同的前提下,计算并比较了三种塞式喷管的推力效率.结果表明,高空下内喷管面积比30.2的8单元接触式的效率最高,而低空下内喷管面积比20.7的8单元间隔式的效率最高.因此在总面积比一定的情况下,若塞式喷管在中低空高度下工作,则单元内喷管面积比可取得小些,以保证中低空下的效率;反之,则单元内喷管面积比尽量取得大些.      

塞式喷管冷流实验和数值模拟研究

开展了塞式喷管冷流实验和数值模拟研究。实验模型采用直排的点膨胀构型，实验测量了两种塞体长度的塞式喷管在不同压比条件下的表面压强分布，拍摄了典型工况下流场的纹影照片。数值模拟采用二阶NND格式求解二维N-S方程，计算得到的壁面压强分布与实验结果吻合较好。通过实验和数值模拟验证了塞式喷管高度补偿特性，得到的结论为塞式喷管设计和性能预报提供了依据。     

固体火箭发动机塞式喷管数值模拟和热态试验

为了考察固体火箭发动机塞式喷管的性能,设计了一套单元内喷管数为16、面积比为25的环簇式固体塞式喷管和一对比用钟型喷管,采用数值模拟方法预估了其推力性能.对塞式喷管进行了地面热态试验,测定了其推力性能.结果表明,环簇式塞式喷管地面燃气流动现象和推力性能的数值模拟结果与试验结果吻合.塞式喷管在地面上的推力效率达92.6%,比同面积比和长度的钟型喷管高9.1%,具有明显的高度补偿效应;但高空效率比用钟型喷管效率降低5.8%.      

水下超声速气流流场非定常特性研究

上浮水雷在工作过程中环境压力大幅度变化,针对上浮水雷火箭发动机具有三种不同工作状态（欠膨胀、完全膨胀、过膨胀）的特点,采用VOF两相流模型,建立了水下火箭发动机在不同工作状态喷射流发展的轴对称计算模型,分别在来流速度0m/s和50m/s两种工况下,研究了水下发动机在不同工作状态时喷射流发展规律及流场脉动特性。结果显示,在静水条件下,欠膨胀工况时颈缩发生位置距喷管出口较远,流场压力和喉部流量没有发生脉动现象,其它工况时颈缩发生位置距喷管出口较近,完全膨胀工况和过膨胀工况喷管出口最大压力振幅分别为5.5MPa,7MPa;喷管喉部流量振动幅度约为5.2%,32.8%;在有流速条件下,三种工作状态发生颈缩、胀鼓和回击现象的位置距离喷管出口较远,完全膨胀喷管出口最大压力振幅为2MPa,喷管喉部流量未发生脉动特性,过膨胀工况喷管出口最大压力振幅6.1MPa,喉部流量振动幅度30.4%。     

爆震发动机超声速进气道的非定常耦合流场数值研究

采用带化学反应的非定常数值仿真方法,对爆震发动机(PDE)内进气道、燃烧室、尾喷管的耦合流场进行了分析,并着重研究了爆震压力波及其反射波对超声速进气道内流动结构的非定常干扰过程,对比了3种进气道/燃烧室耦合方式下进气道流动特性的差异.结果表明:通过在燃烧室内填充以化学恰当比预混的燃料和氧化剂,数值仿真可获得稳定自持的爆震波.在爆震波压强作用下所产生的结尾波系和反射激波会发生耦合干扰作用,它们在流道内的运动传播影响了结尾波系的前传幅度.通过对比不同进气道/燃烧室耦合方式的流场特征发现:采用突扩式进气道/燃烧室连接段可利用燃烧室头部固壁的阻碍以及其对爆震高压波的反射作用将大部分的压缩气体存储于燃烧室内,进而减弱进气道与燃烧室之间的耦合干扰作用,对提高进气道的工作稳定性有利.另外,内凹式燃烧室头部段的引入还为压缩气体提供了额外的存储空间,故可以进一步缓冲爆震高压,提高进气道的稳定工作裕度.      

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

脉冲爆震燃烧室与冲击式涡轮匹配机理及效率的数值研究

针对脉冲爆震燃烧室（PDC）与涡轮共同工作时,涡轮入口为强非稳态气流的特点,结合涡轮损失机理,确定适合脉冲爆震来流的涡轮类型并开展气动设计工作,运用数值方法计算所设计涡轮的效率和探究PDC与涡轮相互作用的机理。研究结果表明:部分进气、小反力度的冲击式涡轮更适合脉冲爆震来流;在喷嘴收缩段、动叶叶片前缘以及动叶压力面生成向上游传播的反射激波,造成能量损失;在设计点,涡轮效率约为75%。上述研究结果可以为脉冲爆震涡轮发动机的涡轮设计提供一定的参考。      

当量比对连续旋转爆轰燃烧室特性的影响

为了探究当量比对甲烷-空气连续旋转爆轰燃烧室（CRDC）特性的影响,利用二维可压缩欧拉方程对CRDC进行了数值研究,分析了爆轰波的发展过程和贫燃熄火过程,对比了不同工况下CRDC特性参数的变化情况。结果表明:CRDC起爆后燃烧场在由不稳定状态到相对稳定状态的过程中发生了2次碰撞,当进气当量比较低时,CRDC未能完全发生2次碰撞过程就已经熄爆。随着进气当量比的降低,爆轰波传播速度、轴向平均速度、出口平均温度、出口平均总压均呈下降趋势;增压比随当量比降低而减小的根本原因在于旋转爆轰燃烧过程和等压燃烧过程的熵增差减小,使吉布斯自由能增量差减小。CRDC的燃料驻留时间处于亚毫秒量级,燃烧热效率保持在99%以上。      

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

圆盘旋转爆震燃烧室中爆震波传播模态

本文的目的在于通过求解14组分19步CH4/O2反应的欧拉方程去研究燃料喷射温度、压力、燃烧室内外直径比值和燃烧室径向长度对圆盘结构下旋转爆震传播模态的影响。研究表明,只有预混气喷射温度范围为500-900K,喷射压力范围为0.5-3.5MPa的条件下能够在固定燃烧室内形成连续的旋转爆震波。受燃烧室结构影响,旋转爆震波的传播模态分为稳定模态和非稳定模态。整个非稳定传播模态根据爆震波对前周期的干涉又可以分为干涉阶段和不干涉阶段。在不干涉阶段,爆震波传播速度略高于Chapman-Jouguet (CJ)速度;在干涉阶段,爆震波传播速度低于CJ速度且爆震波被间断面分割成两部分。非稳定模态爆震波传播速度小于稳定模态,而爆震波夹角、燃烧室出口面积比则反之。     

高总温空气与汽油燃料的旋转爆震验证试验

为验证高总温空气来流条件下汽油燃料旋转爆震的可行性,开展了气液两相旋转爆震发动机试验研究。旋转爆震发动机环形燃烧室外径和内径分别为202mm和166mm,长度为155mm。通过空气加热器模拟高总温空气来流环境,汽油和空气采用分开喷注的方式,分别通过高压喷嘴和环缝进入燃烧室。试验采用垂直安装的预爆震管成功起爆了旋转爆震波,并实现了旋转爆震波的连续稳定传播。试验结果表明：当空气质量流量为1110.0g/s,当量比为0.97,空气总温为713K时,旋转爆震波以双波对撞模态在燃烧室内连续传播,爆震波传播频率为1827.31Hz,与高频压力信号经快速傅里叶变换得到的主频一致,爆震波传播速度为1059.6m/s。在空气质量流量为1110.0g/s,当量比为0.84,空气总温为713K的工况下进行了3s的长程试验,验证了以高总温空气为氧化剂、汽油为燃料的旋转爆震发动机长时间连续稳定工作的可行性,获得的旋转爆震波传播频率为1907.5Hz。     

Experimental research on rotating detonation with liquid hypergolic propellants

This paper describes experimental research into the initiation and propagation of rotating detonation for liquid Nitrogen TetrOxide (NTO) and liquid MonoMethylHydrazine (MMH). An annular rocket-type combustor without nozzle was designed to investigate detonation combustion. The propellants were injected through unlike impingement injectors. The combustion flame fronts and pressure waves were detected using optical diagnostics and dynamic pressure sensors, respectively. The propagation of rotating detonation was established spontaneously by increasing the mass flow rate of propellants. The velocity of propagation of the flame fronts and pressure waves was nearly equal and reaches supersonic speed. Two different detonation combustion patterns are present, single wave mode and double waves mode. And in double waves mode, the two detonation waves are always counter-rotating. The possibility of rotating detonation initiation in a combustor with nozzle was also checked. Stable rotating detonation can be initialized and sustained at similar operating conditions.