Effect of initial disturbance on the detonation front structure of a narrow duct

The effect of an initial disturbance on the detonation front structure in a narrow duct is studied by three-dimensional numerical simulation. The numerical method used includes a high-resolution fifth-order weighted essentially non-oscillatory scheme for spatial discretization, coupled with a third-order total variation diminishing Runge-Kutta time-stepping method. Two types of disturbances are used for the initial perturbation. One is a random disturbance which is imposed on the whole area of the detonation front, and the other is a symmetrical disturbance imposed within a band along the diagonal direction on the front. The results show that the two types of disturbances lead to different processes. For the random disturbance, the detonation front evolves into a stable spinning detonation. For the symmetrical diagonal disturbance, the detonation front displays a diagonal pattern at an early stage, but this pattern is unstable. It breaks down after a short while and it finally evolves into a spinning detonation. The spinning detonation structure ultimately formed due to the two types of disturbances is the same. This means that spinning detonation is the most stable mode for the simulated narrow duct. Therefore, in a narrow duct, triggering a spinning detonation can be an effective way to produce a stable detonation as well as to speed up the deflagration to detonation transition process.     

Pressure profiles in detonation cells with rectangular and diagonal structures

Experimental results presented in this work enable us to classify the three-dimensional structure of the detonation into two fundamental types: a rectangular structure and a diagonal structure. The rectangular structure is well documented in the literature and consists of orthogonal waves travelling independently from each another. The soot record in this case shows the classical diamond detonation cell exhibiting ‘slapping waves’. The experiments indicate that the diagonal structure is a structure with the triple point intersections moving along the diagonal line of the tube cross section. The axes of the transverse waves are canted at 45 degrees to the wall, accounting for the lack of slapping waves. It is possible to reproduce these diagonal structures by appropriately controlling the experimental ignition procedure. The characteristics of the diagonal structure show some similarities with detonation structure in round tube. Pressure measurements recorded along the central axis of the cellular structure show a series of pressure peaks, depending on the type of structure and the position inside the detonation cell. Pressure profiles measured for the whole length of the two types of detonation cells show that the intensity of the shock front is higher and the length of the detonation cell is shorter for the diagonal structures. Received 17 May 2000 / Accepted 29 November 2000     

Numerical simulation of spinning detonation in square tube

A single spinning detonation wave propagating in a square tube is simulated three-dimensionally with the detailed chemical reaction mechanism for hydrogen/air mixture proposed by Petersen and Hanson. The spinning detonation is composed of a transverse detonation rotating around the wall normal to the tube axis, triple lines propagating partially out of phase, and a short pressure trail. The formation of an unburned gas pocket behind the detonation front was not observed in the present simulations because the rotating transverse detonation completely consumed the unburned gas. The calculated profiles of instantaneous OH mass fraction have a keystone shape behind the detonation front. The numerical results for the pitch and track angle on the tube wall agree well with the experimental results. This paper is based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems at Montreal, Canada, from July 31 to August 5, 2005.     

爆轰波在阻尼管道中声吸收的实验研究

实验旨在研究气相爆轰波在阻尼管道 (管壁上衬有吸收材料 )中传播时的衰减现象。先是在光滑管壁的管道中产生稳定的具有胞格结构的爆轰波 ,然后使其通过专门设计的管壁上衬有吸收材料 (金属丝网或不锈钢纤维 )的阻尼段。利用高速摄影、压力传感器和烟迹技术等手段 ,记录和测试了阻尼段对几种混合气体爆轰波的传播速度、爆压及胞格结构产生的影响。实验分别在方管和圆管中进行。发现在某些条件下爆轰波可以被衰减成强爆燃。     

三维气相爆轰动态并行计算程序设计与开发

在三维气相爆轰数值研究中,网格精度和计算域的规模导致网格数占有非常庞大的计算资源,进而给数值模拟带来了极大的挑战。本文针对这一难题,采用5阶WENO格式对带化学反应Euler方程组进行空间离散,基于MPI（MessagePassingInterface）并行模式开发了高精度动态并行代码,并对爆轰波在带有障碍物的三维方形管道中的传播过程进行计算。计算结果表明,高精度动态并行计算能够很好的模拟三维气相爆轰波在大尺寸管道中的传播,不仅提高了计算效率,而且提高了爆轰波阵面的分辨率。与高精度静态并行相比,高精度动态并行计算减少了界面数据通信时间,从而进一步提高了计算效率。因此,高精度动态并行程序为探究三维气相爆轰新的物理机制提供有效的手段。     

