环截面管中的爆燃转爆震实验

为了实现在环截面爆震管封闭端同时点火,采用了先在预燃室中点火,再通过八路支管将火焰引入到爆震管封闭端,各个支管长度一样,因此可看成在环管封闭端多点同时点火.实验结果显示,在爆震管封闭端火焰汇聚形成激波点火,可在环管封闭端附近直接产生爆震波.在爆震管中安放扰流片时,流动阻力使爆震速度偏低,无扰流片时,爆震速度接近理论值.不加扰流片当量比在1.0左右,火焰速度会出现爆燃或爆震两种随机状态,当量比到1.1时可避免这种状况.      

曲面激波诱导斜爆轰的数值模拟

为了研究恰当比预混氢气-空气斜爆轰流场的波系结构和流动特征,基于带化学反应的Navier-Stokes方程,对弹头及楔穿越预混气体时诱导的斜爆轰进行了数值模拟。对流项的离散采用Steger-Warming格式,时间项采用二阶Runge-kutta方法。结果表明,对于弹头：（1）在亚爆轰条件下,能够模拟氢气-空气预混诱导爆轰流场的精细结构;（2）在超爆轰条件下,通过精细调整网格,能够很好地分辨强烈耦合的激波和燃烧波,且与Lehr实验吻合良好。对于楔结构：捕获到了清晰的三波点及其复杂精细的斜爆轰流场结构,预测的诱导燃烧距离、激波角和斜爆轰角与实验吻合良好;通过对流场波系结构变化过程的研究,获得了流场三波点随时间的演化过程。     

旋转爆轰的三维数值模拟

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

旋转爆轰流场的数值模拟

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

水下爆轰燃气泡形态与激波传播过程研究

针对脉冲爆轰发动机在水下工作过程中形成的燃气射流问题,搭建了水下爆轰燃气实验系统,研究了第一个爆轰循环中燃气泡发展变化过程.建立了基于气液两相双流体模型的脉冲爆轰发动机水下喷射模型,采用时-空守恒元和求解元方法,模拟了爆轰波与水相互作用形成的激波的传播及衰减过程.研究结果表明:燃气射流冲击水面时,燃气泡形态在膨胀阶段受...     

预混超声速气流斜激波诱导脱体爆轰研究

为了研究高静温预混超声速气流中爆轰的机理问题,利用自主设计的连续式加热器来产生马赫数为2.6,静温大于700 K的氢/空气预混超声速气流,并采用高速纹影技术来观测斜激波诱导脱体爆轰的直接起爆和发展的动态过程.研究表明,实验得到的起爆前斜激波和起爆后脱体爆轰波的角度与理论分析结果非常一致.并发现由激波诱导爆燃向激波诱导脱体爆轰的转变过程非常快;而且对于当前初温较高当量比较低的混合物,会出现爆轰波突然熄灭与重新起爆的现象.当前研究为爆轰燃烧在高超声速推进中的应用提供了重要参考.      

凹面腔内激波聚焦起爆爆震波过程的数值模拟

运用CFD方法和基元反应机理对环形向心射流产生的激波在凹面腔内的反射聚焦起爆爆震波过程进行了数值模拟,并根据流场分布及变化情况分析了激波聚焦及其在H2-Air混合物中起爆爆震波的全过程。对于本文建立的模型,起爆点在抛物形壁面的底部顶点处,聚焦起爆爆震波后的瞬间压力达到21.3MPa,温度达到4540K。爆震波在凹面腔内向开口端传播过程中仍会出现两次聚焦,压力达到约18MPa,温度达到4000K左右。研究结果表明:利用环形向心射流产生的激波在凹面腔内聚焦可以成功直接起爆爆震波,是一种有效的爆震直接起爆方法。     

激波聚焦火焰聚心起爆器的数值研究

针对短距离内的激波聚焦起爆问题,设计了一种激波聚焦火焰聚心起爆结构,通过数值模拟方法对该结构内部流场变化进行研究。结果表明:加速区是火焰失稳加速的重要区域;聚焦区能量汇聚造成局部爆炸是燃烧转爆轰加速的关键因素,且该区域内由于几何结构形成的卷吸涡能够增强聚焦并阻碍激波回传入加速区;稳燃区主要实现爆轰波从过驱状态到稳定爆轰状态的过渡区。确定中心轴线上温度、压力波动是由于聚焦腔聚焦作用造成的,而爆轰波锋面温度压力变化是由于三波点运动造成。      

气膜冷却对激波与边界层相互作用影响的数值模拟

针对跨声速平面叶栅中气膜冷却对流场的影响,采用数值模拟的方法,分析了激波和边界层的相互作用及引入气膜冷却之后三者之间的影响。结果表明,由于激波形成的逆压力梯度导致边界层出现分离现象,在引入冷却射流以后被部分抑制,流场细节显示在原分离处新形成了两个方向相反的分离旋涡。保持冷却条件不变,随着孔间距的减小,边界层分离现象被抑制的效果更加明显,平面叶栅热力损失系数逐渐减小。当孔径和孔间距之比达到0.67时,相对于没有引入气膜冷却的情况,热力损失系数降低了13%。冷气流量对射流和主流相互作用流场影响显著,冷气出口局部超声速区域显著增大流场损失,降低冷却效果。     

