On the role of turbulence in detonation induced by Mach stem reflection

A series of experiments conducted by Chan has shown that while some shock waves may not be strong enough to induce detonation when they collide with an obstacle the resulting Mach stem will induce detonation if it collides with a subsequent obstruction. A series of numerical simulations, however, failed to demonstrate the expected results if either the Euler or laminar Navier-Stokes equations are solved. On the other hand, calculations using the Favre averaged Navier-Stokes equations with a k--F turbulence model are able to reproduce the experimental results, indicating that turbulent effects may play an important role in the ignition process. A detailed examination of the results shows that turbulence causes the formation of activated kernels in a similar process to that observed in deflagration-detonation transition. The simulations in this paper have been undertaken using a modern high resolution hydrocode and a reduced kinetics mechanism for hydrogen combustion. The paper describes the reduced mechanism, the solution methods employed in the hydrocode and discusses the results of the simulations and their implications. Received 28 October 1997 / Accepted 30 April 1998     

氢氧爆轰波在变截面扩张管道中传播的数值模拟

采用二阶精度频散控制耗散格式(DCD)和8组分20个方程的基元反应模型,对轴对称变截面管道中氢氧爆轰波传播进行数值模拟。结果表明,爆轰波传播至突变截面扩张管道时,由于稀疏波的作用可能会使爆轰波局部熄爆甚至完全熄爆,对于某些敏感度高的反应气体爆轰波可以二次起爆。而在渐变截面扩张管道爆轰波相对不易熄爆。     

流动系统中爆轰波传播特性的数值模拟

基于带化学反应的二维Euler方程,采用氢气/空气的9组分19步基元反应简化模型,对充有当量 比的氢气/空气预混气体的矩形爆轰流场中爆轰波的传播过程进行了数值模拟,讨论了均匀来流对爆轰波传 播的影响。数值结果表明,在均匀来流的影响下,上游方向上的燃烧强度大于C-J爆速,下游方向上的燃烧强 度小于C-J爆速;上游方向传播的爆轰波的阵面压力大于下游方向传播的爆轰波的阵面压力。所以,经典的 C-J爆轰理论并不适用于流动系统中爆轰波传播特性的研究。     

Expansion of the detonation products of lead azide via a supersonic nozzle

The cloud of products formed following the detonation of lead azide (LA) contains gaseous species and solid particles. The dynamics of the detonation products expanding freely or through a supersonic nozzle into vacuum is unraveled via the temporal profiles of the pressure, the emission from Pb atoms and the attenuation of a He-Ne beam. The velocity of the fastest gaseous species is found from the onset of the pressure rise and the emission at a given distance from the LA sample, and the velocity of the fastest solid particles from the attenuation. In the free expansion, the respective velocities are 4.5±0.1 and 3.8±0.2 km/s and in the nozzle expansion 5.1±0.2 and 1.4±0.2 km/s. The expansion into atmospheric pressure air is also monitored and found to be much slower than that into vacuum. The utilization of nozzles as a means for obtaining a particle free, transparent medium of detonation products is stressed in the context of exploiting explosives for achieving chemical lasers in the visible wavelength region.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

单次脉冲爆轰发动机工作过程的数值模拟

通过求解14组分和19个基元反应的CH4-O2-N2详细化学反应动力学机理的二维轴对称Navier-Stokes方程,对爆轰管内半球形高温火团引发的爆轰过程和爆轰波进入外流场后的全流场分布进行数值模拟。模拟结果显示了爆轰波在管内的成长、稳定传播、进入尾喷管后衰减为激波和进入外流场的全过程,以及爆轰管出口端附近区域的复杂涡与激波的相互作用。对轴线上的压力分布和封闭端的压力等进行了讨论,为脉冲爆轰发动机的开发研制提供参考信息。     

On the origin of the double cellular structure of the detonation in gaseous nitromethane and its mixtures with oxygen

An experimental study of the detonation in gaseous nitromethane (NM) and nitromethane-oxygen mixtures has exhibited unambiguously the existence of a double cellular structure in the range of equivalence ratio from 1.3 to 1.75 (NM). Calculations of the reaction zone of the detonation in the same range of equivalence ratio, using a detailed chemical scheme in the ZND model, demonstrate that the chemical energy is released in two main successive distinct exothermic reaction steps characterized by their own induction length which justifies the existence of a two levels detonation cellular structure. This result strengthens the idea that the cellular detonation structure finds its origin in instabilities amplified by delayed local high energy release rate inside the reaction zone.Received: 2 February 2002, Accepted: 27 May 2004, Published online: 7 December 2004[/PUBLISHED]H.-N. Presles: Correspondence to     

Investigation of flows past a plate in the presence of exoergic processes driven by the interaction between a nonequilibrium electric discharge and a propane-air mixture

