旋转爆轰燃烧室高度对流场影响

为了研究旋转爆轰燃烧室几何结构对旋转爆轰波的影响规律,针对旋转爆轰燃烧室高度对爆轰流场结构和爆轰波传播特性问题,采用2维数值模拟方法分析了爆轰波在4种高度的燃烧室内的参数变化和波系特征。结果表明:在4种高度的燃烧室内爆轰波均能自持传播,且爆轰流场结构相似;当燃烧室高度为4 mm时,爆轰波后的反射激波压力较高;在一定范围内,爆轰波传播速度随燃烧室高度的增加而增大,在爆轰波高度为10、16和20 mm时,爆轰波波锋面的压力分别为3.93、4.61、4.69 MPa。     

燃烧室轴向和周向长度对气液两相旋转爆轰特性的影响

为了研究旋转爆轰发动机燃烧室轴向和周向长度对气液两相旋转爆轰特性的影响,采用守恒元和求解元(CE/SE)方法对带化学反应的汽油/富氧空气两相旋转爆轰理论模型进行求解,获得了气液两相旋转爆轰流场结构,并分析了燃烧室轴向和周向长度对燃烧室流场、爆轰波传播特性以及发动机推力性能的影响。计算结果表明:轴向长度对燃烧室上游流场影响甚小,却对下游流场参数影响较显著。随着轴向长度增加,燃烧室出口压力、温度、密度以及周向速度均降低,轴向速度则逐渐增大,同时发动机平均推力密度和燃料比冲先增大后减小。当周向长度过短,燃烧室内难以形成自持传播的爆轰波,随着周向长度增加,上游爆轰波强度增加,对应的流场参数均有所增大,但发动机推力性能略有所降低。     

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段；在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波；第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波；稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s；斜激波的存在使燃烧室出口平面流场产生了较大波动。      

旋转爆轰流场的数值模拟

基于带化学反应的二维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对充有当量比的氢气-空气预混气和空气的环形旋转爆轰流场,从点火燃烧到发展成旋转爆轰的过程进行了数值模拟。根据数值结果分析了波后流场中爆轰产物受激波、高温和离心力等作用而挤向外壁,形成有利于充入燃料,实现持续稳定旋转爆轰的流场特征。还讨论流场中爆轰波、激波与间断面和内外壁面反射或折射,从而形成多个激波相交的波系特征。为认识和理解旋转爆轰流场,开展旋转爆轰的实验研究等具有现实的指导意义。     

连续旋转爆轰燃烧室增压特性的数值研究

为了分析爆轰燃烧室几何尺寸对增压特性的影响,利用二维欧拉方程数值研究了当量H2/Air在连续旋转爆轰燃烧室中的燃烧流场情况。研究表明:在准稳定状态下爆轰燃烧流场特征参数随时间呈周期性振荡,且其振幅保持不变;横波及局部爆炸波的存在是爆轰波后侧附近总压沿x轴呈现波动变化的主要原因;影响连续旋转爆轰燃烧室增压特性的直接因素为斜激波的高度及爆轰波强度,减少轴向尺寸或增加周向尺寸可增强其增压特性。连续旋转爆轰燃烧室的增压比高达2.52,与传统的燃烧室相比,大大提高了燃气的做功能力。     

旋转爆震三维非预混混合特性及流场结构研究

为深入研究旋转爆震三维非预混流场结构,基于7组分8步反应的化学反应模型,开展了H2/Air旋转爆震流场的数值模拟研究,计算中考虑了黏性的影响。计算结果表明:该喷注结构在冷流流场中混合效果较好,能够在较短距离内实现燃料和氧化剂的充分混合,但起爆后燃烧室内形成的高压环境使得外壁面附近来流可燃气难以到达,可燃气体层主要靠近燃烧室内壁面;旋转爆震波沿燃烧室内壁面周向传播,爆震波后沿径向和周向方向形成"双诱导激波"结构,并进一步导致爆震波后出现高温区和高压区"不吻合"的流场现象;在燃烧室入口截面,爆震波后形成"三诱导激波"结构,诱导激波通过空气环缝向上游传播并对来流燃料和氧化剂的喷注产生影响。     

旋转爆轰波与涡轮平面叶栅相互作用数值模拟