Numerical investigation on evolution of cylindrical cellular detonation

Cylindrical cellular detonation is numerically investigated by solving two- dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.     

Cellular detonations in bidispersed gas-particle mixtures

Formation of cellular detonation in bi-fractional stoichiometric mixtures of aluminum particles and oxygen is investigated numerically. The detonation cell size depends on the particle diameters and relative concentration of the fractions. Certain degeneration of cellular detonation is obtained when compared to the monodisperse mixtures. It is characterized by maximal pressure decrease, transverse wave relaxation and detonation front rectification. Complete degeneration of cellular detonations and stable propagation of a plane detonation front is found in some bi-fractional mixtures. The numerical results are confirmed by acoustic analysis of the detonation structures. This paper is based on work that was presented at the 21st International Colloquium on the Dynamics of Explosions and Reactive Systems, Poitiers, France, July 23–27, 2007.     

Propagation mechanism of dust detonations

This paper summarizes the studies on dust detonations at the Stosswellenlabor of RWTH Aachen since 1987. The onset and propagation mechanism of heterogeneous dust detonations are similar to those of marginal gas phase detonations. A self-sustained dust detonation has transverse wave structures that provide the coupling between shock and reaction. Large transition distances and transverse wave spacings require large sized tubes for the propagation of self-sustained dust detonations. The Hugoniot analysis of the Chapman-Jouguet detonation predicts equilibrium detonation states being in reasonable agreement with the self-sustained dust detonations observed. Shock matching calculations at the triple point adequately determine the wave structures of those stable dust detonations.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.     

气相爆轰波反应区结构的平面激光诱导荧光测量

基于平面激光诱导荧光(PLIF)技术对2H2 O2 10Ar的预混气体爆轰波反应区结构进行实验研究.采用高浓度的氩稀释有利于减小爆轰化学反应自发辐射光对OH荧光的影响.合理设置PLIF系统、爆轰波和ICCD之间的同步控制触发延时,得到爆轰波阵面附近的OH荧光分布图像.结果表明:诱导激波后反应阵面不是平面且不稳定.荧光图像上能清晰地看到类似拱顶石的结构,它位于两马赫杆之间,以入射激波、剪切层和反应阵面为边界.无论是在马赫杆后还是在入射激波后,OH浓度分布在诱导区末端急剧增加至最大值.随着离开反应阵面的距离增加,OH浓度快速减小.由于爆轰模式和激光片光方向的影响,从PLIF图像上测得的横波间距值较离散,均小于胞格宽度.     

旋转爆轰的三维结构和侧向稀疏波的影响

基于带化学反应的三维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对圆环形燃 烧室内的旋转爆轰进行了数值模拟,讨论了旋转爆轰波的三维结构及侧向稀疏波对旋转爆轰波阵面的影响。 数值结果表明,爆轰波能够以旋转方式沿预混气层稳定传播。在侧向稀疏波作用下,爆轰波阵面发生变形。 与理想的C-J爆轰相比,爆轰波强度和爆轰参数都有所下降。     

Investigation of the detonation regimes in gaseous mixtures with suspended starch particles

The existence of a secondary discontinuity at the rear of a detonation front shown in experiments by Peraldi and Veyssiere (1986) in stoichiometric hydrogen-oxygen mixtures with suspended 20-m starch particles has not been explained satisfactorily. Recently Veyssiere et al. (1997) analyzed these results using a one-dimensional (1-D) numerical model, and concluded that the heat release rate provided by the burning of starch particles in gaseous detonation products is too weak to support a double-front detonation (DFD), in contrast to the case of hybrid mixtures of hydrogen-air with suspended aluminium particles in which a double-front detonation structure was observed by Veyssiere (1986). A two-dimensional (2-D) numerical model was used in the present work to investigate abovementioned experimental results for hybrid mixtures with starch particles. The formation and propagation of the detonation has been examined in the geometry similar to the experimental tube of Peraldi and Veyssiere (1986), which has an area change after 2 m of propagation from the ignition point from a 69 mm dia. section to a 53 mm 53 mm square cross section corresponding to a 33% area contraction. It is shown that the detonation propagation regime in these experiments has a different nature from the double-front detonation observed in hybrid mixtures with aluminium particles. The detonation propagates as a pseudo-gas detonation (PGD) because starch particles release their heat downstream of the CJ plane giving rise to a non-stationary compression wave. The discontinuity wave at the rear of the detonation front is due to the interaction of the leading detonation front with the tube contraction, and is detected at the farthest pressure gauge location because the tube length is insufficient for the perturbation generated by the tube contraction to decay. Thus, numerical simulations explain experimental observations made by Peraldi and Veyssiere (1986). Received 5 July 1997 / Accepted 13 July 1998     