Supersonic flow (M = 2) past a plate along which propane is injected is investigated within the framework of the solution of problems of combustion initiation and stabilization at low static temperatures and pressures in the presence of a nonequilibrium discharge and metal and dielectric interceptors mounted on the plate surface. The experiments show that two zones with exoergic reactions develop when a metal interceptor is mounted on the plate. One zone is located ahead of the leading separation zone and the other above and behind the interceptor edge, its boundary partially penetrating into the supersonic flow region. Using modern spectroscopic methods, the radiation intensity distributions of a series of plasmochemical reaction products are obtained in the neighborhood of the plate ahead of the interceptor, behind it, and above its edge. It is found that the fuel is intensively converted under the action of the discharge with the occurrence of a series of free radicals, atomic hydrogen and oxygen which are themselves chemically active.     

Initiation and stabilization of combustion at low pressure and temperature by means of a nonequilibrium electric discharge

The pattern of supersonic (M = 2) flow near the surface of a plate at points on which propane jets are injected normal to the air flow direction is investigated. For the initiation and intensification of the chemical reactions a nonequilibrium discharge is used. This is created between an anode oriented along the flow, the plate surface, and a metal interceptor mounted on the plate. The results of a schlieren visualization of the flows developed are presented. Spectroscopic studies show that the distribution of plasmochemical reaction products has a number of fundamental differences as compared with the case of propane injection along the plate surface. A comparative analysis of these distributions for identical gasdynamical experimental conditions is important for testing calculation models of reactive-mixture supersonic flows in which electric discharges are used for ignition and the stabilization of combustion.     

氢燃料缓燃向爆震转捩过程中波与火焰的匹配特性

为了解氢燃料爆震过程中压力波与火焰之间相互匹配的特性,在60mm60mm2000mm 方 爆震管内,用氢气和空气混合物进行了单爆震性能研究。在爆震转捩区内布置压力传感器与离子探针,用来 监控压力波和火焰的信号,同时利用高速摄影仪集中拍摄转捩区域。根据压力波和火焰面在爆震管不同时刻 的强度特性、速度特性及位置特性来分析爆震过程中波与火焰匹配的规律。结果表明:压力波和火焰的强度 呈现为相互正反馈匹配性质;缓燃向爆震转捩(DDT)过程中,压力波和火焰的速度表现为相互交替的变化过 程,且缓燃阶段中火焰速度的增幅大于压力波速度的增幅;当火焰面追赶上激波时,产生过爆,火焰面会临时 位于激波前面;在过爆衰减为正常爆震波的过程中,激波在火焰前面。     

燃烧轻气炮多级渐扩型燃烧室流场特性数值研究

为了分析多级渐扩型燃烧室结构对燃烧轻气炮氢氧燃烧特性的影响, 通过计算流体力学方法, 分别对采用传统圆柱型燃烧室和多级渐扩型燃烧室的燃烧轻气炮氢氧燃烧发射过程进行数值模拟。对比结果表明, 多级渐扩型燃烧室结构能够明显地减小燃烧室压力波动幅度, 提高氢氧燃烧稳定性; 多级渐扩型燃烧室内形成回流区, 可以减小气流轴向运动速度; 火焰扩展形态与渐扩型结构相吻合, 燃烧反应区表面变化平稳; 多级渐扩型燃烧室结构对氢氧火焰传播过程和压力波动现象有着重要影响。     

点火过程和初始条件对燃烧轻气炮内弹道性能的影响

采用计算流体力学方法对燃烧轻气炮膛内燃烧过程进行数值模拟,分析不同的点火点数目和点火能量以及初始温度和压力对燃烧轻气炮内弹道特性的影响。结果表明,采用合理的点火点数目、初始温度和压力条件可以有效控制氢气的燃烧过程,减弱燃烧室内的压力波动。模拟结果对燃烧轻气炮膛内燃烧过程控制具有重要参考价值。     

Numerical study on three-dimensional flow field of continuously rotating detonation in a toroidal chamber

Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonationshock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interesting properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.     

缸内激波对锥顶型燃烧室的活塞破坏机理

以二维数值模拟为基础,研究了锥顶型燃烧室内的冲击波发展的震荡过程,得到作用于活塞不同位置处的超压分布。模拟结果表明:由于燃烧室结构的独特性,导致冲击波能在特定区域进行汇聚,致使该区域超压明显高于其他区域。将该模拟结果与实际破坏失效的活塞进行对比,发现冲击波汇聚区域往往就是活塞被破坏的地方。数值模拟结果与实际破坏结果吻合很好。这为设计燃烧室形状以避免冲击波对活塞造成破坏提供了理论基础。     

C2H4-O2混合气体直接起爆的临界能量

主要基于化学动力学和Ng模型,对C2H4-O2混合气体的爆轰胞格尺寸进行预测;结合Lee表面积能量模型,预测物质在不同初始压力和化学当量比的条件下,直接起爆引起球面爆轰的临界起爆能量。直接起爆实验主要采用高压电点火提供起爆能量,起爆能量通过放电过程中电流的输出信号确定。结果表明,理论预测与实验值较吻合。首先,通过化学动力学计算得出ZND模型的爆轰参数,利用Ng模型得出爆轰胞格尺寸与ZND诱导区长度之间的比例因数A 在不同初始压力与当量比的条件下分别为:A=43.815(1+p1/p0)-0.12371和A=8.531exp(/3.135)+28.644,在此基础上对爆轰胞格尺寸进行定量预测。胞格尺寸的理论预测与实验结果吻合。其次,把爆轰胞格尺寸作为中间特征参数并结合Lee的表面积能量模型,提出可以预测临界起爆能量的定量模型,并得出C2H4-O2混合气体直接起爆的临界起爆能量与初始压力和化学当量比的参量拟合关系分别为Ec=0.332(p1/p0)-2.017和Ec=exp[3.951(-1.401)2-1.9]。     

