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.     

Validation of a new formula for predicting the lower flammability limit of hybrid mixtures

A correlation of the lower flammability limit for hybrid mixtures was recently proposed by us. The experimental conditions including ignition energy and turbulence which play a primary role in a gas or dust explosion were at fixed values. The sensitivity of such experimental conditions to the accuracy of the proposed formula was not thoroughly discussed in the previous work. Therefore, this work studied the effect of varying the ignition energy and turbulence intensity to the formula proposed in our previous paper. For ignition energy effect, results from methane/niacin mixture demonstrated that the MEC and LFL will not be affected by changing ignition energy. There is no distinguishable difference among gas explosion index (K_G) and dust explosion index (K_St) derived from tests with every ignition energy (2.5 kJ, 5 kJ and 10 kJ) in a 36 L vessel. The proposed formula is independent of ignition energy. For turbulence effect, the proposed formula can have a good prediction of the explosion and non-explosion zone if the ignition delay time is within a certain range. The formula prediction is good as the ignition delay time increases up to 100 ms in this work. Propane/niacin and propane/cornstarch mixtures are also tested to validate the proposed formula. It has been confirmed that the proposed formula predicts the explosion and non-explosion zone boundary of such mixtures.     

The effect of vent size and concentration in vented gasoline-air explosions

Experimental data from vented explosion tests using gasoline-air mixtures with concentrations from 0.88 to 2.41% vol. are presented. A 2L vessel was used for the tests with vent sizes of 25 cm², 50 cm² and 100 cm². The tests were focused on the effect of gasoline vapor concentration and vent size on the pressure development and the flame behavior inside and outside the vessel. It was found that the inner flame propagation speed was mainly dependent on the initial concentration, while the maximum flame spreading distance was mainly influenced by the vent size. The external flame speed and duration could be influenced by the combination of the two properties. The internal pressure increases gradually with the flame propagated inside the vessel and decreased sharply when the vent failed. High-pressure durations containing pressure peaks were recorded by transducers in front of the vent and oscillations could be observed besides the vent. At any measure point, the maximum external pressures for A = 25 cm² or 50 cm² were significantly larger than those for A = 100 cm².     

Flame propagation in lycopodium/air mixtures below atmospheric pressure

Industrial processes are often operated at conditions deviating from atmospheric conditions. Safety relevant parameters normally used for hazard evaluation and classification of combustible dusts are only valid within a very narrow range of pressure, temperature and gas composition. The development of dust explosions and flame propagation under reduced pressure conditions is poorly investigated. Standard laboratory equipment like the 20 l Siwek chamber does not allow investigations at very low pressures. Therefore an experimental device was developed for the investigations on flame propagation and ignition under reduced pressure conditions. Flame propagation was analysed by a video analysis system the actual flame speed was measured by optical sensors. Experiments were carried out with lycopodium at dust concentrations of 100 g/m³, 200 g/m³ and 300 g/m³. It was found that both flame shapes and flame speeds were quite different from those obtained at atmospheric pressure. Effects like buoyancy of hot gases during ignition and flame propagation are less strong than at atmospheric conditions. For the investigated dust concentrations the flame reaches speeds that are nearly an order of a magnitude higher than at ambient conditions.     

The acceleration of flames in tube explosions with two obstacles as a function of the obstacle separation distance

The separation distance (or pitch) between two successive obstacles or rows of obstacles is an important parameter in the acceleration of flame propagation and increase in explosion severity. Whilst this is generally recognised, it has received little specific attention by investigators. In this work a vented cylindrical vessel 162 mm in diameter 4.5 m long was used to study the effect of separation distance of two low blockage (30%) obstacles. The set up was demonstrated to produce overpressure through the fast flame speeds generated (i.e. in a similar mechanism to vapour cloud explosions). A worst case separation distance was found to be 1.75 m which produced close to 3 bar overpressure and a flame speed of about 500 m/s. These values were of the order of twice the overpressure and flame speed with a double obstacle separated 2.75 m (83 characteristic obstacle length scales) apart. The profile of effects with separation distance was shown to agree with the cold flow turbulence profile determined in cold flows by other researchers. However, the present results showed that the maximum effect in explosions is experienced further downstream than the position of maximum turbulence determined in the cold flow studies. It is suggested that this may be due to the convection of the turbulence profile by the propagating flame. The present results would suggest that in many previous studies of repeated obstacles the separation distance investigated might not have included the worst case set up, and therefore existing explosion protection guidelines may not be derived from worst case scenarios.     