气相爆轰波近失效状态的传播模式

选用五种碳氢混合气体,采用高压电火花起爆的方法,利用光纤探针测量爆轰波在管道内的传播速度,研究临近失效状态时爆轰波在管道内的传播模式。实验是在自行研制的爆轰管道中进行的,其包括驱动段及内径分别为1.5、3.2、12.7 mm三种规格的测试段。实验结果再次验证了爆轰波在管道内传播时可以有六种不同的传播模式:稳态爆轰、快速波动爆轰、结巴式爆轰、驰振爆轰、低速爆轰和爆轰失效。其中C₂H₂+2.5O₂+70% Ar、C₂H₂+2.5O₂+85% Ar两种组分混合气体（具有较低活化能）,在爆轰波传播过程中只有稳态、快速波动和失效三种模式;而C₃H₈+5O₂、C₂H₂+5N₂O和CH₄+2O₂三种组分混合气体（具有较高活化能）在传播过程中出现六种不同模式。上述结果表明,除气体组分、初始压力等因素外,混合气体的活化能可能对爆轰波在管道内的传播状态也有影响。     

Numerical simulation of H2/air detonation using unstructured mesh

To explore the capability of unstructured mesh to simulate detonation wave propagation phenomena, numerical simulation of H₂/air detonation using unstructured mesh was conducted. The unstructured mesh has several adv- antages such as easy mesh adaptation and flexibility to the complicated configurations. To examine the resolution dependency of the unstructured mesh, several simulations varying the mesh size were conducted and compared with a computed result using a structured mesh. The results show that the unstructured mesh solution captures the detailed structure of detonation wave, as well as the structured mesh solution. To capture the detailed detonation cell structure, the unstructured mesh simulations required at least twice, ideally 5times the resolution of structured mesh solution.      

Behavior of detonation waves at low pressures

With respect to stability of gaseous detonations, unsteady behavior of galloping detonations and re-initiation process of hydrogen-oxygen mixtures are studied using a detonation tube of 14 m in length and 45 mm i.d. The arrival of the shock wave and the reaction front is detected individually by a double probe combining of a pressure and an ion probe. The experimental results show that there are two different types of the re-initiation mechanism. One is essentially the same as that of deflagration to detonation transition in the sense that a shock wave generated by flame acceleration causes a local explosion. From calculated values of ignition delay behind the shock wave decoupled from the reaction front, the other is found to be closely related with spontaneous ignition. In this case, the fundamental propagation mode shows a spinning detonation. Received 10 March 1997 / Accepted 8 June 1997     

周期性非均匀介质中气相爆轰波演变模式研究

气相爆轰波在周期性非均匀介质中的起爆, 稳态传播和失效机制都极为复杂, 很多物理机制尚不明确, 是当前爆轰物理领域研究的热点和难点. 本文使用反应欧拉方程和两步化学反应模型对爆轰波在非均匀介质中的传播机理进行了数值模拟研究, 非均匀性由横向周期性分布的温度扰动体现, 重点分析不同波长、不同幅度的温度扰动对波阵面波系结构的影响. 计算结果表明, ZND爆轰波在温度扰动下向胞格爆轰波的转变主要受制于两种竞争性因素: 一是爆轰波内在的不稳定性; 二是温度扰动的波长和幅度, 前者是内因, 后者是外因. 温度扰动的存在抑制横波的发展, 延迟了ZND爆轰波向胞格爆轰波的演化, 并且内在不稳定性的增加可以减慢这种延迟现象. 这说明, 温度扰动可以在一定的范围内抑制胞格不稳定性的发展, 但是不能够终止这一过程. 温度的不连续性使得爆轰波阵面更为扭曲, 并在横波附近存在较弱的三波点结构, 即温度扰动可增加爆轰波固有的不稳定性, 改变爆轰波阵面的传播机理. 幅值较大的人工温度扰动可抑制爆轰波的传播和爆轰波自身的不稳定性. 爆轰波阵面胞格结构的形成取决于温度扰动与其自身的不稳定性.      

Numerical simulation of detonation reignition in H $_2$-O $_2$ mixtures in area expansions

