细雾水与火焰的相互作用

本文分析了细雾水滴与浮升火焰羽流的相互作用,推导出雾滴在火场下降过程中的运动方程,、汽化速率方程和热量衡算方程。采用龙格－库塔差分方法,数值求解所得到的有关方程。     

Suppression effects of powder suppressants on the explosions of oxyhydrogen gas

Suppression tests of oxyhydrogen gas explosions were performed in an explosion tube with five types of dry powder used as the suppressants. The experimental results showed that the powder with large dust cloud density and small radius has better suppression effect, which agrees well with previous correlative results. Moreover, our results also showed that particles with chemical activity and light material density, their suppression effect are more prominent than that of the inert particles with heavy density. To discover the detailed suppression process of dust powder, governing equations were developed based on the homogeneous reactive two-phase flow. The TVD scheme and the Lax–Wendroff–Rubin scheme were adopted to solve the reactive gas phase and particle phase, respectively. The time splitting technique was employed to handle the stiffness of the coupled equations. Our calculated results showed that the dust cloud has the suppression effect on the explosion of oxyhydrogen gas, and with the increase of dust cloud density or the decrease of particle diameter, its suppression effect become more evident, which is in good agreement with our experimental results, in addition, the numerical results showed that with the same particle diameter, the suppression performance is enhanced with the reduction in particle material density.     

Effects of thermal radiation heat transfer on flame acceleration and transition to detonation in particle-cloud hydrogen flames

The current work examines regimes of the hydrogen–oxygen flame propagation and ignition of mixtures heated by radiation emitted from the flame. The gaseous phase is assumed to be transparent for the radiation, while the suspended particles of the dust cloud ahead of the flame absorb and reemit the radiation. The radiant heat absorbed by the particles is then lost by conduction to the surrounding unreacted gaseous phase so that the gas phase temperature lags that of the particles. The direct numerical simulations solve the full system of two phase gas dynamic time-dependent equations with a detailed chemical kinetics for a plane flames propagating through a dust cloud. It is shown that depending on the spatial distribution of the dispersed particles and on the value of radiation absorption length the consequence of the radiative preheating of the mixture ahead of the flame can be either the increase of the flame velocity for uniformly dispersed particles or ignition either new deflagration or detonation ahead of the original flame via the Zel'dovich gradient mechanism in the case of a layered particle-gas cloud deposits. In the latter case the ignited combustion regime depends on the radiation absorption length and correspondingly on the steepness of the formed temperature gradient in the preignition zone that can be treated independently of the primary flame. The impact of radiation heat transfer in a particle-laden flame is of paramount importance for better risk assessment and represents a route for understanding of dust explosion origin.     

Behavior of flames propagating through lycopodium dust clouds in a vertical duct

The structure of flame propagating through lycopodium dust clouds has been investigated experimentally. Upward propagating laminar flames in a vertical duct of 1800 mm height and 150×150 mm square cross-section are observed, and the leading flame front is also visualized using by a high-speed video camera. Although the dust concentration decreases slightly along the height of duct, the leading flame edge propagates upwards at a constant velocity. The maximum upward propagating velocity is 0.50 m/s at a dust concentration of 170 g/m³. Behind the upward propagating flame, some downward propagating flames are also observed. Despite the employment of nearly equal sized particles and its good dispersability and flowability, the reaction zone in lycopodium particles cloud shows the double flame structure in which isolated individual burning particles (0.5–1.0 mm in diameter) and the ball-shaped flames (2–4 mm in diameter; the combustion time of 4–6 ms) surrounding several particles are included. The ball-shaped flame appears as a faint flame in which several luminous spots are distributed, and then it turns into a luminous flame before disappearance. In order to distinguish these ball-shaped flames from others with some exceptions for merged flames, they are defined as independent flames in this study. The flame thickness in a lycopodium dust flame is observed to be 20 mm, about several orders of magnitude higher than that of a premixed gaseous flame. From the microscopic visualization, it was found that the flame front propagating through lycopodium particles is discontinuous and not smooth.     