CFD simulations of dust dispersion in the 1 m3 explosion vessel

According to standard procedures, flammability and explosion parameters for dusts and dust mixtures are evaluated in 20 L and/or 1 m³ vessels, with equivalent results provided a correct ignition delay time (60 ms in the 20 L vessel; 600 ms in the 1 m³ vessel). In this work, CFD simulations of flow field and dust concentration distribution in the 1 m³ spherical vessel are performed, and the results compared to the data previously obtained for the 20 L. It has been found that in the 1 m³ vessel, the spatial distribution of the turbulent kinetic energy is lower and much more uniform. Concerning the dust distribution, as in the case of the 20 L, dust is mainly concentrated at the outer zones of the vortices generated inside the vessel. Furthermore, an incomplete feeding is attained, with most of the dust trapped in the perforated annular nozzle. Starting from the maps of dust concentration and turbulent kinetic energy, the deflagration index K_St is calculated in both vessels. In the conditions of the present work, the K_St is found to be 2.4 times higher in the 20 L than in the 1 m³ vessel.     

Autoignition temperature data and scaling for amide solvents

Autoignition temperature tests using the ASTM E659 test method have been conducted for N,N-dimethylacetamide (DMAC) and N,N-dimethylformamide (DMF) in test vessels with volumes of 0.5 l, 5 l, and 12 l. Tests were conducted at three different laboratories yielded good agreement (standard deviation with 5 °C) in all cases except for DMAC in the 0.5 l test vessel (standard deviation of 23 °C). Scaling correlations have been developed for the decrease of autoignition temperature with increasing volume and for increasing values of the vessel volume to surface area ratio. The variations for DMAC are steeper than the literature values for almost all other combustible liquids. Cool flames were observed for DMAC at temperatures as much as 44 °C below the autoignition temperature and for DMF at temperatures as much as 171 °C below its autoignition temperature. The DMF cool flame temperatures in the 5-l and 12-l test vessels are approximately equal to the DMF autoignition temperature in a closed 12-l test vessel. Gas samples taken after the cool flame and hot flame tests reveal the presence of high concentrations of diamines and dimethylamino acetonitrile, and small concentrations of many other partial decomposition/oxidation components.     

Configuration predictions of large liquefied petroleum gas (LPG) pool fires using CFD method

Liquefied petroleum gas (LPG) has potential pool fire risks due to its flammability. The configuration of pool fires plays a significant role when applying the solid flame model or point source model to assess the risks from heat radiation. However, no existing correlations can precisely predict the configuration of large LPG (100% propane) pool fires. To enhance the fundamental understanding on how pool diameter and wind velocity can influence the configuration of large LPG pool fires, an experimentally validated Computational Fluid Dynamics (CFD) model is employed to simulate fires using different burning rate models. Fire temperature profiles, flame heights, and flame tilts predicted by the CFD model were compared with empirical models and experimental data. Accordingly, new correlations for flame height and flame tilt as functions of pool diameter D and wind velocity u_w have been developed. The comparisons demonstrate that the new correlations have the best overall accuracy in the prediction of flame height and tilt for large LPG pool fires under different conditions (10 m ≤ D ≤ 20 m, 0 ≤ u_w ≤ 3 m·s⁻¹).     

Efficient simulation of flame acceleration and deflagration-to-detonation transition in smooth pipes

Evaluation of accident scenarios including flame acceleration and deflagration-to-detonation transition (DDT) in chemical plant piping systems increases the need for an efficient numerical simulation tool capable of dealing with this phenomenon. In this work, a hybrid pressure-density-based solver including deflagrative flame propagation as well as detonation propagation is presented. The initial incompressible acceleration stage is covered by the pressure-based solver until the flame velocity reaches the fast flame regime and transition to the density-based solver is done. The deflagration source term is formulated in terms of a turbulent flame speed closure model incorporating various physical effects crucial for flame acceleration at low turbulence conditions (Katzy and Sattelmayer, 2018). Modelling of the detonation source term is based on a quadratic heat release function (Hasslberger, 2017). The presented numerical approach is validated in terms of DDT locations and pressure data from Schildberg (2015) as well as recently completed flame tip position measurements. For this purpose, H₂/O₂/N₂ mixtures ranging from 25.6 vol-% H₂ to 29.56 vol-% H₂ in two different pipe geometries are considered. The focus of the current work is on predicting the DDT location correctly and good agreement is observed for the investigated cases.     

