Investigation of pool spreading and vaporization behavior in medium-scale LNG tests

A failure of a Liquefied Natural Gas (LNG) tanker can occur due to collision or rupture in loading/unloading lines resulting in spillage of LNG on water. Upon release, a spreading liquid can form a pool with rapid vaporization leading to the formation of a flammable vapor cloud. Safety analysis for the protection of public and property involves the determination of consequences of such accidental releases. To address this complex pool spreading and vaporization phenomenon of LNG, an investigation is performed based on the experimental tests that were conducted by the Mary Kay O'Connor Process Safety Center (MKOPSC) in 2007. The 2007 tests are a part of medium-scale experiments carried out at the Brayton Fire Training Field (BFTF), College Station. The dataset represents a semi-continuous spill on water, where LNG is released on a confined area of water for a specified duration of time. The pool spreading and vaporization behavior are validated using empirical models, which involved determination of pool spreading parameters and vaporization rates with respect to time. Knowledge of the pool diameter, pool height and spreading rate are found to be important in calculating the vaporization rates of the liquid pool. The paper also presents a method to determine the vaporization mass flux of LNG using water temperature data that is recorded in the experiment. The vaporization rates are observed to be high initially and tend to decrease once the pool stopped spreading. The results of the analysis indicated that a vaporization mass flux that is varying with time is required for accurate determination of the vaporization rate. Based on the data analysis, sources of uncertainties in the experimental data were identified to arise from ice formation and vapor blocking.     

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.     

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.     

Laboratory scale analysis of the influence of different heat transfer mechanisms on liquid nitrogen vaporization rate

The prediction of the potential hazards associated to accidental liquefied natural gas (LNG) spills has motivated a number of different studies including experimental and numerical approaches. Most of these studies focus on dispersion predictions, however there is limited information regarding source term of it: liquid spill and vaporization. There is a need of further improvements on the understanding of these phenomena and the quantification of the most important parameters that can affect them.The vaporization of cryogenic liquids is governed by the heat transfer phenomena including conduction, convection and thermal radiation mechanisms. The present work investigates the contribution of each of these heat transfer modes to the vaporization rate of cryogenic liquid nitrogen (LN₂) contained in a Dewar flask using well controlled and instrumented laboratory scale experiments. LN₂ vaporization rate was measured with individually controllable contributions from convective (generated by an electric fan) and thermal radiative (generated by light bulb) heat transfer in the presence of a baseline conductive heat transfer rate.In both cases of convection and radiation analysis the experimental study showed that they can play a significant role in the vaporization rate of LN₂. It was observed that the radiative heat absorbed by the LN₂ during the vaporization experiment represents only 50%–65% of the incident radiation that would reach the LN₂ pool surface if no vapour was present. Convective heat transfer generated by the fan was shown to have had the most significant contribution to the total heat transfer. As expected, this contribution was significantly higher than the one from bulb radiation. The experimental data also showed that the liquid level in the Dewar play a key role in the resulting amount of convective heat transfer. This could be attributed to the fact that lower liquid level the side walls of the Dewar were high enough to hold a layer of vapour and limit air motion directly above the liquid surface, thus limiting the heat transfer by convection.     

Development of heat transfer and evaporation model of LNG covered by Hi-Ex foam

Natural gas is a kind of clean, efficient green energy source, which is used widely. Liquefied natural gas (LNG) is produced by cooling natural gas to −161 °C, at which it becomes the liquid. Once LNG was released, fire or explosion would happen when ignition source existed nearby. The high expansion foam (Hi-Ex foam) is believed to quickly blanket on the top of LNG spillage pool and warm the LNG vapor to lower the vapor cloud density at the ground level and raising vapor buoyancy. To identify the physical structure after it contacted with LN₂ and to develop heat transfer model, the small-scale field test with liquid nitrogen (LN₂) was designed. In experiment, three layers including frozen ice layer, frozen Hi-Ex layer and soft layer of Hi-Ex foam were observed at the steady state. By characterizing physical structure of the foam, formulas for calculating the surface of single foam bubble and counting foam film thickness were deduced. The micro heat transfer and evaporation model between cryogenic liquid and Hi-Ex foam was established. Indicating the physical structure of the frozen ice layer, there were a certain number of icicles below it. The heat transfer and evaporation mathematical model between the frozen ice layer and LNG was derived. Combining models above with the heat transfer between LNG, ground and cofferdam, the heat transfer and evaporation mathematical model of LNG covered by Hi-Ex foam was developed eventually. Finally, LN₂ evaporation rate calculated by this model was compared with the measured evaporation rate. The calculated results are 1.2–2.1 times of experimental results, which were acceptable in engineering and proved the model was reliable.     

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.     

