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

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

CE/SE方法在非定常爆轰计算中的应用

建立了非定常爆轰波的理论模型。根据CE/SE方法的基本思想，推导出适合于求解二维N-S方程的长方形网格的CE/SE方法。对于包含刚性反应源项的守恒方程组，将其分为两部分，其中刚性反应源项单独求解。本文以爆轰管内装填C3H8/O2预混气体为例，数值模拟管内爆轰波的传播过程以及爆轰波脱离爆轰管后流场的图谱。计算结果显示，爆轰波在管内传播时，其波系极其复杂。当爆轰波离开出口后，迅速熄灭为非反应的激波。管外不仅有激波、膨胀波的传播，而且还存在激波与涡的相互作用。数值实践表明，该CE／SE方法可以有效捕获爆轰波与激波等强间断。     

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

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

气云爆轰

气云爆轰实验研究在一根专门设计的管中进行，通过对压力和火焰阵面的测试，研究激波的成长与发展过程，气云爆轰的爆速及其影响因素以及激波扣的点火诱导时间。理论研究考虑了液滴在高温，高压和高速气流中的变形，剥离和破碎，液滴的点火和局部爆炸以及气云爆轰的Ｃ－Ｊ判据。     

脉冲爆轰发动机内三维两相爆轰的数值计算

根据脉冲爆轰发动机内脉冲爆轰的特点,建立三维两相爆轰的理论模型。应用国际上最近发展起来的CE/SE方法,数值模拟爆轰管内三维两相燃烧转爆轰的过程。计算结果表明:当点火起爆后,在燃烧转爆轰过程中,三维效应非常明显。当形成稳定爆轰波后,三维效应不明显。计算压力值和实验压力值符合得较好。     

基于符号处理的天平载荷计算方法研究

风洞试验需要使用大量的天平，而经常编写天平公式处理程序给试验处理软件的维护带来了较大的工作量，影响了数据处理的效率和准度。随着试验软件通用性的不断提高，需要研究一种通用的天平公式处理方法。基于符号处理的天平载荷计算方法是一种解决天平公式的自动化处理方法，它运用系统软件的编译原理，通过对公式符号的解释和执行来实现对公式的处理。该方法已运用在风洞试验的数据处理中，可以在风洞试验的数据处理中推广应用。     

脉冲爆轰发动机预爆轰点火数值模拟

本文采用二维数值模拟对脉冲爆轰发动机点火方式中的预爆轰点火进行了分析研究。模拟了预爆轰管在主爆轰管上的两种最基本的布置方式：平行布置方式和正交布置方式。讨论了预爆轰管的布置方式和几何形状对点火的影响。结果表明：正交布置方式优于水平布置方式，因为反射激波在点火过程中起到了非常重要的作用。预爆轰管和主爆轰管的内径比是影响点火的重要参数，而长度的影响则较小。     

高超声速发动机不同燃烧模式的性能比较—斜爆轰发动机性能评价

本文采用热力学分析方法，从熵和Ｙｏｎｇ的概念出发，论证了理想情况下爆轰发动机的有效性，并定量比较了在同样入口压缩条件下，采用爆轰、等压、等Ｍ数等三种燃烧模式的发动机所生能，证明Ｃ－Ｊ斜爆轰是最佳超声速燃烧模式。     

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

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

扩张管和燃料部分填充对PDE效能的影响

采用基于有限速率基元化学反应模型的网格自适应有限体积方法,数值模拟以脉冲爆轰发动机为背景的爆轰波通过扩张喷管的流场,研究了喷管形状及喷管充气状况对脉冲爆轰发动机推进性能的影响。结果表明,氢-氧混合物产生爆轰波的等截面爆轰管冲量在某一比例合适的部分填充时达到最大,而对较大扩张角的喷管来说,其冲量一般随喷管内可燃气体填充比例的增加而增加,比冲通常随可燃气体填充比例的减少而增加;合适的喷管和燃料填充系数能使比冲显著提高。     

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

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

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

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

环形爆震波聚焦起爆

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

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

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

高温区碰并诱导气相爆轰二次起爆的数值研究

针对心形管内气相爆轰二次起爆过程进行了数值研究。采用二阶附加半隐的龙格-库塔法和5阶WENO格式求解二维欧拉方程,采用基元反应来描述爆轰化学反应过程。结果表明:气相爆轰在心形管中传播,由于壁面几何结构的作用,爆轰波熄灭后会发生二次起爆。该二次起爆由高温区碰并所诱导,且在管道不同位置诱导了两处二次起爆。爆轰流场中存在未反应气囊、射流、爆轰波马赫反射、激波与火焰作用及其所诱导的不稳定性等现象。计算得到的二次起爆过程胞格结构演变与实验基本一致。     

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

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

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

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

爆轰波穿越惰性气团时的透射激波参数分析

爆轰波穿越惰性气团再进入可燃混合气体时,透射激波的参数主要受到三个基本过程的影响:头激波两次穿越惰性气团与可燃气体之间接触面导致的激波强度下降、放热反应中止、爆轰波后稀疏波对激波的衰减作用。放热反应中止是透射激波衰减的主要原因,不同种类的惰性气体对透射激波参数影响不大。     

前缘钝化和化学非平衡效应对斜爆震发动机进气道特性的影响

为了探究前缘钝化、化学非平衡效应对斜爆震发动机进气道的工作性能及出口温度边界层分布的影响,采用热完全气体、化学非平衡气体两种模型对不同顶板、唇口前缘钝化半径下斜爆震发动机进气道进行数值模拟,结果表明：相比基准进气道,钝化后进气道上壁面温度边界层厚度较大,下壁面温度边界层厚度较小;在化学非平衡气体模型下,顶板前缘钝化半径(R1)>2mm时进气道的顶板附近和分离区内离解反应较为明显,唇口前缘钝化半径(R2)>2mm时进气道的唇罩、唇口板附近离解反应较为明显;当钝化半径≥4mm时,两种气体模型下进气道出口总压恢复系数和静温的相对变化量绝对值大于0.5%,有必要考虑化学非平衡气体效应对进气道出口性能的影响。     

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

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