Next Article in Journal
Characterizing Spatial and Temporal Variability of Lightning Activity Associated with Wildfire over Tasmania, Australia
Previous Article in Journal
Experimental Study on Fire Sources for Full-Scale Fire Testing of Simple Sprinkler Systems Installed in Multiplexes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of Explosions in Cylindrical Vented Enclosures: Validation of a Computational Model by Experiments

by
Samuel Ogunfuye
1,
Hayri Sezer
2,
Furkan Kodakoglu
1,3,
Hamed Farmahini Farahani
3,
Ali S. Rangwala
3 and
V’yacheslav Akkerman
1,*
1
Center for Innovation in Gas Research and Utilization (CIGRU), Department of Mechanical and Aerospace Engineering, West Virginia University, 1306 Evansdale Drive, Morgantown, WV 26506, USA
2
School of Engineering and Technology, Western Carolina University, 1 University Way, Cullowhee, NC 28723, USA
3
Department of Fire Protection Engineering, Worcester Polytechnic Institute, 100 Institute Road, Worcester, MA 01609, USA
*
Author to whom correspondence should be addressed.
Submission received: 2 December 2020 / Revised: 8 February 2021 / Accepted: 8 February 2021 / Published: 15 February 2021

Abstract

:
Recent explosions with devastating consequences have re-emphasized the relevance of fire safety and explosion research. From earlier works, the severity of the explosion has been said to depend on various factors such as the ignition location, type of a combustible mixture, enclosure configuration, and equivalence ratio. Explosion venting has been proposed as a safety measure in curbing explosion impact, and the design of safety vent requires a deep understanding of the explosion phenomenon. To address this, the Explosion Venting Analyzer (EVA)—a mathematical model predicting the maximum overpressure and characterizing the explosion in an enclosure—has been recently developed and coded (Process Saf. Environ. Prot. 99 (2016) 167). The present work is devoted to methane explosions because the natural gas—a common fossil fuel used for various domestic, commercial, and industrial purposes—has methane as its major constituent. Specifically, the dynamics of methane-air explosion in vented cylindrical enclosures is scrutinized, computationally and experimentally, such that the accuracy of the EVA predictions is validated by the experiments, with the Cantera package integrated into the EVA to identify the flame speeds. The EVA results for the rear-ignited vented methane-air explosion show good agreement with the experimental results.

1. Introduction

