Shock-induced ignition of hydrogen gas during accidental or technical opening of high-pressure tanks

This paper is devoted to the numerical and experimental investigation of hydrogen self-ignition as a result of the formation of a primary shock wave in front of a cold expanding hydrogen gas jet. Temperature increase, as a result of this shock wave, leads to the ignition of the hydrogen–air mixture formed on the contact surface. The required condition for hydrogen self-ignition is to maintain the high temperature in the area for a time long enough for hydrogen and air to mix and inflammation to take place.

Calculations of the self-ignition of a hydrogen jet are based on a physicochemical model involving the gas-dynamic transport of a viscous gas, the kinetics of hydrogen oxidation, the multi-component diffusion, and the heat exchange. We found that the reservoir pressure range, when a shock wave formed in the air during depressurization, has sufficient intensity to produce self-ignition of the hydrogen–air mixture formed at the front of a jet of compressed hydrogen. We present an analysis of the initial conditions (the hydrogen pressure inside the vessel, the temperature of the compressed hydrogen and the surrounding air, and the diameter of the hole through which the jet was emitted), which leads to combustion.     

Numerical study on the spontaneous ignition of pressurized hydrogen release through a tube into air

Spontaneous ignition of pressurized hydrogen release through a tube into air is investigated using a modified version of the KIVA-3V CFD code. A mixture-averaged multi-component approach is used for accurate calculation of molecular transport. Autoignition and combustion chemistry is accounted for using a 21 step kinetic scheme. Ultra fine meshes are employed along with the Arbitrary Lagrangia–Eulerian (ALE) method to reduce false numerical diffusion. The study has demonstrated a possible mechanism for spontaneous ignition through molecular diffusion.

In the simulated scenario, the tube provided additional time to achieve a combustible mixture at the hydrogen–air contact surface. When the tube was sufficiently long under certain release pressure, autoignition would initiate inside the tube at the contact surface due to mass and energy exchange between low temperature hydrogen and shock-heated air through molecular diffusion. Following further development of the hydrogen jet downstream, the contact surface became distorted. Turbulence plays an important role for hydrogen/air mixing in the immediate vicinity of this distorted contact surface and led the initial laminar flame to transit into a stable turbulent flame.     

Numerical study on the spontaneous ignition of pressurized hydrogen during its sudden release into the tube with varying lengths and diameters

Fuel cell vehicles (FCV) and other hydrogen systems with pressurized hydrogen has a safety hazard of spontaneous ignition during its sudden release into the tube. Tube parameter is a key factor affecting the spontaneous ignition of pressurized hydrogen. In this paper, a numerical study on the spontaneous ignition of pressurized hydrogen during its sudden release into the tube with varying lengths and diameters is conducted. The models of Large Eddy Simulation (LES), Eddy Dissipation Concept (EDC), Renormalization Group (RNG), 10-step like opening process of burst disk and 18-step detailed hydrogen combustion mechanism are employed. 6 cases are simulated based on the previous experiments. Numerical results show that the possibility of spontaneous ignition of pressurized hydrogen increases inside the longer and thinner tubes, which agrees with the experimental results. The increasing of tube length has little influence on the shock wave formation and propagation inside the tube. However, there exists critical tube lengths for the generation of Mach disk and the normal shock wave: the maximum and minimum distances for the generation of the Mach disk in 10 mm diameter tube are 7.8 and 6.7 mm, respectively. As for the normal shock wave, these critical values are 22.1 and 19.4 mm, respectively. In addition, the formation times and initial positions of Mach disk and normal shock wave are delayed inside the thicker tube. Due to the shock-affected time increases with the increasing of tube length, the temperature could rise to the critical ignition temperature and triggers the spontaneous ignition due to the sufficient tube length even though the less hydrogen/air mixture and the contact surface with lower temperature is produced inside the thicker tube. Finally, a simple time scale analysis is conducted.     

Numerical analysis of hydrogen deflagration mitigation by venting through a duct

