Simulation of transient convective burning of an n-octane droplet using a four-step reduced mechanism

The transient burning of an n-octane fuel droplet in a hot gas stream is numerically studied using a four-step reduced mechanism, with considerations of droplet surface regression, deceleration due to the drag of the droplet, internal circulation inside the droplet, variable properties, non-uniform surface temperature, and the effect of surface tension. Two different types of the four-step mechanism are examined and found almost identical. The four-step mechanism has earlier instant of the wake-to-envelope transition than the one-step mechanism at low ambient temperature, but this difference between the two mechanisms diminishes when the ambient temperature is increased. The four-step mechanism has smaller mass burning rate for a wake flame but greater mass burning rate for an envelope flame than the one-step mechanism. The two mechanisms have small differences in the critical initial Damkohler number. Lower ambient temperature yields later wake-to-envelope transition and smaller mass burning rate. Higher ambient pressure has greater overall mass burning rate because of greater gas density and thus greater concentrations of reactants for a major part of the lifetime. Greater ambient mass fraction of oxygen yields faster oxidation kinetics and greater Damkohler number. As the ambient mass fraction of oxygen increases, the instant of wake-to-envelope transition advances for an initial wake flame, and finally the initial flame becomes an envelope flame when the ambient mass fraction of oxygen exceeds some critical value. A correlation is developed for the critical initial Damkohler number in terms of the ambient temperature, ambient pressure, and ambient mass fraction of oxygen.     

Transient convective burning of interactive fuel droplets in double-layer arrays

The transient convective burning of n-octane droplets interacting within double-layer arrays in a hot gas flow perpendicular to the layers is studied numerically, with considerations of droplet surface regression, deceleration and relative movement due to the drag of the droplets, internal liquid motion, variable properties, non-uniform liquid temperature and surface tension. Each layer in the double-layer array is a periodic droplet array aligned orthogonal to the free stream direction. The droplets in different layers are arranged either in tandem or staggered. Several different flame structures are found for the double-layer arrays. The transient behaviors of the droplets in both upstream and downstream layers are studied and compared, for various initial relative stream velocity and initial transverse droplet spacing. The average surface temperature and vaporization rate for the front (or upstream) droplets and back (or downstream) droplets are influenced by the flame structure. The front droplets in a double-layer array behave similarly to the droplets in a single-layer array for the streamwise droplet spacing considered in this study. The back droplets approach the front droplets because they generally have lower drag.     

Sooting characteristics of isolated droplet burning in heated ambients under microgravity

This paper presents an investigation into the sooting characteristics of isolated droplets (for fuel n-decane) burning in heated ambients in microgravity. A backlit video view of the droplet was taken to determine the soot shell size and to judge the transient soot generation according to qualitative amount of soot. The independent experiment variables were the ambient temperature and initial droplet diameter. Soot generation was higher for initially larger droplets when compared at the same burning time normalized with the initial droplet diameter squared (called normalized burning time). At the same absolute burning time there existed an obvious initial transient period after ignition in which the stated relationship was not satisfied. This transient time increased with increasing the ambient temperature. There was a peak in the soot generation at about 1000 K throughout the lifetime of the droplet. The soot shell size was generally larger for an initially bigger droplet at the same instantaneous droplet diameter or normalized burning time. At the same absolute burning time, however, an initially smaller droplet exhibited larger relative soot shell sizes (the soot shell size normalized with the initial droplet diameter). The soot shell size increased monotonically with increasing ambient temperature. This is due to the increase in the Stefan flow drag with the larger burning rate at the higher temperature. The consequent result is that the soot shell sizes are much larger for droplets burning in heated ambients than for droplets burning in room-temperature ambients.     

Transport processes and associated irreversibilities in droplet combustion in a convective medium

A mathematical model of the combustion of a droplet surrounded by hot gas with a uniform free stream motion is made from the numerical solution of conservation equations of heat, mass and momentum in both the carrier and the droplet phases. The gas-phase chemical reaction between the fuel vapour and the oxidizer is assumed to be single-step and irreversible. The phenomenon of ignition is recognised by the sudden rise of temperature in the temperature/time histories at different locations in the carrier phase. To ascertain the process irreversibilities, the instantaneous rate of entropy production and its variation with time have been determined from the simultaneous numerical solution of the entropy conservation equations for both the gas and liquid phases. The relative influences of pertinent input parameters, namely the initial Reynolds number Re_i, the ratio of the free stream to initial temperature T_∞ and the ambient pressure on (i) the local and average Nusselt numbers, (ii) the life histories of burning fuel drops, and (iii) the entropy generation rate in the process of droplet combustion have been established.     

Advanced models of fuel droplet heating and evaporation

