Droplets behaviors of flashing jet in accidental release of superheated liquid

Structural failure of an industrial superheated liquid tank or pipe usually results in the flashing jet consisting of a mixture of massive droplets and vapor due to the violent phase transition. In this work, experiments on small-scale releases were carried out with a 20 L storage tank to investigate the droplets behaviors of flashing jets after accidental releases of superheated liquid. Distribution of droplets axial, radial and vertical velocities, as well as droplets diameter along the centreline of flashing jet, were acquired employing a Phase Doppler Anemometry (PDA). The influence of storage pressure, superheat, and nozzle diameter was also studied. Results show that the distribution of actual droplets with different velocities and diameters follows the normal distribution. Droplets mean three-dimensional velocities on the central axis of flashing jet decrease exponentially with the increase of axial distance. The droplets mean diameter first decreases exponentially and finally keeps sable at about 10 μm. Among three-dimensional velocities, the axial velocity is the highest and the vertical velocity is the lowest. Droplets mean three-dimensional velocities increase with the increase of storage pressure, superheat, and nozzle diameter. The droplets mean diameter decreases with the increase of superheat and nozzle diameter but increases with the increase of storage pressure.     

Assessment of the maximum possible liquid superheat during flashing leak flow

During the discharge of flashing liquids through leaks due to abrupt depressurization a transient thermodynamic non-equilibrium in the form of a boiling delay in the superheated liquid flow can occur. As a consequence the actual mass flow quality is smaller than calculated under the assumption of an immediate adjustment of the thermodynamic equilibrium between the phases. For the prediction of the leak mass flow for a given pressure difference the magnitude of this self-adjusting mass flow quality is needed.

Most of the models cited in the literature include only the equilibrium mass quality as limiting quantity and ignore further effects as that of the depressurization velocity or the mean nucleus distance. For the assessment of the maximum possible liquid superheat during flashing only the conduction heat transfer from a stagnant liquid to the bubble surface is used to describe the bubble growth.

The sub-model for the bubble growth due to expansion and mass transfer necessary for the global prediction of the transient thermodynamic non-equilibrium in flashing liquids was validated using bubble radii measured by Hooper et al. [Bubble growth and pressure relationship in the flashing of superheated water. Technical publication 6904, Mechanical Engineering Department, University of Toronto, 1969] for the case of a sudden depressurization of initially saturated water. On this basis the calculated time-dependent temperature field, the actual mass quality, the mean liquid temperature and, in comparison to the corresponding values based on the assumption of immediate thermodynamic equilibrium, the maximum possible liquid superheat are predicted.     

Simulation study on near-field structure and flow characteristics of high-pressure CO2 released from the pipeline

The accidental release of high-pressure carbon dioxide (CO₂) can cause serious damages to both humans and pipeline equipment. Therefore, it is of great significance to have a deeper understanding about the release characteristics of high-pressure CO₂ for improving the safety level of Carbon Capture and Storage (CCS) technologies. Both industrial-scale and laboratory-scale studies have been carried out to predict the release behaviors. In recent years, computational fluid dynamics (CFD) simulation has become a crucial method to study the instantaneous changes and microscopic details of the fluid behaviors. In this paper, the simulation method was employed to study the near-field structure and flow characteristics of high-pressure CO₂ released from pipelines. The Peng-Robinson Equation of State (EOS) was used to compute the thermodynamic properties of high-pressure CO₂, and SST k-ω model was applied to simulate the structure and physical parameters of the under-expanded jet. In addition, the multi-phase mixture model was introduced to study the phase transition. The non-equilibrium liquid/vapor transition is modeled by introducing ‘source terms’ for mass transfer and latent heat. Compared to the experimental results, the simulation results showed good agreement. Furthermore, the influences of operating conditions, including different stagnation pressure, stagnation temperature, and nozzle size, were analyzed.     

Sub-cooled and flashing liquid jets and droplet dispersion I. Overview and model implementation/validation

