Thermo-kinetic modelling of dust explosions

The guidelines for protection and mitigation against hazard coming from dust explosion require the knowledge and then the evaluation either experimentally or theoretically of the thermo-kinetic parameters (i.e. K_St, P_max). We developed a numerical tool for the evaluation of the thermo-kinetic parameters of dust explosion. This model is based on the simulations of the combustion reaction by means of a detailed reaction mechanism assuming that the pyrolysis/devolatilization step is very fast and then gas combustion is controlling dust explosion. The model allows then the determination of the most conservative values of K_St, S_l, P_max. In the present paper we calculated the deflagration index and the laminar burning velocity for dusts utilized in various process industries (i.e. cornstarch, polyethylene, cellulose) as function of dust concentration. The obtained data were successfully compared with the available experimental results.     

Simulation of dust explosions in complex geometries with experimental input from standardized tests

The present paper describes the development of a new CFD-code (DESC) for the assessment of accidental hazards arising from dust explosions in complex geometries. The approach followed entails the estimation of the laminar burning velocity of dust clouds from standardized laboratory-scale tests, and its subsequent use as input to the combustion model incorporated in DESC. The methodology used to obtain the laminar burning velocities is demonstrated by igniting turbulent propane-air mixtures to deflagration in a standard 20-litre USBM-vessel, and extracting the laminar burning velocity from the pressure–time curves; the results are compared with literature data. Laminar burning velocities for clouds of maize starch dust in air were estimated following the same procedure, and the resulting empirical model was used to simulate dust explosions in a 236-m³ silo.     

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.     

Experimental study of the effect of nitrogen addition on gas explosion

Flame propagation and combustion characteristics of methane/air mixed gas in gas explosion were studied in a constant volume combustion bomb. Stretched flame propagation velocity, unstretched laminar flame propagation velocity, unstretched laminar combustion velocity and Markstein length were obtained at various ratios of nitrogen to gas mixture. Combustion stability at various ratios of nitrogen to gas mixture was analyzed by analyzing the pictures of flame propagation. Furthermore, the effect of initial pressure on the flame propagation and combustion characteristics of methane/air mixed gas in gas explosion was analyzed. The results show that the unstretched laminar flame propagation velocity, the unstretched laminar combustion velocity, Markstein length, flame stability, and the maximum combustion pressure decrease distinctly with the increase of nitrogen fraction in the gas mixture. At the same ratios of nitrogen to gas mixture, Markstein length, unstretched laminar flame propagation velocity and unstretched laminar combustion velocity decrease and the maximum combustion pressure increase with the increase of initial pressure of the gas mixture. When nitrogen fraction in the gas mixture is over 20%, the flame will be unstable and is easy to exterminate.     

矿井瓦斯煤尘爆炸传播数值模拟研究

基于连续相、燃烧、颗粒相数理方程建立瓦斯煤尘爆炸传播数理模型,并应用连续相、颗粒相计算方法,依据大型巷道瓦斯爆炸、瓦斯煤尘爆炸传播实验数据,借助普遍应用的流场模拟平台,成功开发瓦斯、煤尘爆炸数值模拟系统。该系统可有效地模拟煤矿瓦斯、煤尘的爆炸事故过程,对瓦斯爆炸的爆燃转爆轰、煤尘是否参与爆炸、爆炸冲击传播速度、衰减规律以及爆炸灾害的波及范围都能进行较准确的模拟。     

Temperature profile across the combustion zone propagating through an iron particle cloud

Knowledge of the mechanism of combustion zone propagation during dust explosion is of great importance to prevent damage caused by accidental dust explosions. In this study, the temperature profile across the combustion zone propagating through an iron particle cloud is measured experimentally by a thermocouple to elucidate the propagation mechanism. The measured temperature starts to increase slowly at a position about 5 mm ahead of the leading edge of the combustion zone, increases quickly at a position about 3 mm ahead of the leading edge, reaches a maximum value near the end of the combustion zone, and then decreases. As the iron particle concentration increases, the maximum temperature increases at lower concentration, takes a maximum value, and then decreases at higher concentration. The relation between the propagation velocity of the combustion zone and the maximum temperature is also examined. It is found that the propagation velocity has a linear relationship with the maximum temperature. This result suggests that the conductive heat transfer is dominant in the propagation process of the combustion zone through an iron particle cloud.     

硬脂酸粉尘爆炸特性试验研究

