改进的蒸气云爆炸模型在LNG储罐爆炸模拟中的应用

针对TNT当量法在LNG储罐蒸气云爆炸模拟中的应用进行了改进,考虑并分析了使用传统TNT模型时所忽略的LNG液池蒸发过程,通过建立LNG与地面的传热模型得出了LNG液池蒸发速率随时间变化的关系,液池的蒸发速率在最初随时间的增长较快,在增至最大值后与时间的平方根成反比逐渐减小。以3万m~3 LNG储罐连续泄漏20 min为例,根据蒸发速率与时间的关系算出了蒸气云团中的燃料量,再结合蒸气云爆炸模型利用Matlab软件进行了事故后果模拟计算,得出发生蒸气云爆炸时的死亡半径为36.629 5 m,重伤半径为83.557 6 m,轻伤半径124.725 m,财产损失半径为109.017 9 m。相较于无蒸发过程的传统模型,此计算结果更加具有参考意义。     

过程液氯贮槽装置封闭厂房及事故氯吸收塔联用安全效用分析

分析、设立液氯贮槽液相泄漏的几种典型事故情景,并利用液相泄漏、液池蒸发、重气扩散和人员中毒死亡概率等模型对比研究封闭厂房及事故氯吸收塔等安全措施对液氯贮槽液相泄漏扩散中毒后果的影响,给出不同事故情景下液氯泄漏速率、液池半径、液池蒸发速率、室外氯气中毒死亡概率等事故后果特征值。对封闭厂房及事故氯吸收塔安全效用进行定量分析和比较研究。结果表明,液氯贮槽的封闭厂房对抑制液氯泄漏扩散中毒事故后果效用明显;事故氯吸收塔能消除液氯贮槽微小孔泄漏所对应的小事故情景,还能对封闭厂房最严重泄漏事故后果起到初期削峰作用。显然,封闭厂房及事故氯吸收塔联用可以降低液氯贮槽事故影响后果,具有良好安全效用。     

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

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

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

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

飞机撞击安全壳燃油泄漏燃烧数值模拟研究

为研究飞机撞击核电厂安全壳后引发燃烧对安全壳的影响,采用Fluent模拟了航空煤油替代燃料在安全壳内部和外部的燃烧过程,分析了安全壳不同位置处的火焰温度变化情况。结果表明,在开始阶段,燃油泄漏时间比较短,燃油缓慢蒸发汽化,燃烧发展过程以扩散蔓延为主,高温区域主要分布在安全壳侧壁面和底面。随着时间推移,液面蒸发速度加快,火焰高度呈现纵向发展的态势,最后在顶部形成较为均匀的温度场。安全壳内部燃烧最高温度可以达到2 229 K,外部燃烧最高温度为1 308 K。此外,撞击位置越高,安全壳顶部会越早出现高温区域,而较低位置处由于液池的堆积缓慢,燃油蒸发速率低,温度上升相对较慢。区别于内部燃烧,外部燃烧火焰不会波及安全壳顶部中心位置。     

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

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

基于Fluent的危化品火灾中毒数值模拟

为了减少事故发生,以苯储罐为例,对危险化学品储罐区的事故后果进行研究.总结了目前一些典型事故后果模型——火灾、泄漏事故后果模型,选定了相适应的应用模型,然后建立适当的几何模型,通过设定各个边界条件,模拟出毒气扩散的范围以及池火灾的热辐射范围,进而采用Statistical软件拟合出泄漏后浓度以及热辐射量与距离之间的关系,为安全间距的设定奠定了良好的基础,对储罐区的本质安全化管理具有很好的指导作用.     

树冠火蔓延模型和数值分析

本文计及森林多孔性燃烧床含水量和有机物含量，以及渗流气体的受热膨胀和质量源，建立了树冠火蔓延的一维物理模型，给出火蔓延的控制方程组和定解条件。通过数值分析，得到树冠火蔓延速度特性曲线，研究水份蒸发和有机物热解对树冠火蔓延速度的影响。     

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

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

有风情况池火灾热辐射下的最小安全距离

防火间距是石油化工企业平面设计中的一个重要参数,开放环境下的火灾热辐射受大气稳定程度的影响.本文从计算流体力学角度出发,应用CFD(Computationsl Fluid Dynamics)软件Fluent,基于SCI爆炸火灾工程试验"Pool Fire A"的大气条件,对"有风情况下,直径为10 m的苯液池火灾"进行数值模拟,得出非绝热条件下,苯燃烧的峰温以及产物组分、池火灾对周围环境热辐射的空间分布.温度最高点在对称面y=0上,最高温度为1 478 K、火焰倾斜角度为32°(与竖直方向的夹角)、火焰高10.2 m.对于锰钢材料、内径为10 m苯储罐,2～3级风力情况下,相邻两储罐间最小安全距离在上风向为20 m,下风向为27 m.最后对模拟结果进行了分析.     

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.     

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.     

Numerical simulation of evaporation of volatile liquids

The paper presents the results of the validation of the developed pool evaporation model using literature and our own experimental data. The proposed model was used to examine the effect of wind velocity and pool sizes on the evaporation rate of volatile liquid (hexane). Contrary to the semi-empirical evaporation model widely used in hazard assessment, stronger dependence of evaporation rate on pool size at low wind speeds is obtained.     

Pool evaporation: Experimental tests at medium-scale with gasoline

Pool evaporation is a major source of flammable vapour clouds. Predicting the evaporation rate of a liquid hydrocarbon pool is therefore a key issue of dispersion modelling for safety concerns. This paper presents small- and medium-scale experiments of pool evaporation carried out with liquid hydrocarbons (pentane, heptane), hydrocarbon “gasoline-like” mixtures and gasoline. Liquid mass loss was measured and the evaporation rate deduced with its evolution in time. Other observations are highlighted, regarding the evolution of liquid temperatures, mixture compositions, and scale effects like the influence of pool length on surface evaporation rate. Comparisons with well-known correlations are then shown. The authors finally suggest a new semi-empirical correlation with a set of parameters fitted on the performed experiments.     

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.     

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.