Effect of scale on freely propagating flames in aluminum dust clouds

The majority of experimental tests done on combustible dusts are performed in constant volume vessels that have limited or no optical access. Over the years, McGill University has been developing alternative experimental techniques based on direct observation of dust flames, yielding reliable fundamental parameters such as flame burning velocity, temperature and structure. The present work describes two new experimental set-ups allowing direct observation of isobaric and freely propagating dust flames at two sufficiently different scales to test the influence of scale on dust flame phenomena. In the laboratory-scale experiments, a few grams of aluminum powder are dispersed in transparent, 30 cm diameter latex balloons that allow for full visualization of the spherical flame propagation. In the field experiments, about 1 kg of aluminum powder is dispersed by a short pulse of air, forming a conical dust cloud with a total volume of about 5 m³. High-speed digital imaging is used to record the particle dispersal and flame propagation in both configurations. In the small-scale laboratory tests, the measured flame speed is found to be about 2.0 ± 0.2 m/s in fuel-rich aluminium clouds. The burning velocity, calculated by dividing the measured flame speed by the expansion factor deduced from thermodynamic equilibrium calculations, correlates well with the previously measured burning velocity of about 22–24 cm/s from Bunsen dust flames. Flame speeds observed in field experiments with large-scale clouds, however, are found to be much higher, in the range of 12 ± 2 m/s. Estimations are presented that show that the presumably greater role of radiative heat transfer in larger-scale aluminium flames is insufficient to explain the six-fold increase in flame speed. The role of residual large-eddy turbulence, as well as the frozen-turbulence effect leading to large-scale dust concentration fluctuations that cause flame folding, are discussed as two possible sources for the greater flame speed.     

Investigations of flame front propagation between interconnected process vessels. Development of a new flame front propagation time prediction model

For the case where a dust or gas explosion can occur in a connected process vessel, it would be useful, for the purpose of designing protection measures and also for assessing the existing protection measures such as the correct placement, to have a tool to estimate the time for flame front propagation along the connecting pipe. Measurements of data from large-scale explosion tests in industrially relevant process vessels are reported. To determine the flame front propagation time, either a 1 m³ or a 4.25 m³ primary process vessel was connected via a pipe to a mechanically or pneumatically fed 9.4 m³ secondary silo. The explosion propagation started after ignition of a maize starch/air mixture in the primary vessel. No additional dust was present along the connecting pipe. Systematic investigations of the explosion data have shown a relationship between the flame front propagating time and the reduced explosion over-pressure of the primary explosion vessel for both vessel volumes. Furthermore, it was possible to validate this theory by using explosion data from previous investigations. Using the data, a flame front propagation time prediction model was developed which is applicable for:

• •gas and dust explosions up to a K value of 100 and 200 bar m s⁻¹, respectively, and a maximum reduced explosion over-pressure of up to 7 bar;
• •explosion vessel volumes of 0.5, 1, 4.25 and 9.4 m³, independent of whether they are closed or vented;
• •connecting pipes of pneumatic systems with diameters of 100–200 mm and an air velocity up to 30 m s⁻¹;
• •open ended pipes and pipes of interconnected vessels with a diameter equal to or greater than 100 mm;
• •lengths of connecting pipe of at least 2.5–7 m.

     

Experiments on the influence of pre-ignition turbulence on vented gas and dust explosions

