A method for simulation of vapour cloud explosions based on computational fluid dynamics (CFD)

The effectiveness of the application of CFD to vapour cloud explosion (VCE) modelling depends on the accuracy with which geometrical details of the obstacles likely to be encountered by the vapour cloud are represented and the correctness with which turbulence is predicted. This is because the severity of a VCE strongly depends on the types of obstacles encountered by the cloud undergoing combustion; the turbulence generated by the obstacles influences flame speed and feeds the process of explosion through enhanced mixing of fuel and oxidant. In this paper a CFD-based method is proposed on the basis of the author’s finding that among the various models available for assessing turbulence, the realizable k-? model yields results closer to experimental findings than the other, more frequently used, turbulence models if used in conjunction with the eddy-dissipation model. The applicability of the method has been demonstrated in simulating the dispersion and ignition of a typical vapour cloud formed as a result of a spill from a liquid petroleum gas (LPG) tank situated in a refinery. The simulation made it possible to assess the overpressures resulting from the combustion of the flammable vapour cloud. The phenomenon of flame acceleration, which is a characteristic of combustion enhanced in the presence of obstacles, was clearly observed. Comparison of the results with an oft-used commercial software reveals that the present CFD-based method achieves a more realistic simulation of the VCE phenomena.     

A novel approach combining bootstrapped non-intrusive reduced order models and unscented transform for the robust and efficient CFD analysis of accidental gas releases in congested plants

The risk assessment for safety-critical, complex systems is a very challenging computational problem when it is performed with high-fidelity models, e.g. CFD, like in the case of accidental gas releases in congested systems. Within this framework, a novel CFD approach, named Source Box Accident Model, has been recently proposed to efficiently model such phenomena by splitting the simulation of the gas release and its subsequent dispersion in the system in two steps. In this view, the present paper proposes a non-intrusive, Proper Orthogonal Decomposition-Radial Basis Functions reduced order model that exploits the two-step nature of the SBAM approach, to mimic the behaviour of the original, long-running CFD model code at a significantly lower computational cost. Moreover, the paper presents a methodology combining the bootstrap and unscented transform approaches to efficiently assess the ROM uncertainty in the safety-critical simulation output quantities of interest, e.g. the flammable volume. The results obtained in a test case involving a high pressure, accidental gas release in an off-shore Oil & Gas plant are in very satisfactory agreement with those produced by CFD, with a relative error smaller than 10% and a reduction in the computational time of about three orders of magnitude.     

LNG accident dynamic simulation: Application for hazardous consequence reduction

The evaluation of exclusion (hazard) zones around the LNG stations is essential for risk assessment in LNG industry. In this study, computational fluid dynamics (CFD) simulations have been conducted for the two potential hazards, LNG flammable vapor dispersion and LNG pool fire radiation, respectively, to evaluate the exclusion zones. The spatial and temporal distribution of hazard in complex spill scenario has been taken into account in the CFD model. Experimental data from Falcon and Montoir field tests have been used to validate the simulation results. With the valid CFD model, the mitigation of the vapor dispersion with spray water curtains and the pool fire with high expansion foam were investigated. The spray water curtains were studied as a shield to prevent LNG vapor dispersing, and two types of water spray curtain, flat and cone, were analyzed to show their performance for reduction and minimization of the hazard influencing distance and area. The high expansion foam firefighting process was studied with dynamic simulation of the foam action, and the characteristics of the foam action on the reduction of LNG vaporization rate, vapor cloud and flame size as well as the thermal radiation hazard were analyzed and discussed.     

Use of water spray curtain to disperse LNG vapor clouds

Effective safety measures to prevent and mitigate the consequences of an accidental release of flammable LNG are critical. Water spray curtain is currently recognized as an effective technique to control and mitigate various hazards in the industries. It has been used to absorb, dilute and disperse both toxic and flammable vapor cloud. It is also used as protection against heat radiation, in case of fighting vapor cloud fire. Water curtain has also been considered as one of the most economic and promising LNG vapor cloud control techniques. Water curtains are expected to enhance LNG vapor cloud dispersion mainly through mechanical effects, dilution, and thermal effects. The actual phenomena involved in LNG vapor and water curtain interaction were not clearly established from previous research. LNG spill experiments have been performed at the Brayton Fire Training Field at Texas A&M University (TAMU) to understand the effect of water curtain in controlling and dispersing LNG vapor cloud. This paper summarizes experimental methodology and presents data from two water curtain tests. The analysis of the test results are also presented to identify the effectiveness of these two types of water spray curtains in enhancing the LNG vapor cloud dispersion.     

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.     