Time-dependent, two-dimensional, numerical simulations of a transmitted detonation show reignition occuring by one of two mechanisms. The first mechanism involves the collision of triple points as they expand along a decaying shock front. In the second mechanism ignition results from the coalescence of a number of small, relatively high pressure regions left over from the decay of weakened transverse waves. The simulations were performed using an improved chemical kinetic model for stoichiometric H-O mixtures. The initial conditions were a propagating, two-dimensional detonation resolved enough to show transverse wave structure. The calculations provide clarification of the reignition mechanism seen in previous H-O-Ar simulations, and again demonstrate that the transverse wave structure of the detonation front is critical to the reignition process. Received 16 July 1998/ Accepted 17 November 1998     

Analysis of the shock structures in a regular detonation

Time-dependent two-dimensional numerical simulations have been used to investigate the detailed shock structures and patterns of energy release in the regions of the triple points and transverse waves in a planar detonation. As the system of shock triple points evolves between collisions, they trace a well shaped cellular pattern characteristic of detonations in argon-diluted, low-pressure mixtures of hydrogen and oxygen. In the region of the triple points, the shock structure evolves continuously from a single Mach structure to a double Mach structure and finally to a complex Mach structure characteristic of spinning detonations. Most of the energy released in the region of the triple points. The amount of energy release increases as the triple point comes closer to a collision with a wall or another triple point. Just before the collision, there is a large region of energy release that covers the length of the interacting transverse waves. The result is a rectangular high-energy region which boosts the propagation of the new detonation cell.     

超音速来流中爆轰波衍射和二次起爆过程研究

预爆管技术被广泛地应用在爆轰波发动机的起爆过程中,但是在超音速来流中基于预爆管技术起始爆轰波的研究并未被广泛地开展。基于此,本文中数值研究了横向超音速来流对半自由空间内爆轰波的衍射和自发二次起爆、及管道内的衍射和壁面反射二次起爆两种现象的影响。数值模拟的控制方程为二维欧拉方程,空间上使用五阶WENO格式进行数值离散,采用带有诱导步的两步链分支化学反应模型。所模拟的爆轰波具有规则的胞格结构,对应于用惰性气体高度稀释过的可爆混合物中形成的爆轰波。结果表明:在半自由空间内,在本文所模拟的几何尺寸下,爆轰波并未成功发生二次起爆现象,但是爆轰波的自持传播距离随着横向超音速来流强度的增强而增加。在核心的三角形流动区域外,波面诱导产生了更多的横波结构;在管道内,横向的超音速来流在逆流侧对出口气流产生了压缩作用,能有效提高波面压力,因此反射后的激波压力也比较高。在同样的几何尺寸下,爆轰波在静止和超音速(Ma=2.0)气流中分别出现了二次起爆失败和成功两种现象,这是由于在超音速来流中化学反应面的褶皱诱导产生了横波结构,横波与管壁以及其他横波之间的碰撞提高了前导激波的强度,并最终促进了爆轰波在超声速流主管道内的成功起始。     

爆轰波强间断问题的伪弧长算法及其人为解验证

针对该问题开展了伪弧长数值算法研究,通过引入弧长参数,使网格按照一定的形式自适应移动,达到在强间断区域自动加密的效果,从而提高网格分辨率。基于伪弧长算法编写了二维程序,并对程序进行人为解方法验证。将伪弧长算法和直接有限体积法的数值结果进行对比,通过误差分析,显示出伪弧长算法能有效提高计算精度。最后将伪弧长算法应用于气相爆轰波在二维管道中的传播问题,研究了波阵面的捕捉效果和爆轰波胞格结构的形成过程。     

矩形凹槽管道中激波传播的数值研究

基于二维Euler方程,结合五阶加权基本无振荡(weighted essentially nonoscillatory,WENO)格式以及自适应网格加密(adaptive mesh refinement,AMR)技术对入射激波在矩形凹槽管道内传播过程进行了数值模拟。数值结果清晰地显示了入射激波传播过程中与多个矩形凹槽作用以及在凹槽内变化的整个过程,且与已有的实验结果吻合较好。另外,结果还揭示了入射激波与单个凹槽作用时,会发生绕射产生膨胀波,还会发生碰撞从而诱导反射激波。膨胀波会导致入射激波压力降低,而反射激波则导致其升高,但膨胀波的影响占主导作用,因而入射激波波阵面强度出现振荡下降。