CFD modeling and simulation of turbulent fluid flow and dust dispersion in the 20-L explosion vessel equipped with the perforated annular nozzle

A three-dimensional CFD model was developed to simulate the turbulent flow field induced by dust feeding and the associated dust dispersion within the 20-L explosion vessel equipped with the perforated annular nozzle. The model was validated against experimental data for pressure and root mean square velocity.Simulation results have shown that the turbulent kinetic energy is rather uniformly distributed and its values are significantly lower than those attained with the rebound nozzle. Furthermore, the perforated annular nozzle is able to generate a uniform dust/air cloud. However, a consistent fraction of the dust remains trapped inside the nozzle and, thus, it does not contribute to the explosion process.     

Effect of separation distance on gas dispersion and vapor cloud explosion in a storage tank farm determined using computational fluid dynamics

Storage tank separation distance, which considerably affects forestalling and mitigating accident consequences, is principally determined by thermal radiation modeling and meeting industry safety requirements. However, little is known about the influence of separation distance on gas dispersion or gas explosion, which are the most destructive types of accidents in industrial settings. This study evaluated the effect of separation distance on gas dispersion and vapor cloud explosion in a storage tank farm. Experiments were conducted using Flame Acceleration Simulator, an advanced computational fluid dynamics software program. Codes governing the design of separation distances in China and the United States were compared. A series of geometrical models of storage tanks with various separation distances were established. Overall, increasing separation distance led to a substantial reduction in vapor cloud volume and size in most cases. Notably, a 1.0 storage diameter separation distance appeared to be optimal. In terms of vapor cloud explosion, a greater separation distance had a marked effect on mitigating overpressure in gas explosions. Therefore, separation distance merited consideration in the design of storage tanks to prevent gas dispersion and explosion.     

Determination of heterogeneous reaction mechanisms: A key milestone in dust explosion modelling

Reaction kinetics is fundamental for modelling the thermal oxidation of a solid phase, in processes such as dust explosions, combustion or gasification. The methodology followed in this study consists in i) the experimental identification of the reaction mechanisms involved in the explosion of organic powders, ii) the proposal of simplified mechanisms of pyrolysis and oxidation, iii) the implementation of the model to assess the explosion severity of organic dusts. Flash pyrolysis and combustion experiments were carried out on starch (22 μm) and cellulose (53 μm) at temperatures ranging from 973 K to 1173 K. The gases generated were collected and analyzed by gas chromatography. In this paper, a semi-global pyrolysis model was developed for reactive systems with low Damköhler number. It is in good agreement with the experimental data and shows that both carbon monoxide and hydrogen are mainly generated during the pyrolysis of the solid, the generation of the latter compound being greatly promoted at high temperature. A simplified combustion model was also proposed by adding two oxidation reactions of the pyrolysis products. In parallel, flame propagation tests were performed in a semi open tube in order to assess the burning velocity of such compounds. The laminar burning velocity of cellulose was determined to be 21 cm s⁻¹. Finally, this model will be integrated to a predictive model of dust explosions and its validation will be based on experimental data obtained using the 20 L explosion sphere. The explosion severity of cellulose was determined and will be used to develop and adjust the predictive model.     

Effects of multiple annular obstacles on flame propagation of local corn starch dust in a vertical pipe

With high-speed camera technology, the propagation behavior of explosion flame for the local dust cloud of corn starch in a semi-open vertical pipe under the action of the annular obstacle was studied experimentally, and the blockage rate and the annular obstacle numbers as well as impact of dust cloud concentration on the flame propagation were investigated. The researches showed that both the blockage rate and the annular obstacle numbers have significant effects on the flame speed and propagation process for the dust cloud explosion of corn starch. The increase of the blockage rate of such annular obstacles will cause that the combustion of dust cloud with high concentration is mainly concentrated in the lower part of the pipe. The increase of the annular obstacle numbers will lead to the acceleration of combustion of the dust cloud. With the increase of the blockage rate and the annular obstacle numbers, the maximum flame speed shows a trend of the first increasing and then decreasing, and the phenomenon of accelerated propagation of the flame becomes more and more obvious, however, the distance of continuous acceleration for the flame is gradually decreased and the maximum flame speed is farther from the outlet of the pipe. Under the action of such annular obstacles, the concentration of dust cloud has a significant effect on the flame speed and shape of the dust cloud of the corn starch. The increase of the concentration of the dust cloud will decrease the acceleration effect of such annular obstacles to result in maximum flame speed showing a trend of the first increasing and then decreasing. However, the acceleration distance of the flame is longer, and the maximum flame speed is closer to the outlet of the pipe. The increasing concentration will make the flame speed develop more slowly, the flame color will be darker, and the flame segmentation phenomenon will be more obvious.     

Numerical simulation study of dust concentration distribution regularity in cavern stope