Validation of geometry modelling approaches for offshore gas dispersion simulations

Computational Fluid Dynamics (CFD) codes are widely used for gas dispersion studies on offshore installations. The majority of these codes use single-block Cartesian grids with the porosity/distributed-resistance (PDR) approach to model small geometric details. Computational cost of this approach is low since small-scale obstacles are not resolved on the computational mesh. However, there are some uncertainties regarding this approach, especially in terms of grid dependency and turbulence generated from complex objects. An alternative approach, which can be implemented in general-purpose CFD codes, is to use body-fitted grids for medium to large-scale objects whilst combining multiple small-scale obstacles in close proximity and using porous media models to represent blockage effects. This approach is validated in this study, by comparing numerical predictions with large-scale gas dispersion experiments carried out in DNV GL's Spadeadam test site. Gas concentrations and gas cloud volumes obtained from simulations are compared with measurements. These simulations are performed using the commercially available ANSYS CFX, which is a general-purpose CFD code. For comparison, further simulations are performed using CFX where small-scale objects are explicitly resolved. The aim of this work is to evaluate the accuracy and efficiency of these different geometry modelling approaches.     

A CFD-based approach for gas detectors allocation

Accidental gas releases are detected by allocating sensors in optimal places to prevent escalation of the incident. Gas release effects are typically assessed based on calculating the dispersion from releasing points. In this work, a CFD-based approach is proposed to estimate gas dispersion and then to obtain optimal gas sensors allocation. The Ansys-Fluent commercial package is used to estimate concentrations in the open air by solving the governing equations of continuity, momentum, energy and species convection-diffusion combined with the realizable κ-ε model for turbulence viscosity effects. CFD dynamic simulations are carried out for potential gas leaks, assuming worst-case scenarios with F-stability and 2 m/s wind speed during a 4 min releasing period and considering 8 wind directions. The result is a scenario-based methodology to allocate gas sensors supported on fluid dynamics models. The three x–y–z geographical coordinates for the sensor allocation are included in this analysis. To highlight the methodology, a case study considers releases from a large container surrounded by different types of geometric units including sections with high obstacles, low obstacles, and no obstacles. A non-redundant set of perfect sensors are firstly allocated to cover completely the detection for all simulations releases. The benefits of redundant detection via a MooN voting arranging scheme is also discussed. Numerical results demonstrate the capabilities of CFD simulations for this application and highlight the dispersion effects through obstacles with different sizes.     

Dispersion of the vapour cloud in the Buncefield Incident

Dispersion of the flammable vapour cloud in the 2005 Buncefield Incident is examined. Footage from security cameras around the site is analysed and the results from Computational Fluid Dynamics (CFD) simulations of the vapour dispersion are presented. It is shown that the shape of the terrain and the presence of obstacles significantly affected the dispersion of vapour from the overflowing tank. The CFD model is shown to produce similar qualitative behaviour to that observed in the incident, both in terms of the arrival time of the vapour cloud and its final depth.     

Numerical analysis of toxic cloud generation and dispersion: A case study of the ethylene oxide spill

The present study examined the accidental spill of ethylene oxide, and a sensitivity analysis of the corresponding consequences was conducted using computational fluid dynamics (CFD). A validation of the gas dispersion CFD model against the experimental data sets included in the model evaluation protocol (MEP) was performed. The effect of the variability of the wind velocity on the extension of the hazardous areas and pool evaporation characteristics was evaluated. Additionally, the mitigation effects of the dike walls surrounding a spill were discussed. CFD simulation results have shown that the mitigation effect of dike walls is determined by their influence on both gas dispersion and pool evaporation and depends strongly on wind velocity in terms of toxic impact distances.     

Numerical study on the initial expansion of two-phase cloud from an instantaneous release

In industries some dangerous liquefied gases may accidentally release and it may form a flammable or toxic mixture after mixing with air. One tool that is being developed in industry for two-phase cloud dispersion modeling is computational fluid dynamics (CFD). In this paper, the dispersion processes of different dangerous materials including liquefied chlorine, liquefied ammonia and liquefied petroleum gas were simulated in the same condition to analyze the characteristics of the initial expansion processes by CFD tool. The heat and mass transfer between droplets and the vapor after an instantaneous release event was calculated by using the Eulerian–Lagrangian method. The results from a number of 3-D CFD based studies were compared with the available small-scale experimental results. The results show that the present model and numerical simulation are reliable.     

Vent Sizing Criteria for Partial Volume Deflagration

