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

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

气相爆轰基元反应模型数值模拟

采用基元反应模型和高精度的ENO格式对气相爆轰进行数值研究，对H2／O2／Ar混合气体起爆和爆轰波传播过程的数值模拟结果表明，计算的爆轰波阵面参数和实验相当符合。数值研究的结果还表明，爆轰波反应区中参与反应的不同组分具有不同类型的变化特征。计算结果的精度随着网格尺寸增加而增加，并能保持较好的收敛性。     

连续爆轰发动机的二维数值模拟研究

介绍了连续爆轰发动机(Continuous detonation engine,简称CDE)的物理模型及爆轰波传播的物理机理,并对CDE流场的波结构进行了分析.使用两步化学反应模型对简化的CDE模型做了二维数值模拟研究.计算了圆柱面展开二维平面的多个爆轰循环过程.模拟结果较好地匹配了相关CDE实验所获得的流场结构.数值计算结果使我们可以更加透彻地理解CDE各方面的参数,对于CDE的进一步研究具有一定的指导意义.      

旋转爆轰的三维数值模拟

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

准爆轰波传播机理的探索研究

本文分析了准爆轰波的传播过程与物理机制,提出了一个全新的准爆轰波物理模型,开展了理论分析与数值模拟研究.在该物理模型中考虑了超声速流动化学反应放热造成的动量损失,计算结果表明准爆轰波发展过程是热壅塞的.计算与实验结果的定性吻合表明了本文提出的理论模型是正确的,获得的计算结果是合理的,对于准爆轰波物理机理的讨论是有重要意义的.     

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

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

与爆轰相关的湍流燃烧

爆轰具有较高的燃烧效率,当应用于推进系统时,有较宽广的前景,故引起国内外学者广泛的关注。与爆轰有关的两个重要现象是——爆轰的生成和传播。爆轰的生成分为直接起爆和间接起爆。直接起爆需要极高的能量,因此间接起爆一直是研究的热点,其中爆燃转爆轰(deflagration-to-detonation transition,DDT)则是关注的重点。然而,经过几十年的实验、理论和数值研究,DDT的机理已被理解,但是准确的预测还无法实现,这其中最核心的因素是湍流燃烧。对于爆轰波的传播,特别是高活化能的可燃气体,其波阵面呈强不稳定结构,存在湍流燃烧,波后出现的未燃气团的燃烧也属于湍流燃烧,这使得爆轰波的精确控制也存在较大的困难,但是爆轰波传播的平均速度仍然可以用CJ速度来预测。本文对DDT和爆轰波结构中的湍流燃烧的研究现状进行了分析和评述,并对未来的研究进行了展望。     

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

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

旋转爆轰流场的数值模拟

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

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

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

共用尾喷管多管脉冲爆震发动机数值模拟研究

采用基元反应模型和显式迎风TVD差分格式数值模拟方法,研究共用尾喷管多管脉冲爆震发动机(Pulse Detonation Engine,简称PDE)各爆震管之间相互影响.采用点隐算法解决化学反应引起的刚性问题,模拟了共用尾喷管内的流场,以及爆震波退化为激波遇到尾喷管收敛斜面反射及反射波向爆震管上游传播的过程.研究结果为多管PDE的共用尾喷管设计提供了理论基础.     

高速飞行弹丸诱导的三维爆轰波数值模拟

采用三维非定常化学反应控制方程组对驻定在带攻角高速飞行弹丸上的斜爆轰波流场进行了数值模拟。对流项和化学反应项采用时间分裂法,分别用TVD格式和基元反应模型进行处理。计算结果表明三维数值模拟可以更准确的描述弹丸飞行攻角对爆轰波的形成及结构的影响,是进一步研究超驱爆轰推进技术的重要手段。     

连续旋转爆震波细致结构及自持机理

基于化学非平衡流解耦方法和详细化学反应机理,在0.1mm的网格尺度下对连续旋转爆震波内流场和爆震波平掠惰性气体界面的传播过程进行了数值模拟。详细分析了连续旋转爆震波的空间结构,由于在出口侧无几何约束,受一系列稀疏波的影响,旋转爆震波强度受到不同程度的削弱,形状也发生了相应的变化。爆震波平掠惰性气体界面时的流场分布与连续旋转爆震波结构极为类似,对比研究了波前可燃气体层高度对爆震波结构和传播过程的影响,当爆震波高度较小时,其抵御惰性气体侧稀疏波衰减的能力弱,波阵面严重变形,爆震波解耦。若爆震波要稳定自持传播,其必须要达到一定的临界高度来抵御稀疏波的影响。     

