A study of the blast wave shape from elongated VCEs

Elongated congestion patterns are common at chemical processing and petroleum refining facilities due to the arrangement of processing units. The accidental vapor cloud explosion (VCE) which occurred at the Buncefield, UK facility involved an elongated congested volume formed by the trees and undergrowth along the site boundary. Although elongated congested volumes are common, there have been few evaluations reported for the blast loads produced by elongated VCEs. Standard VCE blast load prediction techniques do not directly consider the impact of this congested volume geometry versus a more compact geometry.This paper discusses an evaluation performed to characterize the blast loads from elongated VCEs and to identify some significant differences in the resulting blast wave shape versus those predicted by well-known VCE blast load methodologies (e.g., BST and TNO MEM). The standard blast curves are based on an assumption that the portion of the flammable gas cloud participating in the VCE is hemispherical and located at grade level. The results of this evaluation showed that the blast wave shape for an elongated VCE in the near-field along the long-axis direction is similar to that for an acoustic wave generated in hemispherical VCEs with a low flame speed. Like an acoustic wave, an elongated VCE blast wave has a very quick transition from the positive phase peak pressure to the negative phase peak pressure, relative to the positive phase duration. The magnitude of the applied negative pressure on a building face depends strongly on the transition time between the positive and negative phase peak pressures, and this applied negative phase can be important to structural response under certain conditions. The main purpose of this evaluation was to extend previous work in order to investigate how an elongated VCE geometry impacts the resultant blast wave shape in the near-field. The influence of the normalized flame travel distance and the flame speed on the blast wave shape was examined. Deflagration and deflagration-to-detonation transition regimes were also identified for unconfined elongated VCEs as a function of the normalized flame travel distance and flame speed attained at a specified flame travel distance.     

Computer simulation of shock waves transmission in obstructed terrains

Generation and transmission of blast waves in real terrains is of major importance for risk analysis procedures involving accidental explosion scenarios. The problem arises from the impact of overpressure wave on people and structures that may be lethal or catastrophic under certain conditions. In this paper, a CFD simulation of shock wave propagation in obstructed terrain is attempted. Overpressure histories as well as a series of critical parameters, namely the positive and negative peak overpressure, the arrival time, and the positive and negative phase duration at specific points within the domain were obtained during the simulation. Their comparison with experimental measurements from field-scale high explosive blast tests performed by HSE showed a reasonably good agreement indicating that CFD computer programs provide reliable tools for estimating explosive shocks in complex terrains.     

Measurement of gas detonation blast loads in semiconfined geometry

The ignition and explosion of combustible vapor clouds represents a significant hazard across a range of industries. In this work, a new set of gas detonations experiments were performed to provide benchmark blast loading data for non-trivial geometry and explosion cases. The experiments were designed to represent two different accident scenarios: one where ignition of the vapor cloud occurs shortly after release and another where ignition is delayed and a fuel concentration gradient is allowed to develop. The experiments focused on hydrogen-air and methane-oxygen detonations in a semiconfined enclosure with TNT equivalencies ranging from 9 g to 1.81 kg. High-rate pressure transducers were used to record the blast loads imparted on the interior walls of a 1.8 m × 1.8 m × 1.8 m test fixture. Measurements included detonation wave velocity, peak overpressure, impulse, and positive phase duration. A comparison of the pressure and impulse measurements with several VCE models is provided. Results show that even for the most simplified experimental configuration, the simplified VCE models fail to provide predictions of the blast loading on the internal walls of the test fixture. It is shown that the confinement geometry of the experiment resulted in multiple blast wave reflections during the initial positive phase duration portion of the blast loading, and thus, significantly larger blast impulse values were measured than those predicted by analytical models. For the pressure sensors that experienced normally-reflect blast waves for the initial blast impulse, the Baker-Strehlow and TNT equivalency models still struggled to accurately capture the peak overpressure and reflected impulse. The TNO multi-energy model, however, performed better for the case of simple normally-reflected blast waves. The results presented here may be used as validation data for future model or simulation development.     

Using computational fluid dynamics (CFD) for blast wave predictions