为研究硬脂酸粉尘的爆炸特性,采用20 L球型爆炸仪对4个粒径范围的硬脂酸粉尘进行粉尘爆炸试验研究。结果表明:一定浓度范围内增大粉尘浓度能够提升硬脂酸粉尘的爆炸能量和燃烧速率。增大粉尘浓度,爆炸猛烈度先增强后减弱;减小粉尘粒径,能增强爆炸猛烈度和敏感度。粒径小于58 μm粉尘的爆炸猛烈度和敏感度最大,浓度500 g/m3时,该粉尘有最大爆炸压力1.12 MPa和最大升压速率142.00 MPa/s。     

A comparative study between two ignition sources: Electric igniter versus pyrotechnic igniter

The risk assessment of combustible explosive dust is based on the determination of the probability of dust dispersion, the identification of potential ignition sources and the evaluation of explosion severity. It is achieved in most of cases with the two main experimental normalized devices such as the Hartmann tube (spark ignition) and the 20 L spherical bomb (with two 5 kJ pyrotechnic ignitors).Ignition energy of the 5 kJ ignitor is well calibrated and generates a reproducible ignition. But, on the other hand, this ignition is not punctual and the over pressure produced is nearly 2 bar. Moreover, the pyrotechnic igniter accelerates the combustion with multi ignition points in a large volume and that disturbs the flame propagation. In this way, this ignition source does not allow to analyze the combustion products because the composition of the pyrotechnic igniter was found in the combustion products.This paper deals with the comparison of two ignition sources in the 20 L spherical bomb. Different explosive dusts of great industrial interest are studied with electrical and pyrotechnic ignitors, in order to understand, first, the influence of each type of igniter on the explosion behaviour and then to evaluate the possibility of establishing a correspondence between parameters obtained with these two ignition sources.Severity parameters of nicotinic acid, aluminium powder and titanium alloy were measured by using the two types of ignition system in our 20 L spherical bomb equipped with the Kühner dihedral injector. The explosion overpressure $\mathrm{P}$ and the rate of pressure rise $\left(\raisebox{1ex}{$dP$}\!\left/ \!\raisebox{-1ex}{$dt$}\right.\right)$ were measured in a large range of concentration allowing to propose correlations between electrical and pyrotechnic ignition for each parameter and each type of powder. These correlations aim to link the tests used with two different collections of experimental parameters for the same dust. The relevance of these correlations will be discussed.     

Study of dust cloud behaviour in the modified Hartmann tube using the image subtraction method (ISM)

A dispersion of fine particles in the air is needed for a dust explosion to occur since an explosion is the fast combustion of particles in the air. When particles are poorly dispersed, agglomerated, or their concentration is low, the combustion velocity decreases, and deflagration would not occur. The combustion rate is strictly related to dust concentration. Therefore, the maximum explosion pressure rise occurs at dust concentration close to stoichiometric. Conversely, Minimum Explosion Concentration (MEC) is the lower limit at which self-sustained combustion and a pressure rise are possible. Dust explosion tests are designed to reproduce the dispersion and generation of dust clouds in industrial ambiences by using dispersion devices activated by pressurised air pulses. The resulting dust cloud, which has a marked transient character, is considered representative of real clouds by current standards. Over time, several studies have been carried out to optimise these devices (e.g. to reduce the inhomogeneity of the cloud in the 20 L sphere). The Minimum Ignition Energy (MIE) of dust is measured using the Mike3 modified Hartmann tube, where the ignition attempt is made 60–180 ms after dust dispersion regardless of dust characteristics.This work investigates the dust clouds’ actual behaviour inside the modified Hartmann tube before ignition using high-velocity video movies and a new image post-treatment method called Image Subtraction Method (ISM). Movies are recorded with high-speed cameras at a framerate of 2000 fps and elaborated with an on-purpose developed LabVIEW® code. Concentration (mass per volume) and dispersion pressure are varied to evaluate their effect on dust clouds. Maise starch, iron powder and silica powder are chosen to investigate the effect of particle density and size on the cloud structure and turbulence. This approach will help to investigate the structure of the dust cloud, the shape and size of the particle lumps and the change in dust concentration over time. In addition, information on the actual concentration and cloud turbulence at the ignition location and delay time were obtained, which may help identify the local turbulence scale and widen the characterisation of the cloud generated in the Hartmann tube.     

Determination of turbulent burning velocities of dust air mixtures with the open tube method