The accidental spill of volatile solvents or the release of flammable gases within equipment and buildings is likely to form fuel concentration gradients unless efficient mixing is provided. As a consequence, even small amounts of fuel can form flammable clouds, and partial volume deflagrations may occur. Nevertheless, few indications are given in international guidelines for vent sizing and only over-conservative well-mixed stoichiometric assumptions are used. In this paper, we propose a predictive methodology for the evaluation of the dynamics of partial volume deflagration, aiming at defining useful correlations for the design of vent devices, starting from the fundamental equation for the rate of pressure rise and flame propagation in closed vessel. We define a ‘stratified gas deflagration index’ K_G(m), where m is the filling ratio, and use it with the most common design equations for vent sizing. The approach has been validated by means of a CFD code for the simulation of stratified laminar methane–air explosion by varying both filling ratio and volume.     

A CFD based explosion risk analysis methodology using time varying release rates in dispersion simulations

A full probabilistic Explosion Risk Analysis (ERA) is commonly used to establish overpressure exceedance curves for offshore facilities. This involves modelling a large number of gas dispersion and explosion scenarios. Capturing the time dependant build up and decay of a flammable gas cloud size along with its shape and location are important parameters that can govern the results of an ERA. Dispersion simulations using Computational Fluid Dynamics (CFD) are generally carried out in detailed ERA studies to obtain these pieces of information. However, these dispersion simulations are typically modelled with constant release rates leading to steady state results. The basic assumption used here is that the flammable gas cloud build up rate from these constant release rate dispersion simulations would mimic the actual transient cloud build up rate from a time varying release rate. This assumption does not correctly capture the physical phenomena of transient gas releases and their subsequent dispersion and may lead to very conservative results. This in turn results in potential over design of facilities with implications on time, materials and cost of a project.In the current work, an ERA methodology is proposed that uses time varying release rates as an input in the CFD dispersion simulations to obtain the fully transient flammable gas cloud build-up and decay, while ensuring the total time required to perform the ERA study is also reduced. It was found that the proposed ERA methodology leads to improved accuracy in dispersion results, steeper overpressure exceedance curves and a significant reduction in the Design Accidental Load (DAL) values whilst still maintaining some conservatism and also reducing the total time required to perform an ERA study.     

Explosion hazards related to hydrogen releases in nuclear facilities

Hydrogen explosion risk needs to be carefully assessed and evaluated in nuclear facilities because of the potential catastrophic consequences: breakdown of safety equipments, failure of containment, dissemination of radioactive materials in the environment.When studying an indoor release, one possible simplification is to assume a perfect gas mixing inside the room. This assumption is effectively often used to evaluate toxic risks in the environment outside a building (Mastellone, Ponte, & Arena, 2003). However, perfect gas mixing assumption is only a rough approximation, as indoor concentrations can largely differ from mean values, due to buoyancy, recirculation zones or obstacles for example.In order to better evaluate the risk of explosion in case of an accidental release of hydrogen, IRSN conducted a numerical study using FLACS CFD software. Several parameters have been studied to identify dangerous situations and draw a representative picture of the risk: room size, position and direction of hydrogen leak, ventilation characteristics. Hydrogen release flow rates used for numerical simulations have been chosen as the highest leak rate which, by applying the assumption of perfect mixing, produces an average concentration in the room equal to hydrogen lower flammability limit (LFL).Simulation results indicate that in some particular configurations, especially for impinging hydrogen jets, hydrogen concentrations can locally be above LFL and then create explosive atmospheres with significant volumes.     

Comparing tools of varying complexity for calculating the gas dispersion

When handling flammable and/or toxic liquids or gases, the gas dispersion following a release of substance is a scenario to be considered in the risk assessment to determine the lower flammability distance (LFD) and toxicity thresholds. In this work a comparison of different gas dispersion tools of varying complexity ranging from a simple Gaussian model over a boundary layer model (BLM) and a Lagrangian model to CFD (in this case ANSYS CFX v14) is presented. The BLM covers the special case of liquid releases with formation of a pool. It does not only solve the gas dispersion but also calculates the evaporating mass flow out of the pool. The simulation values are compared to each other and to experimental data resulting mainly from our own open air experiments covering the near field and carried out on the Test Site Technical Safety of BAM (BAM-TTS) for different release types (pool evaporation, gas release) and topologies. Other validation data were taken from literature and cover large scale experiments in the range of several 100 m.     

Modified porosity approach and laminar flamelet modelling for advanced simulation of accidental explosions