Recent explosions in Beirut, Lebanon [1], Baltimore, MD, USA [2], and Ajman, UAE [3], which all occurred within one week in August 2020, and claimed numerous lives and expensive properties, reinitialize the importance of fire safety, with a focus on scrutinizing of the nature of the explosions. Specifically, understanding the dynamics of an accidental gaseous explosion would provide a way to develop innovative solutions to prevent subsequent explosions. These novel solutions can help to reduce the frequency of occurrence and mitigate the impacts of an explosion. This is a vital research area, since fossil fuels are still much employed in everyday lives at homes, offices, and industries, which sometimes requires storing them in small quantities.
Explosion venting has been a way of suppressing the drastic effects of an accidental explosion by reducing the maximum overpressure in the enclosure. Earlier studies on this topic have spanned from the works of Bradley and Mitcheson [4,5] to that of Mulpuru and Wilkin [6], along with subsequent collaborative efforts to develop a model predicting the dynamics of explosion for a hydrogen-air mixture. Other empirical models developed to predict the pressure evolution as well as the peak overpressures in an enclosure included (though not limited to) the FM Global models by Bauwens et al. [7,8,9], which were based on the experimental works on hydrogen and hydrocarbon explosions. In particular, it has been suggested that external explosion, flame-acoustic interaction, and the flame wrinkling caused by obstacles are responsible for the multiple peak overpressures noticed in the pressure history. Molkov and Bragin [10] also developed a model predicting peak overpressure for vented hydrogen-air explosion based on the turbulent Bradley number, which was presumed to correlate with the overpressure in an enclosure. Sinha et al. [11] have developed another model predicting vented hydrogen-air explosion, based on the external cloud formation and explosion. In another work of the same team, a simple model was created to predict the peak pressure in the vented explosions of hydrogen-air and hydrocarbon-air mixtures [12]. Here, simplifying assumptions were employed to generate a single equation to predict the overpressure in the vented explosions based on four parameters: two of these parameters depended on the fuel properties, which could be pre-calculated, while the other two were functions of the enclosure geometry. Such a reduced simple model minimizes the numerical efforts needed to estimate the maximum overpressure resulting from the gas explosion in an enclosure. The experimental works [7,8,9,13,14,15,16,17,18] on hydrogen-air, hydrocarbon-air, and natural-gas-air explosions, employing various geometries and conditions, have been used to validate some of the models predicting the dynamics of gaseous explosions. The drawbacks of an explosion experiment, which includes requiring ample time, difficulty in setting-up an experiment, coupled with severity associated with the experiment in the case of an accident, actually makes it a dangerous venture. Computational Fluid Dynamics (CFD) applications were also used in modelling gaseous explosions [19], but they also required a lot of computational time and resources.
To reduce such a demands in the computational time and costs, Ugarte et al. [20] developed and coded a reduced-order transient model named the Explosion Venting Analyzer (EVA), which was built-up on the earlier explosion predicting model equations [6] and is able to promptly predict hydrogen explosions in spheres, cuboids, and cylindrical geometries. A subsequent work by Sezer et al. [21] extended the EVA to predict hydrocarbon (methane and propane) explosions. It should be noted that the EVA has undergone various modifications for better prediction accuracy, including the formulations for the flame shape and the burning velocity as well as the addition of the new enclosure geometry. While Sezer et al. [21] validated and established the EVA capability of predicting hydrocarbon explosion experiments in a cuboid enclosure, later, the experimental work of Kodakoglu et al. [22] analyzed methane-air explosions in a cylindrical geometry. Therefore, the present work extends and employs the EVA to model the conditions and geometry of Ref. [22], which resembles a typical configuration of tubes and pipelines used for the natural gas supply. Specifically, the pressure evolution and peak pressure results from the EVA are tabulated and validated by the experimental data [22] in order to ascertain the accuracy of the EVA prediction in this geometry.
It is noted that the flame velocity model embedded in the EVA is based on the experimental correlations for a specific fuel mixture composition and its range of equivalence ratios. While the EVA can predict the explosion dynamics for various fuel mixtures, experimental correlations for the flame speeds are usually not available for most fuel mixtures, which therefore limits the EVA. Therefore, to fix this uncertainty, the Cantera package [23] is integrated with the EVA to compute the flame speeds, with the results being in good agreement with the available experimental correlations. As a result, by integrating Cantera with the EVA, the earlier limitation of the EVA usability, being restricted to a small range of the fuel mixture compositions and equivalence ratios, is removed.
In the present study, a physics-based phenomenological tool that can potentially simulate the explosion dynamics of any fuel mixture composition and its equivalence ratios is developed. This transient reduced-order model predicting explosions can then be utilized to determine the maximum overpressure when designing safety vents employed to mitigate the consequences of explosions.

2. Model Description

