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

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

过膨胀单边膨胀喷管试验和数值模拟

在低落压比4~10下,对高速二维单边膨胀喷管模型试验件的壁面静压分布进行了试验研究,模拟气体为高压空气;用三维数值计算得到了喷管性能参数,丰富了流场信息。试验和数值模拟结果表明,处于深度过膨胀状态下的单边膨胀喷管内部出现了自由激波、分离激波、两种激波相交及相互影响、自由激波分离和激波诱导分离等复杂现象;上膨胀面大面积的气流分离而形成的低压区使喷管轴向推力系数Cfa仅为0.78~0.87;此外,推力矢量角pδ为18°~25°,这对飞行器的配平和控制面操作带来一定的困难。     

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

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

辅助活门对串联TBCC排气系统的影响研究

在低落压比状态下,对有无辅助活门的两种串联式TBCC组合喷管开展了数值模拟计算,并对包含辅助活门的串联式组合喷管进行了实验研究,通过分析其壁面沿程压力分布与波系结构的变化情况,探究辅助活门对串联式组合喷管气动性能的影响。结果表明：辅助活门结构减少主流喷管外壁面的流动分离区范围,有效的提高了组合喷管内推力。带辅助活门型组合喷管能够引射环境气流,减少引射喷管壁面的回流区范围与底部阻力,从而减少推力损失,提升组合喷管整体性能。提高次流压比或降低主流压比,主流喷管出口处欠膨胀程度降低,声速线往下游方向移动,主流流量减少。提高主流压比或者次流压比都能使辅助活门引射环境气流量增加。辅助活门进气通道距离会改变次流射流与引射喷管的相对位置,距离过小或过大都会造成引射环境流体流量减小。     

单边膨胀喷管试验和数值模拟

在落压比为4～60范围内,对两种不同膨胀型面的高速二元单边膨胀喷管模型试验件(SERN 1和SERN 2)进行了试验研究,获得了壁面静压分布和纹影照片.用三维数值计算得到了喷管性能参数,丰富了流场信息.试验和数值模拟结果表明:下腹板出口气流处于深度过膨胀状态时,喷管内部激波/边界层相互干扰并存在4～5道斜激波,而且SE...     

塞式喷管主喷管角度特性

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

脉冲爆震发动机喷管流动的数值模拟

为了研究喷管形状及充气状况对脉冲爆震发动机推进性能的影响,发展了基于有限速率基元化学反应模型的网格自适应有限体积程序,对影响脉冲爆震发动机推进性能主要过程进行了数值模拟,计算结果和相应的实验结果较为一致。锥形喷管的算例表明,不仅锥角对喷管性能有重要影响,而且同一锥角下喷管内的充气状况对喷管性能也有重要影响。在喷管内充空气情形下,爆轰波进入喷管后迅速衰减为激波,扩张管道的超声速膨胀使得喷管壁面压力显著下降,从而削弱了喷管提高脉冲爆震发动机推进性能的作用。在喷管内充燃气情形下,爆轰波进入喷管后可一直维持到喷管出口,显著提高了脉冲爆震发动机的推动性能。     

激波诱导矢量喷管动态数值模拟

对一二元收扩激波诱导矢量喷管进行了二维动态数值模拟,研究了不同落压比、不同次流加速时间下,喷管气动特性和流场结构随次流增加的变化规律。结果表明:喷管在过膨胀状态下,次流流量达到一定值时,其推力会有较强的振荡;在完全膨胀和欠膨胀状态下,振荡较小。其原因是,次流增加过程中,喷管出口附近有复杂波系存在和回流产生,而回流对提高推力和增加矢量角有益。     

截短单边膨胀喷管的试验和数值模拟

基于最短长度理论喷管设计程序,获得了设计点为马赫数6.5时单边膨胀喷管(Single ExpansionRamp Nozzle,SERN)等熵膨胀型线,并对其截短70%后的喷管进行模型设计和加工。通过试验和数值模拟结合的方法,研究了模型喷管在不同压比下的流场结构和性能。结果表明:喷管长度截短70%对推力性能影响不大,在所有典型压比下,喷管截短70%后推力系数与原喷管相差不超过1%。该喷管截短70%后,其推力系数在落压比大于60时仍能达到0.95以上,随着喷管压比的继续降低,喷管推力系数迅速降低,性能恶化,喷管在过膨胀状态下的性能应该得到重点关注。     

某型喷口加力调节器减压比特性改进设计

为解决某型喷口加力调节器在配装发动机使用过程中存在的涡轮落压比偏离设计值问题,对一维定常可压缩拉瓦尔喷管的气体流动状态进行了理论分析计算,从涡轮落压比的连续性工作要求分析了拉瓦尔喷管正常工作的使用条件,提出了空气减压器二级减压拉瓦尔喷管扩张段流道形状的改进设计方法,并通过了发动机试验验证,结果表明:空气减压器只有在超临界状态工作,并且激波位置在测压点位置后面时,才能保证减压比仅与针塞位置有关,与进气压力大小无关;进气压力较低时,激波位置离测压点较近,会造成减压比相对稳定状态存在偏差;增大喉道面积,可使得相同进气条件下,激波位置后移,远离测压点,有利于提高减压比的稳定性,增大后端角度会导致激波位置前移,不利于减压比的稳定。      

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

上浮水雷在工作过程中环境压力大幅度变化，针对上浮水雷火箭发动机具有三种不同工作状态（欠膨胀、完全膨胀、过膨胀）的特点，采用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以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

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

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

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

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

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

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