Small-scale experimental study of vaporization flux of liquid nitrogen released on ice

One of the LNG accident scenarios is the collision of an LNG carrier on an iceberg during marine transportation. A collision can result in damages to the vessel and lead to the leakage of the contents on ice or an ice-water mixture. When cryogenic liquid comes in contact with ice, it undergoes rapid vaporization due to the difference in temperature between the ice and cryogenic liquid. This process is different from the heat transfer between water and cryogenic liquid as ice is a solid and thus heat transfer to the pool occurs primarily through conduction. In this paper, the heat transfer phenomenon between ice and cryogenic liquid was studied through a small-scale experiment and the resulting vaporization mass fluxes were reported. The experiment involved six spills with varying amount of liquid nitrogen on different ice temperature to determine its effect on vaporization mass flux. The vaporization mass fluxes were determined by direct measurement of the mass loss during the experiment. The results indicated that the vaporization mass flux was a function of release rate and ice temperature. When the release rate and ice temperature was high, the vaporization mass flux follows a decreasing trend. With further reduction in release rate and ice temperature, the vaporization mass flux was found to be independent with time. The one dimensional conduction model was validated against experimental results. The predicted temperatures and heat flux were found to be in good agreement with the experimental data.     

Modeling the vapor source term associated with the spill of LNG into a sump or impoundment area

The recent publication of evaluation protocols for vapor source term models and vapor dispersion models have influenced the modeling approaches that can be used for approval of new and expansion projects at LNG receiving terminals. In the past few years the scientific basis of integral vapor source term models has been questioned with growing concerns regarding their validity. In this paper, the shallow water equations (SWEs) were solved to study the characteristics of the evaporating LNG pool associated with a constant flow rate spill of LNG into a concrete sump. In the early stages of pool spreading, the leading edge thickness profile of the SWE model scales with the square root of the distance from the leading edge as the pool spreads. After the edge of the pool reaches the wall, the reflected wave forms a hydraulic jump that travels back towards the center of the pool at a speed that is considerably slower than the initial spreading of the pool. Once the hydraulic jump reaches the center, the pool assumes a nearly flat free surface for the rest of the spill. The pool spreading and the rate of evaporation from the SWEs were then compared to the solution provided by the integral model, PHAST. The two approaches were found to agree well with one another. The SWE model was also used to demonstrate the influence of an elevated spill source. With an elevated source, the LNG pool spreads faster, significantly increasing the initial rate of vaporization and peak vaporization rate. This increase in the initial rate of vaporization could lead to an increase in the vapor cloud hazard distance. The SWE model was also used to demonstrate the influence of an inclined sump floor in the shape of an inverted cone where the spilling LNG accumulates in the low vertex of the cone. Inclined sump floors can be used to significantly reduce the cumulative evaporation, making them attractive as a possible mitigation approach in cases where a containment sump is located close to a property boundary.     

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.     

水面LNG液池扩展模型的分析与对比研究

为了评价在开阔水面上的液化天然气（LNG）火灾和蒸气云爆炸灾害后果,分析了LNG水面扩展动态过程;对比分析了Fay模型、FERC模型和计算流体力学软件FLACS的计算结果,探讨了LNG液池面积随时间的动态变化过程,分析了泄漏量、泄漏速率等参数对LNG液池扩展半径的影响;根据液池扩展模型的计算结果,确定了LNG液池的最大面积,并以此分析了LNG流淌火灾的辐射危害。研究结果表明:对于相同的泄漏条件,3种方法模拟的泄漏LNG水面扩展动态过程相似,一般情况下,FLACS模型,FERC模型和Fay模型所计算的最大液池半径依次增大;由于FERC模型与FLACS软件的模拟结果接近且偏于保守,故此在一般的工程应用时,采用FERC模型即可方便快捷地获得较为准确的结果。     

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.     

Highlights of LNG risk technology

