Numerical simulation of LNG dispersion under two-phase release conditions

The use of LNG (liquefied natural gas) as fuel brings up issues regarding safety and acceptable risk. The potential hazards associated with an accidental LNG spill should be evaluated, and a useful tool in LNG safety assessment is computational fluid dynamics (CFD) simulation. In this paper, the ADREA-HF code has been applied to simulate LNG dispersion in open-obstructed environment based on Falcon Series Experiments. During these experiments LNG was released and dispersed over water surface. The spill area is confined with a billboard upwind of the water pond. FA1 trial was chosen to be simulated, because its release and weather conditions (high total spill volume and release rate, low wind speed) allow the gravitational force to influence the cold, dense vapor cloud and can be considered as a benchmark for LNG dispersion in fenced area. The source was modeled with two different approaches: as vapor pool and as two phase jet and the predicted methane concentration at sensors' location was compared with the experimental one. It is verified that the source model affect to a great extent the LNG dispersion and the best case was the one modeling the source as two phase jet. However, the numerical results in the case of two phase jet source underestimate the methane concentration for most of the sensors. Finally, the paper discusses the effect of neglecting the ?9.3° experimental wind direction, which leads to the symmetry assumption with respect to wind and therefore less computational costs. It was found that this effect is small in case of a jet source but large in the case of a pool source.     

Numerical simulation of propane detonation in medium and large scale geometries

A modelling strategy has been developed for consequence analysis of medium and large scale gaseous detonation. The model is based on the solution of Euler equations with one-step chemistry. The Van Leer flux limited method which is a total variation diminishing scheme is used for shock capturing. Preliminary calculations were firstly conducted for small domains with fine grids which resolve the wave, relatively coarse grids which have less than 10 grids across the wave and coarse grids in which the minimum grid size is larger than the wave thickness to ensure that the reaction scheme has been properly tuned to capture the correct detonation pressure, temperature and velocity in the resolutions used in the different cases. The model was firstly tested against a medium scale detonation test in a shock tube with U-bends. Reasonably good agreement is achieved on detonation pressure and mean shock wave velocities at different measuring segments of the tube. Following the validation, the detonation of a hypothetical planar propane-air cloud is simulated. The predictions uncovered some interesting features of such large scale detonation phenomena which are of significance in the safety context, especially for accidental investigations. The findings from the present analysis are in line with the forensic evidence on damages in some historic accidents and challenges previous analysis of a major accident in which forensic evidence suggested localised detonation but was considered as the consequence of fire storms by the investigation team.     

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.     

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.     

Overview of a combined computational–experimental evaluation for the assessment of panoramic sunroof impact characteristics for ejection mitigation

Abstract

Objective: To meet increasing customer demand, many vehicle manufacturers are now offering a panoramic sunroof option in their vehicle lineup. Currently, there is no regulatory or consumer test aimed at assessing the potential for ejection mitigation of roof glazing, which leaves manufacturers to develop internal performance standards to guide designs. The goal of this study was to characterize the variety of occupant-to-roof impacts involving unbelted occupants in rollover crashes to determine the ranges of possible effective masses and impact velocities. This information can be used to define occupant retention requirements and performance criteria for roof glazing in occupant ejection protection.

Methods: This study combined computational (MADYMO and LS-Dyna) simulations of occupant kinematics in rollover crashes with laboratory rollover crash tests using the dynamic rollover test system (DRoTS) and linked them through controlled anthropomorphic test device (ATD)-to-roof (“drop”) impact tests. The DRoTS and the ATD drop tests were performed to explore impact scenarios and estimate dummy-to-roof impact impulses. Next, 13 sets of vehicle kinematics and deformation data were extracted from a combination of vehicle dynamics and finite element model simulations that reconstructed variations of rollover crash cases from the field data. Then occupant kinematics data were extracted from a full-factorial sensitivity study that used MADYMO simulations to investigate how changes in anthropometry and seating position would affect occupant–roof impacts across all 13 cases. Finite element (FE) simulations of ATD and Global Human Body Models Consortium (GHBMC) human body model (HBM) roof impacts were performed to investigate the most severe cases from the MADYMO simulations to generate a distribution of head-to-roof impact energies.

Results: From the multiparameter design of experiment and experimental study, kinematics and energy output were extracted and analyzed. Based on dummy-to-roof impact force and dummy-to-roof impact velocity, the most severe rollover scenarios were identified. In the DRoTS experiments followed by the drop tests, the range of identified impact velocities was between 2 and 5.8 m/s. However, computational simulations of the rollover crashes showed higher impact velocities and similar effective masses. The largest dummy-to-roof impact velocity was 11 m/s.

