Numerical simulation of far-field blast loads arising from large TNT equivalent explosives

Large TNT equivalent explosions usually arise from accidents occurring during the transportation, storage, and manufacturing of chemicals relevant to process industries. The blast wave generated by the explosion will spread and interact with the surrounding factories and storehouses, damaging the building structures within several kilometers and causing significant casualties and property losses. This study aims to develop an efficient numerical simulation method to predict blast loads to estimate the consequences of accidents involving far-field free air bursts or surface burst explosions. Before its interaction with the interested target, a blast wave is generated in the numerical model by specifying the initial and boundary conditions of the disturbed air. Based on empirical data of incident overpressure, an explicit formula to calculate the air particle velocity is derived from the governing equations of a perfect inviscid gas. A simplified path line method is proposed to calculate the air density. The proposed method is applied to the LS-DYNA CESE solver to simulate the blast loads on building structures in the far field. Validations against empirical data and experiments indicate that the proposed method is sufficiently accurate for engineering applications and, through a case study, presents a more efficient performance than the LOAD_BLAST_ENHANCED (LBE) and mapping methods.     

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.     

Experimental and numerical investigation of constant volume dust and gas explosions in a 3.6-m flame acceleration tube

This paper describes an experimental investigation of turbulent flame propagation in propane-air mixtures, and in mechanical suspensions of maize starch dispersed in air, in a closed vessel of length 3.6 m and internal cross-section 0.27 m × 0.27 m. The primary motivation for the work is to gain improved understanding of turbulent flame propagation in dust clouds, with a view to develop improved models and methods for assessing explosion risks in the process and mining industries. The study includes computational fluid dynamics (CFD) simulations with FLACS and DESC, for gas and dust explosions respectively. For initially quiescent propane-air mixtures, FLACS over-predicts the rate of combustion for fuel-lean mixtures, and under-predicts for fuel-rich mixtures. The simulations tend to be in better agreement with the experimental results for initially turbulent gaseous mixtures. The experimental results for maize starch vary significantly between repeated tests, but the subset of tests that yields the highest explosion pressures are in reasonable agreement with CFD simulations with DESC.     

A quantitative correlation of evaluating the flame speed for the BST method in vapor cloud explosions

In recent decades, vapor cloud explosions (VCEs) have occurred frequently and resulted in numerous personnel injuries and large property losses. As a main concern in the petrochemical industry, it is of great importance to assess the consequence of VCEs. Currently, the TNT equivalency method (TNT EM), the TNO multi-energy method (TNO MEM), and the Baker-Strehlow-Tang (BST) method are widely used to estimate the blast load from VCEs. The TNO MEM and BST method determine the blast load from blast curves based on the class number and the flame speed, respectively. To quantitatively evaluate the flame speed for the BST method, the experimental data is adopted to validate the confinement specific correlation (CSC) for the determination of the class number in the TNO MEM. As a bridge, a quantitative evaluation correlation (QEC) between CSC correlation and the flame speed is established and the blast wave shapes corresponding to different flame speeds are proposed. CFD software FLACS was used to verify the quantitative correlation with the numerical models of three geometrical scales. It is found that the calculated flame speeds by the QEC are in good agreement with the simulated ones. A petrochemical plant is selected as a realistic scenario to analyze the TNT EM, TNO MEM, BST method and FLACS simulations in terms of the positive-phase side-on overpressure and impulse at different distances. Compared with the flame speed table, the predicted overpressure from BST curves determined by the proposed QEC is closer to that from FLACS and more conservative. Furthermore, the predicted results of different methods are compared with each other. It is found that the estimated positive-phase side-on overpressure and impulse by the TNO MEM are the largest, and the estimated impulse by the TNT EM is the smallest. Moreover, the estimated overpressure and impulse are larger in the higher reactivity gas.     

Evaluating the potential for overpressures from the ignition of an LNG vapor cloud during offloading