The authors have recently undertaken a major review of LNG consequence modeling, compiling a wide range of historical information with more recent experiments and modeling approaches in a book entitled “LNG Risk-Based Safety: Modeling and Consequence Analysis”. All the main consequence routes were reviewed – discharge, evaporation, pool and jet fire, vapor cloud explosions, rollover, and Rapid Phase Transitions (RPT’s). In the book, experimental data bases are assembled for tests on pool spread and evaporation, burn rates, dispersion, fire and radiation and effects on personnel and structures. The current paper presents selected highlights of interest: lessons learned from historical development and experience, comparison of predictions by various models, varying mechanisms for LNG spread of water, a modeling protocol to enable acceptance of newer models, and unresolved technical issues such as cascading failures, fire engulfment of a carrier, the circumstances for a possible LNG BLEVE, and accelerated evaporation by LNG penetration into water.     

Underwater LNG release: Does a pool form on the water surface? What are the characteristics of the vapor released?

One of the scenarios of concern in assessing the safety issues related to transportation of LNG in a marine environment (ship or underwater pipeline) is the release of LNG underwater. This scenario has not been given the same level of scientific attention in the literature compared to surface releases and assessment of consequences therefrom. This paper addresses questions like, (1) does an LNG spill underwater form a pool on the water surface and subsequently evaporate like an LNG spill “on the surface” producing cold, heavier than air vapors?, and (2) what is the range of expected temperatures of the vapor, generated by LNG release due to heat transfer within the water column, when it emanates from the water surface?Very limited data from two field tests of LNG underwater release are reviewed. Also presented are the results from tests conducted in other related industries (metal casting, nuclear fission and fusion, chemical processing, and alternative fuel vehicles) where a hot (or cold) liquid is injected into a bulk cold (or hot) liquid at different depths.A mathematical model is described which calculates the temperature of vapor emanating at the water surface, and the liquid fraction of released LNG that surfaces, if any, to form a pool on the water surface. The model includes such variables as the LNG release rate, diameter of the jet at release, depth of release and water body temperature. Results obtained from the model for postulated release conditions are presented. Comparison of predicted results with available LNG underwater release test data is also provided.     

Evaporating liquid flow in a channel (An integral model based on shallow water flow approximation)

Liquefied Natural Gas (LNG) storage facilities generally include channels to convey potential spills of the liquid to an impoundment. There is increasing concern that dispersion of vapors generated by flow of LNG in a channel may lead to higher than limit vapor concentrations for safety at site boundary from channels that may be close to the dike walls. This issue is of recent concern to regulatory agencies, because the calculation of vapor hazard distance(s) from LNG flow in a channel is not required under existing LNG facility siting standards or regulations.An important parameter that directly affects the calculated LNG vapor dispersion distance is the source strength (i.e., the rate of vaporization of LNG flow from the wetted channel surfaces, as a function of spatial position and time). In this paper a model is presented which considers the variation of the depth of the flowing LNG with spatial location and time, and calculates the spatial and temporal dependence of the mass rate of vapor generation. Self similar profiles for the spatial variation of the thermal boundary layer in the liquid wetted wall and liquid depth variation are assumed. The variation with time of the location of the liquid spread front and the evaporation rate are calculated for the case of a constant LNG spill rate into a rectangular channel. The effects of two different channel slopes are evaluated. Details of the results and their impact on dispersion distances are discussed.     

LNG加气站槽车卸车直供过程泄漏后果分析

为研究LNG加气站槽车直接供液过程泄漏后果严重程度,采用HAZOP辨识槽车供液和储罐供液典型泄漏场景,基于PHAST分析不同泄漏场景下LNG液池半径、蒸汽云扩散距离及积聚时长、爆炸超压和池火热辐射影响范围,定量评价槽车供液可能造成的事故后果扩大程度.结果表明:槽车供液泄漏事故的LNG液池最大半径、蒸汽云最大扩散距离、爆...     

Underwater LNG release test findings: Experimental data and model results