Recent developments in modelling the heating and evaporation of fuel droplets are reviewed, and unsolved problems are identified. It is noted that modelling transient droplet heating using steady-state correlations for the convective heat transfer coefficient can be misleading. At the initial stage of heating stationary droplets, the well known steady-state result Nu=2 leads to under prediction of the rate of heating, while at the final stage the same result leads to over prediction. The numerical analysis of droplet heating using the effective thermal conductivity model can be based on the analytical solution of the heat conduction equation inside the droplet. This approach was shown to have clear advantages compared with the approach based on the numerical solution of the same equation both from the point of view of accuracy and computer efficiency. When highly accurate calculations are not required, but CPU time economy is essential then the effect of finite thermal conductivity and re-circulation in droplets can be taken into account using the so called parabolic model. For practical applications in computation fluid dynamics (CFD) codes the simplified model for radiative heating, describing the average droplet absorption efficiency factor, appears to be the most useful both from the point of view of accuracy and CPU efficiency. Models describing the effects of multi-component droplets need to be considered when modelling realistic fuel droplet heating and evaporation. However, most of these models are still rather complicated, which limits their wide application in CFD codes. The Distillation Curve Model for multi-component droplets seems to be a reasonable compromise between accuracy and CPU efficiency. The systems of equations describing droplet heating and evaporation and autoignition of fuel vapour/air mixture in individual computational cells are stiff. Establishing hierarchy between these equations, and separate analysis of the equations for fast and slow variables may be a constructive way forward in analysing these systems.     

An elementary model for the transition from conductive to penetrative burning in gas-permeable explosives

The experimentally known phenomenon of spontaneous transition from slow conductive to fast penetrative (convective) burning in a confined gas-permeable explosive is discussed. A reduced quasi-linear model, involving only the most essential physical ingredients, is formulated. A good qualitative agreement between theoretical and experimental dependencies is obtained. Similar to the previously studied case of an inert porous matrix filled with an explosive gas, the transition is triggered by a localized autoignition in the extended resistance-induced preheat zone gradually formed ahead of the advancing deflagration. The conventional concept of penetrative burning is re-examined.     

Flame brush characteristics and burning velocities of premixed turbulent methane/air Bunsen flames

The flame brush characteristics and turbulent burning velocities of premixed turbulent methane/air flames stabilized on a Bunsen-type burner were studied. Particle image velocimetry and Rayleigh scattering techniques were used to measure the instantaneous velocity and temperature fields, respectively. Experiments were performed at various equivalence ratios and bulk flow velocities from 0.7 to 1.0, and 7.7 to 17.0 m/s, respectively. The total turbulence intensity and turbulent integral length scale were controlled by the perforated plate mounted at different positions upstream of the burner exit. The normalized characteristic flame height and centerline flame brush thickness decreased with increasing equivalence ratio, total turbulence intensity, and longitudinal integral length scale, whereas they increased with increasing bulk flow velocity. The normalized horizontal flame brush thickness increased with increasing axial distance from the burner exit and increasing equivalence ratio. The non-dimensional leading edge and half-burning surface turbulent burning velocities increased with increasing non-dimensional turbulence intensity, and they decreased with increasing non-dimensional bulk flow velocity when other turbulence statistics were kept constant. Results show that the non-dimensional leading edge and half-burning surface turbulent burning velocities increased with increasing non-dimensional longitudinal integral length scale. Two correlations to represent the leading edge and half-burning surface turbulent burning velocities were presented as a function of the equivalence ratio, non-dimensional turbulence intensity, non-dimensional bulk flow velocity, and non-dimensional longitudinal integral length scale. Results show that the half-burning surface turbulent burning velocity normalized by the bulk flow velocity decreased as the normalized characteristic flame height increased.     

Laminar burning velocities and flame stability analysis of syngas mixtures at sub-atmospheric pressures

The effects of low pressure on the laminar burning velocity and flame stability of H₂/CO mixtures and equimolar H₂/CO mixtures diluted with N₂ and CO₂ were studied experimentally and theoretically. Experiments were conducted at real sub-atmospheric conditions in three places located at high altitudes 500 m.a.s.l. (0.947 atm), 1550 m.a.s.l. (0.838 atm), and 2300 m.a.s.l. (0.767 atm). Flames were generated using contoured slot-type nozzle burners and Schlieren images were used to determine the laminar burning velocity with the angle method. The behavior of the laminar burning velocity at low pressures depends on the equivalence ratio considered; it decreases at lean and very rich equivalence ratios when pressure is increased. However, a contrary behavior was obtained at equivalence ratios corresponding to the highest values of the laminar burning velocity, where it increases as pressure increases. Numerical calculations were also conducted using a detailed reaction mechanism, and these do not reproduce the behavior obtained experimentally; a sensitivity analysis was carried out to examine the differences found. At lean equivalence ratios, flame instabilities were observed for all the syngas mixtures. The range of equivalence ratios where flames are stable increases at lower pressures. This behavior is due to the increase of the flame thickness, which considerably reduces the hydrodynamic instabilities in the flame front.     

