Thermal explosions in spherical vessels at large Rayleigh numbers. (December 2017)

• Record Type: Journal Article
• Authors: Iglesias, I.; Moreno-Boza, D.; Sánchez, A.L.; Liñán, A.; Williams, F.A.
• Is Part Of: International journal of heat and mass transfer. Volume 115(2017)Part B
• Page Start: 1042
• Page End: 1053
• Publication Date: 2017-12
• Subjects: Thermal explosions -- Laminar reacting flows -- Natural convection -- Boundary-layer theory; 00-01 -- 99-00; Heat -- Transmission -- Periodicals; Mass transfer -- Periodicals; Chaleur -- Transmission -- Périodiques; Transfert de masse -- Périodiques; Electronic journals; 621.4022
• DOI: 10.1016/j.ijheatmasstransfer.2017.08.109
• Languages: English
• ISSNs: 0017-9310

Abstract

Highlights: Numerical integrations are used to assess buoyancy effects on thermal explosions in spherical vessels. Frank-Kamenetskii explosion curves are determined over a wide range of Rayleigh numbers Ra . The boundary-layer structure of the flow in the limit Ra ≫ 1 is investigated by asymptotic methods. The asymptotic analysis leads to accurate predictions of explosion limits. Abstract: This paper investigates effects of buoyancy-driven motion on the "slowly reacting" mode of combustion, and its thermal-explosion limits, of an initially cold gaseous mixture enclosed in a spherical vessel with a constant wall temperature. As in Frank-Kamenetskii's seminal analysis, the strong temperature dependence of the effective overall reaction is modeled with a single irreversible reaction with an Arrhenius rate having a large activation energy. Besides the classical Damköhler number Da, measuring the ratio of the heat-release rate by chemical reaction evaluated at the wall temperature to the rate of heat removal by heat conduction to the wall, the solution is seen to depend on the Rayleigh number Ra, measuring the effect of buoyancy-induced motion on the heat-transport rate. For values of Da below a critical value Da c the system evolves in a slowly reacting mode where the heat losses to the wall limit the temperature increase associated with the chemical reaction, whereas for Da > Da c the initial stage of slow reaction ends abruptly at a well-defined ignition time, at which a thermal runaway occurs. Transient numerical integrations of the initial stage of slow reaction, formulated in the distinguished limit Da ∼ 1 and Ra ∼ 1 with account taken of the effects of the temporal pressure variation, are used to investigate influences of natural convection on thermal-explosion development, including changes in ignition times for Da > Da c and modified explosion curves. Our analysis reveals that Frank-Kamenetskii's criterion for the determination of critical explosion conditions, based on the investigation of existence of steady solutions, provides values of Da c ( Ra ) that are identical to those extracted from the transient computations. Specific consideration is given to the structure of the steady solution in the asymptotic limit Ra ≫ 1 in which the flow includes a thin chemically frozen near-wall boundary layer of downward moving cold gas bounding a central inviscid region of slowly rising reacting flow driven by the boundary-layer entrainment, with the critical explosion conditions predicted to occur for Da c ∼ Ra 1 / 4 . The mathematical structure of the resulting boundary-layer problem is fundamentally similar to that found in the unrelated problem of flow in curved pipes at large Dean numbers. As in that similar problem, the boundary layer exhibits a region of recirculating flow, so that the problem must be formulated as a boundary-value problem accounting for the self-similar local solutions that exist near the upper and lower stagnation points. The problem is solved with use made of an approximate integral method. The resulting asymptotic prediction for the critical Damköhler number Da c = 0.655 Ra 1 / 4 is found to be in excellent agreement with the results obtained by integration of the complete conservation equations for Ra ≫ 1 .