Many accidents involve two-phase releases of hazardous chemicals into the atmosphere. This paper describes the results of a third phase of a Joint Industry Project (JIP) on liquid jets and two-phase droplet dispersion. The aim of the project is to increase the understanding of the behaviour of sub-cooled non-flashing and superheated flashing liquid jets, and to improve the prediction of initial droplet size, droplet dispersion and rainout.Phase III of the JIP first included scaled experiments for materials with a range of volatilities (water, cyclohexane, butane, propane and gasoline). These experiments were carried out by Cardiff University including measurements of flow rate and initial droplet size across the full relevant range of superheats. See the companion paper II for further details of these experiments and the derivation of new refined correlations for droplet size distribution and Sauter Mean Diameter. Furthermore large-scale butane experiments were carried out by INERIS (France) to ensure that for more realistic scenarios the derived droplet size correlations are accurate.Model validation and model improvements were carried out by DNV Software, including validation of release rate and initial droplet size against the above scaled and large-scale experiments. New correlations for droplet size distribution and Sauter Mean Diameter (SMD) were implemented into the Phast discharge model. These were compared against a range of other droplet size and rainout correlations published in the literature, in conjunction with validation against an extensive set of experiments. It was shown that the new droplet size correlation agrees better against experimental data than the existing Phast correlation. To further improve the rainout prediction, the Phast dispersion model (UDM) was also extended to allow simultaneous modelling of a range of droplet sizes and distributed rainout (rather than rainout at one point).     

Phast validation of discharge and atmospheric dispersion for pressurised carbon dioxide releases

The consequence modelling package Phast examines the progress of a potential incident from the initial release to the far-field dispersion including the modelling of rainout and subsequent vaporisation. The original Phast discharge and dispersion models allow the released substance to occur only in the vapour and liquid phases. The latest versions of Phast include extended models which also allow for the occurrence of fluid to solid transition for carbon dioxide (CO₂) releases.As part of two projects funded by BP and Shell (made publicly available via CO2PIPETRANS JIP), experimental work on CO₂ releases was carried out at the Spadeadam site (UK) by GL Noble Denton. These experiments included both high-pressure steady-state and time-varying cold releases (liquid storage) and high-pressure time-varying supercritical hot releases (vapour storage). The CO₂ was stored in a vessel with attached pipework. At the end of the pipework a nozzle was attached, where the nozzle diameter was varied.This paper discusses the validation of Phast against the above experiments. The flow rate was predicted accurately by the Phast discharge models (within 10%; considered within the accuracy at which the BP experimental data were measured), and the concentrations were found to be predicted accurately (well within a factor of two) by the Phast dispersion model (UDM). This validation was carried out with no fitting whatsoever of the Phast extended discharge and dispersion models.     

Influence of the elastic behaviour of liquids on transient flow from pressurised containments

Release of liquids from vessels and pipelines is a very important scenario in consequence assessment. When the operating pressure is significant, its evolution and the corresponding discharge rate are affected by the elastic behaviour of the liquid. Bulk modulus is the key-driver of the related fluid-dynamics and is expected to govern the elastic-to-atmospheric pressure transition. Despite the well known physical background, the availability of release models relevant to the elastic phase is rather poor, compared to those relevant to vapours, gases and atmospheric liquids. On the other hand, if the operating pressure is high, release times and dynamics are expected to be strongly affected by the elastic behaviour and prediction of release time and flow rate is fundamental. This article carries outs a general analysis of the elastic behaviour of liquids, including water and hydrocarbons. Lumped (vessel-type) and distributed (pipeline) parameters systems have been modelled with the aim to characterize the pressure evolution during the elastic phase and to evaluate the significance of the associated duration. The findings have shown that, depending on pressure and geometric features, the importance of the elastic phase is high and that the availability of reliable and relatively simple models could be beneficial in consequence assessment. The approach is strictly valid for non-boiling liquid, but accounting for vaporisation of flashing or boiling liquids such as LPG and LNG would not imply significant complications.     

Experimental study on the characteristics of the strong boiling induced by pressure relief at the top of vertical vessels

In this study, the strong boiling process of a superheated liquid induced by pressure relief at the top of a vertical vessel was experimentally investigated. Through monitoring the response of the pressure at the top of the vessel and the full field morphology of the two-phase flow, the correlations between the pressure rise and the two-phase flow as well as the trend of the characteristic pressure values under different discharge areas and filling rates were analyzed. The results indicated that the expansion of the mist-like two-phase flow, which was generated due to the strong boiling and bubble collapse, was the direct reason for the two pressure rebounds at the top of the vessel. And under the effect of the intermittent expansion of the two-phase flow, the pressure rising rate in the second rebound stage fluctuated. When the discharge area changed, the characteristic pressure values presented different relativeness under different filling rates. In addition, the average pressure rising rate during the second pressure rebound process presented an exponential growth with the increase of the discharge area, and the exponent coefficient decreased with the increase of the liquid filling rate.     

