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.     

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.     

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.     

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.     

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.     

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

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

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.     

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.     

Evaluation of CFD simulations of transient pool fire burning rates

Fire is the most commonly occurring major accident hazard in the chemical and process industries, with industry accident statistics highlighting the liquid pool fire as the most frequent fire event. Modelling of such phenomena feeds heavily into industry risk assessment and consequence analyses. Traditional simple empirical equations cannot account for the full range of factors influencing pool fire behaviour or increasingly complex plant design. The use of Computational Fluid Dynamics (CFD) modelling enables a greater understanding of pool fire behaviour to be gained numerically and provides the capability to deal with complex scenarios.This paper presents an evaluation of the Fire Dynamics Simulator (FDS) for predictive modelling of liquid pool fire burning rates. Specifically, the work examines the ability of the model to predict temporal variations in the burning rate of open atmosphere pool fires. Fires ranging from 0.4 to 4 m in diameter, involving ethanol and a range of liquid hydrocarbons as fuels, are considered and comparisons of predicted fuel mass loss rates are compared to experimental measurements.The results show that the liquid pyrolysis sub-model in FDS gives consistent model performance for fully predictive modelling of liquid pool fire burning rates, particularly during quasi-steady burning. However, the model falls short of predicting the subtleties associated with each phase of the transient burning process, failing to reliably predict fuel mass loss rates during fire growth and extinction. The results suggest a range of model modifications which could lead to improved prediction of the transient fire growth and extinction phases of burning for liquid pool fires, specifically, investigation of: ignition modelling techniques for high boiling temperature liquid fuels; a combustion regime combining both infinite and finite-rate chemistry; a solution method which accounts for two- or three-dimensional heat conduction effects in the liquid-phase; alternative surrogate fuel compositions for multi-component hydrocarbon fuels; and modification of the solution procedure used at the liquid-gas interface during fire extinction.     

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.     

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.     

高温热表面油液蒸发的时变性热质传递模型与实验研究

针对高温热表面油液蒸发热质传递过程的时变性,考虑这一过程中的对流传质传热,建立了热环境作用下油液蒸发的热质传递模型方程,通过无量纲变换,求得空间浓度分布和温度场随时间的变化规律。以庚烷为试验对象,对高温热表面油液蒸发过程进行了实验研究。理论分析与实验表明:庚烷蒸发过程中,刘易斯数大于1,传热速率大于传质速率;蒸发导致的质量损失与时间平方根的成正比,与液面的面积成正比,且与质量扩散系数的平方根成正比,饱和蒸气浓度越大,蒸发速率也越大。油液蒸发计算结果与试验结果基本一致,表明了模型的有效性。     

喷雾碰壁液膜蒸发三维数值模拟

针对喷雾过程喷注碰壁形成的液膜蒸发现象进行研究,通过理论分析和实验为数据建立了相应的数值计算模型。在典型浴盆型燃烧室内,采用一个四孔喷嘴进行高压喷雾,给定初始涡流比和燃烧室壁四温度,对包含液膜蒸发的喷雾过程进行详细的三维数值模拟,获得了液滴轨迹图以及空气运动速度的空间分布和时间进程,为研究空气运动和壁温等参数对喷雾过程和液膜蒸发过程的影响以及它们之间的相互耦合作用打下了坚实的基础。     

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.     

氨分解装置中液氨储罐火灾爆炸危险性分析

以某金属处理企业氨分解装置中液氨储罐罐区为例,对液氨泄漏后火灾爆炸事故及其伤害范围进行了研究,用池火、蒸气云爆炸和沸腾液体扩展蒸气爆炸模型进行计算分析,给出火灾、爆炸事故的人员伤害和财产损失范围。结果表明:围堤堤内池火或罐内池火时,罐区建构筑物内的汽化器、管道等设备会因直接过火或热辐射导致损坏,建筑内人员死亡,但难以波及罐区之外;蒸气云爆炸产生相当于1192.72kgTNT爆炸的当量,爆炸的后果严重,应重点防范,防范的重点为液氨泄漏、点火源;沸腾液体扩展蒸气爆炸的火球半径56.1m,持续时间8.7s,死亡半径27.2m,其源于储罐受热或系统突然失效,液体瞬时泄漏汽化并遇点火源而发生,具有突发性且后果严重,企业应高度重视并严格储罐及系统的定期检验与校验、密切关注系统的有效运行。     

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.     

烃类池火灾热辐射量化分析模型探讨

为澄清有关烃类池火灾热辐射量化分析模型选择中存在的问题,针对"池火计算方法"模型进行量纲分析和物理意义方面的讨论,提出该模型在概念使用、量纲关系、热释放速率的计算以及池火焰模化等方面存在的问题。系统阐述点源模型、Shokri-Beyler模型、Mudan模型等常用烃类池火灾热辐射通量分析模型的适用条件和应用范围。点源模型适用于被辐射目标物从池火焰接受的热辐射通量小于5 kW.m-2情况下池火灾热辐射的量化分析;Shokri-Beyler模型主要应用于估算被辐射目标物从池火焰接受的热辐射通量大于5 kW.m-2的情况;Mudan模型可用于估算无风或有风条件下被辐射目标物从池火焰接受的热辐射通量。     

The modelling of pool vaporization

This paper describes a mathematical model which calculates the time dependencies of the flow rate and composition of the vapour emerging from a pool. A large variety of accidental cases can be covered: continuous or instantaneous spills, on confined or unconfined ground, ideal or non-ideal liquid mixtures in boiling or evaporating conditions. The boiling, when present, is modelled through an equation system comprising the Rachford-Rice relation and the energy balance of the pool, which is assumed to be well-mixed. In the case of a volatile pool, interfacial mass rates are determined taking into account the Stefan flux, and the thermal resistance inside the liquid phase is also considered. In all situations, the energy balance includes the contribution of ground, sun and air. Known experimental data have been used to validate the model.     

太阳辐射对纯液体蒸发行为影响的实验研究

通过蒸发实验,考察了太阳辐射对液体苯、甲苯和乙醇蒸发行为的影响.实验结果表明,苯和甲苯表面温度随时间呈现出先降低后升高的趋势,而乙醇则呈现出先降低后恒定的趋势,且环境风速越大,这种变化趋势就越明显.液体表面温度的变化不但与太阳辐射强度以及环境风速的大小有关,还与液体本身的物理性质(如热容、蒸发潜热和饱和蒸气压)有关.因此,在预测液体的蒸发速率时,必须考虑太阳辐射或周围热辐射源对液体表面温度的影响.