An underwater LNG release test was conducted to understand the phenomena that occur when LNG is released underwater and to determine the characteristic of the vapor emanating from the water surface. Another objective of the test was to determine if an LNG liquid pool formed on the water surface, spread and evaporated in a manner similar to that from an on-the-surface release of LNG.A pit of dimensions 10.06 m × 6.4 m and 1.22 m depth filled with water to 1.14 m depth was used. A vertically upward shooting LNG jet was released from a pipe of 2.54 cm diameter at a depth of 0.71 m below the water surface. LNG was released over 5.5-min duration, with a flow rate of 0.675 ± 0.223 L/s. The wind speed varied between 2 m/s and 4 m/s during the test.Data were collected as a function of time at a number of locations. These data included LNG flow rate, meteorological conditions, temperatures at a number of locations within the water column, and vapor temperatures and concentrations in air at different downwind locations and heights. Concentration measurements were made with instruments on poles located at 3.05 m, 6.1 m and 9.14 m from the downwind edge of the pit and at heights 0.46 m, 1.22 m, and 2.13 m. The phenomena occurring underwater were recorded with an underwater video camera. Water surface and in-air phenomena including the dispersion of the vapor emanating from the water surface were captured on three land-based video cameras.The lowest temperature recorded for the vapor emanating from the water surface was −1 °C indicating that the vapor emitted into air was buoyant. In general the maximum concentration observed at each instrument pole was progressively at higher and higher elevations as one traveled downwind, indicating that the vapor cloud was rising. These findings from the instrument recorded data were supported by the visual record showing the “white” cloud rising, more or less vertically, in air. No LNG pool was observed on the surface of water. Discussions are provided on the test findings and comparison with predictions from a previously published theoretical model.     

Medium scale LNG-related experiments and CFD simulation of water curtain

This work is a continuation of the experimental research on LNG releases and their consequence mitigation methods, which has been carried on by the Mary Kay O’Connor Process Safety Center (MKOPSC) at the Texas A&M University since 2004.A series of medium scale experiments to test the ability of the water curtains to hold up and disperse a vapor has been performed. Colored smoke has been used as an analog of the LNG vapor for easier tracking of the vapor path through the water curtain. The results and some analysis of the experimental data are presented. The CFD software FLACS (GexCon AS) was used to simulate the effects of the water curtains on vapor dispersion. The results of the simulations were compared with experimental data.     

BLEVE: The case of water and a historical survey

After a short update of the current more accepted definition of BLEVE, the special features of water BLEVEs are analyzed. The stronger overpressure wave generated in the case of water as compared to that of other substances is justified in terms of volume change. Through a comparison with liquefied pressurized propane, three possibilities are analyzed: the simultaneous contribution of both the liquid and the preexisting vapor, the contribution of the liquid flash vaporization, and the contribution of the pre-existing vapor. Also a historical survey on a set of 202 BLEVE accidents –the largest sample of BLEVE accidents surveyed until now– is presented. LPG was the most common substances in this set of accidents. However, water and LNG (11% of water and 4% of LNG in the studied cases) have also been involved. Impact failure (44.8%) and human factor (30.3%) were the most common causes of BLEVEs. Transport, storage, process plants, and transfer were the activities in which more accidents occurred.     

Quantification of turbulence in cryogenic liquid using high speed flow visualization

A high speed flow visualization experiment was conducted to characterize the boiling induced turbulence when a cryogenic liquid is released on water. The advective transport of turbulent structures traversing through the liquid was captured and reconstructed using image processing software to obtain information on velocity components. The numerical results obtained from image processing were used to determine turbulence parameters like turbulent intensity, turbulent kinetic energy and eddy dissipation rate. An interesting aspect of the study was the formation of wavy structures called ‘thermals’ which were characteristic of turbulent convection. The thermals were found to act as a catalyst in increasing heat transfer and turbulence between water and cryogenic pool. The turbulent intensity was influenced by the turbulent velocity and had direct effects on the vaporization flux. Among the turbulence parameters, increase in turbulent kinetic energy resulted in faster vaporization of cryogenic liquid through enhanced mixing, whereas variations in the eddy dissipation rate had weak dependence on vaporization. Additionally, the initial height of cryogenic liquid was also found to strongly affect the vaporization mass flux.     