由于旋转爆轰燃烧室具有自增压特性，可提高热力循环效率，因此将旋转爆轰燃烧室应用于燃气轮机可进一步提高系统的性能。基于非稳态雷诺时均Navier-Stokes方法，采用剪切应力输运k-ω湍流模型，建立旋转爆轰燃烧室与涡轮平面叶栅耦合计算模型，研究旋转爆轰燃烧室内的复杂波系与涡轮叶片的相互作用，分析涡轮叶栅对高频爆轰压力振荡的抑制作用。结果表明：旋转爆轰燃烧室内的燃气在涡轮叶栅内加速，并且在斜激波后的局部区域马赫数的增加更为明显。斜激波与涡轮静转子叶片的前缘、压力面、吸力面以及尾缘相互作用，由于旋转爆轰波不同的传播方向，使得斜激波与静子叶片呈相互垂直或平行，进而形成两种不同的波系结构。涡轮叶栅对高频压力振荡存在明显的抑制作用，涡轮叶栅上下游高频压力振荡幅值的衰减率达到80%以上。研究结果展示了旋转爆轰波作用下涡轮叶栅内复杂波系结构特征，并对基于爆轰燃烧推进技术的应用提供了一定的理论基础。     

旋转爆轰波中多波流动模式的数值研究

为了研究旋转爆轰燃烧室内复杂的波系结构,特别是爆轰波多波传播模式的影响因素,采用两步诱导-放热总包反应模型,对简化的二维旋转爆轰波进行数值模拟。在不考虑环形爆轰燃烧室曲率的情况下,研究了进口总温和周向尺寸对单波及多波流场结构的影响。结果表明,一定范围内进口总温的增加会使得爆轰波的波头数目增加,且双波流场结构中存在双波对撞和双波同向两种传播形式,双波同向传播和三波同向传播之间则稳定存在着两组双波对撞流场,且有对应的温度范围;对于双波对撞模式,持续减小周向尺寸能够使得流场结构转变成单波模式;周向尺寸的增加则会使得双波对撞向双波同向传播模式转变。     

旋转爆轰的三维数值模拟

旋转爆轰利用爆轰波在燃烧室的供气端传播,而爆轰产物从另一端排出。为了揭示其宏观特性,采用带化学反应的三维Euler方程,在贴体坐标系中,利用带限制函数的波传播算法,对圆环形燃烧室内的旋转爆轰进行了数值研究,讨论了爆轰传播过程中的波系结构。数值结果表明,即使仅在燃烧室内壁面的部分区域充入可燃气,爆轰波也能够以旋转方式在该区域稳定传播。     

连续旋转爆轰波在无内柱圆筒内的数值模拟

在对同轴圆管内的连续旋转爆轰的已有研究成果进行总结分析后,提出无内柱的连续旋转爆轰的燃烧室模型.利用结合了一步化学反应模型的圆柱坐标系下的Euler方程三维数值模拟了在这种燃烧室内爆轰波的传播行为.数值结果表明:爆轰波由初始起爆的单向传播最终自动收敛稳定到若干个典型爆轰波头的流场结构.在当前进气条件下,流量可以达到450kg/m2s,其基于反应物的比冲约为1900s.与同以H₂/Air作为燃料时同轴连续旋转爆轰发动机(RDE)燃烧室模型的2000s相比,基于反应物的比冲十分接近,但由于内部有可燃气以爆燃形式燃烧等原因,有5%的损失.      

燃烧室宽度对液态燃料旋转爆轰发动机影响实验研究

为了研究燃烧室宽度对液态燃料旋转爆轰发动机工作特性的影响,搭建了气液两相旋转爆轰实验系统,以汽油/富氧空气为工质,氢气/氧气预爆轰管作为点火装置,在不同燃烧室宽度下开展了一系列实验研究,分析了爆轰波的起爆过程,以及燃烧室宽度对爆轰波传播特性与发动机推力性能的影响。实验结果表明：点火后,燃烧室内需要经过一个爆燃转爆轰过程才能形成自持传播的爆轰波;爆轰波在不同燃烧室宽度下均以双波对撞模态传播,对应的波速分布在850~1025m/s内,随着当量比增加,波速整体呈增加趋势;当燃烧室宽度减小,波速整体有所降低;不同燃烧室宽度下推力性能存在显著差异,其中燃烧室宽度在16.5mm下,发动机的推力和燃料比冲要明显低于11.5mm和9mm的;随着燃烧室宽度减小,内外壁面边界层在流场中的作用更为突出,降低了发动机推力的稳定性。     

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

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