The production of oil and gas is an inherently hazardous task. Therefore it is crucial to provide reliable estimates of the risks involved. The major contributors to the risk level of an offshore installation, for example, arise from accidents involving explosion and fire. Computational Fluid Dynamics (CFD) can be a powerful tool to help with the calculation of accidental explosion scenarios. In this context, the present work suggests a novel implementation of a model based on a modified Porosity Distributed Resistance (MPDR) approach within an unstructured 3D Navier–Stokes solver. The model operates by representing parts of the filtered geometry from the original model through porosity values attributed to an unstructured tetrahedral mesh. Extra resistance terms are added in the momentum equation as well as extra sources of turbulence. Two extra sources of turbulence are modelled. The first of these is due to the shear layers of the non-resolved obstacles, whilst the second is due to the presence of wakes behind the non-resolved obstacles. Results for 2D and 3D test cases are compared against both experimental data and simulations with fully-resolved geometry and good agreement is observed.     

Dynamic simulation of hazard analysis of radiations from LNG pool fire

Because of its highly flammable nature, any accidental release of liquefied natural gas (LNG) could possibly pose significant fire hazard. In this study, a computational fluid dynamics (CFD) model was used to analyze this hazard around an existing LNG station. By assuming an LNG pool fire occurring in an impoundment area, dynamic simulations of flame development have been carried out. In order to provide more reliable simulation results, a study was first conducted to determine the mesh independence and suitable time step. The results of CFD simulations were also compared with those using the commonly-used phenomenological model. The simulation results showed that LNG tanks in the neighbor dike area could withstand the received radiant heat flux, and the areas involving human activities, such as security office and public area, were also secure enough for people to escape from the hazards. LNG vaporizers, which are often located close to tank area, could possibly receive relatively higher radiant heat flux. High temperature achieved on vaporizers could cause material failure. CFD calculations have also indicated that increasing the spacing distance or using flowing water curtain could reduce this temperature. It is concluded that CFD method is significantly more effective to account for LNG hazard analysis and provide realistic results for complicated scenarios, thus providing meaningful information for safety consideration.     

穿越城市生活区的天然气管道泄漏连锁爆燃后果评估研究*

为评估城市天然气管道泄漏连锁爆燃事故后果,基于计算流体力学(CFD)方法构建穿越城市区域的天然气管道泄漏连锁爆燃后果预测与评估模型,以某城市生活区域为例,在城市生活区域建筑物内风场流动计算的基础上,模拟风场作用下可燃气体在城市建筑物空间内的运移规律,预测可燃气云的积聚区域;考虑意外点火的情况,计算城市生活区域内可燃气云爆燃灾害特征,预测爆燃超压、热辐射和高温的影响。研究结果表明:由于建筑物之间的阻挡与反射作用,建筑物下风向有明显的低风速区域,并在一定时间段后扩散过程趋于稳定;在爆燃火焰作用下,高温和热辐射会造成建筑物部分钢结构发生失效变形。     

大尺度障碍物与泄爆面对天然气内爆炸的协同作用规律研究

为阐释民用建筑内部大尺度物品与门窗等泄爆面对天然气爆炸灾害的协同作用机制，基于典型厨房空间布局及内部物品特征，借助计算流体动力学技术研究了不同泄爆面开启压力和不同大尺度障碍物体积阻塞率条件下天然气内爆炸火焰速度、爆炸超压的分布规律。研究结果表明:大尺度障碍物与泄爆面对室内天然气爆炸过程具有显著的协同作用，共同促进火焰速度与爆炸超压的显著增长，并缩短峰值超压到达时间；大尺度障碍物的存在虽然显著降低了室内天然气的体积，但从增加房间内湍流源和相对长径比的角度进一步促进了泄爆效应；大尺度障碍物与泄爆面协同作用下，室内火焰速度呈现明显的阶段性特征，并在泄爆面附近发生波动。研究结论可为民用建筑物内气体爆炸事故调查分析和灾害评估提供科学依据。     

Un-ignited and ignited high pressure hydrogen releases: Concentration – turbulence mapping and overpressure effects

Safety studies for production and use of hydrogen reveal the importance of accurate prediction of the overpressure effects generated by delayed explosions of accidental high pressure hydrogen releases. Analysis of previous experimental work demonstrates the lack of measurements of turbulent intensities and lengthscales in the flammable envelope as well as the scarceness of accurate experimental data for explosion overpressures and flame speeds. AIR LIQUIDE, AREVA STOCKAGE ENERGIE and INERIS join in a collaborative project to study un-ignited and ignited high pressure releases of hydrogen.The purpose of this work is to map hydrogen flammable envelopes in terms of concentration, velocity and turbulence, and to characterize the flame behaviour and the associated overpressure. These experimental results (dispersion and explosion) are also compared with blind FLACS modelling.