Visualization of turbulent combustion of TNT detonation products in a steel vessel

Mixing and afterburning of TNT detonation products in a steel vessel are recorded by the use of the Schlieren visualization system and high speed photography. The vessel is filled with air or 50% oxygen enriched air. Overpressure histories at the vessel wall are also recorded by using pressure transducers. In these experiments nitrogen, air or 50% oxygen enriched air are used as vessel fillers. The Oppenheim-Kuhl theory of thermodynamics of closed systems is applied to estimate the released energy on the basis of pressure histories. Received 29 August 1999 / Accepted 21 January 2000     

铝粉/空气二维黏性两相爆轰的数值模拟

为了深入研究爆轰波形成和传播的机理与特性,建立了管内铝粉/空气二维黏性两相爆轰过程的数学模型,采用守恒元与求解元方法进行数值计算,并对其物理参数的分布进行了分析。结果表明:管内燃烧转爆轰的初期,压力沿径向变化明显,与壁面碰撞有明显的反射波;燃烧转爆轰的中后期,压力沿轴向变化明显,但径向效应仍不能忽视,碰撞形成的反射波对最终稳定爆轰波的形成影响较大。研究结果同时表明:不仅铝粉颗粒初始半径对爆轰波的形成与传播有一定的影响;而且气体的黏性作用在研究爆轰管内近壁面处流场时不容忽视。研究结果有利于进一步揭示铝粉燃烧转爆轰的机理。     

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.      

An exact Riemann solver for detonation products

Numerical methods based upon the Riemann Problem are considered for solving the general initial-value problem for the Euler equations applied to real gases. Most of such methods use an approximate solution of the Riemann problem when real gases are involved. These approximate Riemann solvers do not yield always a good resolution of the flow field, especially for contact surfaces and expansion waves. Moreover, approximate Riemann solvers cannot produce exact solutions for the boundary points. In order to overcome these shortcomings, an exact solution of the Riemann problem is developed, valid for real gases. The method is applied to detonation products obeying a 5th order virial equation of state, in the shock-tube test case. Comparisons between our solver, as implemented in Random Choice Method, and finite difference methods, which do not employ a Riemann solver, are given.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.     

基于CE/SE方法模拟空气中RDX-Al悬浮粉尘的两相爆轰

通过CE/SE方法模拟了空气中炸药-铝粉尘的两相爆轰过程, 研究了双粉尘爆轰过程中粒子不同密度对爆轰波速度、压力的影响, 得到密度与波速、爆压间的线性关系。模拟得到悬浮粉尘在复杂通道中的爆轰波传播过程, 研究了双粉尘爆轰的流场演化过程, 选取流场中的一些点对该处流场的压力及温度随时间的变化进行重点研究, 对比了单铝粉尘在同种条件下的爆轰过程, 发现双粉尘爆轰明显提高了爆轰波波速和流场的压力及温度。模拟结果表明CE/SE方法可以成功模拟双粉尘的爆轰过程, 可为多粉尘爆轰的研究提供参考。     

爆轰产物驱动飞片运动数值模拟研究

利用数值模拟方法,研究了主装药与飞片之间的空气隙厚度变化对飞片中载荷幅度和载荷上升时间、飞片速度和飞片变形量的影响。数值模拟与验证实验结果吻合得较好。结果表明,带空气隙的化爆加载装置可以在飞片中获得较小的载荷幅度和较长的载荷上升时间,并且可以降低飞片击靶速度,从而实现较低的冲击压力,但飞片变形量随着空气隙厚度增大而增大。     

PBX9502炸药爆轰产物的状态方程

采用vanderWaals等效单组分流体模型和Ross硬球微扰理论软球修正模型,计算了爆轰气相产 物的状态方程;用石墨相、金刚石相、类石墨液相和类金刚石液相4种相态描述凝聚成分,由Gibbs自由能最 小确定了不同状态下的凝聚产物相态。对爆轰产物混合系统采用自由能最小原理,通过化学平衡方程组求解 了炸药爆轰产物系统的平衡组分。使用该理论计算了高含碳炸药PBX9502Chapman-Jouguet(CJ)点的爆轰 参数,计算值与实验值符合很好;同时计算了3条等温线,并与Sesame库比较,发现温度超过1000 ℃时,计 算值与Sesame库的计算结果比较接近。在计算的5802K 等温线上发现了一个拐点,分析发现是由于在此 处游离态的碳发生了相变。     

Numerical modelling of detonation structure in two-phase flows

A one-dimensional physical model and a numerical method for the simulation of heterogeneous detonation were proposed based on an Eulerian approach for heterogeneous flows. The combination of modern shock-capturing schemes in combination with a dynamically moving, adaptive grid ensure the properresolution of both reaction zones and flow discontinuities. Numerical examples illustrate the effect of the heat release due to heterogeneous combustion. Received August 4, 1995 / Accepted December 12, 1995     

Study of detonation initiation in hydrogen/air flow