激波诱导矢量喷管流场的数值模拟及试验

为研究固体火箭发动机激波诱导矢量控制效率的影响因素,及验证数值模拟方法的准确性,对激波诱导轴对称推力矢量喷管模型进行了壁面测压试验,采用二阶精度Roe格式和k-ωSST两方程湍流模型求解强守恒型Navier-Stokes方程对矢量喷管复杂干扰内流场进行数值模拟。根据试验和数值模拟结果分析了喷管内主流和次流相互作用产生的复杂流场结构,比较了在不同喷管落压比NPR和次主流压力比SPR下喷管壁面静压的分布情况。结果表明,数值计算和试验结果基本吻合,验证了计算方法的准确性;在一定范围内减小喷管落压比,增大次主流压力比可以增大喷管周向壁面静压差,提高喷管的推力矢量偏角。     

环形火焰引发爆震的数值研究

采用改进的适于多组分反应流的波传播算法,对直管道中甲烷．空气预混气的爆震引发过程进行数值研究,考虑了两种不同构型的初始环形火焰对引发爆震的影响。对这两种构型下爆震引发过程中的压力场、温度场、反应放热率以及组分浓度进行了计算,由此探讨了“热点”形成机制及其对爆震引发的作用。结果表明,火焰发展产生的燃烧压缩波在管壁或管轴心发生汇聚、碰撞所形成的复杂波系与火焰相互作用,导致火焰前锋与压缩波之间的局部区域反应剧烈,反应放热率明显增加,从而促使热点形成。形成的热点迅速放大成为过驱爆震波,并逐渐衰减为稳定爆震波。稳定爆震波阵面参数与CJ理论结果进行了对比,两者符合得较好。     

障碍物对激波触发爆震波影响的数值研究

为了研究不同几何形状障碍物对激波触发爆震波的影响,对氢气/空气混合物在不同马赫数的激波作用下的爆震波触发过程进行了数值模拟与分析。研究发现,对于不同几何结构的障碍物,触发爆震波的机理不同,热点的产生可以是激波与激波的相互作用,也可以是激波与壁面的相互作用。对于堵塞比较小的障碍物,合理的结构设计可以使成功触发爆震波的激波马赫数进一步减小。     

二维爆震波数值模拟

针对爆震发动机的多维效应,利用CFD计算软件对单管脉冲爆震发动机进行了二维的数值模拟。以氢气和氧气混合物为例,仿真了爆震波在管内的传播过程以及传出管外后的流场情况。模拟所得的爆震波参数与文献中的计算结果和实验值基本一致,验证了模拟的正确性。爆震波传出爆震管后退化成一道激波,对该激波的流动情况进行了分析。     

跨声速和超声速流中激波/边界层干扰数值模拟

对SST湍流模型中的Bradshaw常数a1进行了修正,并对跨声速和超声速流中激波/边界层干扰进行了数值模拟研究,空间离散采用二阶精度差值的低耗散通量分裂格式(LDFSS),时间离散采用对称高斯-赛德尔(SGS)算法。结果表明:在跨声速流动中,计算得到的壁面压力分布、分离区长度和速度剖面都与实验值吻合较好,而且很好地模拟了典型的λ激波结构;在超声速流动中,修正后模型的计算精度较原始模型有了较大改善,计算得到壁面压力分布和分离点的位置都和实验值吻合较好。     

环形射流喷口位置对激波聚焦起爆的影响分析

为了研究环形射流喷口位置对激波聚焦起爆爆震的影响,以氢气和空气混合物为例,对不同环形射流喷口位置条件下的激波聚焦起爆过程进行了数值模拟。结果表明,随着喷口距凹壁面距离L增大,直接聚焦起爆和凹壁面反射起爆所需的临界入射压力逐渐增大,相应的聚焦时刻会滞后,形成焦点的位置距凹壁面越远。     

爆震波点火器管路传播过程数值模拟

为了进一步研究液体火箭发动机爆震波点火器的工作过程，结合实际试验发动机用爆震波点火器的工作条件，在混合比为3，初始氢氧混气压力为0．1MPa情况下，采用常用的两步显式MacCormak差分格式和有限速率化学反应动力模型，对爆震波点火器管路中爆震波的传播过程进行了数值模拟，研究了爆震波形成后在一维管路中传播特性及性能参数。从数值模拟结果来看，稳定传播的爆震波温度可达3500K，压力达到2MPa，其传播速度可达3430m／s，从爆震波形成到其传播到2m距离的出口，用时不到0．6ms。并与试验结果进行了对比分析，数值模拟结果和试验结果符合很好，为爆震波点火器的工程设计提供了理论依据。     

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.