Structure of flames propagating through aluminum particles cloud and combustion process of particles

Structure of flames propagating through aluminum particles clouds and combustion processes of the particles have been examined experimentally to understand the fundamental behavior of a metal dust explosion. The combustion process of individual aluminum particles in a flame propagating through the aluminum particles cloud has been recorded by using a high-speed video camera with a microscopic optical system, and analyzed. The flame is shown to be consisted of a preheat zone of about 3 mm thick, followed by a combustion zone of 5–7 mm thick. In the combustion zone, discrete gas phase flames are observed around each aluminum particle. Also an asymmetric flame around a particle is observed, which might be caused by an ejection of aluminum vapor from a crack of oxide shell surrounding the particle.     

Numerical analysis of phase behavior during rapid decompression of rich natural gases

The effect of the condensation process on the gas and liquid phase behavior during rapid decompression of rich natural gases is studied in the paper numerically. A one-dimensional mathematical model of transient thermal two-phase flow of compressible multi-component natural gas mixture and liquid phase in a shock tube is developed. The set of mass, momentum and enthalpy conservation equations are solved for the gas and liquid phases. The approach to model a liquid condensation process during rapid decompression of rich natural gas mixture is proposed. The mass transfer between the gas and the liquid is taken into account by introducing the appropriate terms into the governing equations. Thermo-physical properties of multi-component natural gas mixture are calculated by solving the Equation of State (EOS) in the form of the Soave–Redlich–Kwong (SRK-EOS) model. The proposed liquid condensation model is integrated into the GDP model. A simple case of GDP model, where the liquid was not considered, was extensively validated on base and dry natural gases. The proposed two-phase model is validated against the experiments where the decompression wave speed was measured in rich natural gases at low temperature. It shows a good agreement with the experimental data.     

Experimental study on water sealing of fire barriers and explosion venting of large-scale pipeline with low-concentration gas

Low-concentration gas transported in pipelines may lead to explosion accidents because gas with a concentration of less than 30% is prone to explode. To reduce the incidence of gas explosions, water sealing of fire barriers is implemented, and explosion venting devices are installed along the pipeline. To investigate their suppression effect on low-concentration gas explosion, experiments using methane–air premixed gas under different conditions were implemented on a DN500 pipeline test system. The effects of three types of explosion venting forms (rupture disc, asbestos board, and plastic film) on explosion overpressure and flame were compared and analysed. Results show that the rupture disc, asbestos board, and plastic film can achieve adequate explosion venting, causing the peak decay rates of explosion overpressure to reach 82.37%, 81.72%, and 90.79%, respectively. The foregoing indicates that the greater the static activation pressure of the explosion venting form, the higher the peak explosion overpressure at each measurement point. Moreover, the shorter the explosion flame duration, the greater the flame propagation velocity. The research results provide an essential theoretical foundation for the effective suppression of gas explosion accidents in the process of low-concentration gas transportation.     

水平管分层流下微孔泄漏特性数值模拟和实验分析

在油气田开发过程中，通常采用气液相混输模式，管道受腐蚀等因素影响容易出现穿孔而发生两相流泄漏。为分析两相流泄漏特性，对管内常见流型分层流下的微孔泄漏特性进行数值和实验分析；采用VOF耦合Level set算法分析了不同影响因素下的气液两相泄漏特性，设计了1种管道泄漏收集装置，进行室内两相流泄漏实验，并验证了数值预测模型的准确性。研究结果表明:气液两相流经过管壁泄漏口时会发生相分离，泄漏特性受小孔方位、管路内外压差、气液相流速影响较大；泄漏口位于管路侧壁时的泄漏特性与其他角度下的泄漏特性有所不同，可用泄漏影响区内的气液分布进行解释；当泄漏口位于管路底部时，存在临界液相分流系数，当液相分流比小于此临界值时，泄漏流体为单相液体。VOF耦合Level set算法的数值方法可为管路泄漏量预测和相分离特性分析提供参考。     

