A computational study of low oxygen flammability limit for thick solid slabs
Introduction
The limiting oxygen index (LOI) is commonly used as one of the numerical measure of the flammability of solids. It is defined as the minimum level of oxygen (expressed in volume percentage) in a flowing homogeneous mixture of oxygen and nitrogen that will barely support combustion of the material in downward spreading mode (candle-like). Experimentally, the LOI is determined by following specific test procedures described in various equivalent test standards (e.g., ASTM D2863, BS2782 PART 141, etc.).
The specimens used in the tests are generally thick (thin film-type specimen can also be used) in the form of rods and bars. A specimen is mounted vertically on a specimen holder and ignited at the top in normal gravity and in an imposed-forced flow (3–5 cm/s, in ASTM D2863) of known oxygen and nitrogen composition. The top of the specimen is engulfed in the flame while the flame front spreads downward around the exterior. Sibulkin and Little [1] and Halli and T'ien [2] have reported their observations of burning phenomenon for downward flame spread over cylindrical PMMA rods where a conical fuel shape formed underneath the flame. As the oxygen level in the flowing gas mixture was reduced, the flame front, which normally anchors on the sides of the fuel (side-stabilized flame), moved back (downstream) relative to the pyrolysis front. The flame-spread rate decreased and the conical surface flattened. The retreated flame remained in the wake region (i.e., a wake flame) behind the flattened top end of the specimen (if given sufficient experimental time) and eventually extinguished as a wake flame upon further reduction in oxygen level. Fig. 1A shows the flame stabilization and the solid specimen shape as described above in a typical LOI testing procedure.
Early LOI studies by Finemore and Martin [3] showed that an LOI would be fairly constant within the velocity range 3–12 cm/s, a conclusion which was later questioned by several investigators [4], [5], [6], [7]. Wharton [5] reports several previous investigations on the subject: while some found LOI values to increase nonlinearly with an increase in gas flow velocity, others reported little or no dependence of LOI on gas flow velocity. Wharton [5], [6] tried to resolve this discrepancy. He reported that ambient air could enter the column under certain conditions of small gas inlet velocity and alter the inlet gas composition and hence the LOI value, but at high velocity where there is no entrained ambient air, the LOI is independent of velocity. McIlhagger and Hill [8] conducted a study on thin polypropylene films where LOI showed an upward trend with increase in volumetric flow rate of the gas. The LOI first increased and then remained approximately constant at higher flow rates. Further, Zhevlakov et al. [7] described their LOI-like device for PMMA and fiberglass-reinforced polymer specimens of cross section and length 15–20 cm. Their observation of flame stabilization was similar to that of Refs. [1], [2]. In their work they defined the LOI as the oxygen level where the flame stopped propagating along the side surface (i.e., the transition limit to wake flame). From 0 to 6 m/s their LOI value increased with increasing forced-flow velocity.
Some of the above-noted observations suggest that flame extinction mechanism for thick fuel specimens may be different from that for thin film-type specimens. For the thicker fuel specimen the flame can stabilize in the wake region; however, for thin specimens the flame is always anchored on the sides. Therefore, one of the objectives of this work is to study the effect of imposed-flow velocity on flame stabilization and extinction of the two modes of flames, namely the side-stabilized flame and the wake flame. The other objective of this study is to investigate the effect of gravity on flame stabilization and extinction phenomena. All previous LOI experiments have been done in normal earth gravity, but the LOI values and trends may not be the same in a zero-gravity environment. Therefore, an assessment of the effect of gravity on LOI is necessary to understand the effectiveness/limitation of these earth-based tests for determining the flammability of material for space application.
These objectives are accomplished by studying a numerical model where a nonspreading flame is stabilized over a two-dimensional (thick), nonregressing solid slab in a mixed-flow environment of gravity and forced flow. Fig. 1B shows the schematic of the model problem showing that the flame is anchored at top of the rectangular fuel specimen. The gravity vector acts downward and the forced flow with upstream velocity is imposed from the bottom and is directed upward. This configuration approximates a realistic situation with small burning rate, near the extinction limit. The specimen considered here is a 1.0 cm (long) slab. In this model fuel sample, the upper 0.5 cm is the fuel vapor source (including both side surfaces and the top surface) and the lower 10.0 cm side surfaces are inert. Note that in the model computation the fuel surfaces do not regress but are geometrically fixed vapor sources. In this configuration both side surface flame stabilization and wake flame stabilization can be achieved. The choice of 0.5 cm fuel surface length in the model is based on the combined considerations of physics and computational cost. This length is larger than thermal-diffusion length at the flame base so that an unambiguous side-stabilized flame can be identified and at the same time this length is not too large so that there is saving in the cost of computation on the combined computation domain for both side-stabilized flame and wake flame. Fig. 1C shows the model configuration, the computational domain, and the grid distribution. The specimen is positioned vertically along with gravitational vector pointing downward in the upstream direction. Because of symmetry about , only half of the domain needs to be considered. marks the beginning of the fuel surface which extends up to . The top surface of the solid, which is also a fuel surface, extends from to (since the width is 1.0 cm). The domain extends to 25 cm in the Y direction and 40 cm in the downstream X direction measured from the beginning of the pyrolysis front. The upstream extent is prescribed at 10 cm .
Section snippets
Numerical model and numerical scheme
A previously established numerical model [9], [10] of opposed-flow flame spread over a thin solid is adopted and modified for the present configuration. The opposed-flow model has been successful in predicting the observed flame-spread behavior (flame-spread rate, flame shape and size, and extinction characteristics) over thin solids [9], [10]. The model consists of two-dimensional, steady, laminar, full Navier–Stokes equations along with conservation equations of mass, momentum, energy, and
Results
As noted, earlier previous experiments [1], [2] in normal gravity indicate the presence of two modes of flame stabilization: side flame and wake flame. The transition from the former to the latter occurs when the oxygen level is sufficiently reduced. In the first part of the results, the flame stabilization characteristics of the two modes are studied in a normal-gravity LOI device (i.e., mixed buoyant-forced-flow conditions). The flame transition oxygen limits (side to wake and vice versa) and
Conclusions
A detailed two-dimensional numerical model was used to study the flame extinction phenomena over a model solid fuel slab in a flow configuration similar to that in the limiting oxygen index test. The emphases of the work are on the effects of imposed-forced oxygen velocity, the gravity level, and the wake flame on the limiting oxygen index. The major findings are:
- (1)
In a high oxygen environment, the diffusion flame in normal gravity is stabilized on the sides of the thick solid sample. As oxygen
Acknowledgements
We thank Ioan Feier for reviewing the manuscript and for his kind comments. This research was supported by NASA Grants NCC-633 and NCC-669 under the technical monitoring of Dr. Kurt Sacksteder at NASA Glenn Research Center and Dr. Mickey King at NASA Headquarters.
References (16)
- et al.
Combust. Flame
(1978) - et al.
Combust. Flame
(2003) - et al.
Combust. Flame
(1995) - Y. Halli, J.S. T'ien, National Bureau of Standards, NBS-GCR-86-507,...
- et al.
J. Mod. Plastics
(1966) - et al.
SPE J.
(1971) J. Fire Sci.
(1983)J. Fire Mater.
(1981)
Cited by (18)
Three-dimensional image reconstruction of a thermoplastic rod burning with and without melting
2022, Combustion and FlameA review of near-limit opposed fire spread
2021, Fire Safety JournalExperimental study on flammability limits of electrolyte solvents in lithium-ion batteries using a wick combustion method
2019, Experimental Thermal and Fluid ScienceCitation Excerpt :The flow velocity is one of the environmental variables that may affect the material flammability and fire dynamics in an opposed flow diffusion flame configuration [30,42]. According to the blow-off mechanism of a wire or candle-like flame [43,44], the increased axial flow velocity can shorten the fluidic residence time and lead to an earlier extinction (corresponding to a higher LOC). The axial flow velocity of the wick-LOC method should also be specified properly for the flammability comparison of solvents.
Experimental study on flame stability limits of lithium ion battery electrolyte solvents with organophosphorus compounds addition using a candle-like wick combustion system
2019, Combustion and FlameCitation Excerpt :The total heat release of the flame is generally proportional to the fuel consumption rate and reflected by the flame volume. While, the quenching mechanism of flame extinction always caused by the excessive heat conduction from the flame reacting zone to the outside leading to the temperature drop, especially for a small-scale flame [19,20]. Thus, the ratio of heat loss to heat generation is dominated by the surface volume ratio (S/V) of the flame.
Transition from opposed flame spread to fuel regression and blow off: Effect of flow, atmosphere, and microgravity
2019, Proceedings of the Combustion Institute
- 1
Present address: Department of Aerospace Engineering, Indian Institute of Technology Madras, Chennai 600036, India.