Correlating turbulent burning velocity to turbulence intensity and basic flame parameters-like laminar burning velocity for dust air mixtures is not only a scientific challenge but also of practical importance for the modelling of dust flame propagation in industrial facilities and choice of adequate safety strategy. The open tube method has been implemented to measure laminar and turbulent burning velocities at laboratory scale for turbulence intensities in the range of a few m/s. Special care has been given to the experimental technique so that a direct access to the desired parameters was possible minimising interpretation difficulties. In particular, the flame is propagating freely, the flame velocity is directly accessible by visualisation and the turbulence intensity is measured at the flame front during flame propagation with special aerodynamic probes. In the present paper, those achievements are briefly recalled. In addition, a complete set of experiments for diametrically opposed dusts, starch and aluminium, has been performed and is presented hereafter. The experimental data, measured for potato dust air mixtures seem to be in accordance with the Bray Gülder model in the range of 1.5 m/s<u′<3.5 m/s. For a further confirmation, the measurement range has been extended to lower levels of turbulence of u′<1.5 m/s. This could be achieved by changing the mode of preparation of the dust air mixture. In former tests, the particles have been injected into the tube from a pressurised dust reservoir; for the lower turbulence range, the particles have been inserted into the tube from above by means of a sieve–riddler system, and the turbulence generated from the pressurised gas reservoir as before. For higher levels of turbulence, aluminium air mixtures have been investigated using the particle injection mode with pressurised dust reservoir. Due to high burning rates much higher flame speeds than for potato dusts of up to 23 m/s have been obtained.     

Effect of initial temperature on explosion characteristics of 2, 3, 3, 3–Tetrafluoropropene

The flammability of refrigerants is a major cause of refrigerant explosion incidents. Studying the explosion characteristics of refrigerants at different initial temperatures can provide significant benefits for solving the safety problems of refrigerants under actual working conditions. This paper studied the effects of the initial temperature and refrigerant concentration on the explosion characteristics of refrigerant 2, 3, 3, 3-tetrafluoropropene (R1234yf) at 0.1 MPa. The curves of explosion characteristics with different initial temperature revealed the same variation trend ranged from 25 °C to 115 °C. Specifically, as the refrigerant concentration was raised, the peak overpressure, the maximum rate of pressure rise, and laminar burning velocity increased initially and decreased afterwards, along with maximum values at the refrigerant concentration of 7.6%. When the refrigerant concentration was 7.6%, the peak overpressure declined exponentially with the initial temperature rise, while the maximum rate of pressure rise increased linearly. The laminar burning velocity calculated from the spherical expansion method indicated that the flame propagation was gradually accelerated by the increase of initial temperature, which coincided with the change of the maximum rate of pressure rise. Meanwhile, experiments and CHEMKIN simulation results demonstrated the effects of elevated temperature from 20 °C to 50 °C on the explosion limits of R1234yf. The lower explosion limit reduced and the upper explosion limit increased with rising initial temperature. In general, R1234yf exhibited moderate combustion and lower explosion risk, compared with traditional refrigerants.     

Study of the severity of hybrid mixture explosions and comparison to pure dust-air and vapour-air explosions

The explosion behaviour of heterogeneous/homogeneous fuel-air (hybrid) mixtures is here analysed and compared to the explosion features of heterogeneous fuel-air and homogeneous fuel-air mixtures separately.Experiments are performed to measure the pressure history, deflagration index and flammability limits of nicotinic acid/acetone-air mixtures in a standard 20 L Siwek bomb adapted to vapour-air mixtures. Literature data are also used for comparison.The explosion tests performed on gas-air mixtures in the same conditions as explosion tests of dust-air mixtures, show that the increase in explosion severity of dust/gas-air mixtures has to be addressed to the role of initial level of turbulence prior to ignition.At a fixed value of the equivalence ratio, by substituting the dust to the flammable gas in a dust/gas-air mixture the explosion severity decreases. Furthermore, the most severe conditions of dust-gas/air mixtures is found during explosion of gas-air mixture at stoichiometric concentration.     

On the determination of the laminar burning velocity from closed vessel gas explosions