Numerical study on the spontaneous ignition of pressurized hydrogen release through a tube into air

Spontaneous ignition of pressurized hydrogen release through a tube into air is investigated using a modified version of the KIVA-3V CFD code. A mixture-averaged multi-component approach is used for accurate calculation of molecular transport. Autoignition and combustion chemistry is accounted for using a 21 step kinetic scheme. Ultra fine meshes are employed along with the Arbitrary Lagrangia–Eulerian (ALE) method to reduce false numerical diffusion. The study has demonstrated a possible mechanism for spontaneous ignition through molecular diffusion.

In the simulated scenario, the tube provided additional time to achieve a combustible mixture at the hydrogen–air contact surface. When the tube was sufficiently long under certain release pressure, autoignition would initiate inside the tube at the contact surface due to mass and energy exchange between low temperature hydrogen and shock-heated air through molecular diffusion. Following further development of the hydrogen jet downstream, the contact surface became distorted. Turbulence plays an important role for hydrogen/air mixing in the immediate vicinity of this distorted contact surface and led the initial laminar flame to transit into a stable turbulent flame.     

Influence of water storage on deflagration characteristics of methane in confined space

Comparative experiments were conducted under different water level heights and methane concentration conditions using a self-designed explosion experiment pipeline. The results showed that, in comparison with the scenario without water storage, when the water level was 2 cm and methane concentration was 6.5–12.5%, there was a dual effect including a pressure decrease caused by the endothermic cooling of liquid water and a pressure increase caused by the expansion of water vapour. These effects caused the pressure time history curve to exhibit a double-peak or multi-peak structure, and the average decrease in the peak deflagration pressure was 23.76%. The heat of vaporisation absorbed by the stored water and barrier effect of water vapour on the transfer of the heat slowed down the increase in the deflagration temperature. The average decrease in the peak deflagration temperature of methane was 13.82%, and the time to reach the peak deflagration temperature was extended as a whole, with an average delay of 0.22 s. Water storage also changed the shape of the deflagration flame front, which exhibited ‘knife’, ‘V’, and ‘crescent’ structures. Moreover, the flame propagation speed was significantly reduced, with the peak and average flame propagation speeds decreased by 83.3% and 83.6%, respectively. The research results can provide a certain reference for preventing gas explosions in typical confined spaces, and also help to explore new anti-explosion methods, which can be applied to marine equipment such as ships.     

大型相连容器中火焰传播的研究

为了进一步了解相连装置中粉尘爆炸的火焰传播行为和压力发展,为该结构的安全防护设计提供有价值的信息,采用大型实验装置对相连容器中玉米淀粉/空气混合物爆炸时的火焰传播行为进行了实验研究,同时采用已开发的数值模型对实验进行仿真计算。实验表明:粉尘浓度的变化对粉尘爆炸的火焰传播行为有重要影响;在粉尘浓度很低的情况下,火焰仍然能够在管道中加速传播且爆炸发展的最终结果相当猛烈。数值模型采用欧拉-拉格朗日方法模拟两相流现象,通过求解非稳态的湍流两相反应流守恒方程对实验进行二维仿真,计算结果与实验结果符合性较好,表明该模型可以很好地应用于粉尘爆炸火焰传播的研究。     

Modelling bund overtopping using a shallow water CFD model

This paper describes the development of a model for assessing the hazards arising from liquid spills from storage tanks containing flammable liquids. The model is based on computational fluid dynamics (CFD) solutions to the shallow water equations and is able to account for the interaction of a liquid spill with sloping ground and walls/bunds. These features mean that it has a much wider range of applicability than simple models and yet, as it is based on the solution to the two-dimensional shallow water equations, is relatively quick to run. The key parameter of interest in risk assessments is the fraction of a tank contents that would overtop the surrounding wall/bund that has been designed to retain it. A simple ‘sub-model’ has been developed to take account of this three-dimensional phenomenon within the two-dimensional model. The numerical solution of the shallow water equations is obtained through an adaptive method using two shock-capturing numerical techniques; the weighted average flux (WAF) and random choice methods (RCM). These provide a robust solution algorithm. A simple test case has been used to demonstrate the applicability of the model.     