This paper presents a model and simulation results for the mitigation of a hydrogen–air deflagration by venting through a duct. A large eddy simulation (LES) model, applied previously to study both closed-vessel, and open atmosphere hydrogen–air deflagrations, was developed further to model a hydrogen–air explosion vented through a duct. Sub-grid scale (SGS) flame wrinkling factors were introduced to model major phenomena which contribute to the increase of flame surface area in vented deflagrations. Simulations were conducted to validate the model against 20% hydrogen–air mixture deflagrations (vent diameters 25 and 45 cm) and 10% hydrogen–air mixture deflagration (vent diameter 25 cm). There was reasonable correlation between the simulations and the experimental data. The comparative importance of different physical phenomena contributing to the flame wrinkling is discussed.     

Experimental investigation on shock waves generated by pressurized gas release through a tube

An experimental investigation on the flow structures and the strength of shock waves generated by high-pressure gas release through a tube into air was conducted. The results demonstrated that a leading shock wave was generated in front of the compressed gas jet and the shock wave speed increased firstly, then decreased and finally kept constant with an increase of the propagation distance in the tube. The experimentally measured Mach numbers of shock waves were close to those calculated from the theory of ideal shock tube flow. After spouting out of the tube, the normal shock quickly developed into a hemispherical shape. The Mach disk was observed in the under-expanded jet. For high-pressure combustible gas release, the concept of theoretical critical pressure of ignition was introduced and several theoretical critical pressures of common gaseous fuels were obtained.     

Self-ignition and explosion during discharge of high-pressure hydrogen

The phenomenon of self-ignition and explosion during discharge of high-pressure hydrogen was investigated. To clarify the ignition conditions of high-pressure hydrogen jets, rapid discharge of the high-pressure hydrogen was examined experimentally. A diaphragm was used to allow rapid discharge of the high-pressure hydrogen. The burst pressure was varied from 4 to 30 MPa. The downstream geometry of the diaphragm was a flange and extension pipes, with the pipe length varying from 3 to 300 mm. The diameter of the nozzle was 5 or 10 mm. When short pipes were used, the hydrogen jet did not ignite. However, the hydrogen jet showed an increasing tendency to ignite in the pipe as the length of the pipe became longer. At higher burst pressures, a diffusion jet flame was formed from the pipe. The blast wave from the fireball formed on self-ignition of the hydrogen jet resulted in an extremely rapid pressure rise.     

A methodology to predict shock overpressure decay in a tunnel produced by a premixed methane/air explosion

A methodology for estimating the blast wave overpressure decay in air produced by a gas explosion in a closed-ended tunnel is proposed based on numerical simulations. The influence of the tunnel wall roughness is taken into account in studying a methane/air mixture explosion and the subsequent propagation of the resulting shock wave in air. The pressure time-history is obtained at different axial locations in the tunnel outside the methane/air mixture. If the shock overpressure at two, or more locations, is known, the value at other locations can be determined according to a simple power law. The study demonstrates the accuracy of the proposed methodology to estimate the overpressure change with distance for shock waves in air produced by methane/air mixture explosions. The methodology is applied to experimental data in order to validate the approach.     

Effects of premixed methane concentration on distribution of flame region and hazard effects in a tube and a tunnel gas explosion

Study of flame distribution laws and the hazard effects in a tunnel gas explosion accident is of great importance for safety issue. However, it has not yet been fully explored. The object of present work is mainly to study the effects of premixed gas concentration on the distribution law of the flame region and the hazard effects involving methane-air explosion in a tube and a tunnel based on experimental and numerical results. The experiments were conducted in a tube with one end closed and the other open. The tube was partially filled with premixed methane-air mixture with six different premixed methane concentrations. Major simulation works were performed in a full-scale tunnel with a length of 1000 m. The first 56 m of the tunnel were occupied by methane–air mixture. Results show that the flame region is always longer than the original gas region in any case. Concentration has significant effects on the flame region distribution and the explosion behaviors. In the tube, peak overpressures and maximum rates of overpressure rise (dp/dt)_max for mixtures with lower and higher concentrations are great lower than that for mixtures close to stoichiometric concentration. Due to the gas diffusion effect, not the stoichiometric mixture but the mixture with a slightly higher concentration of 11% gets the highest peak overpressure and the shock wave speed along the tube. In the full-scale tunnel, for fuel lean and stoichiometric mixture, the maximum peak combustion rates is achieved before arriving at the boundary of the original methane accumulation region, while for fuel rich mixture, the maximum value appears beyond the region. It is also found that the flame region for the case of stoichiometric mixture is the shortest as 72 m since the higher explosion intensity shortens the gas diffusion time. The case for concentration of 13% can reach up to a longest value of 128 m for longer diffusion time and the abundant fuel. The “serious injury and death” zone caused by shock wave may reach up to 3–8 times of the length of the original methane occupied region, which is the widest damage region.     