A methodology to determine the laminar burning velocity from closed vessel gas explosions is explored. Unlike other methods which have been used to measure burning velocities from closed vessel explosions, this approach belongs to the category which does not involve observation of a rapidly moving flame front. Only the pressure–time curve is required as experimental input. To verify the methodology, initially quiescent methane–air mixtures were ignited in a 20-l explosion sphere and the equivalence ratio was varied from 0.67 to 1.36. The behavior of the pressure in the vessel was measured as a function of time and two integral balance models, namely, the thin-flame and the three-zone model, were fitted to determine the laminar burning velocity. Data on the laminar burning velocity as a function of equivalence ratio, pressure and temperature, measured by a variety of other methods have been collected from the literature to enable a comparison. Empirical correlations for the effect of pressure and temperature on the laminar burning velocity have been reviewed and two were selected to be used in conjunction with the thin-flame model. For the three-zone model, a set of coupled correlations has been derived to describe the effect of pressure and temperature on the laminar burning velocity and the laminar flame thickness. Our laminar burning velocities are seen to fall within the band of data from the period 1953–2003. A comparison with recent data from the period 1994–2003 shows that our results are 5–10% higher than the laminar burning velocities which are currently believed to be the correct ones for methane–air mixtures. Based on this observation it is concluded that the methodology described in this work should only be used under circumstances where more accurate methods can not be applied.     

Effects of ignitors and turbulence on dust explosions

The aim of this work is in an attempt to increase the understanding of the acting behaviour of pyrotechnic ignitors and their effects on confined dust explosions. Flame visualization has shown that pyrotechnic ignitors can initiate an explosion by instantaneous jet-like volumetric and/or multipoint ignition. Hence, the rate of pressure rise and also the apparent burning velocity will be increased to some extent, depending on the ignitor energy and the reactivity of the mixtures. The ignitor effect is more important for the early stages of flame propagation and would be more significant in small explosion chambers. Thus, for dust explosion tests with various purposes, use of pyrotechnic ignitors should be made carefully, and the ignitor effect must be accounted for in the data interpretation. Turbulence induced by dust dispersion is a dominant factor in affecting dust explosions. At different ignition delays, however, the turbulence influence will be coupled with that of ignitors. This complicates further the interpretation of explosion data measured under turbulent conditions.     

Laminar burning velocities of hydrogen–air mixtures from closed vessel gas explosions

The laminar burning velocity of hydrogen–air mixtures was determined from pressure variations in a windowless explosion vessel. Initially, quiescent hydrogen–air mixtures of an equivalence ratio of 0.5–3.0 were ignited to deflagration in a 169 ml cylindrical vessel at initial conditions of 1 bar and 293 K. The behavior of the pressure was measured as a function of time and this information was subsequently exploited by fitting an integral balance model to it. The resulting laminar burning velocities are seen to fall within the band of experimental data reported by previous researchers and to be close to values computed with a detailed kinetics model. With mixtures of an equivalence ratio larger than 0.75, it was observed that more advanced methods that take flame stretch effects into account have no significant advantage over the methodology followed in the present work. At an equivalence ratio of less than 0.75, the laminar burning velocity obtained by the latter was found to be higher than that produced by the former, but at the same time close enough to the unstretched laminar burning velocity to be considered as an acceptable conservative estimate for purposes related to fire and explosion safety. It was furthermore observed that the experimental pressure–time curves of deflagrating hydrogen–air mixtures contained pressure oscillations of a magnitude in the order of 0.25 bar. This phenomenon is explained by considering the velocity of the burnt mixture induced by the expansion of combusting fluid layers adjacent to the wall.     

Fast and tiny: A model for the flame propagation of nanopowders

To avoid the influence of external parameters, such as the vessel volume or the initial turbulence, the explosion severity should be determined from intrinsic properties of the fuel-air mixture. Therefore, the flame propagation of gaseous mixtures is often studied in order to estimate their laminar burning velocity, which is both independent of external factors and a useful input for CFD simulation. Experimentally, this parameter is difficult to evaluate when it comes to dust explosion, due to the inherent turbulence during the dispersion of the cloud. However, the low inertia of nanoparticles allows performing tests at very low turbulence without sedimentation. Knowledge on flame propagation concerning nanoparticles may then be modelled and, under certain conditions, extrapolated to microparticles, for which an experimental measurement is a delicate task. This work focuses on a nanocellulose with primary fiber dimensions of 3 nm width and 70 nm length. A one-dimensional model was developed to estimate the flame velocity of a nanocellulose explosion, based on an existing model already validated for hybrid mixtures of gas and carbonaceous nanopowders similar to soot. Assuming the fast devolatilization of organic nanopowders, the chemical reactions considered are limited to the combustion of the pyrolysis gases. The finite volume method was used to solve the mass and energy balances equations and mass reactions rates constituting the numerical system. Finally, the radiative heat transfer was also considered, highlighting the influence of the total surface area of the particles on the thermal radiation. Flame velocities of nanocellulose from 17.5 to 20.8 cm/s were obtained numerically depending on the radiative heat transfer, which proves a good agreement with the values around 21 cm/s measured experimentally by flame visualization and allows the validation of the model for nanoparticles.     