Explosions will, in most cases, generate blast waves. While simple models (e.g., Multi Energy Method) are useful for simple explosion geometries, most practical explosions are far from trivial and require detailed analyses. For a reliable estimate of the blast from a gas explosion it is necessary to know the explosion strength. The source explosion may not be symmetric; the pressure waves will be reflected or deflected when hitting objects, or even worse, the blast waves may propagate inside buildings or tunnels with a very low rate of decay. The use of computational fluid dynamics (CFD) explosion models for near and far field blast wave predictions has many advantages. These include more precise estimates of the energy and resulting pressure of the blast wave, as well as the ability to evaluate non-symmetrical effects caused by realistic geometries, gas cloud variations and ignition locations. This is essential when evaluating the likelihood of a given leak source as cause of an explosion or equally when evaluating the potential risk associated with a given leak source for a consequence analysis.In addition, unlike simple methods, CFD explosion models can also evaluate detailed dynamic effects in the near and far field, which include time dependent pressure loads as well as reflection and focusing of the blast waves. This is particularly valuable when assessing actual near-field blast damage during an explosion investigation or potential near-field damage during a risk analysis for a facility. One main challenge in applying CFD, however, is that these models require more information about the actual facility, including geometry details and process information. Collecting the necessary geometry and process data may be quite time consuming. This paper will show some blast prediction validation examples for the CFD model FLACS. It will also provide examples of how directional effects or interaction with objects can significantly influence the dynamics of the blast wave. Finally, the challenge of obtaining useful predictions with insufficient details regarding the geometry will also be addressed.     

Pressure characteristics and dynamic response of coal mine refuge chamber with underground gas explosion

Coal mine refuge chambers are new devices for coal mine safety which can provide basic survival conditions after gas explosion. In order to simulate the propagation of underground methane/air mixture blast wave, and check structural safety of coal mine mobile refuge chamber, an underground tunnel model and a refuge chamber model have been established based on explicit nonlinear dynamic ANSYS/LS-DYNA 970 program. Results show that the reflected wave pressure on the impact surface was about two times higher than that on the incident one. The relationship between the pressure fields of the chamber was analyzed. The maximum pressure of gas explosion reached about 0.71 MPa, and the pulse width was 360 ms. The maximum absolute displacement and stress occurs at the main door center and the connection of stiffeners and the front plate, respectively. The entire coal mine mobile refuge chamber was in elastic state and its strength and stiffness meet the safety requirements. The cabin door, the front plate and the connecting flange at cabin back as well as the stiffeners on each side were the most critical components. Suggestions were put forward for the refuge chamber.     

Blast phenomena in urban tunnel systems

Traffic infrastructure in urbanised areas is increasingly projected in tunnels underground or covered over, these days. A consequence is that in case of an incident with hazardous materials the safety level for fellow road users in tunnels is considerably less than it is in surface infrastructure. To reduce the consequences of incidents for fellow tunnel users, urban tunnels are sometimes interrupted by open spaces of limited length. Open spaces allow, for instance, the release of smoke in case of a fire. In this way, possible lethal effects are limited to the tunnel section in which the incident occurred. To what extent an open space may also be effective in the mitigation of blast effects from an explosion in a tunnel system is subject of this paper.

To this end, the blast effects originating from the rupture of a 50 m³ LPG pressure vessel in an urban tunnel system have been computed by numerical simulation. The results show that an open space in a tunnel system has a significant mitigating effect on the blast effects indeed. However, as a consequence of the ingress of a high-velocity jet flow that follows on a primary blast wave, a second blast wave develops in the tunnel section following on an open space. The strength of this second blast wave is not very dependent on the length of the open space. It shows that an open space in a tunnel system may not always limit the lethal effects of explosion incidents in tunnels to the tube in which the incident occurred. The second blast wave in the tunnel section following on an open space may have lethal consequences for fellow tunnel users by car window failure.     

A fuzzy set analysis to estimate loss intensity following blast wave interaction with process equipment

Blast waves are able to produce structural damage to process equipment even at great distances from the source point of an explosion. A loss of containment may follow and, if hazardous substances are released, relevant secondary scenarios may be triggered, resulting in domino effects.The present study was focused on the assessment of the expected structural damage and of the associated intensity of loss of containment of process vessels loaded by blast waves. Hence, a knowledge-based fuzzy set analysis was used to assess the expected overall probability of occurrence of different damage states defined for several categories of process equipment items. The fuzzy approach was also used to obtain specific threshold values for the escalation sequences (domino effects), taking into account the hazard due to the expected secondary scenarios caused by the loss of containment following blast wave impact.     