Combustion modeling in large scale volumes using EUROPLEXUS code

Most of the numerical benchmarks on combustion in large scale volumes for hydrogen safety, which were performed up until today have demonstrated, that current numerical codes and physical models experience poor predictive capabilities at the industrial scale, both due to under-resolution and deficiencies in combustion modeling. This paper describes a validation of the EUROPLEXUS code against several large scale experimental data sets in order to improve its hydrogen combustion modeling capabilities in industrial settings (e.g. reactor buildings). The code is based on the Euler equations and employs an algorithm for the propagation of reactive interfaces, RDEM, which includes a combustion wave, as an integrable part of the Reactive Riemann problem, propagating with a fundamental flame speed (being a function of initial mixture properties as well as gas dynamics parameters). Validation of the first combustion model implemented in the code is based on obstacle-laden channels, interconnected reactor-type compartments, vented enclosures and covers all major premixed flame combustion regimes (slow, fast and detonation) with an aim to obtain conservative results. An improvement of this model is found in a direction of transient interaction of flame fronts with regions of elevated integral length scales presented in the velocity gradient field due to e.g. interactions with geometrical non-uniformities and pressure waves.     

水雾抑制气体爆炸火焰传播的实验研究

利用自行设计的全程透明的火焰加速管系统和细水雾实验系统 ,对不同水雾条件下的气体火焰传播现象进行了实验研究。运用光电传感器与CCD摄像技术 ,笔者分析了不同水雾条件下的甲烷预混气体火焰传播速度、传播火焰阵面轨迹 ;探讨了水雾抑制气体火焰传播的机理及条件。实验发现了在一定条件水雾作用下的气体传播火焰阵面拉伸与火焰驻留的现象与条件 ,实验结果表明 :水雾对气体爆炸火焰传播的抑制是由于水雾作用于火焰阵面反应区 ,降低了反应区内火焰温度和气体燃烧速度 ,减缓了火焰阵面传热与传质的进行 ,从而使传播火焰得以抑制 ;而水雾对气体爆炸火焰传播的抑制效果与水雾通量、雾区浓度、水雾区长度以及火焰到达水雾区的火焰传播速度有关     

基于实验与仿真的三维水下气体泄漏预测与安全评估研究

水下气体羽流特性是海底气体泄漏风险评估的重要基础。为准确预测水下气体羽流行为，基于计算流体动力学(CFD)方法，建立1种考虑气体卷吸湍流特性的三维水下气体羽流数值预测模型；采用欧拉-欧拉流体体积模型捕捉气液作用界面，以大涡模拟(LES)方式预测羽流上浮及卷吸过程中的湍流特性，从而实现对水下气体羽流行为的预测；搭建小尺度实验平台，对比仿真与实验条件下的气体羽流形态，验证数值模型的可行性及预测精度；应用建立的数值模型对工程条件下的水下气体泄漏事故进行预测和评估，以某浅层气井喷事故为例，评估水下气体羽流上浮时间、海面影响范围和涌流高度。结果表明:基于欧拉多相流与大涡模拟的数值模型对水下气体羽流预测结果与实验具有较好的吻合度，该模型能够较好捕捉羽流的湍流特性，可为水下气体泄漏羽流行为评估提供参考。     

Flame propagation through aluminum particle cloud in a combustion tube