Conclusions: This study combined computational and experimental analyses to determine a range of possible unbelted occupant-to-roof impact energies. These results can be used to determine design parameters for an impactor for the assessment of the risk of roof glazing ejection for unbelted occupants in rollover crashes.     

A theoretical and methodological framework for studying and modelling drivers’ mental representations

This article provides a synthetic overview of the research programme carried out at INRETS-LESCOT over the last 10 years, in view to studying and modelling the mental representations of car drivers. Theoretically, this research is in line with two complementary scientific traditions: Human Information Processing theories on the one hand, and theories of Operative Activity on the other. As discussed in Section 1 of the article, attention is given to the functional representations of drivers, constructed “by” and “for” the action, such as they are implemented in a driving situation, and taking into account both implicit and explicit dimensions. This aim directly impacts on the methodological approach implemented, in so far as it entails defining an “experimental continuum” ranging from naturalistic observations of the driving activity (on the open road), to setting up more controlled experimental protocols in order to permit in-depth, systematic and reproducible scientific investigations of drivers’ cognition. Section 2 of the article presents a synthetic view of these methods, while Section 3 presents several significant results obtained with them. Lastly, the final part of the article focuses on the computational formalism defined at INRETS-LESCOT (i.e. the driving schemas) designed to model driver knowledge and mental representations, and developed in a COgnitive Simulation MOdel of the DRIVEr called COSMODRIVE.     

Seismic fragility analysis of a coupled tank-piping system based on artificial ground motions and surrogate modeling

The catastrophic consequences of recent NaTech events triggered by earthquakes highlighted the inadequacy of standard approaches to seismic risk assessment of chemical process plants. To date, the risk assessment of such facilities mainly relies on historical data and focuses on uncoupled process components. As a consequence, the dynamic interaction between process equipment is neglected. In response to this gap, researchers started a progressive integration of the Pacific Earthquake Engineering Research Center (PEER) Performance-Based Earthquake Engineering (PBEE) risk assessment framework. However, a few limitations still prevent a systematic implementation of this framework to chemical process plants. The most significant are: (i) the computational cost of system-level simulations accounting for coupling between process equipment; (ii) the experimental cost for component-level model validation; (iii) a reduced number of hazard-consistent site-specific ground motion records for time history analyses.In response to these challenges, this paper proposes a recently developed uncertainty quantification-based framework to perform seismic fragility assessments of chemical process plants. The framework employs three key elements: (i) a stochastic ground-motion model to supplement scarcity of real records; (ii) surrogate modeling to reduce the computational cost of system-level simulations; (iii) a component-level model validation based on cost-effective hybrid simulation tests. In order to demonstrate the potential of the framework, two fragility functions are computed for a pipe elbow of a coupled tank-piping system.     

Experimental and numerical study of CO2 plume diffusion in a confined space

In this paper, heavy gas diffusion in a confined space has been investigated. The effects of barrier and source intensity on CO₂ diffusion are explored by the small-scale experiments and computational fluid dynamics (CFD) methods. Six different turbulence models are selected to predict the gas concentrations. By comparing these experimental values with the simulated ones, it is found that all models can effectively predict the concentration variation with time, and SST k-ω model is most close to the ideal model compared with others. Three source-barrier distances and three CO₂ flow rates have been set up for the study. In this confined space, the main flow is concentrated in the region near the ground. The existence of barriers in the space will have a dilution effect on the high-concentration plume near the ground of the near-source area and a barrier effect on the low-concentration plume in the far-source area. The changes in source intensity have notable impact on the gas concentrations. This study can provide an experimental basis for the risk assessment in the confined spaces, as well as an experimental and data reference for large-scale CFD simulations.     

Advanced multi-perspective computer simulation as a tool for reliable consequence analysis

Major accidents involving hazardous materials are a crucial issue for the chemical and process industries. Many accidental events taken place in the past showed that dangerous substances may pose a severe threat for people and property. Aiming at loss prevention, a series of actions have been instituted through international regulations concerning hazardous installations safety preparedness. These actions involve efficient land-use planning, safety studies execution, as well as emergency response planning drawing up. A key factor for the substantial consideration of the above is the effective prediction of possible accident forms and their consequences, for the estimation of which, a number of empirical models have been developed so far. However, (semi-)empirical models present certain deficiencies and obey to certain assumptions, thus leading to results of reduced accuracy. Another approach that could be used for this purpose and it is discussed in this work, is the utilization of advanced computational fluid dynamics (CFD) techniques in certain accident forms modeling. In particular, composite CFD-based models were developed for the simulation of several characteristic accident forms involving isothermal and non-isothermal heavy gas dispersion, confined and unconfined explosion in environment of complex geometry, as well as flammable cloud fire. The simulation cases were referred to real-scale trials allowing us to conclude about the validity of the quantitative results. Comparisons of the computational predictions with the experimental observations showed that obtained results were in good agreement with the experimental ones, whereas the evaluation of statistical performance measures proved the simulations to be statistically valid.     