Experiments were performed on the influence of pre-ignition turbulence on the course of vented gas and dust explosions. A vertical cylindrical explosion chamber of approximately 100 l volume and a length-to-diameter ratio (l/d) of 4.7 consisting of a steel bottom segment and three glass sections connected by steel flanges was used to perform the experiments. Sixteen small fans evenly distributed within the chamber produced turbulent fluctuations from 0 to 0.45 m/s. A Laser-Doppler-anemometer (LDA) was used to measure the flow and turbulence fields. During the experiments the pressure and in the case of dust explosions the dust concentration were measured. In addition, the flame propagation was observed by a high-speed video camera. A propane/nitrogen/oxygen mixture was used for the gas explosion experiments, while the dust explosions were produced by a cornstarch/air mixture.It turned out that the reduced explosion pressure increased with increasing turbulence intensity. This effect was most pronounced for small vents with low activation pressures, e.g. for bursting disks made from polyethylene foil. In this case, the overpressure at an initial turbulence of 0.45 m/s was twice that for zero initial turbulence.     

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.     

Propane-air mixture deflagration hazard with different ignition positions and restraints in large-scale venting chamber

In order to better assess the hazards of explosion accidents, propane-air mixture deflagrations were conducted in a large-scale straight rectangular chamber (with a cross-section of 1.5 m × 1.5 m, length of 10 m, and total volume of 22.5 m³). The effect of initial volume, ignition position, and initial restraints on the explosion characteristics of the propane-air mixtures was investigated. The explosion overpressure, flame propagation, and flame speed were obtained and the computational fluid dynamics (CFD) software was used to simulate the flame-propagation process and field flow for auxiliary analysis. The hazards of large-scale propagation explosion under weak and strong constraints were evaluated and the different phases of flame propagation under weak and strong constraints were discriminated. Results indicate that the hazards caused by propane deflagration under weak constraint are mainly caused by flame spread. And the maximum overpressure under strong constraint appeared at the front part of the chamber under the large-scale condition, which is consistent with the previous small-scale test. Moreover, the simulations of flame structures under weak and strong constraint are in good agreement with experimental results, which furthers the understanding of large-scale propane deflagration under different initial conditions in large-scale spaces and provides basic data for three-dimensional CFD model improvement.     

Gas flame acceleration in long ducts

In many practical situations, a flame may propagate along a pipe, accelerate and perhaps transform into a devastating detonation. This phenomenology has been known, more or less qualitatively, for a long time and mitigation techniques were proposed to try and avoid this occurrence (flame arresters, vents,...). A number of parameters need to be known and in particular the “distance to detonation” and more generally the flame acceleration characteristic scales. Very often, the ratio between the detonation run-up distance and the pipe diameter is used without any strong justification other that using a non-dimensional parameter (L/D). In this paper, novel experimental evidence is presented on the basis of relatively large scale experiments using 10 cm and 25 cm inner diameter duct with a length between 7 and 40 m. Homogeneous C₂H₄-air, CH₄-air, C₃H₈-air and H₂-air mixtures were used and different ignition sources. The interpretation suggests that the self-acceleration mechanism of the flame may be much better represented by flame instabilities than by turbulence build-up. One consequence would be that the maximum flame velocity and, following, the maximum explosion overpressure, would be rather linked with the run-up distance than with the L/D ratio.     

Effect of concentration and ignition position on vented methane–air explosions

Experiments were conducted in a 1 m³ vessel with a top vent to investigate the effect of methane concentration and ignition position on pressure buildup and flame behavior. Three pressure peaks (p₁, p₂, and P_ext) and two types of pressure oscillations (Helmholtz and acoustic oscillations) were observed. The rupture of vent cover results in p₁ that is insensitive to methane concentration and ignition position. Owing to the interaction between acoustic wave and the flame, p₂ forms in the central and top ignition explosions when the methane–air mixture is near–stoichiometric. When the methane–air mixture is centrally ignited, p₂ first increases and then decreases with an increase in the methane concentration. The external explosion-induced P_ext is observed only in the bottom ignition explosions with an amplitude of several kilopascals. Under the current experimental conditions, flame–acoustic interaction leads to the most serious explosions in central ignition tests. Methane concentration and ignition position have little effect on the frequency of Helmholtz and acoustic oscillations; however, the Helmholtz oscillation lasts longer and first decreases and then increases as the methane concentration increases for top ignition cases. The ignition position significantly affects the Taylor instability of the flame front resulting from the Helmholtz oscillation.     

Structural response for vented methane–air deflagrations: Effects of volumetric blockage ratio