While extensive studies have been conducted concerning the formation of detonation waves in various combustible gaseous mixtures under static conditions since the 1950s, there is very little experimental work on simple flowing systems. In this study, experiments on the deflagration to detonation transition (DDT) of a hydrogen–air flow system were carried out to see the effects of tube diameter, equivalence ratio, and flow types in a premixed and non-premixed flow. Tube diameters used were 25, 50, and 100 mm. The premixed experiments show that the larger tube diameter provides a wider range in run-up distance, reduction of L _DDT/D (ratio of the run-up distance, L _DDT to tube diameter), and expansion of the detonable concentration limit by spreading the cell width. The result of the non-premixed experiments show that similar values of the run-up distance to the premixed experiments are obtained at an equivalence ratio of about 1.0, however, fluctuations of DDT occur near the DDT concentration limit. Under laminar flow conditions at a Reynolds number of less than 2,300, the difference between the two systems could not be observed. However, when the Reynolds number increases towards turbulent conditions, the DDT run-up distance decreases compared to that of static flow conditions. 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.     

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.     

正向爆轰驱动高焓激波风洞的数值模拟

对充满氢氧可燃气体、带扩容腔的正向爆轰驱动的激波风洞进行了数值模拟。计算采用了欧拉方程,频散可控耗散差分格式（DCD）和改进的二阶段化学反应模型。在扩容腔附近采用二维轴对称计算模型,而在驱动段和被驱动段的直管道部分则采用一维计算模型。本文分析了爆轰波在管道中的传播、反射和绕射过程。计算结果表明扩容腔的尺寸对爆轰波的传播、反射、汇聚等起着决定性的作用;带扩容腔的正向爆轰驱动的激波风洞能够得到平稳的持续时间较长的气流,提高了实验的精确度和可重复性。     

爆轰产物组成和爆轰参数计算方法的理论研究

为实现爆轰产物组成和爆轰参数的计算,采用拉格朗日乘数法和牛顿迭代的方法预测爆轰产物组成,利用BKW状态方程预测爆轰参数,在0~600 GPa和300~15 000 K压力温度范围内选取金刚石作为碳的生成相;对爆轰产物系统采用最小自由能原理,结合牛顿迭代法求解爆轰产物的化学平衡方程组;对BKW状态方程参数提出修订,取α=0.5,β=0.298,θ=6 620,κ=9.50;采用自编程序实现计算过程。使用此方法和Hugoniot关系计算密度为1.77 g/cm3的PETN爆轰CJ点爆轰参数验证计算精度,结果显示计算与实验结果的误差均小于1%。利用此方法结合Hugoniot关系预测出爆轰CJ点的产物密度为2.43 g/cm3。     

多循环脉冲爆震发动机流场数值研究

针对二维带有收敛扩张喷管的脉冲爆震发动机模型，采用带基元化学反应的Euler方程组和H2、空气的9组分20基元反应，对发动机在前六个工作循环的流场进行了数值模拟。通过对前几个循环流场进行比较，发现脉冲爆震发动机在第五个循环后流场就基本稳定，单循环得到的流场和多循环稳定后的流场有很大的差别，同时喷管对发动机内流场影响特别大。     

Numerical predictions of oblique detonation stability boundaries

Oblique detonation stability was studied by numerically integrating the two-dimensional, one-step reactive Euler equations in a generalized, curvilinear coordinate system. The integration was accomplished using the Roe scheme combined with fractional stepping; nonlinear flux limiting was used to prevent unphysical solution oscillations near discontinuities. The method was verified on one- and two-dimensional flows with exact solutions, and its ability to correctly predict one-dimensional detonation stability boundaries was demonstrated. The behavior of straight oblique detonations attached to curved walls was then considered. Using the exact, steady oblique detonation solution as an initial condition, the numerical simulation predicted both steady and unsteady oblique detonation solutions when a detonation parameter known as the normal overdrive, , was varied. For a standard test case with a specific heat ratio of , a dimensionless activation energy of , and dimensionless heat release of , an oblique detonation with a constant dimensionless component of velocity tangent to the lead shock, , underwent transition to unstable behavior at . This is slightly higher than the threshold of predicted by one-dimensional theory; thus, the two-dimensionality renders the flow slightly more susceptible to instability. Received 4 August 1996 / Accepted 5 March 1996     

Numerical solution of the two-dimensional nonstationary problem of the motion of a shell under the action of detonation products

The problem of the motion of an incompressible cylindrical shell with an explosive charge is solved numerically for the propagation of a plane detonation wave from the end of the charge. The strength of the shell is not taken into account. A three-term equation of state [1] is assumed for the detonation products. A comparison is made with the one-dimensional case.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 76–79, July–August, 1972.The authors thank G. S. Roslyakov and V. M. Paskonov for assistance in the work and for helpful advice.     

爆轰波在楔面上反射数值分析

应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波在楔面反射进行了数值模拟,计算中氢氧混合物的化学反应采用了8种组分20个反应方程式,在处理化学反应引起的刚性问题时采用了时间算子分裂的方法,模拟了爆轰波在楔面反射由马赫反射向规则反射转变的过程,得到了反射转变临界角,同时考虑了初始压力和组分的影响,并和实验及理论分析结果进行了比较,结果是令人满意的。     

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.