Experimental study and prediction of liquid nitrogen pool evaporation process on concrete surface

The vaporization rate of pool boiling process of a liquid nitrogen spilled on concrete surface was investigated by a visual experiment platform. The boiling curve for liquid nitrogen on concrete cooling process was obtained. The shapes of bubbles in three typical boiling regimes were observed. Based on the experimental results, the coefficients of the empirical formulas for nuclear boiling and film boiling are modified, and the empirical formulas for boiling of liquid nitrogen-concrete surfaces are obtained. Combined with the calculation formula of the non-steady-state semi-infinite one-dimensional heat conduction temperature, a coupled calculation model for the heat transfer and vaporization process of liquid pool on the liquid nitrogen-concrete surface is proposed. Application of this model can better predict the quality of liquid nitrogen vaporization.     

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.     

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

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

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.     

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.     

Experimental research of LNG accidental underwater release and combustion behavior

Evaluating potential hazards caused by accidental LNG release from underwater pipelines or vessels is a significant consideration in marine transportation safety. The aim of this study was to capture the dynamic behavior of LNG jet released under water and to analyze its vapor dispersion characteristics and combustion characteristics on the water surface during different release scenarios. Controlled experiments were conducted where LNG was jet released from a cryogenic storage tank. The dynamic process of LNG being jet released from orifices of different sizes and shapes, as well as the rising plume structure, were captured by a high-speed camera. The leakage flow rate and pipeline pressure were recorded by a flow meter and pressure gauge, respectively. The concentration distribution that emanated from the water surface was measured utilizing methane sensors in different positions with various wind speeds. The flame combustion characteristics of LNG vapor clouds, which immediately ignited upon the enclosed water tank, were also recorded. Additionally, the mass burning rate of the flame on the water surface was evaluated, and a new correlation between the ratio of flame length and width was established. The results indicated a large dimensionless heat release rate (Q*) and a continuous release flow rate in a limited burning area. This study could provide greater understanding of the mechanisms of LNG release and combustion behavior under water.     

On the effectiveness of mitigation strategies for cryogenic applications

The need for sustainable energy sources, as well as the current energetic crisis involving the majority of markets, has promoted the use of cryogenic liquefaction for the transportation and storage of natural gas (i.e., LNG). To guarantee the development of a robust and safe infrastructure, a complete understanding of the main phenomena occurring at low temperatures is paramount. In this sense, the largest grey areas are the characterization of the combustion at low-initial temperature and the interactions between water and cryogenic liquid. For these reasons, this work presents an experimental campaign on the possible mitigation strategies for the mitigation of consequences related to the accidental release of LNG. Particular emphasis was posed on the direct and indirect effects of water on cryogenic pool fire. The former resulted in a significant increase in the dimensions of fire (∼+50%) and burning rate (∼300%) with respect to the case with no direct contact between water and LNG, whereas the latter generated an abrupt decrease in the measured temperatures (<100 °C). The use of an emergency flare to empty an LNG tank was tested, as well. The spatial distribution of temperature was monitored along with the time to guarantee the safe operability of this equipment in the case of LNG combustion. The explanations for the observed phenomena and trends were provided, allowing for the development of safe procedures for the emergency response related to cryogenic fuels.     

冷冻堵漏法在危险化学品泄漏事故中的应用