Laminar burning characteristics of premixed methane-dissociated methanol-diluent mixtures

The influence of dissociated methanol (DM) and diluent (CO₂ and N₂) addition on methane was investigated in a constant volume chamber under initial conditions of 3 bar and 343 K. CO was also added in separate proportions instead of DM under the same conditions to assess its effect. The laminar burning velocity, Markstein length and flame instability were analyzed systematically under various equivalence ratios (0.8–1.4), dissociated methanol gas ratios (40 and 80%), CO ratios (40 and 80%) and dilution ratios (0–15%). Furthermore, the flame speed of the fuel mixture and the production rate of key reactants were analyzed based on the calculation results of the Aramco Mech 2.0 mechanism to determine the influence principles of dilution. The results show that dissociated methanol gas increases the flame speed of the mixtures and promotes instability of the flame, and H₂ is the dominant component in enhancing the combustion process. Within the dilution ratio range of this study, the diluents decrease the laminar burning velocity of the mixtures since the addition of diluent gas decreases the concentration of key reactants, such as H and OH. The addition of diluent gas can inhibit the flame instability, but the effect is not clear. Compared with N₂, the effect of CO₂ is more significant.     

Measurements of laminar burning velocities and Markstein lengths for methanol-air-nitrogen mixtures at elevated pressures and temperatures

The laminar burning velocities and Markstein lengths for the methanol-air mixtures were measured at different equivalence ratios, elevated initial pressures and temperatures, and dilution ratios by using a constant volume combustion chamber and high-speed schlieren photography system. The influences of these parameters on the laminar burning velocity and Markstein length were analyzed. The results show that the laminar burning velocity of the methanol-air mixture decreases with an increase in initial pressure and increases with an increase in initial temperature. The Markstein length decreases with an increase in initial pressure and initial temperature, and increases with an increase in the dilution ratio. A cellular flame structure is observed at an early stage of flame propagation. The transition point is identified on the curve of flame propagation speed against stretch rate. The reasons for the cellular structure development are also analyzed.     

Laminar burning velocities and flame stability analysis of hydrogen/air premixed flames at low pressure

An experimental and numerical study on laminar burning velocities of hydrogen/air flames was performed at low pressure, room temperature, and different equivalence ratios. Flames were generated using a small contoured slot-type nozzle burner (5 mm × 13.8  mm). Measurements of laminar burning velocity were conducted using the angle method combined with Schlieren photography. Numerical calculations were also conducted using existing detailed reaction mechanisms and transport properties. Additionally, an analysis of the intrinsic flame instabilities of hydrogen/air flames at low pressure was performed. Results show that the behavior of the laminar burning velocity is not regular when decreasing pressure and that it depends on the equivalence ratio range. The behavior of the laminar burning velocity with decreasing pressure can be reasonably predicted using existing reaction mechanisms; however changes in the magnitude of the laminar burning velocity are underestimated. Finally, it has been found experimentally and proved analytically that the intrinsic flame instabilities are reduced when decreasing the pressure at sub-atmospheric conditions.     

Turbulent burning velocity, burned gas distribution, and associated flame surface definition

Experimental studies of premixed, turbulent, gaseous explosion flames in a fan-stirred bomb are reported. The turbulence was uniform and isotropic, while changes in the rms turbulent velocity were achieved by changes in the speed of the fans. Central spark ignitions created mean spherical flame propagation. The spatial distributions of burned and unburned gases during the propagation were measured from the Mie scattering of tobacco smoke in a thin planar laser sheet. The plane was located just in front of the central spark gap and was generated by a copper vapor laser operating at a pulse rate of 4.5 kHz. High-speed schlieren images also were captured simultaneously.The distributions of the proportions of burned and unburned gases around circumferences were found for all radii at all stages of the explosion, and mean values of these proportions were derived as a function of the mean flame radius. The flame brush thickness increased with flame radius. The way the turbulent burning velocity is defined depends on the chosen associated flame radius. Various definitions are scrutinized and different flame radii presented, along with the associated turbulent burning velocities. Engulfment and mass turbulent burning velocities are compared. It is shown how the latter might conveniently be obtained from schlieren cine images. In a given explosion, the burning velocity increased with time and radius, as a consequence of the continual broadening of the effective spectrum of turbulence to which the flame was subjected. A decrease in the Markstein number of the mixture increased the turbulent burning velocity.     

Structures and burning velocity of biomass derived gas flames