Developing a CFD heat transfer model for applying high expansion foam in an LNG spill

Liquefied natural gas (LNG) is widely used to cost-effectively store and transport natural gas. However, a spill of LNG can create a vapor cloud, which can potentially cause fire and explosion. High expansion (HEX) foam is recommended by the NFPA 11 to mitigate the vapor hazard and control LNG pool fire. In this study, the parameters that affect HEX foam performance were examined using lab-scale testing of foam temperature profile and computational fluid dynamics (CFD) modeling of heat transfer in vapor channels. A heat transfer model using ANSYS Fluent® was developed to estimate the minimum HEX foam height that allows the vapors from LNG spillage to disperse rapidly. We also performed a sensitivity analysis on the effect of the vaporization rate, the diameter of the vapor channel, and the heat transfer coefficient on the required minimum height of the HEX foam. It can be observed that at least 1.2 m of HEX foam in height are needed to achieve risk mitigation in a typical situation. The simulation results can be used not only for understanding the heat transfer mechanisms when applying HEX foam but also for suggesting to the LNG facility operator how much HEX foam they need for effective risk mitigation under different conditions.     

Mutual effects between dynamic leakage behavior and the pressure/temperature in a LNG tank with external heat fluxes

The objective of this work is to investigate and model the mutual effects between the dynamic pressure/temperature in the LNG tank and the leakage behavior with external heat fluxes. The results suggest that the pressure and temperature in tank during leakage change with the comparison results between the heat flux consumed in liquid boil-off and the external heat flux supplied. At the liquid leakage stage, when the external heat flux is not very high, the pressure in tank tends to increase significantly, even results in tank explosion. It increases more and more heavily with higher and higher external heat fluxes. At the vapor leakage stage, large amount of vapor spray out, which results in a high generation rate of vapor by the liquid boil-off. The pressure in tank normally decreases to be low, which is unfavorable for the LNG tank explosion. Therefore, at this vapor leakage stage, blocking the leakage hole as soon as possible is not always a right choice for fire fighters. Finally, it is suggested that reducing the heat flux into the tank, either at the liquid leakage stage or in vapor one, is key to the tank safety.     

Application of fire suppression materials on suppression of LNG pool fires

An LNG pool fire is considered one of the main hazards of LNG, together with LNG vapor dispersion. Suppression methods are designed to reduce the hazard exclusion zones, distance to reach radiant heat of 5 kW/m², when an LNG pool fire is considered. For LNG vapor dispersion, the hazard exclusion zone is the distance travelled by the LNG vapor to reach a concentration of 2.5% v/v (half of the LNG lower flammability limit).Warming the LNG vapor to reach positive buoyancy faster is one way to suppress LNG vapor dispersion and reduce evaporation rate (thus fire size and its associated radiant heat) and that is the main objective in LNG pool fire suppression. Based on previous research, the use of high expansion foam has been regarded as the primary method in suppressing LNG pool fires. However, in 1980, another method was introduced as an alternative pool fire suppression system, Foamglas^®. The research concluded that 90% of the radiant heat was successfully reduced. Currently-called Foamglas^® pool fire suppression (Foamglas^® PFS) is a passive mitigation system and is deployed after the leak occurs. Foamglas^® PFS is non-flammable, and has a density one-third of the density of LNG, thus floats when an LNG pool is formed.This paper describes the study and confirmation of Foamglas^®PFS effectiveness in suppressing LNG pool fires. In addition, while Foamglas^® PFS is not expected to suppress LNG vapor dispersion, further investigation was conducted to study the effect of Foamglas^®PFS on LNG vapor dispersion. An LNG field experiment was conducted at Brayton Fire Field. The experimental development, procedures, results and findings are detailed in this paper.     

Modeling of a cryogenic liquid pool boiling by CFD simulation