Ignition of natural gas (composed primarily of methane) is generally not considered to pose explosion hazards when in unconfined and low- or medium-congested areas, as most of the areas within LNG regasification facilities can typically be classified. However, as the degrees of confinement and/or congestion increase, the potential exists for the ignition of a methane cloud to result in damaging overpressures (as demonstrated by the recurring residential explosions due to natural gas leaks). Therefore, it is prudent to examine a proposed facility’s design to identify areas where vapor cloud explosions (VCEs) may cause damage, particularly if the damage may extend off site.An area of potential interest for VCEs is the dock, while an LNG carrier is being offloaded: the vessel hull provides one degree of confinement and the shoreline may provide another; some degree of congestion is provided by the dock and associated equipment.In this paper, the computational fluid dynamics (CFD) software FLACS is used to evaluate the consequences of the ignition of a flammable vapor cloud from an LNG spill during the LNG carrier offloading process. The simulations will demonstrate different approaches that can be taken to evaluate a vapor cloud explosion scenario in a partially confined and partially congested geometry.     

Review of the DESC project

Dust Explosion Simulation Code (DESC) was a project supported by the European Commission under the Fifth Framework Programme. The main purpose of the project was to develop a simulation tool based on computational fluid dynamics (CFD) that could predict the potential consequences of industrial dust explosions in complex geometries. Partners in the DESC consortium performed experimental work on a wide range of topics related to dust explosions, including dust lifting by flow or shock waves, flame propagation in vertical pipes, dispersion-induced turbulence and flame propagation in closed vessels, dust explosions in closed and vented interconnected vessel systems, and measurements in real process plants. The new CFD code DESC is based on the existing CFD code FLame ACceleration Simulator (FLACS) for gas explosions. The modelling approach adopted in the first version entails the extraction of combustion parameters from pressure–time histories measured in standardized 20-l explosion vessels. The present paper summarizes the main experimental results obtained during the DESC project, with a view to their relevance regarding dust explosion modelling, and describes the modelling of flow and combustion in the first version of the DESC code. Capabilities and limitations of the code are discussed, both in light of its ability to reproduce experimental results, and as a practical tool in the field of dust explosion safety.     

Lithium-ion energy storage battery explosion incidents

Utility-scale lithium-ion energy storage batteries are being installed at an accelerating rate in many parts of the world. Some of these batteries have experienced troubling fires and explosions. There have been two types of explosions; flammable gas explosions due to gases generated in battery thermal runaways, and electrical arc explosions leading to structural failure of battery electrical enclosures. The thermal runaway gas explosion scenarios, which can be initiated by various electrical faults, can be either prompt ignitions soon after a large flammable gas mixture is formed, or delayed ignitions associated with late entry of air and/or loss of gaseous fire suppression agent. The electrical explosions have entailed inadequate electrical protection to prevent high energy arcs within electrical boxes vulnerable to arc induced high pressures and thermal loads. Estimates of both deflagration pressures and arc explosion pressures are described along with their incident implications.     

Modelling large-scale vented gas explosions in a twin-compartment enclosure

Natural gas and LPG are common fuels that have been used relatively safely in the home for many decades. However, when there is a release of gas within a dwelling, or gas from a leaking external pipeline migrates into a building, an explosion may occur. Most of the experimental research into vented gas explosions has been conducted in single enclosure, cuboid or spherical geometries which are not representative of accidental explosions in dwellings or process industries. This paper discusses the findings of a comprehensive large-scale experimental programme undertaken by British Gas Research and Development and also compares FLACS CFD (Computational Fluid Dynamics) simulations against a number of these experiments. The results suggest that the software is useful in gaining a greater understanding of the dynamics of explosion development in dwellings. The paper highlights areas of good performance of the software as well as areas of shortcomings where further understanding and modelling effort is needed.     

New investigation findings on the 2006 Danvers,MA explosion

On November 22, 2006 the largest explosion in the history of Massachusetts occurred in Danvers, MA at approximately 2:46 am. This paper presents a detailed analysis into the potential causes and lessons learned from the Danvers explosion. Other investigative groups concluded that the cause of the explosion was an overheated production tank. However, the analyses presented here demonstrate that their proposed scenario could not have occurred and that other potential causes are more likely.Using the computational fluid dynamics tool FLACS, it was possible to investigate the chain of events leading to the explosion, including: (1) evaluating various leak scenarios by modeling the dispersion and mixing of gases and vapors within the facility, (2) evaluating potential ignition sources within the facility of the flammable fuel–air mixture, and (3) evaluating the explosion itself by comparing the resulting overpressures of the exploding fuel–air cloud with the structural response of the facility and the observed near-field and far-field blast damage. These results, along with key witness statements and other analyses, provide valuable insight into the likely cause of this incident. Based on the results of our detailed analysis, lessons learned regarding the investigative procedure and methods for mitigating this and future explosions are discussed.     