气/固/液三相混合物燃烧转爆轰过程实验研究

利用多相燃烧爆炸实验系统,通过高压喷粉/喷雾以及高能点火等过程,对化学当量比条件下3种 典型燃料空气炸药,即硝基甲烷/铝粉/空气、硝酸异丙酯/铝粉/空气、乙醚/铝粉/空气三相混合物的燃烧转爆 轰过程进行了实验研究,同时根据实验结果对比了3种三相混合物的燃爆性能。得到了三相悬浮混合物燃烧 转爆轰过程的宏观规律以及三相混合物燃爆性能随质量浓度变化的规律。     

直管内胞格爆轰的基元反应数值研究

基于基元反应和二维欧拉方程,对直管内胞格爆轰进行了数值模拟。采用5阶WENO(weighted essentially nonoscillatory scheme)求解对流项,采用2阶附加半隐的龙格-库塔法处理化学反应源相引起的刚性。获得了密度、压力、温度和典型组元质量分数流场及数值胞格结构等。结果表明：网格精度的差异明显影响胞格的规则性和爆轰的平衡模数,随着网格尺度的减小,胞格由不规则变为规则。预混气组成、初压、初温及管道宽度给定,三波点数收敛为确定值。足够强度的初始扰动可再现胞格爆轰,最终形成的自持胞格爆轰模数与初始扰动的形状、大小、位置均无关。沿胞格中心线,爆轰波速度变化范围为0.88DCJ~1.5DCJ,爆轰波平均速度与CJ爆轰速度仅偏差0.88％。峰值压力与初压之比为14~50。计算爆轰波平均速度、胞格宽长比与实验值基本一致,但计算胞格宽度比实验值略小。数值模拟加深了对横波的产生和发展、未反应气囊、爆轰胞格的二次起爆等胞格爆轰特性的认识。     

T型管内流动气体中爆轰绕射过程的数值模拟

基于带化学反应的二维Euler方程,对H₂、O₂、Ar体积比为2:1:1的混合气体系统在T型管内的爆轰绕射进行了数值模拟。用二阶附加半隐的Runge-Kutta法和五阶WENO格式分别离散欧拉方程的时间和空间导数项,采用9组分48步基元反应简化模型描述爆轰波在静止系统和流动系统中的传播过程,得到了温度、压力、典型组元H质量分数的分布及数值胞格结构。结果表明:在流动系统中,迎风面上波阵面为斜爆轰结构,静止系统两侧和顺风面上的波阵面为完全解耦的前导激波;在水平管中,波阵面与上下壁面经历一系列马赫碰撞后,最终形成正爆轰;在流动系统中,胞格结构明显向下游偏移;横向爆轰波的产生对爆轰波的再生起到了关键作用。     

动态存储方法在气相爆轰波数值模拟中的应用

为减少反应流计算花费的时间,采用动态存储建立热力学数据表（In situ adaptive tabulation, ISAT）的方法来取代反应流计算中化学反应的直接积分（Direct integral, DI）,并针对气相H2/O2爆轰波的传播过程进行了数值模拟。采用两种不同的ISAT误差放大判据对一维爆轰波的压力、温度及组分变化进行了计算,并和DI的结果进行了比较以考察其计算精度。此外,还探讨了ISAT方法的效率。计算结果表明,和DI方法相比,在爆轰波传播计算中ISAT方法计算误差小于3％,化学反应计算速度提高了8倍。     

爆轰波通过扩张喷管的双曝光全息实验和数值研究

结合实验和数值模拟方法,对以脉冲爆轰发动机为背景的爆轰波通过扩张喷管的流动进行了系列研究。实验采用双曝光全息干涉方法对爆轰波绕射流场进行测量,得到了比传统的纹影法更清晰和可定量化的照片。发展了基于非结构四边形网格自适应有限体积程序,结合基元化学反应模型对扩张喷管中爆轰化学反应流场进行了数值模拟,模拟结果与实验照片吻合较好。实验和数值模拟结果表明,爆轰波绕射具有许多和激波绕射不同的流场特征,其中包括二次起爆现象、化学反应面与前导激波相脱离而引起的复杂流场等,同时初始压力和扩张角度变化也对爆轰波绕射过程产生较大影响,初始压力越低,化学反应区和前导激波分离现象越明显,且前导激波的曲率越大。     

Large-Scale Streamwise Vortices in the Supersonic Part of a Permeable Nozzle

Large-scale streamwise vortices in the vicinity of a perforated wall in the supersonic part of the nozzle are studied. The governing effect of gas inflow through a perforated wall on origination and parameters of streamwise vortices is experimentally established.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 5, pp. 68–75, September–October, 2005.