The propagation of a flame is investigated experimentally and theoretically for a large, horizontal combustion tube containing a mixture of air and aluminum powder with pre-existing turbulence. One end of the tube is closed and the other is connected to a large dump-tank. Twenty dispersion systems are used on the tube to produce a uniform suspension of aluminum dust in the tube with a mean diameter of 6 μm. The characteristics of a flame front from the ignitors at the closed end are measured using photodiodes and the development of pressure is monitored by transducers. Experimental results revealed the entire process of an accelerating flame and the development of shock waves. A set of conservation equations for two-phase turbulent combustion flow is derived, using the two-fluid model, k–ε model, Hinze–Tchen model and EBU-Arrhenius model for turbulent combustion. The SIMPLE scheme usually applied to the homogeneous turbulent combustion is extended to fit this two-phase, reactive behavior. The results of calculations show the positive feedback coupling among combustion, expansion and turbulence during flame propagation. Computed and measured results are generally in good agreement.     

A numerical method to predict flame fractal dimension during gas explosion

Accurate prediction of the flame propagation velocity during a gas explosion is essential to assess its consequences and to evaluate the risk level. The propagating premixed flame is self-turbulized due to the hydrodynamic instability, resulting in a fractal flame structure. It is therefore important for accurate prediction of flame speed to understand the flame’s fractal structure in detail and to predict its fractal dimension in particular. Numerical simulations of spherically-propagating flames have been previously attempted for such purposes. There are, however, difficulties to accurately predict the fractal dimension from the result of the numerical simulation of a spherically-propagating flame. In this study, we propose a method to easily predict the fractal dimension based on the numerical simulation of a planar flame. Planar flame propagation is simulated for different sizes of computational domain. The fractal dimension can be determined from the dependence of flame speed on computational domain size. The determined fractal dimension is favorably compared with previous experimental results.     

火焰沿斜面可燃固体板传播的数值模拟

本文采用数值模拟方法，详细研究了倾角对火焰沿斜面可燃固体板传播的影响；计算结果表明：随着可燃固体板的倾角增加，火焰在其上面的传播速率不断增大，这一结论对建筑物的火灾安全设计具有指导意义。     

LES – DFSD modelling of vented hydrogen explosions in a small-scale combustion chamber

Accidental explosions are a plausible danger to the chemical process industries. In the event of a gas explosion, any obstacles placed within the path of the flame generate turbulence, which accelerates the transient flame and raises explosion overpressure, posing a safety hazard. This paper presents numerical studies using an in-house computational fluid dynamics (CFD) model for lean premixed hydrogen/air flame propagations with an equivalence ratio of 0.7. A laboratory-scale combustion chamber is used with repeated solid obstacles. The transient compressible large eddy simulation (LES) modelling technique combined with a dynamic flame surface density (DFSD) combustion model is used to carry out the numerical simulations in three-dimensional space. The study presented uses eight different baffle configurations with two solid obstructions, which have area blockage ratios of 0.24 and 0.5. The flame speed, maximum rate of pressure-rise as well as peak overpressure magnitude and timing are presented and discussed. Numerical results are validated against available published experimental data. It is concluded that, increasing the solid obstacle area blockage ratio and the number of consecutive baffles results in a raised maximum rate of pressure rise, higher peak explosion overpressure and faster flame propagation. Future model development would require more experimental data, probably in a more congested configuration.     

水雾作用下气体爆炸火焰传播的实验研究

利用自行设计的全程透明的火焰加速管和水喷雾系统，对不同水雾条件下的气体火焰传播现象进行了实验研究。运用光电传感器与数字摄像技术分析了不同浓度的甲烷在不同水雾条件下的火焰传播速度、火焰阵面轨迹以及火焰结构特性，并通过对传播火焰反应区温度的测量，探讨了水雾抑制气体爆炸火焰传播的内在机理及所需的条件。实验结果表明：由于水雾作用于火焰反应区，降低了火焰反应区内的温度和气体燃烧速度，延长火焰阵面的预热区，减缓火焰阵面传热与传质的进行，从而使传播火焰得以抑制。水雾对气体爆炸火焰传播的抑制效果与水雾流量速度、雾区浓度以及火焰到达水雾区的火焰传播速度有关。