To investigate the effects of cylinders placed parallel to the venting direction on the structural response of the vessel walls to an explosion, 25 batches of vented explosion tests were conducted in a 1 m³ rectangular vessel. Two types of structural response with different amplitudes and frequency distributions were observed and evaluated by comparing the vibration data with both the pressure data and high-speed videos. A low-amplitude structural response of approximately 150–250 m/s², which increased slightly as VBR increased, was triggered by a combination of the initial flame propagation, external explosion, Helmholtz oscillations, and the Taylor instability. A high-amplitude structural response of approximately 9500 m/s² was also observed, which decreased sharply as VBR increased. Additionally, the high amplitude response was never observed when more than two cylinders were present in the vessel. The high amplitude response was triggered due to the coupling between the acoustic wave, the flame, and the resonance of the vessel. The presence of obstacles did not increase the severity of the structural responses under the current experimental conditions. To the contrary, the presence of obstacles in the container attenuated or even inhibited the high-amplitude vibration of the container caused by the explosion.     

Direction of pressure wave propagation during combustion-induced rapid phase transition

In the work presented in this paper, the propagation direction of the pressure waves generated during combustion-induced Rapid Phase Transition (cRPT) was investigated. To this end, explosion tests were performed for CH₄/O₂/N₂/CO₂ mixtures in a tubular reactor. Ignition was provided at the top or at the bottom of the vessel. Pressure time histories were recorded by two transducers positioned one at the top and the other one at the bottom.Results have shown that the preferential direction for the pressure waves is that of the flame propagation. When the cRPT phenomenon is weak, an over-adiabatic pressure peak (of around 10–20 bar) can be measured by only one transducer and, in particular, by the transducer far away from the ignition point. Conversely, when the cRPT phenomenon becomes severe, over-adiabatic peaks (as high as 250 bar) can be detected even by the other transducer. Such peaks are the result of separate cRPT events that occur very close to the transducers and, thus, are not damped along the vessel length. In spite of the fact that the upward flame propagation is faster, the cRPT phenomenon is more severe in the case of downward flame propagation.     

Large-scale dust explosions in vessel-pipe systems

This study investigates dust explosions in vessel-pipe systems to develop a better understanding of dust flame propagation between interconnected vessels and implications for the proper application of explosion isolation systems. Cornstarch dust explosions were conducted in a large-scale setup consisting of a vented 8-m³ vessel and an attached pipe with a diameter of 0.4 m and a length of 9.8 m. The ignition location and effective dust reactivity were varied between experiments. The experimental results are compared against previous experiments with initially quiescent propane-air mixtures, demonstrating a significantly higher reactivity of the dust explosions due to elevated initial turbulence, leading to higher peak pressures and faster flame propagation. In addition, a physics-based model developed previously to predict gas explosion dynamics in vessel-pipe systems was extended for dust combustion. The model successfully predicts the pressure transients and flame progress recorded in the experiments and captures the effects of ignition location and effective dust reactivity.     

Experimental study on unconfined methane explosion: Explosion characteristics and overpressure prediction method

Explosion accidents have become the main threat for the high-efficiency use of cleaner gas energy sources, such as natural gas. During an explosion, obstacle causing flame acceleration is the main reason for the increase of the explosion overpressure, which still remains to be fully understood. In this research, field experiments were conducted in a 1 m³ cubic frame apparatus to investigate the effect of built-in obstacles on unconfined methane explosion. Cage-like obstacles were constructed using square steel rods with different cross section size. The results demonstrated that the flame could get accelerated due to the hydrodynamic instability and obstacle-induced turbulence, which enhanced the explosion overpressure. In the near field, the overpressure wave travelled slower and the maximum overpressure could almost keep constant. Reducing the cross section size, or increasing the obstacle height or the obstacle number per layer could determine the rise of the maximum overpressure, the maximum pressure rising rate and the overpressure impulse. For uniformly constructed obstacles, self-similar theory was chosen to measure the influence of the hydrodynamic instability, and a parameter β was adopted to measure the flame acceleration caused by obstacle-induced turbulence, the value of which was 2 in this research. Based on the acoustic theory, an overpressure prediction model was proposed and the predicted results agreed with the measured values better than previous models, such as TNT equivalency model and TNO multi-energy model.