化学品突发性事故预测的不确定性分析

介绍化学品突发性事故预测的不确定性概率统计和模糊数学分析方法。对参数和模式的不确定性分别进行了讨论,计算机系统包括化学品常用数据库、事故识别模块,事故后果分析模块和不确定性分析模块,对液氯泄漏排放事故、液化气（丙烷）两相闪蒸爆炸和管道煤气（CO）泄漏伤害模式实例进行了分析。对连续源和瞬时源的讨论说明泄漏时间的确定对事故分析极为重要。     

Experimental investigation on spontaneous ignition caused by pressurized hydrogen suddenly release into an S-shaped tube

An experimental investigation on the effects of continuous semicircular curved structure on spontaneous ignition during pressurized hydrogen suddenly release was conducted. An S-shaped tube with 700 mm in length and 10 mm in diameter was used in our experiments, and a straight tube with the same configuration was adopted for comparison. The results show that the continuously generated rarefaction waves and reflected shock waves make the pressure curves in the S-shaped tube more complicated. Meanwhile, the mean velocity and intensity of the leading shock wave undergo considerable attenuation when it propagates in the S-shaped structure. By comparing with the straight tube, the minimum critical pressure condition for spontaneous ignition in the S-shaped tube is slightly difficult to reach, but the difference is not huge. Nevertheless, the S-shaped structure can effectively promote hydrogen-air mixing and make combustion more intense. A secondary overpressure peak detected by the pressure transducer near the nozzle occurs in the spontaneous ignition cases and no such pressure increase is caught in the non-ignition cases. The transition from spontaneous combustion flame to a jet flame at the nozzle and the complete out-tube jet flame development process are captured and discussed.     

Expansion-controlled evaporation: a safe approach to BLEVE blast

A new method is presented to calculate the blast effects originating from an exploding vessel of liquefied gas. Adequate blast calculation requires full knowledge of the blast source characteristics, i.e., the release and consequent evaporation rate of the flashing liquid. As the conditions that allow explosive evaporation are not entirely clear and the evaporation rate of a flashing liquid is unknown, safe assumptions are the starting point in the modelling. The blast effects from a BLEVE are numerically computed by imposing the vapour pressure of a flashing liquid as boundary condition for the gas dynamics of expansion. The numerical modelling is quantitatively explored just for liquefied propane. In addition, it is demonstrated that often an estimate of BLEVE blast effects is possible with very simple acoustic volume source expressions.

The modelling shows that the rupture of a pressure vessel containing a liquefied gas in free space only develops a blast of significant strength if the vessel nearly instantaneously disintegrates. Even if a rupture and the consequent release and evaporation of a liquefied gas extend over just a short period of time, the blast effects are minor.     

A simple,explicit formula for the critical pressure of homogenous two-phase nozzle flows

The critical pressure ratio of the homogeneous two-phase nozzle flow model known as the Omega method is expressed in function of the Omega Parameter as the exact numerical solution of a transcendental equation. A well fitting, easy to use, explicit approximation for flashing and non-flashing flows is presented here. The validation against the exact numerical solution proves that this new formula is better fitting than the other ones in the technical literature for both single and two-component flows.     

Sub-cooled and flashing liquid jets and droplet dispersion II. Scaled experiments and derivation of droplet size correlations

This paper describes the results of the first stage of Phase III of a Joint Industry Project (JIP) on liquid jets and two-phase droplet dispersion. This stage included scaled experiments for water, gasoline, and cyclohexane for a range of superheats and nozzles with different aspect ratios. Additional experiments for butane and propane were conducted as a validation exercise and are discussed in the companion paper. Moreover this paper provides recommendations for atomisation correlations in the regimes of mechanical break-up, transition to flashing, and fully flashing. The objectives of this scaled experimental programme are to : (i) provide confidence in the previously proposed modelling methodology (Phase II) across a broad range of initial conditions (ii) update the models’ correlations to generalise further its applicability (iii) recommend further model improvements. Development of new correlations for Sauter Mean Droplet diameter (SMD) and droplet size distribution is based on a best fit of the current scaled experimental data. The new data endorses the previous tri-functional Phase II approach including regimes for mechanical break-up, transition to flashing, and fully flashing, which is hence updated in the new Phase III SMD model. Considerable effort is devoted to capturing the full droplet size range under low-superheat conditions. Also, new enhancements in PDA technology were adopted to enable better quality data under high-superheat conditions. The priority recommendation for further model improvement is better characterisation of the poor quality releases under low-superheat conditions, where this work indicates that extremely broad droplet size distributions are likely. A companion paper (Part I) presents a more general overview of the dispersion problem, implementation of the correlations and subsequent performance against both the current scaled experiments and additional large-scale butane experiments.     

