圆柱坐标系下连续旋转爆轰发动机的数值模拟

使用两步化学反应模型对连续旋转爆轰发动机(Continuous rotating detonation engine,CRDE)进行了二维数值模拟研究。数值计算获得了同轴圆管腔中间层曲面上连续爆轰的多个循环过程。分析了燃料入射、提前燃烧、爆轰波结构和波传播速度等几个关键问题。计算获得燃烧室内流场结构与之前实验研究结果定性符合。计算中燃料以最大1400m/s沿轴向入射,爆轰波可约以2400m/s沿周向连续循环传播。计算的燃烧室内爆轰波循环频率约为2300Hz。     

连续旋转爆轰发动机流场三维数值模拟

对连续旋转爆轰发动机(continuous rotating detonation engine,简称CRDE)流场进行了三维数值模拟.数值计算获得同轴圆管腔内CRDE的多个爆轰循环过程.依据计算获得流场,分析了可燃气入射、提前燃烧、爆轰波结构等实现连续爆轰的几个关键机理问题.通过多个算例对不同来流总压时的旋转爆轰的推进性能进行了比较与分析,最后计算获得基于燃料的比冲约为2 200 s,燃料流量随总压变化呈线性增长.      

基元化学反应一维爆轰波的数值模拟

应用基元化学反应模型对一维氢氧混合气体爆轰波进行数值模拟。采用一端高温高压点火起爆的方式,模拟爆轰波的传播。对一维爆轰波的结构、传播特征以及各个化学组分的变化趋势进行了研究。通过将模拟的结果同国内外的研究结果进行对比,发现了一些规律,为多维爆轰波的数值模拟和脉冲爆轰发动机的研究奠定一定的基础。     

旋转爆轰流场的数值模拟

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

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

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

连续爆轰发动机热力学性质的数值模拟

在二维圆柱腔中连续爆轰发动机数值模拟的基础上,对喷入流场中的粒子进行跟踪,得到经历不同燃烧过程的粒子轨迹图.分析影响粒子轨迹及其上物理参量变化的主要因素.数值模拟得到相应的示功图和示热图,计算连续爆轰发动机的机械功及热效率,将数值模拟结果与理想ZND(Zeldovich-Neumann-Doring)模型进行定性和定量的比较,所得数值模拟结果与理论结果符合 得很好,验证了连续爆轰发动机的优势.      

斜爆轰发动机燃烧机理试验研究

为了研究斜爆轰发动机的稳定燃烧机理,本文开展了飞行马赫数9的斜爆轰发动机的数值模拟研究和试验研究。首先,设计了全尺度斜爆轰发动机模型,发动机的总长度为2.8m。采用两级进气道压缩,每级压缩角度均为15°。利用三个小支板在进气道前缘主流核心区中进行氢气的喷射和混合。其次,对氢气混合过程和发动机燃烧过程进行了数值模拟研究。控制方程采用带化学反应的雷诺平均的N-S方程,湍流模型采用SST k-ω模型,化学反应模型采用9组份19步反应的基元反应模型。数值结果表明,氢气在进气道内混合得比较均匀,在燃烧室内获得了稳定的斜爆轰流场和正爆轰流场。接下来,在激波风洞中开展了马赫数9状态下的斜爆轰发动机稳定燃烧机理试验研究,在50 ms的风洞有效试验时间内获得了持续稳定的斜爆轰流场,试验结果与数值模拟结果吻合较好,表明在试验中形成了斜爆轰波。研究结果证明了斜爆轰发动机的技术可行性。     

双管爆轰单引射增推机理的数值研究

基于二维Euler方程,结合高精度Roe/HLL(Harten-Lax-Van Leer)混合格式与自适应网格加密技术,对引射模态下的双管爆轰非定常单次引射复杂流场进行了数值模拟.结果分别描述了方管爆轰波喷出及引射管内二次爆轰波的传播过程.双爆轰波喷出爆轰管后,以欠膨胀射流形式向外传播,同时,它们相互作用延长了引射管尾...     

气相爆轰波传播特性的数值模拟及实验对照

本文应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波进行了二维数值模拟.氢氧混合物的化学反应模型考虑了8种组分20个反应方程式.在处理化学反应引起的刚性问题时采用了时间算子分裂的方法.本文首先对爆轰波数值结果和实验结果进行了对照验证,然后对爆轰波在楔面反射由马赫反射向规则反射转变的过程进行了数值分析,得到了反射转变临界角,并和实验结果及理论分析进行了比较,结果是令人满意的;本文还对爆轰波的多波结构进行了初步的数值分析.     

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

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

预爆管爆震波衍射特性数值研究

为了实现预爆管出口爆震波向主爆管的成功传播,针对预爆管出口不同几何形状,开展了数值模拟研究,获得了爆震波从预爆管向主爆管传播的衍射特性.研究表明:当预爆管出口为直喷管时,不能实现向预爆管出口截面上游传播的爆震燃烧;当预爆管出口扩张角为45°时,可以实现爆震波在主爆管内向衍射面上、下游传播,研究结果获得的爆震波衍射特性得到了试验验证.      

斜向环形激波聚焦诱导直接起爆的大涡模拟

基于大涡模拟(LES)数值方法,利用基元反应模型对斜向环形激波聚焦诱爆氢气-空气可燃混合气直接起爆进行了详细的数值模拟.研究结果表明:在进气马赫数为2.5的条件下,斜向环形激波聚焦产生的高能区可以在可燃气体中实现直接起爆形成爆震波,并且燃料能够高效完全反应;而垂直环形激波则没能诱导出爆震燃烧过程.可见,斜向环形入射比垂直环形入射更加有利于燃料混合气聚焦诱导起爆成爆震燃烧.      