Comments on explosion problems for hydrogen safety

The future widespread use of hydrogen as an energy carrier brings in safety issues that have to be addressed before public acceptance can be achieved. The prediction of the consequences of a major accident release of hydrogen into the atmosphere or the contamination of high-pressure hydrogen storage facilities by air entrainment requires a good knowledge of the explosion parameters of hydrogen–air mixtures. The present paper reviews and comments on the current knowledge of dynamic parameters of hydrogen detonation for hazard assessment. The major problem that remains to be resolved involves the understanding of the effect of turbulence on the cellular detonation structure, the propagation of high-speed deflagrations and the transition from deflagration to detonations. It is recommended that future research should be aimed towards experiments that permit the quantitative understanding of the mechanisms of high-speed turbulent combustion rather towards large-scale tests in complex geometries where minimal quantitative information of fundamental significance could be extracted. In spite of its wide flammability and sensitivity to ignition and detonation initiation, it is felt that hydrogen can be produced, stored and handled safely with the appropriate considerations in the design of the hydrogen facilities.     

Scale effects on hydrogen–air fast deflagrations and detonations in small obstructed channels

An experimental study of flame propagation, acceleration and transition to detonation in hydrogen–air mixture in 2-m-long rectangular cross-section channel filled with obstacles located at the bottom wall was performed. The initial conditions of the hydrogen–air mixture were 0.1 MPa and 293 K and stoichiometric composition (29.6% H₂ in air). The channel width was 0.11 m and blockage ratio was 0.5 in all experiments. The effect of channel geometrical scale on flame propagation was studied by using four channel heights H of 0.01, 0.02, 0.04, and 0.08 m. In each case, the obstacle height was equal to H/2 and the obstacle spacing was 2H.

The propagation of flame and pressure waves was monitored by four pressure transducers and four ion probes. The pairs of transducers and probes were placed at various locations along the channel in order to get information about the progress of the phenomena along the channel.

As a result of the experiments, the deflagration and detonation regimes and velocities of flame propagation in the obstructed channel were established.     

Experimental investigation on spontaneous ignition caused by pressurized hydrogen suddenly release into an S-shaped tube

An experimental investigation on the effects of continuous semicircular curved structure on spontaneous ignition during pressurized hydrogen suddenly release was conducted. An S-shaped tube with 700 mm in length and 10 mm in diameter was used in our experiments, and a straight tube with the same configuration was adopted for comparison. The results show that the continuously generated rarefaction waves and reflected shock waves make the pressure curves in the S-shaped tube more complicated. Meanwhile, the mean velocity and intensity of the leading shock wave undergo considerable attenuation when it propagates in the S-shaped structure. By comparing with the straight tube, the minimum critical pressure condition for spontaneous ignition in the S-shaped tube is slightly difficult to reach, but the difference is not huge. Nevertheless, the S-shaped structure can effectively promote hydrogen-air mixing and make combustion more intense. A secondary overpressure peak detected by the pressure transducer near the nozzle occurs in the spontaneous ignition cases and no such pressure increase is caught in the non-ignition cases. The transition from spontaneous combustion flame to a jet flame at the nozzle and the complete out-tube jet flame development process are captured and discussed.     

Explosions and deflagration-to-detonation transitions in epoxy propane/air mixtures

The explosion and deflagration-to-detonation transition (DDT) in epoxy propane (E.P.) vapor/air mixture clouds under weak ignition conditions has been studied in an experimental tube of diameter 199 mm and length 29.6 m. E.P. vapor clouds were formed by injecting liquid E.P. into the experimental tube and evaporating of the fine E.P. droplets. The dimension and the evaporating process of the E.P. droplet were measured and analyzed. The E.P. vapor/air mixture clouds were ignited by an electric spark with an ignition energy of 40 J. The characteristics and the stages of the DDT process in the E.P. vapor/air mixtures have been studied and analyzed. A self-sustained detonation wave formed, as was evident from the existence of a transverse wave and a cellular structure. Moreover, a retonation wave formed during the DDT process in the E.P. vapor/air mixture. The influence of the E.P. vapor concentration on the DDT process has been studied. The minimum E.P. vapor concentration for the occurrence of the DDT in the E.P. vapor/air mixture has been evaluated and the variation of DDT distance with E.P. vapor concentration has been analyzed.     