Experimental and computational fluid dynamic (CFD) simulation of leak shapes and sizes for gas pipeline

In this paper, a parametric study has been carried out to predict the exit velocity of air through the leakage in the pipe with the help of CFD software ANSYS Fluent. The effect of air pressure in the pipe and the shape of leakage have been studied. Further experiments were also carried out to determine the exit velocity for the defined shape of leakage by varying the air pressure in the pipe. Experimentally, the velocity at a distance of 8 cm from the location of a leak in the horizontal plane was obtained with the help of differential pressure transducers. Using the experimental results, the computational results were validated. The results of the parametric simulation study showed that even for a pressure of 2 bars the velocity profile at the leak location indicates the supersonic state where the Mach number is greater than 1. The study is useful because it may be used as a foundation for risk assessment and safety management in the case of flammable gas leaks through gas pipes.     

Computer-aided modeling of the protective effect of explosion relief vents in tunnel structures

This paper aims at contributing to the efficient design of explosion protection systems against confined explosions. The issue addressed concerns the quantitative estimation of the protective effect of explosion relief vents in the case of confined explosions inside tunnels. A series of virtual experiments performed by computer simulation, revealed how the number of vents, their diameter, as well as the angle between the vents and the tunnel, influences the blast wave attenuation. The computational study was performed considering a complicated large-scale tunnel configuration with branches on its half portion. The purpose was the calculation of the attenuation effect due to the presence of vents by comparing the total explosion-specific impulse developing at antidiametric positions inside the tunnel. Simulations were carried out via a three-dimensional numerical model built in the computational fluid dynamics code CFX 5.7.1, which has been validated in previous papers against experimental overpressure histories data demonstrating reasonable performance. Computer results showed that the use of branch vents provides an effective method for shock wave attenuation following an explosion, whereas their statistical elaboration revealed that the attenuation is mainly affected by the number of vents and their diameter. In contrast, the angle between the side vents and the main tunnel appeared to slightly affect the pressure wave weakening. Eventually, the quantitative influence of the above parameters was effectively illustrated in functional diagrams, so that the total attenuation effect may be promptly estimated, if the design variables are known. In addition, two statistical models with reasonable fitting to the calculated data are proposed, which express the attenuation effect as a dependent variable of the design variables including their interactions.     

Numerical simulation of water curtain application for ammonia release dispersion

Using water curtain system to forced mitigate ammonia vapor cloud has been proven to be an effective measure. Currently, no engineering guidelines for designing an effective water curtain system are available, due to lack of understanding of complex interactions between ammonia vapor cloud and water droplets, especially the understanding of ammonia absorption into water droplets. This paper presents numerical calculations to reproduce the continuous ammonia release dispersion with and without the mitigating influence of a downwind water curtain using computational fluid dynamic (CFD) software ANSYS Fluent 14.0. The turbulence models k–ɛ and RNG were used to simulate the ammonia cloud dispersion without downwind water curtain. The simulated results were compared with literature using the statistical performance indicators. The RNG model represents better agreement with the experimental data and the k–ɛ model generates a slightly lesser result. The RNG model coupled with Lagrangian discrete phase model (DPM) was used to simulate the dilution effectiveness of the water curtain system. The ammonia absorption was taken into account by means of user-defined functions (UDF). The simulated effectiveness of water curtains has good agreements with the experimental results. The effectiveness of water mitigation system with and without the ammonia absorption was compared. The results display that the effectiveness mainly depends on the strong air entrainment enhanced by water droplets movement and the ammonia absorption also enhances the effectiveness of water curtain mitigation system. The study indicates that the CFD code can be satisfactorily applied in design criteria for an effective mitigation system.     

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.     

Premixed methane-air deflagrations in a completely adiabatic pipe and the effect of the condition of the pipe wall