Based on the theory of gas-solid two-phase flow and the characteristics of cavern stope a model of dust migration was established. The dust concentration changing of cavern stope by ventilation in 20 min after blasting and the dust trajectory in different wind speed were simulated by Fluent Software. The results show that distribution of dust concentration is significantly affected by flow field of airway in cavern, and the dust concentration of inlet is higher than that of outlet and the highest one on the corner of inlet’s side. In the stope, the smaller the wind speed of inlet is, the shorter of dust can be captured, settled and discharged, the more obviously affected by the trajectory of gas flow field. It goes into the stage of clean cycle emissions after 60 s, the speed of dust concentration dropped is the biggest between 0 s and 70 s, the main dust in stope is respirable dust after 70 s, it needs much time to settlement.According to the measured data of metal mining, approximately 87% of dust was generated during the drilling and blasting in the mine (Wang, 1979). A lot of dust with high concentrations was produced during the cavern stope blasting and it was difficult to be discharged. It can help choose the right speed to rule out the dust quickly which produced during cavern blasting, if the dust concentration distribution and the dust migration law of different inlet velocity in the cavern can be verified, what’s more, the labor productivity can be increased. It has great significance for choosing reasonable ventilation parameters, reducing dust hazards of stope to researching the dust concentration distribution regularity in the stope.     

CFD simulations of dust dispersion in the 20 L vessel: Effect of nominal dust concentration

Measurements of flammability and explosion parameters for dust/air mixtures require uniform dispersion of the dust cloud inside the test vessel. In a previous work, we showed that, in the standard 20 L sphere, the dust injection system does not allow generation of a uniform cloud, but rather high gradients of dust concentration are established. In this work, we used a previously validated three-dimensional CFD model to simulate the dust dispersion inside the 20 L sphere at different dust nominal concentrations (and fixed dust diameter). Results of numerical simulations have shown that, as the dust nominal concentration is increased, sedimentation prevails and, thus, when ignition is provided, the dust is mainly concentrated at the vessel walls.     

Effects of premixed methane concentration on distribution of flame region and hazard effects in a tube and a tunnel gas explosion

Study of flame distribution laws and the hazard effects in a tunnel gas explosion accident is of great importance for safety issue. However, it has not yet been fully explored. The object of present work is mainly to study the effects of premixed gas concentration on the distribution law of the flame region and the hazard effects involving methane-air explosion in a tube and a tunnel based on experimental and numerical results. The experiments were conducted in a tube with one end closed and the other open. The tube was partially filled with premixed methane-air mixture with six different premixed methane concentrations. Major simulation works were performed in a full-scale tunnel with a length of 1000 m. The first 56 m of the tunnel were occupied by methane–air mixture. Results show that the flame region is always longer than the original gas region in any case. Concentration has significant effects on the flame region distribution and the explosion behaviors. In the tube, peak overpressures and maximum rates of overpressure rise (dp/dt)_max for mixtures with lower and higher concentrations are great lower than that for mixtures close to stoichiometric concentration. Due to the gas diffusion effect, not the stoichiometric mixture but the mixture with a slightly higher concentration of 11% gets the highest peak overpressure and the shock wave speed along the tube. In the full-scale tunnel, for fuel lean and stoichiometric mixture, the maximum peak combustion rates is achieved before arriving at the boundary of the original methane accumulation region, while for fuel rich mixture, the maximum value appears beyond the region. It is also found that the flame region for the case of stoichiometric mixture is the shortest as 72 m since the higher explosion intensity shortens the gas diffusion time. The case for concentration of 13% can reach up to a longest value of 128 m for longer diffusion time and the abundant fuel. The “serious injury and death” zone caused by shock wave may reach up to 3–8 times of the length of the original methane occupied region, which is the widest damage region.     

Assessment of four turbulence models in simulation of large-scale pool fires in the presence of wind using computational fluid dynamics (CFD)

Pool fires are the most common of all process industry accidents. Pool fires often trigger explosions which may result in more fires, causing huge losses of life and property. Since both the risk and the frequency of occurrence of pool fires are high, it is necessary to model the risks associated with pool fires so as to correctly predict the behavior of such fires.Among the parameters which determine the overall structure of a pool fire, the most important is turbulence. It determines the extent of interaction of various parameters, including combustion, wind velocity, and entrainment of the ambient air. Of the various approaches capable of modeling the turbulence associated with pool fires, computational fluid dynamics (CFD) has emerged as the most preferred due to its ability to enable closer approximation of the underlying physical phenomena.A review of the state of the art reveals that although various turbulence models exist for the simulation of pool fire no single study has compared the performance of various turbulence models in modeling pool fires. To cover this knowledge-gap an attempt has been made to employ CFD in the assessment of pool fires and find the turbulence model which is able to simulate pool fires most faithfully. The performance of the standard k–? model, renormalization group (RNG) k–? model, realizable k–? model and standard k–ω model were studied for simulating the experiments conducted earlier by Chatris et al. (2001) and Casal (2013). The results reveal that the standard k–? model enabled the closest CFD simulation of the experimental results.     

