Numerical simulation for explosion wave propagation of combustible mixture gas

Abstract

A two-dimensional multi-material code was indigenously developed to investigate the effects of duct boundary conditions and ignition positions on the propagation law of explosion wave for hydrogen and methane-based combustible mixture gas. In the code, Young’s technique was employed to track the interface between the explosion products and air, and combustible function model was adopted to simulate ignition process. The code was employed to study explosion flow field inside and outside the duct and to obtain peak pressures in different boundary conditions and ignition positions. Numerical results suggest that during the propagation in a duct, for point initiation, the curvature of spherical wave front gradually decreases and evolves into plane wave. Due to the multiple reflections on the duct wall, multi-peak values appear on pressure—time curve, and peak pressure strongly relies on the duct boundary conditions and ignition position. When explosive wave reaches the exit of the duct, explosion products expand outward and forms shock wave in air. Multiple rarefaction waves also occur and propagate upstream along the duct to decrease the pressure in the duct. The results are in agreement with one-dimensional isentropic gas flow theory of the explosion products, and indicate that the ignition model and multi-material interface treatment method are feasible.

Keywords

combustible mixture gas / explosion wave / interface treatment / combustion function mode

References

[1] LinB.-q., ZhangR.-g., LüH.-hong.. Research on acceleration mechanism and flame transmission in gas explosion [J]. Journal of China Coal Society, 1999, 24(1): 164-167

[2] LuJ., NingJ.-g., WangC., LinB.-quan.. Study on flame propagation and acceleration mechanism of city coal gas [J]. Explosion and Shock Waves, 2004, 24(4): 305-311

[3] VaagsaetherK., KnudsenV., BjerketvedtD.. Simulation of flame acceleration and DDT in H₂-air mixture with a flux limiter centered method [J]. International Journal of Hydrogen Energy, 2007, 32(13): 2186-2191

[4] OranE. S., GamezoV. N.. Origins of the deflagration-to-detonation transition in gas-phase combustion [J]. Combustion and Flame, 2007, 148(1/2): 4-47

[5] FairweatherM., HargraveG. K., IbrahimS. S.. Studies of premixed flame propagation in explosion tubes [J]. Combustion and Flame, 1999, 116(4): 504-518

[6] DunnR. D., MccannM. A.. Overpressure from nondetonating baffle-accelerated turbulent flame in tubes [J]. Combustion and Flame, 2000, 120(4): 504-514

[7] GamezoV. N., OgawaT., OranE. S.. Numerical simulations of flame propagation and DDT in obstructed channels filled with hydrogen-air mixture [C]. Proceedings of the Combustion Institute, 2007, 31(2): 2463-2471

[8] NaamansenP., BaraldiD., HjertagerB. H., SolbergT., CantS.. Solution adaptive CFD simulation of premixed flame propagation over various solid obstructions [J]. Journal of Loss Prevention in the Process Industries, 2002, 15(3): 189-197

[9] NingJ.-g., WangC., LuJ.. Explosion characteristics of coal gas under various initial temperatures and pressures [J]. Shock Waves, 2006, 15(6): 461-472

[10] WangC., NingJ.-g., LeiJ.. Numerical study on propagation of explosion wave in H₂-O₂ mixtures [J]. Lecture Notes in Computer Science, 2006, 4222: 816-819

[11] OranE. S., WeberJ. E., StefaniwE. I., LefebvreM. H., AndersonJ. D.. A numerical study of two-dimensional H₂-O₂-Ar detonation using a detailed chemical reaction model [J]. Combustion and Flame, 1998, 113(1/2): 147-163

[12] BielertU., SicheM.. Numerical simulation of premixed combustion processes in closed tubes [J]. Combustion and Flame, 1998, 114(3/4): 397-419

[13] WangC.-j., XuS.-li.. Numerical study on re-initiation phenomenon of gaseous detonation induced by the coalescence of high temperature regions [J]. Acta Aerodynamica Sinica, 2007, 25(1): 70-75

[14] ZhangD.-l., XieW., GuoC.-m., HuX.-yu.. Numerical simulation of cellar structures and mach reflection of gaseous detonation waves [J]. Explosion and Shock Waves, 2001, 21(3): 161-163

[15] NingJ.-g., ChenL.-wei.. Fuzzy interface treatment in Eulerian method [J]. Science in China (Series E): Engineering & Materials Science, 2004, 47(5): 550-568

[16] OuJ.-p., MaA.-c., ZhanS.-h., ZhouJ.-m., XiaoZ.-qiang.. Dynamic simulation on effect of flame arrangement on thermal process of regenerative reheating furnace [J]. Journal of Central South University of Technology, 2007, 14(2): 243-247

[17] WangG., LiuS.-j., LiL.. FEM modeling for 3D dynamic analysis of deep-ocean mining pipeline and its experimental verification [J]. Journal of Central South University of Technology, 2007, 14(6): 808-813

[18] ZhangB.-p., ZhangQ.-m., HuangF.-lei.Detonation physics [M], 2001, Beijing, Ordnance Industry Press

[19] PilliodJ. E., PuckettE. G.. Second-order accurate volume-of-fluid algorithms for tracking material interfaces [J]. Journal of Computational Physics, 2004, 199(2): 465-502

[20] FAN Xi-sheng. The actuality and study orientation of the combustible gas explosion in pipeline [J]. Industrial Safety and Dust Control, 1998(6): 10–15. (in Chinese)

[21] XuG.-y., GuD.-sheng.. Experimental investigation on generation mechanism of explosive stress wavels [J]. Transactions of Nonferrous Metals Society of China, 1999, 9(2): 433-436

[22] XieF.-h., LuoG.-w., LiaoJ.-sheng.. Numerical analysis of explosive welding interface layer based on δ function [J]. The Chinese Journal of Nonferrous Metals, 2007, 17(12): 2029-2033