A boiling model is developed by Computational Fluid Dynamics (CFD) code to calculate the source term of a cryogenic liquid spill. The model includes the effect of the changing ground temperature on the vaporization rate of the cryogenic liquid. Simulations are performed for liquid nitrogen. The model can describe different boiling regimes (film, transition and nucleate). The heat flux calculated for each boiling regimes are compared to the experimental data from literature. The developed numerical model seems to have a good ability to predict the heat flux for the film boiling stage. Model development is still necessary to improve the prediction of the nucleate boiling regime. Overall, the approach shows very promising results to model the complex physical phenomena involved in in the vaporization of cryogenic liquid pool spilled on ground.     

A mechanistic model for hydrocarbon plumes rising through water

For releases of hydrocarbons from a subsea pipeline, riser, or production facility, the shape of the plume rising through the water must be predicted prior to any assessment of gas dispersion, liquid pools, or fire above the water surface. The location and size of the plume at the water surface are key parameters for subsequent consequence modeling. A mechanistic model has been developed to predict the plume trajectory and size, based on mass and momentum balances and an empirical water entrainment ratio from the literature. With suitable physical property values available, the model is applicable to releases of gas and/or liquid hydrocarbons, predicting the vaporization and vapor expansion due to decreasing hydrostatic pressure as the plume rises through the water. Some validation of the model was obtained with 16 tests in a small-scale transparent tank. The data cover a wide range of flow rates, including both choked and unchoked flow. The predicted and measured trajectories (centerline displacement) agreed reasonably well. Predictions of the model are presented for three fluids. The model is valuable for assessing the consequences of underwater hydrocarbon releases, providing input for subsequent modeling of gas dispersion or liquid pools and pool fires.     

高倍泡沫加速泄漏LNG扩散的有效性模拟试验

泡沫灭火剂在扑灭液体火灾中起到重要作用,关于低温液体蒸气云扩散控制的研究也逐渐得到应用。通过小尺寸模拟试验验证高倍泡沫加速泄漏LNG扩散的有效性,设计并进行了低温液体自然蒸发和高倍泡沫覆盖低温液体两个对照试验,测量了竖直方向上10个高度处的温度及装置整体质量,从而获取了低温液体蒸气到达泡沫层顶端时温度及蒸发速率的变化情况。结果表明,与未添加泡沫的情况对比,高倍泡沫的覆盖使泄漏低温液体在1 800 s内的蒸发量减少了6.4%,如果时间更长则减少的比例更多,且蒸发出的低温液体穿过泡沫层后蒸气温度可达0℃左右,而未添加泡沫时同等高度处蒸气温度为-75℃左右。0℃时,LNG蒸气密度已明显小于空气密度,此温度下LNG蒸气会迅速向上扩散,而不至在地表积聚,由此证明高倍泡沫能够加速泄漏低温液体蒸气向上扩散,减小了低温液体蒸气在地面积聚并引发火灾爆炸事故的可能性,从而证实了高倍泡沫加速泄漏LNG扩散的有效性。     

Validation of liquid nitrogen vaporisation rate by small scale experiments and analysis of the conductive heat flux from the concrete

The vaporisation of a liquid nitrogen pool spilled on concrete ground was investigated in small scale field experiments. The pool vaporisation rate and the heat transfer from the concrete ground were measured using a balance and a set of embedded heat flux sensors and thermocouples. The ability to predict the concrete's thermal properties based on these measurements was investigated. This work showed that a simple, one-dimensional theoretical model, assuming heat conduction through a semi-infinite ground with ideal contact between the cryogenic liquid and the ground, commonly used to describe the heat transfer from a ground to the LNG, can be used to match the observed vaporisation rate. Though estimated parameters, thermal conductivity and thermal diffusivity, do not necessary represent real values. Although the observed vaporization rate follows a linear trend, and thus can be well represented by the model, the overall model prediction seems to be overestimated. The temperature profile inside the concrete is slightly over-predicted at the beginning and under-predicted at later stage of the spill. This might be an effect of the dependence of the concrete's thermal properties on the temperature or may indicate an incorrect modelling and a varying temperature of the ground surface.