External pressures generated by vented gas and dust explosions

Several different data correlations have been developed for the external pressures associated with vented gas explosions and dust explosions. These correlations, which are applicable to external locations in the direct line-of-sight of the enclosure vent, are reviewed here. In addition, the application of spherically symmetric and of ellipsoidal blast wave models is explored as a possible means of calculating external pressures over a wider range of conditions than is possible with the existing data correlations. Results indicate that the spherically symmetric blast wave model can obtain a comparable accuracy (8–9 kPa standard deviation) for line-of-sight locations as the more recent data correlations. In the case of the lower blast pressures at locations perpendicular to the vent line-of-sight, the ellipsoidal blast wave provides significantly better agreement with data (to within 1 kPa standard deviation for the one set of available test data) than the spherically symmetric model.     

Validation of FLACS against experimental data sets from the model evaluation database for LNG vapor dispersion

The siting of facilities handling liquefied natural gas (LNG), whether for liquefaction, storage or regasification purposes, requires the hazards from potential releases to be evaluated. One of the consequences of an LNG release is the creation of a flammable vapor cloud, that may be pushed beyond the facility boundaries by the wind and thus present a hazard to the public. Therefore, numerical models are required to determine the footprint that may be covered by a flammable vapor cloud as a result of an LNG release. Several new models have been used in recent years for this type of simulations. This prompted the development of the “Model evaluation protocol for LNG vapor dispersion models” (MEP): a procedure aimed at evaluating quantitatively the ability of a model to accurately predict the dispersion of an LNG vapor cloud.This paper summarizes the MEP requirements and presents the results obtained from the application of the MEP to a computational fluid dynamics (CFD) model – FLACS. The entire set of 33 experiments included in the model validation database were simulated using FLACS. The simulation results are reported and compared with the experimental data. A set of statistical performance measures are calculated based on the FLACS simulation results and compared with the acceptability criteria established in the MEP. The results of the evaluation demonstrate that FLACS can be considered a suitable model to accurately simulate the dispersion of vapor from an LNG release.     

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.     

Computational modeling of vapor cloud explosions in off-shore rigs using a flame-speed based combustion model

Computational modeling is a useful tool in determining the consequences from vapor cloud explosions. Here an approach that uses a flame-speed based combustion model is evaluated. Various scenarios of explosions in full-scale off-shore modules are simulated and compared to available test data. The ignition location of the cloud and available venting paths are found to affect the overpressure field in and outside the module. For end ignition cases, the combustion of gas pushed out of the module is found to play a key role. Using the flame-speed based model with appropriate effective flame speeds is found to provide accurate simulations.     

Validation of CFD-model for hydrogen dispersion

To be able to perform proper consequence modelling as a part of a risk assessment, it is essential to be able to model the physical processes well. Simplified tools for dispersion and explosion predictions are generally not very useful. CFD tools have the potential to model the relevant physics and predict well, but without proper user guidelines based on extensive validation work, very mixed prediction capability can be expected. In this article, recent dispersion validation effort for the CFD tool FLACS–HYDROGEN is presented. A range of different experiments is simulated, including low-momentum releases in a garage, subsonic jets in a garage with stratification effects and subsequent slow diffusion, low momentum and subsonic horizontal jets influenced by buoyancy, and free jets from high-pressure vessels. LH₂ releases are also considered. Some of the simulations are performed as blind predictions.     

Consequence prediction for dust explosions involving interconnected vessels using computational fluid dynamics modeling

Combustible dust explosions continue to present a significant threat toward operating personnel and pneumatic conveyance equipment in a wide variety of processing industries. Following ignition of suspended fuel within a primary enclosure volume, propagation of flame and pressure fronts toward upstream or downstream interconnected enclosures can result in devastating secondary explosions if not impeded through an appropriate isolation mechanism. In such occurrences, an accelerated flame front may result in flame jet ignition within the secondary vessel, greatly increasing the overall explosion severity. Unlike an isolated deflagration event with quantifiable reduced pressures (vent sizing according to NFPA 68 guidance), oscillation of pressure between primary and secondary process vessels leads to uncertain overpressure effects. Dependent on details of the application such as relative enclosure volumes, relief area, fuel type, suspended concentration, duct size, and duct length, the maximum system pressure in both interconnected vessels can be unpredictable. This study proposes the use of FLame ACceleration Simulator (FLACS) computational fluid dynamics (CFD) modeling to provide reliable consequence predictions for specific case scenarios of dust deflagrations involving interconnected equipment. Required minimum supplement to the originally calculated relief area (A_v) was determined through iterative simulation, allowing for reduced explosion pressures (P_red) to be maintained below theoretical enclosure design strengths (P_es).     