过渡段对爆震波衍射过程的影响分析

为了研究过渡段对爆震波衍射过程的影响,以氢气和氧气混合物为例,对不同过渡段张角下爆震波的衍射过程进行了数值模拟,结果表明,爆震波会在过渡段壁面产生反射激波,其强度随张角的减小而增大,同时,反射激波能将传入过渡段的爆震波分为3段,并在分界处产生高温、高压点,对爆震波从预爆管向主爆管的传播具有加强作用。模拟结果对预爆管起爆主爆震管实验中过渡段的选型有一定的参考意义。     

激波聚焦诱导气液两相爆震燃烧的数值模拟

对以激波聚焦和增加障碍物方式诱导煤油-空气气液两相爆震燃烧的过程进行了数值模拟.采用欧拉-拉格朗日方法建立了脉冲爆震发动机(PDE)中气液两相流的喷射、雾化、掺混过程.研究发现环形爆震波在爆震管凹腔内经过反射、汇聚后能够引燃可燃混合物.而在障碍物处,激波的反射和再反射聚焦能够形成高温高压点(2700K,25MPa),产生局部爆炸,有助于形成稳定的脉冲爆震燃烧(波面速度为1900m/s,温度为2 950K),有效地缩短由缓燃向爆震转变(DDT)距离至0.45m.     

障碍物对脉冲爆震发动机性能影响的数值模拟

运用计算流体力学方法计算了理想脉冲爆震发动机（PDE）单次爆震的性能参数,并和Wintenberger半分析模型进行了对比,研究了圆环形障碍物阻塞比和节距内径比对空气流动损失和脉冲爆震发动机单次爆震性能的影响,并与连续的多脉冲爆震实验的相关结果做了比较。研究结果表明,数值模拟计算得到的冲量和Wintenberger半分析模型得到的冲量非常接近,佐证了数值模拟方法和结果的正确性;障碍物阻塞比增大、节距内径比减小或进口速度增大,空气流过障碍物时的总压恢复系数减小,流动损失增大;障碍物阻塞比增大,单次爆震的体积比冲减小,单位燃油消耗率增大,性能降低;当阻塞比为41％时,体积比冲为理想爆震管体积比冲的92％,损失了8％;障碍物节距内径比增大,单次爆震的体积比冲先减小后增大,而单位燃油消耗率先增大后减小,不同节距内径比的体积比冲为理想爆震管体积比冲的94％左右,不同节距内径比对性能影响规律与多次脉冲爆震实验结果趋势符合。     

弯曲方管内爆燃向爆震转变特性

为缩短爆震管轴向尺寸,采用曲管爆震管替代直管爆震管.基于弯曲方管实验器,进行了二维数值仿真计算,研究了点火位置、点火能量及弯曲段曲率半径对爆震波形成及传播的影响规律.数值模拟结果表明:爆燃向爆震转变(DDT)时间随弯曲段曲率半径增大而减小,DDT距离在弯曲段曲率半径最大时达到最小值;点火能量越大,爆震管内激波强度越大,DDT时间越小;在距离爆震管进口300mm处点火时,爆震管内爆震波压力峰值最大,DDT时间最短.     

六管气动阀PDE协调工作研究

为了得到六管气动阀脉冲爆震发动机(PDE)的协调工作,在共用进气道条件下,主要采用离子探针技术进行了试验研究,同时对爆震波激波后燃烧波的传播速度进行了定量测量,并用压力传感器进行了试验验证.研究结果表明:在单管频率f＝30Hz(总频率f＝90Hz)时,六管PDE能够协调、稳定的工作;采用局部蒸发器改善燃油蒸发后,爆震管...     

特征线法在脉冲爆震发动机性能计算中的应用

在x-t平面建立爆震波在爆震室内传播过程的物理模型,应用特征线算法对脉冲爆震发动机(PDE)一维非定常流进行了数值计算,在此基础上给出了脉冲爆震发动机的性能参数的计算方法。并举算例进行了计算和分析,结果表明应用特征线算法计算PDE性能是可行而且准确的。      

球形点火腔对脉冲爆震DDT过程影响的数值研究

为了减小脉冲爆震发动机PDE（pulse detonation engine）单次循环时间，缩短爆燃向爆震转变DDT（deflagration to detonation transition）距离是关键.为此，提出了一种在头部添加球形点火腔的新型爆震管道结构设计方案，并采用丙烷和空气为可爆混合物，通过对新型爆震管DDT过程的二维数值模拟，研究其对DDT距离以及DDT时间的影响.数值模拟结果表明，头部采用球形点火腔后，点火腔中的压缩波经过多次反射后，能够在爆震管中更快地促使爆震波形成，DDT距离大大缩短；当球形点火腔直径为1.5倍等直爆震管直径时，相对于常规等直管爆震管结构，其DDT距离和过程时间分别减少了14%和16.26%.      

脉冲爆震发动机的飞行性能计算与分析

为了获得脉冲爆震发动机的飞行性能参数随飞行高度和飞行速度的变化规律,在应用最小G ibbs自由能理论建立求解爆震燃烧的平衡成分及其浓度、平衡参数的基本关系式,据此对爆震波特性参数进行数值模拟的基础上,通过引入进气道和尾喷管建立了一定飞行状态下的多管自吸气式脉冲爆震发动机理论性能计算模型,对脉冲爆震发动机的性能随飞行高度、速度的变化趋势进行研究。结果表明,在不同飞行高度和飞行马赫数的状态下,脉冲爆震发动机性能有一定的变化规律,同时模型也考虑了加入尾喷管对性能的影响,为脉冲爆震发动机的性能分析与理论设计提供了依据。