Biomass derived gases produced via gasification, pyrolysis, and fermentation are carbon neutral alternative fuels that can be used in gas turbines, furnaces, and piston engines. To make use of these environmentally friendly but energy density low fuels the combustion characteristics of these fuels have to be fully understood. In this study the structure and laminar burning velocity of biomass derived gas flames are investigated using detailed chemical kinetic simulations. The studied gaseous fuels are the air-blown gasification gas, co-firing of gasification gas with methane, pyrolysis gases, landfill gases, and syngas, a mixture of carbon monoxide and hydrogen. The simulated burning velocities of reference fuel mixtures using two widely used chemical kinetic mechanisms, GRI Mech 3.0 and the San Diego mechanism, are compared with the experimental data to explore the uncertainties and scattering of the simulation data. The different chemical kinetic mechanisms are shown to give a reasonable agreement with each other and with experimental data, with a discrepancy within 7% over most of the conditions. The results show that the structures of typical landfill gas flames and co-firing of methane/gasification gas flames share essential similarity with methane flames. The reaction zones of these flames consist of a thin inner layer and a relatively thick CO/H₂ oxidation layer. In the inner layer hydrocarbon fuel (methane) is converted through chain reactions to intermediates such as CH₃, CH₂O, CO, H₂, etc. The structures of gasification gas flames, pyrolysis gas flames, syngas flames share similarity with the oxidization layer of the methane/air flames. Overall, the chemical reactions of all biomass derived gas flames occur in thin zones of the order of less than 1 mm. The thickness of all BDG gas flames is inversely proportional to their respective laminar burning velocity. The laminar burning velocities of landfill gases are found to increase linearly with the mole fraction of methane in the mixtures, whereas for gasification gas, syngas and pyrolysis gas where hydrogen is present, the laminar burning velocities scale linearly with the mole fraction of hydrogen.     

Laminar burning velocities and flame instabilities of 2,5-dimethylfuran–air mixtures at elevated pressures

An experimental investigation on laminar burning velocities and onset of flame instabilities on spherically expanding flames in 2,5-dimethylfuran–air mixtures at elevated pressures was conducted over a wide range of equivalence ratios. The laminar burning velocities, laminar burning fluxes and Markstein lengths at different equivalence ratios and initial pressures were obtained. Furthermore, the diffusional–thermal and hydrodynamic effects on flame front instabilities were specified, and the onset of cellularity was reported. Results show that laminar burning velocities are decreased with increasing initial pressure due to the increase of the free-stream density and the progressively more important three-body termination reactions. With increasing initial pressure, Markstein length decrease, while the laminar burning flux increases. Onsets of flame instabilities, expressed in terms of critical radius or Peclet number, were found to be promoted with increasing equivalence ratio and initial pressures, due to the combined effects of diffusional–thermal and hydrodynamic instabilities.     

Transient flame spread during convective ignition of solid fuel in a sudden-expansion combustor

The convective ignition of solid fuel (PMMA) in a sudden-expansion combustor (Re_h = 6200, u₀ = 22 m/s, [O₂] ～ 11.7%, T₀ = 810 °C) is investigated from the perspective of flame–vortex interactions. Three phases of the transient flame spread are identified via the diagnostics of flow visualization and particle image velocimetry (PIV). The dominance of small/large vortices is revealed, respectively, in the pre-/post-ignition regimes, which demonstrates the small-to-large vortex transformation due to heat release. Attributed to the decreased characteristic reaction time and enhanced mixing, the first ignition is observed at the downstream end of the fuel, after which a primitive flame is formed and initiates the opposed flame spread. During the spread, the rolling behavior of flame kernels are considered to be dominated by the small eddies. The combined effects of broken vortices and continuing pyrolysis introduce the periodical extinction–reignition around the reattachment region. At the final phase, the entrainment of flame kernels into the shear layer is facilitated by the large shedding vortices, and a sustained diffusion flame is established. The study not only provides novel insights into the convective ignition of solid fuel in the separation–reattachment flow, but also serves as a basis for the advancement of ignition control.     

Effect of preferential diffusion and flame stretch on flame structure and laminar burning velocity of syngas Bunsen flame using OH-PLIF

By using OH-PLIF technique, experiments were conducted for laminar Bunsen flame of premixed CO/H₂/air mixtures with equivalence ratio ranging from 0.5 to 1.8. Reynolds number was varied from 800 to 2200, X_H2 = H₂/(H₂+CO) in the mixture was varied from 20% to 100% to study the effects of both preferential diffusion and flame curvature on flame structures and laminar flame burning velocities. Results showed that the combined effects of preferential diffusion and curvature gave an interesting phenomenon of the flame OH radical distributions on high hydrogen content flames. Furthermore, with the increase of H₂ fraction in fuel mixture, the effects of both preferential diffusion and flame curvature were increased. Interpretation of flame stretch effect on laminar burning velocity is also provided in this paper.