甲烷爆炸传播过程膜状障碍物的激励效应研究

针对管状空间内膜状障碍物对甲烷爆炸传播的激励效应现象,基于机理分析进行了数值模拟和实验研究,计算分析薄膜附近爆炸冲击波压力峰值大小与火焰速度变化,同时运用激波管道进行相同工况条件下的实验,并对两者结果对比分析,发现有无膜状障碍物的压力峰值相差6倍以上。研究表明,膜状障碍物的激励效应是破膜以后形成的带压燃烧,提高了燃烧速率,导致甲烷爆炸的火焰传播速度剧增。实验结果一定意义诠释了同样数量的甲烷气体爆炸在不同环境内后果上的巨大差异,研究结果对矿井瓦斯爆炸事故调查及防治具有指导意义。     

Analysis of the experimental results of the initiation of detonation in short tubes with kerosene–oxidizer mixtures

The paper describes the experimental investigation of detonation initiation in a mixture of kerosene–oxidant in a short test tube. Various mixtures of oxygen and nitrogen were used as an oxidant, from pure oxygen to the composition of air. The goal of the study was to determine the minimum diameter of the tube and the minimum level of energy needed for the direct initiation of detonation. As a result of the measurements the pressure courses were obtained for two kinds of cases: with and without (only shock waves) of fuel injection. The results of both kinds of measurements were compared, providing information about the initiation of detonation in a fuel–oxidizer mixture. Brief analyses of the results for different initiators and different oxidizers were performed and compared with the shock wave and Chapman–Jouget velocity.     

200 MW燃气流风洞高压氧气系统安全操作分析

为保证200 MW燃气流风洞高压氧气系统安全运行,从初始能量出发,对高压氧气系统充气、供气、排气时管道内的激波管流动、绝热压缩等过程进行安全分析,并提出针对性安全措施。结果表明:对于充气管道内存在的激波管流动,当驱动气体压力为20 MPa、被驱动气体压力为0.1 MPa时,激波反射后末端气体温度远远高于200 ℃,通过减小阀门开启速度,对阀前管道进行充气以减小上下游压差,可避免因绝热压缩产生的高温;供气管道充填时,管道内最高温度为73 ℃,通过控制充填速度,可进一步降低管道内氧气温度;通过高压排气、低压排气2种模式,可满足国标中对氧气流速的要求。研究结果可为氧气管道远程安全操作提供参考。     

尿素合成塔爆炸事故机理的研究

尿素合成塔工艺条件复杂,体积庞大,内部物料状态复杂,需要在高温、高压下操作,国内外曾发生过几起著名的尿素合成塔爆炸事故,事故后果非常严重。工艺条件的复杂性加大了事故分析的难度,目前对事故原因众说纷纭。利用数值模拟软件,对尿素合成塔爆炸过程进行研究,对爆炸机理进行分析,对事故发生的原因进行探讨,为预防同类事故提供一定的依据。分析了顶部气相空间爆炸性混合体系的形成和爆炸过程,结合数值模拟得出关键点的压力变化和气相冲击波最终的压力,依据能量守恒推导出液相中冲击波的压力。利用数值模拟得出冲击波在液相中的传播速度和叠加的情况,综合反射、直接透射两部分冲击波进入液相的时间差,找到位于筒体底部的叠加点,冲击波在该处叠加造成筒体破裂,过热液体瞬间气化喷出,引发破坏更大的二次蒸汽大爆炸。     

Numerical simulations on the attenuation effect of a barrier material on a blast wave