Numerical study of dynamic response and failure analysis of spherical storage tanks under external blast loading

The performance of energy infrastructures under extreme loading conditions, especially for blast and impact conditions, is of great importance despite the low probability for such events to occur. Due to catastrophic consequences of structural failure, it is crucial to improve the resistance of energy infrastructures against the impact of blasts. A TNT equivalent method is used to simulate a petroleum gas vapor cloud explosion when analyzing the dynamic responses of a spherical tank under external blast loads. The pressure distribution on the surface of a 1000 m³ spherical storage tank is investigated. The dynamic responses of the tank, such as the distribution of effective stress, structural displacement, failure mode and energy distribution under the blast loads are studied and the simulation results reveal that the reflected pressure on the spherical tank decreases gradually from the equator to the poles of the sphere. However, the effects of the shock wave reflection are not so evident on the pillars. The structural damage of the tank subjected to blast loads included partial pillar failure from bending deformation and significant stress concentration, which can be observed in the joint between the pillar and the bottom of the spherical shell. The main reason for the remarkable deformation and structural damage is because of the initial internal energy that the tank obtained from the blast shock wave. The liquid in the tank absorbs the energy of impact loads and reduces the response at the initial stage of damage after the impact of the blast.     

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.     

Impact of a shock wave on a structure on explosion at altitude

The number of explosive attacks on civilian buildings has recently increased and the pattern of damage inflicted on structures when an explosion takes place at altitude remains quite difficult to predict. The primary aim of the work reported here was to enhance the understanding of how blast waves from an explosion at altitude interact with the ground and with a structure. Small-scale experiments were conducted using a propane–oxygen stoichiometric mixture as explosive. This approach is original because it models high-explosive detonation in terms of gaseous charge explosion using TNT equivalents. Several non-dimensional laws are expressed and validated by experiments. These relationships allow determination of the propagation of a blast wave and its interaction with a structure as a function of the position of the explosive charge when the explosion occurs at altitude. Then, from knowledge of the blast loading, using Hopkinson's scaling law and TNT equivalents, we can predict the interaction of blast waves with the ground and a structure on a real scale. Simulations were performed using the Autodyn code, and good correlation with the experimental results was obtained.     

Numerical and reduced-scale experimental investigation of blast wave shape in underground transportation infrastructure

When an explosion occurs in a tunnel, the study of the blast wave quickly becomes complicated, owing to the multiple propagation patterns of the blast wave (incident wave, regular and Mach reflections) and to the geometrical conditions. Considering this problem, two patterns can be revealed. Near the explosive, the well-known free-field pressure wave can be observed. After multiple reflections on the tunnel's walls, this overpressure behaves like a one-dimensional (1D) wave. One aim of this paper is to determine the position of this transition spherical-to-planar wave propagation in a tunnel using both numerical and reduced-scale experiments, and thereby validate the dedicated law established in a previous work.For this purpose, a detonation of TNT in a tunnel with a cross-section of up to 55 m² is considered. Results show good agreement between the numerical simulations and experiments. The transition zone between the three-dimensional (3D) and the 1D wave is well detected. An application to a simplified subway station is also investigated which shows that significant planar waves can be transmitted to the neighboring stations via the junction tunnels.     

基于风险的液化烃罐区建筑物抗爆荷载评估方法

为合理确定液化烃罐区周边建筑物的抗爆设防荷载，有效进行抗爆设计和防护，建立1套系统的抗爆设防荷载定量评估方法。以某液化烃罐区建筑物为例，计算172个爆炸场景，获得4组累积爆炸频率曲线，基于风险控制标准确定抗爆设防荷载。结果表明:爆炸场景发生频率应包括初始泄漏频率、气象概率、泄漏方向概率和延迟爆炸概率；获得的爆炸超压-累积频率曲线是确定抗爆设防荷载的基础，在爆炸超压较低时，与爆炸源中心距离不同的4面墙体的超压累积频率曲线极为接近；随着爆炸超压的继续增大，累积发生频率的差异逐渐明显；液化烃罐区建筑物的抗爆设防荷载应同时满足2个准则，即万年1次的风险可接受准则和风险可接受范围内爆炸超压最大化准则；根据该准则确定的液化烃罐区附近建筑物东墙的爆炸冲击波峰值入射超压为44.6 kPa，正压作用时间为89.3 ms。     