Predictive correlations for leaking heat transfer fluid aerosols in air

Mist or aerosol explosions present a serious hazard to process industries. Heat transfer fluids are widely used in the chemical process industry, are flammable above their flash points, and can cause aerosol explosions. Though the possibility of aerosol explosions has been widely documented, knowledge about their explosive potential is limited. Studying the formation of such aerosols by emulating leaks in process equipment will help define a source term for aerosol dispersions and aid in characterizing their explosion hazards.Current research conducted at the Mary Kay O’Connor Process Safety Center involves the non-intrusive measurement of heat transfer fluid aerosol sprays using a Malvern Diffraction Particle Analyzer. A predictive correlation relating aerosol droplet diameters to bulk liquid pressures, temperatures, thermal and fluid properties, leak sizes, and ambient conditions is presented. This correlation can be used to predict the conditions under which leaks will result in the formation of aerosols and ultimately help in estimating the explosion hazards of heat transfer fluid aerosols. The goal is to provide information that will help improve safety in process industries.     

Effect of ignition position on the run-up distance to DDT for hydrogen-air explosions

The method described in this paper enabled reliable and accurate positioning of an overdriven detonation by calculation of shock wave velocities (detonation and retonation) for hydrogen explosions in a closed 18 m long horizontal DN150 pipe. This enabled an empirical correlation between the ignition position and the run-up distance to DDT to be determined. It was shown that the initial ability of the flame to expand unobstructed and the piston-like effect of burnt gas expanding against the closed end of the tube contributed to initial flame acceleration and hence were able to affect the run-up distance to overdriven detonation. Flame speeds and rates of initial pressure rise were also used to explain how these two competing effects were able to produce a minimum in the run-up distance to DDT. The shortest run-up distance to DDT, relative to the ignition position, for this pipe and gas configuration was found when the ignition position was placed 5.6 pipe diameters (or 0.9 m) from the closed pipe end. The shortest run-up distance to DDT relative to the end of the pipe was recorded when the ignition source was placed 4.4 pipe diameters or 0.7 m from the pipe end.     

Large-scale dust explosions in vessel-pipe systems

This study investigates dust explosions in vessel-pipe systems to develop a better understanding of dust flame propagation between interconnected vessels and implications for the proper application of explosion isolation systems. Cornstarch dust explosions were conducted in a large-scale setup consisting of a vented 8-m³ vessel and an attached pipe with a diameter of 0.4 m and a length of 9.8 m. The ignition location and effective dust reactivity were varied between experiments. The experimental results are compared against previous experiments with initially quiescent propane-air mixtures, demonstrating a significantly higher reactivity of the dust explosions due to elevated initial turbulence, leading to higher peak pressures and faster flame propagation. In addition, a physics-based model developed previously to predict gas explosion dynamics in vessel-pipe systems was extended for dust combustion. The model successfully predicts the pressure transients and flame progress recorded in the experiments and captures the effects of ignition location and effective dust reactivity.     

An alternative CFD tool for gas dispersion modelling of heavy gas

The numerical simulation of gas dispersion is of great importance in various areas of engineering such as optimisation, synthesis of chemical process, petroleum industry and process safety. The OpenFOAM (Open Field Operation and Manipulation) code is a free and open source computational fluid dynamics (CFD) program. The current research is focused on the development and customisation of a computational tool for handling gas dispersion of heavy gases, such a LNG and CO₂. The novel CFD tool relies on OpenFOAM framework. The core of the work is based on the OpenFOAM solver rhoReactingBuoyantFoam to handle gas dispersion. A series of CFD simulations has been performed for methane and CO₂. The source term of the former is modelled by HSM (Hybrid Switch Model). The model comprises contribution from HEM (Homogeneous Equilibrium Model) approach, frozen model and non-equilibrium model for CO₂ leak. The novel approach switches between equilibrium and non-equilibrium conditions based on the meta-stable parameter on the grounds of thermodynamics and experimental observations. Good agreement with experimental data is observed. Numerical findings for methane leakage from the proposed CFD tool are compared with experimental data and FLACS. Good agreement is observed.