Effect of particle morphology on metal dust deflagration sensitivity and severity

Combustible dust explosions continue to present a significant threat toward industries processing, storing, or pneumatically conveying metal dust hazards. Through recent years, investigations have observed the influence of particle size, polydispersity, and chemical composition on dust explosion sensitivity and severity. However, studies characterizing the effect of particle shape (or morphology) on metal dust explosibility are limited and merit further consideration. In this work, high-purity aluminum dust samples of three unique particle morphologies were examined (spherical granular, irregular granular, and dry flake). To maintain consistency in results obtained, all samples were procured with similar particle size distribution and polydispersity, as verified by laser diffraction particle size analysis. Scanning electron microscopy (SEM) imaging and Brunauer-Emmett-Teller (BET) experiments were executed to confirm supplier claims on morphology and to quantify the effective surface area associated with each sample, respectively. Investigations performed in a Kühner MIKE3 minimum ignition energy apparatus and a Siwek 20 L sphere combustion chamber resulted in the direct characterization of explosion sensitivity and severity, respectively, as a function of suspended fuel concentration and variable particle morphology. Recommendations to standard risk/hazard analysis procedures and to existing design guidance for the mitigation of deflagrations that originate from ignition of distinctively processed metal dust fuels have been provided.     

超细聚苯乙烯微球粉体的燃爆特性研究

为研究超细聚苯乙烯微球粉体的燃爆特性,通过粉尘层最低着火温度测试装置、MIE-D1.2最小点火能测试装置、20 L球形爆炸测试装置,对其最低着火温度、最大爆炸压力、最小点火能量(MIE)等爆炸特性参数进行测定,探讨了加热温度、点火延滞时间、粉尘质量浓度、粉尘粒径对粉体燃爆特性的影响。结果表明:超细聚苯乙烯微球粉尘层在350℃左右时会发生无焰燃烧,且加热温度越高,粉体粒径越小,粉尘层发生着火时所需的时间越短;当粉体质量浓度为250 g/m3时,最大爆炸压力达到0.65 MPa,质量浓度为500 g/m3时,最大爆炸压力的上升速率达90 MPa/s以上;随点火延滞时间增加,最小点火能表现出先缓慢减小再急剧增大的规律;随粉尘质量浓度增加,最小点火能逐渐降低,当粉尘质量浓度超过500g/m3后逐渐趋于稳定。     

不同粒径镁铝合金粉尘爆炸与抑爆特性研究

为了研究镁铝合金粉爆炸危险特性,利用20L球形爆炸容器进行测试,结果表明:180目 (80 μm)、 120目(125 μm) 和60目(250 μm)3种粒径下的金属粉尘爆炸下限浓度分别为45 g/m3,55 g/m3和95 g/m3。相同浓度下最大爆炸压力随粒径增大的而减小。以碳化硅和石墨为代表的研究中,60目,120目和180目的镁铝合金粉以10%的浓度梯度加入碳化硅浓度分别至50%,70%和80%,石墨浓度至30%,50%和60%时,镁铝合金粉不会发生爆炸。表明碳化硅及石墨等惰性粉尘都能对粉尘爆炸有抑制作用,其中石墨对镁铝合金粉的抑爆作用明显优于碳化硅。     

2,5-二甲基呋喃-空气燃烧过程研究

利用球形发展火焰研究了常温常压下不同当量比,不同相态时2,5-二甲基呋喃-空气的层流燃烧速度和马克斯坦长度,分析了火焰拉伸对火焰传播速度的影响。研究结果表明:随着当量比的增加,2,5-二甲基呋喃-空气混合气的马克斯坦长度减少,火焰的稳定性减弱。并且分别计算出当量比为1.25和1.5的层流燃烧速度,分别为:1.189m/s,1.135m/s.。对于同一当量比1.5的情况下,不同相态的2,5-二甲基呋喃-空气混合物,在相同时刻的气液两相混合物的火焰半径已经拉伸火焰传播速度远远大于纯气相的混合物。