时空守恒元和解元方法的爆震波一维数值模拟

运用基元化学反应模型和流体力学的基本原理，建立了脉冲爆震发动机一维非定常反应流数学模型。应用一种改进的时空守恒元和解元(CESE)方法建立控制方程的差分格式，针对化学反应源项的刚性问题，在对非线性源项进行线性化处理的基础上，应用隐式梯形方法实现源项的积分。对氢氧爆震波的形成和传播过程进行数值计算，结果表明，应用CESE方法和基元化学反应模型能够成功地模拟氢氧爆震波的精细结构。     

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

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

脉冲爆震发动机喷管性能数值分析

采用详细的化学反应机理,对脉冲爆震发动机的不同结构喷管进行了二维轴对称单循环数值模拟.分析了爆震波在爆震室内的传播过程以及传出后的流场情况,模拟所得的爆震波参数与STANJAN值较为一致.对4种单喷管和4种组合喷管的性能模拟结果表明:扩张喷管可以产生最大的瞬时推力峰值和比冲;带引射的组合喷管产生的平均推力比单喷管产生的平均推力要大.最后对非稳态引射过程进行了分析,为提高脉冲爆震发动机的性能提供参考.      

进气温度和当量比对脉冲爆轰发动机工作过程影响

为了研究进气温度和当量比对脉冲爆轰发动机工作过程的影响,建立了带简单化学反应的气液两相爆轰欧拉-拉格朗日模型,使用二维时空守恒元与求解元(CE/SE)方法和变步长4阶龙格-库塔法分别求解气液两相爆轰方程.计算结果表明:提高进气温度,能加速液滴雾化、蒸发,缩短燃烧转爆轰距离和时间,但是会降低爆轰波的峰值压力;当量比小于1.1时,增加当量比,能缩短燃烧转爆轰距离和时间,提高爆轰波的峰值压力,加快爆轰波传播速度;当量比大于等于1.1时,增加当量比,能小幅提高爆轰波速传播度和缩短燃烧转爆轰距离和时间,但对爆轰波的峰值压力影响很小.      

斜爆轰波并行数值模拟的化学加速算法性能研究

考虑基元反应条件下的爆轰波精细结构的数值模拟计算量巨大,发展高精度和高效率的计算方法十分必要。本文以Ma马赫数为7的H2/O2/N2预混气来流形成的斜爆轰波为数值模拟对象,研究了一种基于并行计算架构的用于加速化学反应计算的建表算法的计算性能,考察了不同建表策略和不同并行分区数量对算法性能的影响。研究结果显示,本文采用的建表算法能够很好地再现斜爆轰结构,其计算精度不受建表策略和并行计算分区数量的影响；而算法的计算效率则取决于不同分区对应的数据表之间操作的同步性,其中,数据表中节点数据的取回率和数据表设定的尺寸上限都会影响数据表操作的同步性。本文采用的TP和PLP两种建表策略的计算结果表明,TP策略数据表中节点取回率高于PLP策略,故计算效率更高；而计算分区数量越少,则分区对应的数据表尺寸上限越大,数据表的同步性就越好,计算效率也越高。上述研究为发展适于气相斜爆轰波数值模拟的高保真和高效率的计算方法提供了有益参考。     

环形爆震波聚焦起爆数值模拟

对三种带不同结构形式的环形预爆管的爆震发生器进行了数值模拟,来研究环形爆震波聚焦起爆现象及其气动特性.数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动的NavierStokes方程来模拟化学动力学和流体动力学过程.研究发现用低的点火能量对环形预爆管中的燃料和氧化剂点火产生层流火焰,层流火焰在狭窄管壁的...     

环形爆震波聚焦起爆

本文对三种带不同结构形式的环形预爆管的爆震发生器进行了数值模拟,来研究环形爆震波聚焦起爆现象及其气动特性。数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动的Navier-Stokes方程来模拟化学动力学和流体动力学过程。研究发现用低的点火能量对环形预爆管中的燃料和氧化剂点火产生层流火焰,层流火焰在狭窄管壁的作用下完成爆燃向爆震转捩,形成环形爆震波,环形爆震波在聚焦腔经过反射、汇聚,最终引爆主爆震室中的可燃混合物;并对稳定爆震的距离和强度进行了对比研究。通过分析对环形爆震波衍射、反射、聚焦的气体动力学特性及复杂波系发展规律有了更深的认识,为进一步试验提供参考。