Effects of release pressure on the transportation characteristics in pipeline and the diffusion behaviors in enclosure space of typical gas extinguishing agent

Reducing the release pressure of fire extinguishing systems can decrease potential safety hazards in large transport airplanes. To explore whether reducing the release pressure can achieve the release effect required by the airworthiness standards or not, the transportation characteristics in the pipeline and diffusion behaviors in the enclosure space of a typical fire extinguishing agent (Halon 1301) were investigated under five release pressures in the present study. The effects of the release pressure on the degree of superheat, injection duration, jet structure, and concentration distribution of Halon 1301 were analyzed. The results show that both of the degree of superheat and the injection duration decrease with an increase in the release pressure. The supplement of bubble expansion in the fire extinguishing agent can slow down the pressure decrease in the vessel. Both of the maximum and mean values of the pipeline differential pressure increase with an increase in release pressure. The maximum value of the jet angle decreases linearly with the increase in release pressure, and the jet deflects upward owing to the effects of buoyancy. The maximum concentration value decreases with an increase in the distance from the nozzle. The maximum concentration values in the near field from the nozzle increase with an increase in the release pressure. Under five release pressures, the concentration and holding time (duration above 6% volume concentration) of Halon 1301 on the centerline of the jet meet the requirements of airworthiness provisions.     

Two-phase jet releases,droplet dispersion and rainout I. Overview and model validation

Many accidents involve two-phase releases of hazardous chemicals into the atmosphere. This paper describes the results of the fourth phase of a Joint Industry Project (JIP) on liquid jets and two-phase droplet dispersion. The objective of Phase IV of the JIP was to generate experimental rainout data for non-flashing experiments, and to develop recommendations for the best methodology to predict rainout [total rainout mass and its spatial distribution (‘distributed’ rainout)].Phase IV of the JIP first included rainout experiments by the UK Health and Safety Laboratory (HSL) for sub-cooled releases of water and xylene with a range of orifice sizes and stagnation pressures. See the companion paper II for further details. Secondly model validation was carried out by DNV Software for these experiments using different correlations for the initial droplet size (Sauter Mean Diameter, SMD), i.e. the CCPS SMD correlation and the Phase III JIP SMD correlation. The validation includes flow rates, droplet size, distributed rainout and cloud temperature drop. Subsequently validation was considered for a wider range of experiments from the literature (sub-cooled and superheated releases) for both SMD and total rainout. Adopted rainout methods comprised both methods including explicit modelling of the droplets (using an extended version of Phast dispersion model UDM), as well as more simple methods based on rainout correlations without droplet modelling. Recommendations are made for the most accurate droplet size and rainout modelling. A modified CCPS UDM droplet size correlation has been shown to agree best against experimental rainout data.     

Modelling of discharge and atmospheric dispersion for carbon dioxide releases

This paper discusses the modelling of the discharge and subsequent atmospheric dispersion for carbon dioxide releases using extensions of models in the consequence modelling package Phast. Phast examines the progress of a potential incident from the initial release to the far-field dispersion including the modelling of rainout and subsequent vaporisation. The original Phast discharge and dispersion models allow the released chemical to occur only in the vapour and liquid phases. As part of the current work these models have been extended to also allow for the occurrence of liquid to solid transition or vapour to solid transition. This applies both for the post-expansion state in the discharge model, as well as for the thermodynamic calculations by the dispersion model. Solid property calculations have been added where necessary. The above extensions are generally valid for fluid releases including CO₂. Using the extended dispersion formulation, a sensitivity study has been carried out for mixing of solid CO₂ with air, and it is demonstrated that solid effects may significantly affect the predicted concentrations.     

液氮泡沫发生器内流动特性研究