A completely adiabatic pipe that is similar to a coal-mine coal or rock roadway was simulated using the computational software AutoReaGas. A partially adiabatic pipe was established using an experimental steel pipe with heat-insulating material installed in the inner wall, and a non-adiabatic pipe was also established using the experimental steel pipe without the heat-insulating material. Premixed methane/air deflagrations were studied in the three types of pipe to reveal the influence of the condition of the pipe wall on gas explosions. The results showed that in the completely adiabatic pipe, the maximum explosion overpressure was dynamic and decreased and increased with increasing distance; however, the flame-propagation speed increased gradually. In the partially adiabatic pipe and the non-adiabatic pipe, the maximum explosion overpressure and flame-propagation speed increased initially and then gradually decreased with increasing distance. The majority of explosion overpressure and flame-propagation speed values at each gauge in the completely adiabatic pipe were larger than those of the partially adiabatic pipe. Both measurements at each gauge in the partially adiabatic pipe were much greater than those of the non-adiabatic pipe. The condition of the pipe wall has a large influence on the maximum explosion overpressure and the flame-propagation speed. In future explosion experiments, heat insulating materials should be installed in the inner wall of steel pipes to obtain data for application to the prevention and control of gas explosions in underground coal mines.     

Propane-air mixture deflagration hazard with different ignition positions and restraints in large-scale venting chamber

In order to better assess the hazards of explosion accidents, propane-air mixture deflagrations were conducted in a large-scale straight rectangular chamber (with a cross-section of 1.5 m × 1.5 m, length of 10 m, and total volume of 22.5 m³). The effect of initial volume, ignition position, and initial restraints on the explosion characteristics of the propane-air mixtures was investigated. The explosion overpressure, flame propagation, and flame speed were obtained and the computational fluid dynamics (CFD) software was used to simulate the flame-propagation process and field flow for auxiliary analysis. The hazards of large-scale propagation explosion under weak and strong constraints were evaluated and the different phases of flame propagation under weak and strong constraints were discriminated. Results indicate that the hazards caused by propane deflagration under weak constraint are mainly caused by flame spread. And the maximum overpressure under strong constraint appeared at the front part of the chamber under the large-scale condition, which is consistent with the previous small-scale test. Moreover, the simulations of flame structures under weak and strong constraint are in good agreement with experimental results, which furthers the understanding of large-scale propane deflagration under different initial conditions in large-scale spaces and provides basic data for three-dimensional CFD model improvement.     

断面突然扩大巷道流场特征及局部阻力特性研究

为研究突扩巷道流场特征和局部阻力特性随壁面粗糙度的变化规律,采用计算流体力学(CFD)方法对实际巷道进行数值模拟,并基于流体相似理论,搭建相似比为1∶20的实验模型,使用粒子图像测速(PIV)系统等实验装置和仪器进行突扩巷道流场测试实验与局部阻力测定实验,将实验结果对模拟进行验证,根据模拟结果,对不同壁面粗糙度的突扩巷道流场特征和局部阻力特性进行分析。研究结果表明：在流场测试方面,随风速增大,突扩后涡流区长度先增大后保持不变;在局部阻力测定方面,数值模拟结果与实验测定结果相差在10%以内,随着巷道粗糙度增加,突扩巷道局部阻力系数呈非线性增大,且当突扩比为1∶2时,粗糙度分别为0.02,0.04,0.06,0.08 m的突扩巷道局部阻力系数ξ₁分别为0.373 6,0.386 3,0.395 0,0.401 6。研究结果对于新掘巷道的局部通风阻力预测工作,以及为矿井智能通风提供准确风阻参数具有重要意义。     

Vented confined explosions in Stramberk experimental mine and AutoReaGas simulation

Research Mining Institute, Inc., Ostrava-Radvanice, in cooperation with Dept. of Theory And Technology of Explosives of University of Pardubice and Klokner Institute of CTU in Prague, has performed three series of experiments examining methane–air mixture explosions and their impact on 14 and 29 cm thick wall. The project was named ‘Modeling Pressure Fields Effects on Engineering Structures During Accidental Explosions of Gases in Buildings’ and was sponsored by Grant Agency of Czech Republic (project No. 103/01/0039). The project is aimed at deeper understanding of pressure field effect upon the structures. Methane-air mixture explosion was used to generate the blast wave. The geometrical configuration of the environment resembled a room of an average size, such as larger kitchen. Preliminary simulations were made by AutoReaGas code (Century Dynamics and TNO). The design phase was followed by tests in an experimental mine in Stramberk. Two masonry dams were build in the mine, with cross-section areas of 10.2 m² and longitudinal distance of 5.7 m, creating an explosion chamber with a volume of 58 m³. Two vent openings with an adjustable free cross-section were used to control the maximum overpressure inside the chamber. The concentration of methane-air mixture was approximately 9.5% (vol.) and the volumes of the clouds were 5.25, 10.2 and 15.3 m³ respectively. The generated blast wave overpressures inside the chamber ranged between 1 and 150 kPa. According to experimental results a calibration of the code was performed. After the calibration it is possible to make relatively accurate simulations in similar geometry and to calculate the pressure loading of the structure at any spot in the simulated space. This paper describes the experiments performed and compares experimental and computational results.     