Effect of turbulence spatial distribution on the deflagration index: Comparison between 20 L and 1 m3 vessels

In this work, the effect of spatial distribution and values of the turbulent kinetic energy on the pressure-time history and then on the explosion parameters (deflagration index and maximum pressure) was quantified in both the standard vessels (20 L and 1 m³).The turbulent kinetic energy maps were computed in both 20 L and 1 m³ vessels by means of CFD simulations with validated models. Starting from these maps, the turbulent flame propagation of cornstarch was calculated, by means of the software CHEMKIN. Then, the pressure-time history was evaluated and from this, the explosion parameters.Calculations were performed for three cases: not uniform turbulence level as computed from CFD simulations, uniform turbulence level and equal to the maximum value, uniform profile and equal to the minimum value. It was found that the cornstarch in the 20 L vessel get variable classes (St-1, St-2, St-3) with respect to the 1 m³ (St-1). However, simulations performed on increasing the ignition delay time, shown that the same results can be attained only using 260 ms as ignition delay time in the 20 L vessel.     

基于Factsage的金属硫化矿尘爆炸过程产物分析*

为研究金属硫化矿尘爆炸的过程与产物,预防与控制其粉尘发生爆炸,在20 L爆炸球中开展金属硫化矿尘爆炸试验,发现存在磁黄铁矿含量越高爆炸产物颜色越红的现象。为揭示这一现象的本质,应用XRD对爆炸产物进行分析,基于Factsage软件模拟计算黄铁矿矿尘与磁黄铁矿矿尘爆炸过程。结果表明:产物中含有的Fe2O3是主要致色原因;最终产物除固体Fe2O3外还存在气体SO2,SO3,反应中间过程主要固体产物为Fe2(SO4)3;在磁黄铁矿与黄铁矿含硫量相同或质量相同时,Fe2O3的生成量磁黄铁矿较黄铁矿多,且磁黄铁矿含硫量越多Fe2O3生成量越多,模拟计算结果与试验现象一致,研究结果可为分析金属硫化矿尘爆炸过程提供理论与数据支撑。     

抽出式通风煤巷掘进过程中粉尘浓度分布规律的数值模拟

根据气固两相流理论，针对矿井掘进工作面的特点，采用计算流体力学的离散相模型（DPM）对掘进工作面通风过程中粉尘浓度进行数值模拟，总结抽出式通风掘进巷道中粉尘浓度的沿程分布及变化规律。     

Hydrogen inhibition in wet dust removal systems by using L-aspartic: A feasible way of hydrogen explosion control measures

Alloy dust generated from automobile wheel hub grinding, after entering the wet dust collector, will react with water to produce hydrogen, thus exposing the entire ventilation and dust removal system to potential hydrogen explosion. In this paper, the inhibition mechanism and kinetic characteristics of different concentrations of L-Aspartic acid (L-Asp) on the reaction of Al_0.9Mg_0.1 alloy with water were studied with respect to adsorption morphology, chemical kinetic modeling and molecular dynamics (MD), using L-Asp as the environmentally-safe hydrogen inhibitor. The results show that within a given temperature interval, the hydrogen production rate of Al_0.9Mg_0.1 alloy dust decreases with increasing L-Asp concentration. When the L-Asp concentration exceeds 1.0g/L, the hydrogen evolution rate is almost zero. The calculated results of chemical kinetics agree with the Langmuir adsorption model, confirming that L-Asp is an ideal monolayer physical adsorption system on the surface of alloy particles. The FTIR and MD simulation results show that –NH₂ and –COOH groups in L-Asp molecules contribute greatly to the adsorption. The research results of this paper can help fundamentally avoid hydrogen generation in wet dust collectors and guarantee intrinsic safety.     

An optimal level of dust explosion risk management: Framework and application

The current research provides guidance on the prevention and mitigation of dust explosion using a Quantitative Risk Management Framework (QRMF). Using concepts drawn from previous studies, the framework consists of three main steps: (i) a new combined safety management protocol, (ii) the use of DESC (Dust Explosion Simulation Code) and FTA (Fault Tree Analysis) to assess explosion consequences and likelihood, respectively, and (iii) application of the hierarchy of controls (inherent, engineered and procedural safety). QRMF assessment of an industrial case study showed that the original process was at high risk. DESC simulations and Probit equations determined the destructive percentages. FTAs revealed high probabilities of explosion occurrence; in addition, detailed individual and societal risks calculations were made, before and after the framework was applied. Based on the hierarchy of controls technique, the framework showed significant risk reduction to the point where the residual risk was acceptable for the process.