低温液氮与泡沫混合液直接接触产生氮气泡沫是一种新型的掺混形式,利用液氮高汽化比的特点,搭建液氮泡沫可视化实验装置,进行氮气-水两相流及液氮泡沫流动特性的研究。结果表明,液氮相变产生大量氮气,其与泡沫液混合产生泡沫,温度有所回升,最终趋于泡沫混合液温度;管路沿程压降较小;液氮射流破碎及流动过程可分为6个区域:低温液氮区、向上循环翻滚区、滞留区、泡沫与泡沫混合液混合区、致密泡沫区、泡沫混合液区。流体向下游流动过程中持续发泡;为防止管路结冰,需合理控制泡沫混合液与液氮流量。     

前混合式磨料射流喷嘴的基础研究

本文对前混合式磨料射流喷嘴进行了较为系统的研究。利用液固两相流的理论进行了磨料射流喷嘴内的流动分析，采用无量纲分析法，列出磨料粒子在喷嘴内的质量方程、运动方程和能量方程，编写了计算程序，通过计算机进行数值求解，发观了喷嘴收缩段锥角、磨料浓度对磨料粒子的加速、压力损失以及能量传递效率的影响，确定出了最优参数。根据理论分析，对几种不同几何形状、不同几何参数的喷嘴，通过改变压力、磨料浓度等参数进行寿命磨损试验，定期测量其磨损，观察喷嘴内部的磨损状况，研究了喷嘴材料、几何参数及动力参数对喷嘴磨损的影响，发现了一些新现象，通过对试验数据的分析，得出了一些有价值的结论。     

控制排气放空噪声的消声器之研究与应用

排气放空噪声是工业上常遇到的噪声源。本文介绍了根据喷注噪声发射理论研制成的小孔喷注与节流降压相复合的新型消声器，适用于高压气源排放，取得了较好的消声效果。鉴于喷注发射噪声的强度取决于喷注的驻压与喷口直径等参量，对节流降压层的设计与控制小孔喷注层的驻压大小，提出了新的见解，可供设计研究喷注消声器参考使用。     

自进式旋转钻头钻孔修复理论及关键参数研究

瓦斯预抽钻孔在抽放过程中经常出现塌孔、堵孔等现象。针对这一问题,提出采用自进式旋转钻头修复失效钻孔的新方法。其利用自进式旋转钻头后置喷嘴的喷射反冲力作为动力,自行至钻孔堵塞段处,前置喷嘴形成的高压水射流对堵孔煤渣进行破碎,将其与水混合返出孔外,实现清孔排渣、钻孔修复的目的。对自进式旋转钻头钻孔破煤清渣的临界射流压力及流量进行分析,基于摩擦动力学计算钻头自驱的临界射流压力,在此基础上设计钻头结构及优化喷嘴布置,数值分析不同喷嘴组合情况下射流速度分布,得出最优喷嘴布置方式。并应用于鹤煤八矿-655轨道石门揭煤瓦斯预抽,实验结果表明:该方法能够实现失效钻孔的修复,并有效的提高了瓦斯抽采效果,缩短了瓦斯抽放时间。     

Blast overpressures from medium scale BLEVE tests

The measured blast overpressures from recent tests involving boiling liquid expanding vapour explosions (BLEVE) has been studied. The blast data came from tests where 0.4 and 2 m³ ASME code propane tanks were exposed to torch and pool fires. In total almost 60 tanks were tested, and of these nearly 20 resulted in catastrophic failures and BLEVEs. Both single and two-step BLEVEs were observed in these tests. This paper presents an analysis of the blast overpressures created by these BLEVEs. In addition, the blast overpressures from a recent full scale fire test of a rail tank car is included in the analysis.The results suggest that the liquid energy content did not contribute to the shock overpressures in the near or far field. The liquid flashing and expansion does produce a local overpressure by dynamic pressure effects but it does not appear to produce a shock wave. The shock overpressures could be estimated from the vapour energy alone for all the tests considered. This was true for liquid temperatures at failure that were below, at and above the atmospheric superheat limit for propane. Data suggests that the two step type BLEVE produces the strongest overpressure. The authors give their ideas for this observation.The results shown here add some limited evidence to support previous researchers claims that the liquid flashing process is too slow to generate a shock. It suggests that liquid temperatures at or above the Tsl do not change this. The expansion of the flashing liquid contributes to other hazards such as projectiles, and close in dynamic pressure effects. Of course BLEVE releases in enclosed spaces such as tunnels or buildings have different hazards.