ABSTRACT
Theoretical predictions of dust-explosion characteristics, such as the minimum explosive concentration (MEC), can help improve the accuracy of quantitative risk assessment. This paper presents a continuum model to simulate one-dimensional dust-cloud combustion. It is found that the dilution of dust by thermal expansion reduces the propagation speed of a combustion wave through a dust cloud. At the same time, the delay in the increase of particle velocity due to relaxation time counteracts the dilution effect. Therefore, neglecting the relaxation time results in the underestimation of propagation speed. In the absence of heat loss, extinction occurs when the adiabatic combustion temperature decreases to the ignition temperature; MEC is independent of particle size in such cases. With the presence of heat loss, on the other hand, extinction occurs at a higher dust concentration; MEC increases with particle size. Accurate evaluation of ignition temperature is key to the quantitative prediction of MEC. A particular focus is given to gravity effects, reflecting recent experimental efforts to obtain fundamental insights into propagation mechanisms without being affected by gravity. In upward propagation, particles fall toward the reaction front, enhancing combustion and increasing the propagation speed. The situation in downward propagation is the opposite, and extinction occurs when the gravity level is beyond a critical value. In upward propagation, MEC can decrease with an increase in particle size, owing to the influence of gravity.
Acknowledgements
This work was partly supported by JSPS KAKENHI Grant Numbers JP18H03822, JP19H01807, JP21H04593, and JP21K14379 and partly by the front-loading research project of the Japan Aerospace Exploration Agency (JAXA).
Disclosure statement
No potential conflict of interest was reported by the authors.
Nomenclature
= | dust concentration [kg·m−3] | |
= | specific heat [J·kg−1·K−1] | |
= | particle diameter [m] | |
= | acceleration of gravity [m·s−2] | |
= | Heaviside step function | |
= | heat-loss coefficient [W·m−3·K−1] | |
= | reaction rate constant [m2·s−1] | |
= | inter-particle spacing [m] | |
= | particle number density [m−3] | |
= | heat of combustion [J·kg−1] | |
= | particle Reynolds number [–] | |
= | temperature [K] | |
= | ignition temperature [K] | |
= | velocity [m·s−1] | |
= | moving coordinate [m] | |
Greek symbols | = | |
= | emissivity [–] | |
= | thermal conductivity [W·m−1·K−1] | |
= | viscosity [Pa·s] | |
= | density [kg·m−3] | |
= | thermal expansion ratio, , [–] | |
= | Stefan-Boltzmann constant [W·m−2·K−4] | |
= | relaxation time [s] | |
Subscripts | = | |
0 | = | unburnt state |
c | = | convection |
g | = | gas phase |
p | = | particle |
r | = | radiation |