Numerical study of medium to large scale BLEVE for blast wave prediction

So far, the prediction of blast wave generated from the Boiling Liquid Expanding Vapour Explosion (BLEVE) has been already broadly investigated. However, only a few validations of these blast wave prediction models have been made, and some well-established methods are available to predict BLEVE overpressure in the open space only. This paper presents numerical study on the estimation of the near-field and far-field blast waves from BLEVEs. The scale effect is taken into account by conducting two different scale BLEVE simulations. The expansion of pressurized vapour and evaporation of liquid in BLEVE are both modelled by using CFD method. Two approaches are proposed to determine the initial pressure of BLEVE source. The vapour evaporation and liquid flashing are simulated separately in these two approaches. Satisfactory agreement between the CFD simulation results and experimental data is achieved. With the validated CFD model, the results predicted by the proposed approaches can be used to predict explosion loads for better assessment of explosion effects on structures.     

A numerical study of the evolution of the blast wave shape in rectangular tunnels

When the explosion of condensed materials occurs in square or circular cross-section tunnel, the subsequent blast wave reveals two patterns: three-dimensional close to the explosive charge and one-dimensional far from the explosion. Pressure decays for these two patterns have been thoroughly studied. However, when the explosion occurs in rectangular cross-section tunnel, which is the most regular geometry for underground networks, the blast wave exhibits a third, two-dimensional, patterns. In order to assess the range of these three patterns, several numerical simulation of blast waves were carried out varying the width and the height of the rectangular cross-section as well as the mass of the charge. Laws are presented to localize the transition zones between the 3D and the 2D patterns, and between the 2D and the 1D patterns, as functions of non-dimensional width and height. The numerical results of the overpressure are compared to existing 3D and 1D laws. An overpressure decay law is proposed to represent the 2D pattern. Knowing the two transition zones and the overpressure decays within these zones, an algorithm is presented to efficiently predict an overpressure map. This algorithm is validated by comparison with experimental data.     

爆炸冲击波在建筑群中传播规律的数值模拟研究

为了研究建筑群的布局特性对爆炸冲击波传播的干扰作用及爆炸对周围环境的影响,笔者基于任意拉格朗日多物质流固耦合算法,对空气采用ALE网格,爆轰产物采用JWL状态方程,利用LS-DYNA程序对TNT理想爆炸源的爆炸进行了数值模拟,并与经验公式的计算结果进行了对比,得出了数值模拟方法可行的结论。在该基础上,以一个建筑小区为例,对其遭受爆炸后的响应作了进一步模拟研究。结果表明,建筑物的密度和布局方式对爆炸场的数值与形态都有敏感作用,因此,在城市规划建筑布局阶段,应考虑建筑密度和布局对爆炸冲击波的影响。     

On the mechanism of a BLEVE occurrence due to fire engulfment of tanks with overheated liquids

The BLEVE (boiling liquid expanding vapor explosion) effect that involves the formation of a fireball occurs at the engulfment by fire of a tank with a highly flammable liquid or a liquid gas. Heating of the tank causes elevation of the liquid phase temperature and pressure inside the tank. A partial rupture of the dry tank walls is possible, with the formation of a rarefaction wave propagating into the liquid phase. An evaporation wave moves after the rarefaction wave and cause a rapid increase of pressure, exceeding the initial pressure before depressurization. Rapid violent destruction of the tank occurs. The mechanism of a BLEVE initiation is considered using Van der Waals isotherms. The following criterion for the possibility of a BLEVE was formulated. If the final state is located on an unstable part of the Van der Waals isotherm, a BLEVE takes place. Limiting values of the temperatures for overheating of certain highly flammable liquids and liquid gases (propane, n-butane, n-pentane, isopentane) were calculated using the proposed method, and were found to be in good agreement with experimental data available in the literature.     