为了研究液化气体泄漏冷冻堵漏的堵漏机制,运用流体力学、传热学等知识对液化石油气（LPG）储罐（槽罐）泄漏时泄漏口处产生局部低温的现象进行了研究,探讨了LPG液相泄漏和气相泄漏2种不同泄漏形式的低温效应。结果表明:液相泄漏时,泄漏口处温度下降程度与泄漏口面积成正比,且随着罐体内部压力的减小而减弱,推导出喷水冷冻堵漏的成冰时间公式;气相泄漏时,对罐内压力与温度的平衡关系进行模拟并建立了数学模型;发现由于LPG气、液相之间对流换热和汽化吸热效应的差异,导致液相与气相之间的温度差,此温度差是罐体外壁产生结霜分层现象的主要原因。     

Experimental study of the evaporation of spreading liquid nitrogen

The investigation of cryogenic liquid pool spreading is an essential procedure to assess the hazard of cryogenic liquid usage. There is a wide range of models used to describe the spreading of a cryogenic liquid pool. Many of these models require the evaporation velocity, which has to be determined experimentally because the heat transfer process between the liquid pool and the surroundings is too complicated to be modeled. In this experimental study, to measure the evaporation velocity when the pool is spreading, liquid nitrogen was continuously released onto unconfined concrete ground. Almost all of the reported results are based on a non-spreading pool in which cryogenic liquid is instantaneously poured onto bounded ground for a very short period of time. For the precise measurement of pool spreading and evaporation weight with time, a cone-type funnel was designed to achieve a nearly constant liquid nitrogen release rate during discharge. Specifically, three nozzles with nominal flow rates of 3.4 × 10⁻² kg/s, 5.6 × 10⁻² kg/s and 9.0 × 10⁻² kg/s were used to investigate the effect of the release rate on the evaporation velocity. It is noted that information about the release rate is not necessary to measure the evaporation velocity in case of the non-spreading pool. A simultaneous measurement of the pool location using thermocouples and of the pool mass using a digital balance was carried out to measure the evaporation velocity and the pool radius. A greater release flow rate was found to result in a greater average evaporation velocity, and the evaporation velocity decreased with the spreading time and the pool radius.     

低温液体在水中快速相变过程流动与传热数值模拟

用试验和理论分析的方法对快速相变爆炸强度的预测缺乏定量模型,因此建立了一种欧拉-欧拉双流体多相流模型与传热模型相互耦合的数值模型,并通过与Clarke H将液氮喷射入水的快速相变试验数据对比来验证模型的可靠性和正确性。通过数值计算得出快速相变过程中流场、压力场、温度场随时间变化的情况,探讨了快速相变的传热机理。结果表明:快速相变是强制对流、膜态沸腾、爆发沸腾和核态沸腾之间的转换过程;相间换热系数随时间的无量纲变化关系可以用3个高斯分布的叠加来描述。     

吨量煤体的自燃过程实验模拟研究

为更好地弄清矿井实际的煤自然发火规律,利用装煤量达5吨的大型实验台对两种烟煤分别进行了自燃模拟实验。大煤量的实验能够很好地模拟煤矿中煤低温氧化和传热传质共同作用导致的发火过程,实验得到的自然发火期与煤矿实际发火期也是一致的。实验中煤样从缓慢氧化变为快速氧化的临界温度为100～110℃。当煤温低于,临界温度时,煤样的升温受到空气流带走热量和向外界散热的影响很大,因此夹层水的保温作用就很关键。当煤温超过临界温度后,反应加快,温度急剧上升,散热的影响明显降低,反应主要受限于氧气的供应情况。     

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.     

Modeling of pool spreading of LNG on land

This paper presents a source term model for estimating the rate of spreading of LNG and other cryogenic mixtures on unconfined land. The model takes into account the composition changes of a boiling mixture, the varying thermodynamic properties due to preferential boiling within the mixture and the effect of the various boiling regimes on conductive heat transfer. A sensitivity analysis is conducted to determine the relative effect of each of these phenomena on pool spread. The model is applied to continuous and instantaneous spills. The model is compared to literature experimental data on cryogenic pool spreading.     

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.