ABSTRACT
Deflagration of fire and explosion caused by ignition of a methane-air mixture is a major safety concern in industrial settings, such as the mining and petrochemical industries. The aim of this study is to investigate numerically the effect of tube size on the flame and pressure wave propagations for methane-air deflagration in a tube closed at one end. A partially premixed combustion model that avoids the need to specify the flame speed is deployed based upon the Flamelet Generated Manifold (FGM) model. Good agreement is achieved between the predicted results and the benchmark experimental data collected using a large-scale detonation tube (L/D = 66). Subsequently, the explosion behaviour immediately after the ignition of methane-air mixtures and the propagation characteristics of the flame front and pressure wave through the tube are examined, covering a broad range of L/D, that is from 26 to 526, in which the diameter is changed and the length is kept fixed. The results show that the pressure wave propagates significantly faster in narrower tubes and hence decouples from the flame front shortly after ignition, which in turn results in a low overpressure at the flame front. Moreover, abrupt changes in gas properties are observed in narrow tubes with L/D ≥ 132. The peak overpressure increases as the tube diameter increases; however, the local maximum pressure decreases substantially in large tubes when approaching the tube vent but remains almost constant throughout in the narrow tubes. Similarly, the flame propagates faster in narrower tubes. A correlation that estimates the distance the flame propagates in the exponential acceleration stage is proposed as a function of tube size and time. Deviations of less than 7% are obtained when comparing the predicted results using the correlation against the experimental data. The results provide local information that aids the theoretical interpretation of experimental observations and the understanding of the fuel combustion and explosion phenomena in different sized tubes.
Acknowledgments
The authors wish to acknowledge the financial support by Australian Coal Association Low Emission Technologies Ltd (ACALET), Australian Department of Resources, Energy and Tourism, The University of Newcastle, Australia. This research/project was undertaken with the assistance of resources and services from the National Computational Infrastructure (NCI), which is supported by the Australian Government.
A combustion model that avoids the need to specifying the flame speed is validated.
Global maximum pressure is found at the flame front in narrow tubes.
Local maximum pressure decreases in large tubes when approaching the tube vent.
Abrupt changes in gas properties are observed in narrow tubes.
A correlation predicting the flame propagation distance after ignition is proposed.
Notation
Symbols
= | Constant, m/s | |
= | Constant, 1/s | |
= | Reaction progress variable, - | |
= | Mixture specific heat, J/(kg⋅K) | |
= | Local sound speed, m/s | |
= | Effective diffusion coefficient, kg/(m⋅s) | |
= | Mixture fraction, - | |
= | Maximum slope of temperature profile across the flame front, K/m | |
= | Enthalpy, J | |
= | Thermal conductivity, W/(m⋅K) | |
= | Tube length, m | |
= | Flame thickness, m | |
= | Reactant mass at the stoichiometric condition, kg | |
= | Pressure, Pa | |
= | Joint PDF of reaction-progress (c) and mixture fraction (f), - | |
= | Prandtl number, - | |
= | Energy generation rate from methane-air combustion, W | |
= | Total energy liberated from the reaction per unit volume, Pa | |
= | Radial position inside DT, m | |
= | Radius of DT, m | |
= | Reynolds number, m | |
= | Local position along the DT centreline, m | |
= | Flamelet source term, 1/s | |
= | Finite-Rate flamelet source term from the flamelet library, 1/s | |
= | Time, s | |
= | Temperature, K | |
= | Auto-ignition temperature of methane, K | |
= | Spatial density weighted velocity in the ith direction, m/s | |
= | Velocity, m/s | |
= | Mass fraction, - | |
= | Reynolds averaged un-normalized progress variable, - |
Greek letters
= | Constants, zero for reactants and unity for a few product species, - | |
= | Ratio of specific heat, - | |
= | Pressure drop, Pa | |
= | Distance between burnt and unburnt gases boundaries, m | |
= | Turbulence dissipation rate, m2/s3 | |
= | Gas thermal expansion factor, - | |
= | Turbulence kinetic energy, m2/s2 | |
= | Flame thickness, m | |
= | Viscous diffusion rate (kinematic viscosity), m2/s | |
= | Density, kg/m3 | |
= | Scaled growth rate of flame acceleration, - | |
= | Relative humidity, - | |
= | Scalar dissipation rate, 1/s | |
= | Prescribed maximum scalar dissipation within the premixed flamelet, 1/s | |
= | Reaction rate of the progress variable CC, mol/(s⋅m3) | |
= | Molar production rate of species k, mol/(s⋅m3) |
Subscripts or superscripts
= | Chemical equilibrium | |
= | Flame | |
= | Flame front | |
= | Laminar | |
= | The kth species | |
= | Stoichiometric | |
= | Unburnt mixture |