旋转爆震发动机喷管非定常流动特性的数值模拟研究

为分析旋转爆震喷管内非定常流动特性,本文基于特征线方法和旋转爆震燃烧室出口时均参数设计了一种塞式喷管模型,通过改变塞锥式喷管的进出口压比（压比分别为15、30和45）,分析了喷管内非定常流场结构与喷管工作性能之间的相互关系,并探究了喷管进出口压比变化对喷管工作性能的影响。研究结果表明：旋转爆震燃烧室下游斜激波进入喷管后,在喷管内形成一道绕喷管壁面螺旋分布的激波,激波结构主要由波前气流角和激波前气流速度决定;爆震喷管局部工作状态依据相对螺旋激波位置分别存在自由膨胀状态、欠膨胀状态、欠膨胀与过膨胀过渡状态以及过膨胀状态四种情况,喷管处于第二种工作状态（欠膨胀）时,喷管内存在激波,并随着激波的传播不断增强,在第四种工作状态达到最大（喷管处于过膨胀状态）,且喷管在第二到第四种工作状态内工作性能变差;在低压比（压比15）工作条件下喷管性能优于高压比（压比30）,其最佳推力和燃料比冲分别为992.05N和1769.52s,因而旋转爆震喷管设计应选择较高的进出口压比作为设计点,使旋转爆震喷管更多的工作在低压比工作状态,以提高旋转爆震喷管工作效率。     

旋转爆震波发动机二维数值模拟

采用一种改进的化学非平衡流解耦方法对Euler反应流方程解耦处理,对流项采用五阶WENO格式离散,化学反应源项的刚性采用简化的隐式方法处理,时间步进采用二阶精度的Runge-Kutta方法,对H2/Air预混气旋转爆震发动机内流场进行了二维数值模拟。模拟结果给出了不同发动机尺寸条件下的详细起爆过程,结果表明当发动机尺寸小于临界直径时无法成功起爆;详细分析了流场结构和爆震波形状,旋转爆震波的传播速度与理论预测值吻合;性能分析结果表明在喷注总压低于燃烧室平均压力时仍可实现推进剂喷注,没装尾喷管的情况下发动机比冲达到176.5s。     

旋转爆震发动机燃烧室壁面烧蚀热防护研究

为了研究碳化烧蚀材料对旋转爆震燃烧室的被动热防护作用效果,建立了旋转爆震燃烧室热环境计算模型和多层复合壁一维烧蚀热响应分析模型。基于旋转爆震燃烧过程的数值模拟,确定了沿燃烧室轴向的壁面平均温度与热流的分布特征;针对典型的壁面热负荷输入条件,采用动边界隐式差分格式,对碳化型烧蚀材料的烧蚀过程和传热过程进行了耦合计算,分析了碳化层剥蚀、热解热、热解气体质量流率、材料厚度等因素对燃烧室壁面热防护效果的影响。研究结果表明,旋转爆震燃烧室温度和热流密度分布不均匀,斜激波扫掠区域热负荷更为剧烈;碳化层剥蚀后会使烧蚀材料的烧蚀率增加,壁面温度迅速上升;增大热解热及热解气体质量流率都有利于延长工作时间、降低壁面温度; 旋转爆震燃烧室不同区域的烧蚀有一定差异,增加烧蚀材料厚度无疑有利于壁面热防护,但是在实际应用中应综合考虑壁面重量和空间尺寸。     

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.     

超燃冲压发动机内外流场三维燃烧数值模拟

以某型煤油燃料超燃冲压发动机整机模型为研究对象,对发动机外流场的流动与发动机内的流动燃烧过程开展了耦合三维数值模拟研究.燃烧室壁面静压计算结果与实验数据吻合良好,平均相对误差为4.9%,验证了数值方法的有效性.结果表明:内流场的流动燃烧状况既受到外流场捕获空气流量的限制,又会反过来影响着外流场的激波结构;泄压孔对保持发动机中各部件流量匹配起到关键作用;由于进气道内激波与边界层的相互作用以及泄压孔诱发的斜激波的影响,流场在第1处燃油喷射孔附近呈现出燃烧流场不对称特征.      

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

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