Effect of ignition position on vented hydrogen–air deflagration in a 1 m3 vessel

This study investigates the effect of the ignition position on vented hydrogen-air deflagration in a 1 m³ vessel and evaluates the performance of the commercial computational fluid dynamics (CFD) code FLACS in simulating the vented explosion of hydrogen-air mixtures. First, the differences in the measured pressure-time histories for various ignition locations are presented, and the mechanisms responsible for the generation of different pressure peaks are explained, along with the flame behavior. Secondly, the CFD software FLACS is assessed against the experimental data. The characteristic phenomena of vented explosion are observed for hydrogen-air mixtures ignited at different ignition positions, such as Helmholtz oscillation for front ignition, the interaction between external explosion and combustion inside the vessel for central ignition, and the wall effect for back-wall ignition. Flame-acoustic interaction are observed in all cases, particularly in those of front ignition and very lean hydrogen-air mixtures. The predicted flame behavior agree well with the experimental data in general while the simulated maximum overpressures are larger than the experimental values by a factor of 1.5–2, which is conservative then would lead to a safe design of explosion panels for instance. Not only the flame development during the deflagration was well-simulated for the different ignition locations, but also the correspondence between the pressure transients and flame behavior was also accurately calculated. The comparison of the predicted results with the experimental data shows the performance of FLACS to model vented mixtures of hydrogen with air ignited in a lab scale vessel. However, the experimental scale is often smaller than that used in practical scenarios, such as hydrogen refueling installations. Thus, future large-scale experiments are necessary to assess the performance of FLACS in practical use.     

LES – DFSD modelling of vented hydrogen explosions in a small-scale combustion chamber

Accidental explosions are a plausible danger to the chemical process industries. In the event of a gas explosion, any obstacles placed within the path of the flame generate turbulence, which accelerates the transient flame and raises explosion overpressure, posing a safety hazard. This paper presents numerical studies using an in-house computational fluid dynamics (CFD) model for lean premixed hydrogen/air flame propagations with an equivalence ratio of 0.7. A laboratory-scale combustion chamber is used with repeated solid obstacles. The transient compressible large eddy simulation (LES) modelling technique combined with a dynamic flame surface density (DFSD) combustion model is used to carry out the numerical simulations in three-dimensional space. The study presented uses eight different baffle configurations with two solid obstructions, which have area blockage ratios of 0.24 and 0.5. The flame speed, maximum rate of pressure-rise as well as peak overpressure magnitude and timing are presented and discussed. Numerical results are validated against available published experimental data. It is concluded that, increasing the solid obstacle area blockage ratio and the number of consecutive baffles results in a raised maximum rate of pressure rise, higher peak explosion overpressure and faster flame propagation. Future model development would require more experimental data, probably in a more congested configuration.     

Quantification of source-level turbulence during LNG spills onto a water pond

The use of computational fluid dynamics (CFD) models to simulate LNG vapor dispersion scenarios has been growing steadily over the last few years, with applications to LNG spills on land as well as on water. Before a CFD model may be used to predict the vapor dispersion hazard distances for a hypothetical LNG spill scenario, it is necessary for the model to be validated with respect to relevant experimental data. As part of a joint-industry project aimed at validating the CFD methodology, the LNG vapor source term, including the turbulence level associated with the evaporation process vapors was quantified for one of the Falcon tests.This paper presents the method that was used to quantify the turbulent intensity of evaporating LNG, by analyzing the video images of one of the Falcon tests, which involved LNG spills onto a water pond. The measured rate of LNG pool growth and spreading and the quantified turbulence intensity that were obtained from the image analysis were used as the LNG vapor source term in the CFD model to simulate the Falcon-1 LNG spill test. Several CFD simulations were performed, using a vaporization flux of 0.127 kg/m² s, radial and outward spreading velocities of 1.53 and 0.55 m/s respectively, and a range of turbulence kinetic energy values between 2.9 and 28.8 m²/s². The resulting growth and spread of the vapor cloud within the impounded area and outside of it were found to match the observed behavior and the experimental measured data.The results of the analysis presented in this paper demonstrate that a detailed and accurate definition of the LNG vapor source term is critical in order for any vapor cloud dispersion simulation to provide useful and reliable results.