This paper numerically modeled previous experimental results and quantitatively revealed the attenuation effect of a barrier material on a blast wave. Four fluids were considered in the present study: the detonation products, water, foamed polystyrene, and air. These fluids were modeled by Jones-Wilkins-Lee (JWL), stiffened gas, and ideal gas equations of state. A mixture of water and foamed polystyrene was used as a barrier to encircle a 0.1 kg mass of spherical pentolite, and the interface problem between the barrier and the blast wave was investigated. The simulation parameters were the radius and the water volume fraction of the barrier. To elucidate the effect of the barrier, we conducted two series of numerical simulations; one without a barrier, and another with a barrier of 50 or 100 mm in outer radius and 0–1 in the water volume fraction. Peak overpressure, positive impulse, and pressure history all agreed well with the previous experimental results. We focused on the energy transfer from high-pressure detonation products to other fluids. The sum of the kinetic energies of the detonation products and the barrier induced by the blast wave could quantitatively estimate the attenuation effect of the blast wave and was minimized when the water volume fraction was 0.5, as was the case in the previous experiment.     

Methane/coal dust/air explosions and their suppression by solid particle suppressing agents in a large-scale experimental tube

Methane/coal dust/air explosions under strong ignition conditions have been studied in a 199 mm inner diameter and 30.8 m long horizontal tube. A fuel gas/air manifold assembly was used to introduce methane and air into the experimental tube, and an array of 44 equally spaced dust dispersion units was used to disperse coal dust particles into the tube. The methane/coal dust/air mixture was ignited by a 7 m long epoxypropane mist cloud explosion. A deflagration-to-detonation transition (DDT) was observed, and a self-sustained detonation wave characterized by the existence of a transverse wave was propagated in the methane/coal dust/air mixtures.The suppressing effects on methane/coal dust/air mixture explosions of three solid particle suppressing agents have been studied. Coal dust and the suppressing agent were injected into the experimental tube by the dust dispersion units. The length of the suppression was 14 m. The suppression agents examined in this study comprised ABC powder, SiO₂ powder, and rock dust powder (CaCO₃). Methane/coal dust/air explosions can be efficiently suppressed by the suppression agents characterized by the rapid decrease in overpressure and propagating velocity of the explosion waves.     

Real-time numerical simulations and experimental research for the propagation characteristics of shock waves and gas flow during coal and gas outburst

When coal and gas outburst occurs, high-speed gas flow and air shock wave with high kinetic energy could be created. In this paper, the formation process of outburst shock waves and gas flow has been analyzed firstly. Afterwards, the numerical simulation models of the roadways with right-angled intersection have been established, by which real-time simulation of the propagation of outburst gas flow and the process of gas transport has been conducted. Gas pressure, gas velocity and gas concentration can be simulated and shown. From analyzing the simulation results, qualitative and quantitative conclusions that the characteristics and patterns of the propagation and attenuation of outburst shock waves and gas flow can be arrived at. Finally, experimental models have been carried out to investigate the outburst shock waves and gas flow at the roadways with the similar shapes as the simulated ones. The results indicate that when shock wave and gas flow passes the intersection, most of the shock wave and gas flow will flow into the roadway of section opposite the intersection, and a little of it would flow into the roadway below the intersection. And turbulence will appear, shock wave reflects and diffracts at branches with more influence on the roadway below the intersection.     

Influence of laneway support spacing on methane/air explosion shock wave

A laneway support system provides an available way to solve problems related to ground movements in underground coal mines, but also poses another potential hazard. Once a methane/air explosion occurs in a laneway, inappropriate design parameters of the support system, especially the support spacing, likely have a negative influence on explosion disaster effects. The commercial software package AutoReaGas, a computational fluid dynamics code suitable for gas explosions, was used to carry out the numerical investigation for the methane/air explosion and blast process in a straight laneway with different support spacing. The validity of the numerical method was verified by the methane/air explosion experiment in a steel tube. Laneway supports can promote the development of turbulence and explosion, and also inhibit the propagation of flame and shock wave. For the design parameters in actual laneway projects, the fluid dynamic drag due to the laneway support plays a predominant role in a methane/air explosion. There is an uneven distribution of the peak overpressure on the same cross section in the laneway, and the largest overpressure is near the laneway walls. Different support spacing can cause obvious differences for the distributions of the shock wave overpressure and impulse. Under comparable conditions, the greater destructive effects of explosion shock wave are seen for the laneway support system with larger spacing. The results presented in this work provide a theoretical basis for the optimized design of the support system in coal laneways and the related safety assessments.