The mathematical model employed in this work, the EVA, was developed by Ugarte et al. [20] and it has undergone several modifications as mentioned in Ref [21]. The EVA is a transient, reduced-order model used to predict the maximum (over)pressure, the mass transfer, and the flame speed during explosions in vented and unvented enclosures. The EVA is capable to estimate the explosion characteristics in various geometries by solving the mass and energy balance equations, along with the burning rate formulations, as summarized below.
The mass balance:
d d t ( m u m i ) + d d t ( n ) + d d t ( m v m i ) = 0 ,
m u m i = P ¯ 1 γ u ( 1 V ¯ ) ,       n = m b m i .  
Here, t and m stand for the time and mass, respectively, the subscripts u , b , v , and i designate the unburned, burned, vented, and initial conditions, P ¯ = P / P i is instantaneous-to-initial pressures ratio, V ¯ is the initial volume occupied by the burnt gas, and γ u = C p / C v represents the specific heat ratio of the unburnt gas, employed from the NASA-CEA solver [24] embedded in the EVA to calculate the thermal-physical properties of the respective fuel-air mixtures.
The burning rate formulation:
d n d t = A V i S T   P ¯ 1 / γ u K d d t ( m v m i ) b .  
Equation (3) determines the rate at which the burnt matter is generated, as a function of the burning velocity with respect to the fuel mixture denoted as S T , the surface area of the flame front A , and the initial volume of the fuel mixture denoted as V i . The term K in Equation (3) is used to specify the type of gas vented, with K = 0 if the unburned gas is vented, and K = 1 if the burnt gas is vented.
The energy balance:
d d t [ P ¯ ( 1 V ¯ ] + γ u 1 γ b 1 d d t ( P ¯ V ¯ ) = b d n   d t P ¯ ( 1 1 γ u ) d d t ( m v m i ) ,
b = ( e u o e b o + T 0 [ C V b C V u   ] ) / C V u T i ,           V ¯ = V b / V i .
Here, e o is specific energy of formation, along with the reference temperature T 0 , the specific heat at constant volume C V , and the specific heat ratio of the burnt gas γ b . Further details of the original EVA model are given in Refs. [20,21].

2.1. Models for the Flame Shape

The flame shape is one of the major parameters imposed into the EVA solver, being either spherical or ellipsoidal. The latter is one of the modifications made to the original EVA platform. It is noted that the pioneering Mulpuru–Wilkin predictive explosion model [6] considered a spherical flame shape. Such a choice was suitable to predict the explosion dynamics in the configuration used in the experimental work, as there was good match between the experimental work and the model. However, for most geometries, the spherical flame assumption is not accurate, leading to a notable discrepancy between the results of the predictive explosion model as compared to the experimental results.
According to Ref. [21], the flame shape influences the maximum peak pressure values as well as the explosion behavior. In particular, the ellipsoidal shape of the flame front in Ref. [21] provided better predictions of the explosion characteristics in the cylindrical geometry as compared to the spherical flame shape. Therefore, in the present work, where a cylindrical geometry is being considered, the ellipsoidal flame shape, illustrated in Figure 1, is used. The fuel-air mixture for this work was ignited at the rear end far away from the vent area. Figure 1b therefore depicts the flame shape behavior in this study better.

2.2. Calculation of the Flame Speed

The laminar burning velocity S L is a key parameter used in predictive models characterizing explosions. Consequently, accurate determination of the value S L yields a substantial effect on the accuracy of the results of the predictive model. When a particular S L -value is not available, it is computed in the EVA by using the correlation formula [25,26]
S L = S L , 0 ( T u T u , 0 ) α ( P u P u , 0 ) β
as a function of unburnt gas temperature, T u , and pressure, P u , scaled by their initial values T u , 0 and P u , 0 . Here S L , 0 is the laminar burning velocity at the initial temperature and pressure, S L , 0 = S L ( T u , 0 , P u , 0 ) , while the exponents α and β are usually functions of the fuel-to-oxidizer equivalence ratio ϕ . Stone et al. [25] suggested the following relations for the methane-air mixtures:
S L , 0 = 37.6 + 15.1 ( ϕ 1 ) 221 ( ϕ 1 ) 2 458 ( ϕ 1 ) 3 358 ( ϕ 1 ) 4 , α = 1.42 + 1.98   ( ϕ 1 ) ,     β = 0.314 + 608   ( ϕ 1 ) .
As a result, in the present work, we use Equations (6) and (7) to calculate S L .
Another approach used in earlier works [27,28] to calculate the burning velocity is to employ a combustion modeling software such as Cantera—an open-source chemical kinetics software having various available chemical kinetics model, embedded into it, to compute S L . Specifically, Cantera identifies the laminar burning velocity through modelling a one-dimensional (1D) freely propagating laminar premixed flame front. Then the laminar flame speed is calculated as [29,30]
  S L =   d 2 m d t 2   1 ρ u ,
where ρ u is the density of the unburnt gas, while the quantity d 2 m / d t 2 for the mass flux through a freely propagating flame is calculated according to Ref. [31]. Further details regarding how Cantera calculates the flame speed for various fuel-oxidizer mixtures can be found in Refs. [29,30,31].

2.3. Description of the Experimental Set-Up

We next briefly describe the experimental setup and conditions of Kodakoglu et al. [22], which are simulated in this present study using the EVA solver. Specifically, the EVA results for fuel-lean ( ϕ = 0.8 ), stoichiometric ( ϕ = 1 ), and fuel-rich ( ϕ = 1.2 ) methane-air mixtures were compared with the experiments for rear ignition occurring in a cylindrical chamber having the diameter of 19 cm and the length of 30 cm such that the internal volume was 8,505.6 cm3. Three various circular venting openings have been considered, namely: a small vent (SV) of the venting area 67.9 cm2, a medium vent (MV) of area 88.6 cm2, and a large vent (LV) of area 132.7 cm2. The experiments [22] are illustrated in Figure 2. Specifically, Figure 2a shows a schematic of the experimental setup. Here, the front panel accommodates the vent while the rear panel is used for injection of the gas mixture and ignition probes. The vent was covered with a thin foil layer and was cut prior to ignition. Figure 2b presents evolutions of the elliptic flames for all three vents considered; see the top, middle, and bottom rows for the small, middle, and large vents, respectively.

3. Results and Discussion

In the past, Dong et al. [19] reported a combined experimental and computational investigation of vented explosions of a hydrogen-air mixture in a cylinder of length 100 cm and diameter 18 cm. Consequently, we started the present work with employing the EVA to predict the pressure of hydrogen-air explosions in such a cylindrical enclosure. The initial conditions were taken as P u , 0 = 1 atm and T u , 0 = 25   ° C . Figure 3 compares the EVA results with the experimental data and the CFD modelling of [19]. Specifically, time evolution of the scaled pressure, P 0 / P i , in the case of rear ignition (RI) is depicted in Figure 3, with good agreement between the EVA and Ref. [19] seen.
In addition, we have tested how the flame shape and the ignition location influence the accuracy of the EVA predictions. This is shown in Figure 4, where the maximum peak pressures are presented versus the vent area, with two options for the flame morphology—the spherical flame (SF) shape and the ellipsoidal flame (EF) shape—considered for both central ignition (CI) and rear ignition (RI). Here, the EVA results are compared with the experimental data, shown by the markers in this figure. Again, we see good agreement between the EVA and the experiments in Figure 4. More details about these simulations can be found in Ref. [21].
Unlike Figure 3 and Figure 4, devoted to hydrogen-air explosions, in the present work we mainly focus on vented explosions of a stoichiometric methane-air mixture, with various vent areas considered. The model configuration imitates that of Ref. [22], with RI employed as the location of ignition. Below, the corresponding pressure evolutions will be shown and validated by the experiments [22].

3.1. Pressure Evolution

Figure 5 shows the overpressure versus time for various vent areas, namely, A v = 67.9 cm2, 86.6 cm2, and 132.7 cm2 in Figure 5a–c, respectively. In all three cases, the EVA predictions are validated by the experiments [22]. Specifically, in Figure 5a, the transient pressure, recorded from the experiment, and the EVA simulation results matched initially at the onset. Quantitatively, the peak pressure were 0.036 bar-g and 0.031 bar-g for the EVA and the experiments, respectively. While the maximum overpressure from the EVA appeared 14% higher than that reported in the experiments; in fact, such an over-prediction is good for safety purposes (the key point of the EVA is to predict the worst scenario and to not underestimate the risks). The EVA simulations attained the maximum overpressure quicker as compared to the experiments, which shows some timing discrepancy. It is also noted that after attaining the maximum overpressure, the pressure history results from the EVA stopped matching the experimental results; such a deviation can be due to combustion instabilities or other factors from the experiments not accounted in the EVA.
Figure 5b is a counterpart of Figure 5a for the medium vent area, A v = 86.6 cm2. Both figures demonstrate qualitatively similar trends. The maximum overpressures for the middle vent case, Figure 5b, are 0.023 and 0.015 bar-g from the EVA and the experiments, respectively. Here, the EVA over-predicted the experimental values peak pressure by 53%—much stronger as compared to the 14% in the small vent case, Figure 5a. Regarding Figure 5b, it is also noted that the time the EVA attained the peak pressure fits the experiment peak pressure time well. We also observe that the EVA transient pressure results of Figure 5b match the experimental values until the peak pressure is attained; thereafter, the EVA pressure results dropped and transitioned to an almost constant value.
Finally, when the vent area has been further increased, to A v = 132.7 cm2 (the large vent case, Figure 5c), the EVA pressure history and timing generally match the experimental results before attaining the peak pressure. Quantitatively, the EVA predicted the overpressure of 0.010 bar-g in this case, while the experimental value was 0.016 bar-g.

3.2. Maximum Overpressure

It is recalled that the Cantera software has been integrated with the EVA—in order to be able to compute the laminar flame speeds by means of Cantera and, thereby, provide an alternative to the laminar flame velocity correlations originally used in the EVA. This integration was implemented at the very beginning of the present study in order to enhance the versatility of the EVA code, making it able to be used for a wider range of fuels.
Figure 6 compares the peak pressures obtained from the original EVA (blue) as well as from the EVA with Cantera (black) with those from the experiments [22] (red) for various vent areas. The results calculated by means of the NFPA 68 standards [32] are shown by the green circles. The results are presented in both the normal, Figure 6a, and the semi-logarithmic, Figure 6b, scales. It is visually seen from Figure 6 that the EVA predictions generally agree with the experiments [22]. Moreover, it is clearly seen that both EVA cases (with and without Cantera) show much better agreement with the experiments [22] than the NFPA 68 standards. The corresponding peak pressure values as well as the errors of the approaches are tabulated in Table 1. It is seen that the EVA-Cantera model yields higher overpressures than the original EVA. Obviously, this is because Cantera provided the higher laminar burning velocity value as compared to the original EVA correlations. It is seen that it was particularly beneficial to use the integrated EVA-Cantera model for the large vent area case, A v = 132.7 cm2: the EVA-Cantera model predicted the experimental peak pressure better than the original EVA model in this case. In contrast, for other two cases of the medium and small vent areas, 86.6 cm2 and 67.9 cm2, the original EVA model predicted the experiments [22] better than the integrated EVA-Cantera.

4. Conclusions

The present work showcased the capability of using the Explosion Venting Analyzer (EVA) model [20,21] to predict the pressure evolution in the process of methane-air mixture explosions in vented cylindrical vessels. The EVA results are validated by the experimental measurements [19,22]. It is shown that the EVA over-predicts the peak pressure for small and medium vent area of 67.9 cm2 and 86.6 cm2 but under-predicts the peak pressure when the vent area is as large as 132.7 cm2. This might have resulted from the combustion instabilities or other factors not accounted by the EVA. Regarding the pressure evolution with time, for a small vent, 67.9 cm2, the EVA predicted faster pressure raise than that in the experiments [22], while for the 86.6 cm2 and 132.7 cm2 vent areas, the EVA results generally agree with the experimental data. Also, Cantera was integrated with the EVA to compute the laminar flame velocity. As a result, the EVA-Cantera model predicts higher overpressure as compared to the original EVA. Overall, the EVA is proven to be an acceptable model as its results have good match with the experimental values. This therefore demonstrate the usability of the EVA in determining the peak pressures of gas explosions in cylindrical enclosure.

Author Contributions

Conceptualization, S.O., H.S. and V.A.; methodology, H.S., H.F.F., A.S.R. and V.A.; software, S.O. and H.S.; validation, F.K.; investigation, S.O., H.S., F.K. and H.F.F.; resources, F.K. and H.F.F.; data curation, S.O., F.K. and H.F.F.; writing—original draft preparation, S.O.; writing—review and editing, H.S. and V.A.; supervision, H.F.F., A.S.R. and V.A.; project administration, A.S.R., H.S. and V.A.; funding acquisition, V.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the Alpha Foundation for the Improvement of Mine Safety and Health, Inc., through Award #ASTI14FO-82 (V.A.), prior to which it was supported by the U.S. National Science Foundation (NSF) through CAREER Award No. 1554254 (V.A.) as well as by the West Virginia Higher Education Policy Commission through Grant #HEPC.dsr.18.7 (V.A.). Any views, opinions, findings, conclusions or recommendations expressed in this material are solely those of the authors and do not necessarily reflect the views of the NSF, and do not imply any endorsement by the Alpha Foundation, its directors and staff.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beirut Explosion: What We Know So Far—BBC News. Available online: https://www.bbc.com/news/world-middle-east-53668493 (accessed on 7 September 2020).
  2. Baltimore Explosion: 1 Dead and 7 Hospitalized after Homes Explode—CNN. Available online: https://www.cnn.com/2020/08/10/us/baltimore-maryland-house-explosion/index.html (accessed on 7 September 2020).
  3. Massive Fire Breaks Out at Ajman Market in UAE | Fox News. Available online: https://www.foxnews.com/world/massive-fire-market-ajman-dubai (accessed on 7 September 2020).
  4. Bradley, D.; Mitcheson, A. The venting of gaseous explosions in spherical vessels. I.-Theory. Combust. Flame 1978, 32, 221–236. [Google Scholar] [CrossRef]
  5. Bradley, D.; Mitcheson, A. The venting of gaseous explosions in spherical vessels. II-Theory and experiment. Combust. Flame 1978, 32, 237–255. [Google Scholar] [CrossRef]
  6. Mulpuru, S.; Wilkin, G. A Model for Vented Deflagration of Hydrogen in a Volume; Atomic Energy of Canada Ltd.: Chalk River, ON, Canada, 1982. [Google Scholar]
  7. Bauwens, C.R.; Chaffee, J.; Dorofeev, S. Effect of ignition location, vent size, and obstacles on vented explosion overpressures in propane-air mixtures. Combust. Sci. Technol. 2010, 182, 1915–1932. [Google Scholar] [CrossRef]
  8. Bauwens, C.; Chaffee, J.; Dorofeev, S. Vented explosion overpressures from combustion of hydrogen and hydrocarbon mixtures. Int. J. Hydrogen Energy 2011, 36, 2329–2336. [Google Scholar] [CrossRef]
  9. Bauwens, C.; Chao, J.; Dorofeev, S. Effect of hydrogen concentration on vented explosion overpressures from lean hydrogen–air deflagrations. Int. J. Hydrogen Energy 2012, 37, 17599–17605. [Google Scholar] [CrossRef]
  10. Molkov, V.; Bragin, M. Hydrogen–air deflagrations: Vent sizing correlation for low-strength equipment and buildings. Int. J. Hydrogen Energy 2015, 40, 1256–1266. [Google Scholar] [CrossRef]
  11. Sinha, A.; Rao, V.C.M.; Wen, J.X. Performance evaluation of empirical models for vented lean hydrogen explosions. Int. J. Hydrogen Energy 2019, 44, 8711–8726. [Google Scholar] [CrossRef] [Green Version]
  12. Sinha, A.; Wen, J.X. A simple model for calculating peak pressure in vented explosions of hydrogen and hydrocarbons. Int. J. Hydrogen Energy 2019, 44, 22719–22732. [Google Scholar] [CrossRef]
  13. Kumar, R.K.; Dewit, W.A.; Greig, D.R. Vented explosion of hydrogen-air mixtures in a large volume. Combust. Sci. Technol. 1989, 66, 251–266. [Google Scholar] [CrossRef]
  14. Daubech, J.; Leprette, E.; Duclos, A.; Proust, C. Accounting for turbulence in gas explosion venting design. In 12. International Symposium on Hazards, Prevention, and Mitigationof Industrial Explosions (ISHPMIE); HAL: Kansas City, MO, USA, 2018. [Google Scholar]
  15. Skjold, T.; Hisken, H.; Lakshmipathy, S.; Atanga, G.; Bernard, L.; Van Wingerden, M.; Olsen, K.; Holme, M.; Turøy, N.; Mykleby, M. Vented hydrogen deflagrations in containers: Effect of congestion for homogeneous and inhomogeneous mixtures. Int. J. Hydrogen Energy 2019, 44, 8819–8832. [Google Scholar] [CrossRef]
  16. Chao, J.; Bauwens, C.; Dorofeev, S. An analysis of peak overpressures in vented gaseous explosions. Proc. Combust. Inst. 2011, 33, 2367–2374. [Google Scholar] [CrossRef]
  17. Harrison, A.J.; Eyre, J.A. External explosions as a result of explosion venting. Combust. Sci. Technol. 1987, 52, 91–106. [Google Scholar] [CrossRef]
  18. Tomlin, G.; Johnson, D.; Cronin, P.; Phylaktou, H.; Andrews, G. The effect of vent size and congestion in large-scale vented natural gas/air explosions. J. Loss Prev. Process. Ind. 2015, 35, 169–181. [Google Scholar] [CrossRef]
  19. Fan, B.; Jiang, X.; Chen, Z.; Ye, J.; Dong, G. Investigation on external explosions during venting. Fire Saf. Sci. 2005, 8, 1365–1374. [Google Scholar] [CrossRef]
  20. Ugarte, O.J.; Akkerman, V.; Rangwala, A.S. A computational platform for gas explosion venting. Process. Saf. Environ. Prot. 2016, 99, 167–174. [Google Scholar] [CrossRef]
  21. Sezer, H.; Kronz, F.; Akkerman, V.; Rangwala, A.S. Methane-induced explosions in vented enclosures. J. Loss Prev. Process. Ind. 2017, 48, 199–206. [Google Scholar] [CrossRef]
  22. Kodakoglu, F.; Farahani, H.F.; Rangwala, A.S.; Akkerman, V. Dynamics of explosion venting in a compartment with methane-air mixtures. J. Loss Prev. Process. Ind. 2020, 67, 104230. [Google Scholar] [CrossRef]
  23. Cantera: An Object-Oriented Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport Processes. Version 2.2.0. Available online: https://www.researchgate.net/publication/308982625_Cantera_An_Object-oriented_Software_Toolkit_for_Chemical_Kinetics_Thermodynamics_and_Transport_Processes_Version_220 (accessed on 7 September 2020).
  24. Computer Program for Calculation of Complex Chemical Equilibrium Compositions and Applications, I. Analysis. Available online: https://www.grc.nasa.gov/www/CEAWeb/RP-1311.htm (accessed on 7 September 2020).
  25. Stone, R.; Clarke, A.; Beckwith, P. Correlations for the Laminar-Burning Velocity of Methane/Diluent/Air Mixtures Obtained in Free-Fall Experiments. Combust. Flame 1998, 114, 546–555. [Google Scholar] [CrossRef]
  26. Liao, S.; Jiang, D.; Cheng, Q. Determination of laminar burning velocities for natural gas. Fuel 2004, 83, 1247–1250. [Google Scholar] [CrossRef]
  27. Strakey, P. Oxy-combustion Fundamentals for Direct Fired Cycles; National Energy Technology Laboratory (NETL): Pittsburgh, PA, USA, 2017.
  28. Le Dortz, R.; Bellenoue, M.; Bonneau, L. Laminar burning velocities and Markstein lengths of jet fuel surrogate/air mixtures in a spherical chamber. In Proceedings of the 8th European Combustion Meeting, Dubrovnik, Croatia, 18–21 April 2017. [Google Scholar]
  29. Kee, R.J.; Coltrin, M.E.; Glarborg, P.; Zhu, H. Chemically Reacting Flow; Wiley: Hoboken, NJ, USA, 2017. [Google Scholar]
  30. One-Dimensional Flames, Cantera. Available online: https://cantera.org/science/flames.html#kee2017 (accessed on 9 September 2020).
  31. Smooke, D.M.; Miller, J.A.; Kee, R.J. Determination of adiabatic flame speeds by boundary value methods. Combust. Sci. Technol. 1983, 34, 79–90. [Google Scholar] [CrossRef]
  32. National Fire Protection Association. NFPA 68: Standard on Explosion Protection by Deflagration Venting; National Fire Protection Association: Quincy, MA, USA, 2013. [Google Scholar]
Figure 1. The ellipsoidal flame geometry for (a) central ignition (CI) and (b) rear ignition (RI) [21].
Figure 1. The ellipsoidal flame geometry for (a) central ignition (CI) and (b) rear ignition (RI) [21].
Fire 04 00009 g001
Figure 2. Experiments for vented methane-air explosions with rear ignition (RI) [22]: A sketch of the experimental setup (a) and consecutive photos for the evolution of stoichiometric ( ϕ = 1 ) methane-air, elliptic-like flames inside a cylinder with vent areas: 67.9 cm2 (top), 88.6 cm2 (middle), and 132.7 cm2 (bottom) (b).
Figure 2. Experiments for vented methane-air explosions with rear ignition (RI) [22]: A sketch of the experimental setup (a) and consecutive photos for the evolution of stoichiometric ( ϕ = 1 ) methane-air, elliptic-like flames inside a cylinder with vent areas: 67.9 cm2 (top), 88.6 cm2 (middle), and 132.7 cm2 (bottom) (b).
Fire 04 00009 g002
Figure 3. Validation of the EVA by Ref. [19]: The scaled pressure P 0 / P i versus time for a hydrogen-air explosion in a vented cylindrical enclosure of length 100 cm and diameter 18 cm.
Figure 3. Validation of the EVA by Ref. [19]: The scaled pressure P 0 / P i versus time for a hydrogen-air explosion in a vented cylindrical enclosure of length 100 cm and diameter 18 cm.
Fire 04 00009 g003
Figure 4. The maximum scaled peak pressure versus the vent area for a hydrogen-air explosion in a vented cylindrical enclosure of length 77.5 cm and diameter 59.4 cm.
Figure 4. The maximum scaled peak pressure versus the vent area for a hydrogen-air explosion in a vented cylindrical enclosure of length 77.5 cm and diameter 59.4 cm.
Fire 04 00009 g004
Figure 5. Time evolution of the overpressure for stoichiometric methane-air explosions in a cylindrical enclosure with the vent areas 67.9 cm2 (a), 86.6 cm2 (b), and 132.7 cm2 (c).
Figure 5. Time evolution of the overpressure for stoichiometric methane-air explosions in a cylindrical enclosure with the vent areas 67.9 cm2 (a), 86.6 cm2 (b), and 132.7 cm2 (c).
Fire 04 00009 g005
Figure 6. The peak pressures versus the vent area in the normal (a) and semi-logarithmic (b) scales: Validation of the EVA without (blue diamonds) and with (black circles) Cantera by the experiments [22] (red squares). The green squares show the results calculated through the NFPA 68 standards [32].
Figure 6. The peak pressures versus the vent area in the normal (a) and semi-logarithmic (b) scales: Validation of the EVA without (blue diamonds) and with (black circles) Cantera by the experiments [22] (red squares). The green squares show the results calculated through the NFPA 68 standards [32].
Fire 04 00009 g006
Table 1. The peak pressures as well as the errors of the approaches for the rear-ignited (RI) stoichiometric ( ϕ = 1 ) methane-air explosions in the cylindrical enclosures.
Table 1. The peak pressures as well as the errors of the approaches for the rear-ignited (RI) stoichiometric ( ϕ = 1 ) methane-air explosions in the cylindrical enclosures.
Vent Area (cm2)Peak Pressure (bar-g)Error (%)
ExperimentEVAEVA +
Cantera
EVAEVA +
Cantera
Av = 67.9—the small vent (SV)0.0310.0360.04616.1248.39
Av = 88.6—the medium vent (MV)0.0150.0230.02953.3393.33
Av = 132.7—the large vent (LV)0.0160.0100.01337.518.75
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ogunfuye, S.; Sezer, H.; Kodakoglu, F.; Farahani, H.F.; Rangwala, A.S.; Akkerman, V. Dynamics of Explosions in Cylindrical Vented Enclosures: Validation of a Computational Model by Experiments. Fire 2021, 4, 9. https://0-doi-org.brum.beds.ac.uk/10.3390/fire4010009

AMA Style

Ogunfuye S, Sezer H, Kodakoglu F, Farahani HF, Rangwala AS, Akkerman V. Dynamics of Explosions in Cylindrical Vented Enclosures: Validation of a Computational Model by Experiments. Fire. 2021; 4(1):9. https://0-doi-org.brum.beds.ac.uk/10.3390/fire4010009

Chicago/Turabian Style

Ogunfuye, Samuel, Hayri Sezer, Furkan Kodakoglu, Hamed Farmahini Farahani, Ali S. Rangwala, and V’yacheslav Akkerman. 2021. "Dynamics of Explosions in Cylindrical Vented Enclosures: Validation of a Computational Model by Experiments" Fire 4, no. 1: 9. https://0-doi-org.brum.beds.ac.uk/10.3390/fire4010009

Article Metrics

Back to TopTop