水不耦合装药对深孔预裂爆破应力波能量的影响

从理论上探讨了在深孔预裂爆破过程中,水不耦合装药对爆轰冲击波的形成和应力波传播规律的影响,利用Matlab编程求解了冲击波初始参数和沿爆破孔径向方向上传播到爆破孔孔壁处的参数;并根据弹性理论提出并求解了正入射情况下爆破孔孔壁表面上的初始冲击压力。在淮南矿业集团丁集矿进行了水不耦合装药深孔预裂爆破实验,并对实验获得的应力波信号进行归一化频谱分析研究。数值计算和现场实验均表明:通过在煤层中进行深孔预裂爆破实验研究发现,当ξ的值为1.5左右时,爆生应力波在10~50Hz内的分能量最强,同时具有最高的信噪比,这与理论计算结果相吻合。因而,在利用煤层爆破强化增透技术时应将不耦合装药系数设计为1.5左右,同时适当加大爆破药柱的用量,以获得最佳的增透效果。     

城市天然气工程环境风险评价

结合广州、深圳、东莞、佛山等城市的天然气工程情况,对城市天然气工程的环境风险评价进行了讨论,并采用穆尔哈斯(Moorhowse)和普里恰特(Prichard)提出的热辐射预测模式和爆炸冲击波预测模式对城市天然气工程进行了风险评价.结果表明,如果发生天然气泄漏并引起火灾,假设在10 min以内,城市门站、调压站和城市高中压管道的火球对建筑物和设备的严重损害范围( A 级)最远距离分别为35.8 m、18.0 m和28.7 m,爆炸冲击波严重损害范围( C1 级)分别为距事故处54.8 m、27.8 m和44.4 m.最后从城市门站、调压站选址及输气管道选线、安全防范距离、作业过程中的风险控制与管理以及事故应急对策四方面提出了风险事故的防范措施与对策.     

瓦斯爆炸传播过程中矿车激励效应的实验研究

为了分析矿车在矿井瓦斯爆炸传播过程中的激励效应,模拟井下工作环境,构建了实验管道,按照比例建造了矿车模型,分别设定一个矿车和三个矿车链接两种实验工况进行试验。根据冲击波传播过程中激励效应的形成机制,以爆炸传播过程中火焰和冲击波波阵面变化测试为主要测试目标,采用激光纹影的测试方法,分别观测了两种工况条件下下矿车区域附近爆炸冲击波畸变情况,捕获了冲击波经过矿车附近时,冲击波变化特征。根据试验观测,当冲击波经过障碍物附近时,波阵面明显发生畸变,而当三个矿车存在时,火焰出现了回流和绕流现象,燃烧强度和速度都明显增加,这一实验结果一方面证实了矿车对瓦斯爆炸存在明显的激励效应,同时指出矿车数量多,激励效应会更强。本文研究结果对瓦斯爆炸事故阻爆、隔爆技术开发和事故调查有一定的指导意义。     

A new standard for predicting lung injury inflicted by Friedlander blast waves

An important blast injury mechanism is the rupture of the lungs and the gastrointestinal tract. In explosives safety studies and threat analysis the empirical model of Bowen is often used to quantify this mechanism. The original model predicts the lethality for a person in front of a reflecting surface caused by simple Friedlander blast waves. Bowen extended the applicability to persons in prone position and standing in the free field by making assumptions about the pressure dose at these positions. Based on new experimental data, some authors recently concluded that the lethality for a person standing in the free field is the same as for a person in front of a reflecting surface, contrary to Bowen's assumptions.In this article, we show that only for a short duration blast wave, the load on a person standing in the free field is comparable to that on a person in front of a reflecting surface. For long positive phase durations, a safe and conservative assumption is that the load on a person standing in the free field is the sum of the side-on overpressure and the dynamic pressure. This hypothesis is supported by common knowledge about blast waves and is illustrated with numerical blast simulations.In a step by step derivation we present a new standard for the prediction of lethality caused by Friedlander blast waves, which will be included in the NATO Explosives Safety Manual AASTP-4. The result is a comprehensive engineering model that can be easily applied in calculations.