Flammability conditions for ultra-lean hydrogen premixed combustion based on flame-ball analyses

It has been reasoned that the structures of strongly cellular flames in very lean mixtures approach an array of flame balls, each burning as if it were isolated, thereby indicating a connection between the critical conditions required for existence of steady flame balls and those necessary for occurrence of self-sustained premixed combustion. This is the starting assumption of the present study, in which structures of near-limit steady sphericosym-metrical flame balls are investigated with the objective of providing analytic expressions for critical combustion conditions in ultra-lean hydrogen-oxygen mixtures diluted with N2 and water vapor. If attention were restricted to planar premixed flames, then the lean-limit mole fraction of H2 would be found to be roughly ten percent, more than twice the observed flammability limits, thereby emphasizing the relevance of the flame-ball phenomena. Numerical integrations using detailed models for chemistry and radiation show that a onestep chemical-kinetic reduced mechanism based on steady-state assumptions for all chemical intermediates, together with a simple, optically thin approximation for water-vapor radiation, can be used to compute near-limit fuel-lean flame balls with excellent accuracy. The previously developed one-step reaction rate includes a crossover temperature that determines in the first approximation a chemical-kinetic lean limit below which combustión cannot occur, with critical conditions achieved when the diffusion-controlled radiation-free peak temperature, computed with account taken of hydrogen Soret diffusion, is equal to the crossover temperature. First-order corrections are found by activation-energy asymptotics in a solution that involves a near-field radiation-free zone surrounding a spherical flame sheet, together with a far-field radiation-conduction balance for the temperature profile. Different scalings are found depending on whether or not the surrounding atmosphere contains wáter vapor, leading to different analytic expressions for the critical conditions for flame-ball existence, which give results in very good agreement with those obtained by detailed numerical computations.

Soot measurements in steady and pulsed ethylene/air diffusion flames using Laser-Induced Incandescence

Combustion-generated carbon black nano particles, or soot, have both positive and negative effects depending on the application. From a positive point of view, it is used as a reinforcing agent in tires, black pigment in inks, and surface coatings. From a negative point of view, it affects performance and durability of many combustion systems, it is a major contributor of global warming, and it is linked to respiratory illness and cancer. Laser-Induced Incandescence (LII) was used in this study to measure soot volume fractions in four steady and twenty-eight pulsed ethylene diffusion flames burning at atmospheric pressure. A laminar coflow diffusion burner combined with a very-high-speed solenoid valve and control circuit provided unsteady flows by forcing the fuel flow with frequencies between 10 Hz and 200 Hz. Periodic flame oscillations were captured by two-dimensional phase-locked LII images and broadband luminosity images for eight phases (0° – 360°) covering each period. A comparison between the steady and pulsed flames and the effect of the pulsation frequency on soot volume fraction in the flame region and the post flame region are presented. The most significant effect of pulsing frequency was observed at 10 Hz. At this frequency, the flame with the lowest mean flow rate had 1.77 times enhancement in peak soot volume fraction and 1.2 times enhancement in total soot volume fraction; whereas the flame with the highest mean flow rate had no significant change in the peak soot volume fraction and 1.4 times reduction in the total soot volume fraction. A correlation (fvRe-1 = a + b·Str) for the total soot volume fraction in the flame region for the unsteady laminar ethylene flames was obtained for the pulsation frequency between 10 Hz and 200 Hz, and the Reynolds number between 37 and 55. The soot primary particle size in steady and unsteady flames was measured using the Time-Resolved Laser-Induced Incandescence (TIRE-LII) and the double-exponential fit method. At maximum frequency (200 Hz), the soot particles were smaller in size by 15% compared to the steady case in the flame with the highest mean flow rate.

Asymmetric hydrogen flame in a heated micro-channel: role of Darrieus-Landau and thermal-diffusive instabilities

Present work examines numerically the asymmetric behavior of hydrogen/air flame in a micro-channel subjected to a non-uniform wall temperature distribution. A high resolution (with cell size of 25 μm × 25 μm) of two-dimensional transient Navier–Stokes simulation is conducted in the low-Mach number formulation using detailed chemistry evolving 9 chemical species and 21 elementary reactions. Firstly, effects of hydrodynamic and diffusive-thermal instabilities are studied by performing the computations for different Lewis numbers. Then, the effects of preferential diffusion of heat and mass transfer on the asymmetric behavior of the hydrogen flame are analyzed for different inlet velocities and equivalence ratios. Results show that for the flames in micro-channels, interactions between thermal diffusion and molecular diffusion play major role in evolution of a symmetric flame into an asymmetric one. Furthermore, the role of Darrieus–Landau instability found to be minor. It is also found that in symmetric flames, the Lewis number decreases behind the flame front. This is related to the curvature of flame which leads to the inclination of thermal and mass fluxes. The mass diffusion vectors point toward the walls and the thermal diffusion vectors point toward the centerline. Asymmetric flame is observed when the length of flame front is about 1.1–1.15 times of the channel width.

Modeling of Turbulent Natural Gas and Biogas Flames of the Delft Jet-in-Hot-Coflow Burner: Effects of Coflow Temperature, Fuel Temperature and Fuel Composition on the Flame Lift-Off Height

The structure of a turbulent non-premixed flame of a biogas fuel in a hot and diluted coflow mimicking moderate and intense low dilution (MILD) combustion is studied numerically. Biogas fuel is obtained by dilution of Dutch natural gas (DNG) with CO2. The results of biogas combustion are compared with those of DNG combustion in the Delft Jet-in-Hot-Coflow (DJHC) burner. New experimental measurements of lift-off height and of velocity and temperature statistics have been made to provide a database for evaluating the capability of numerical methods in predicting the flame structure. Compared to the lift-off height of the DNG flame, addition of 30 % carbon dioxide to the fuel increases the lift-off height by less than 15 %. Numerical simulations are conducted by solving the RANS equations using Reynolds stress model (RSM) as turbulence model in combination with EDC (Eddy Dissipation Concept) and transported probability density function (PDF) as turbulence-chemistry interaction models. The DRM19 reduced mechanism is used as chemical kinetics with the EDC model. A tabulated chemistry model based on the Flamelet Generated Manifold (FGM) is adopted in the PDF method. The table describes a non-adiabatic three stream mixing problem between fuel, coflow and ambient air based on igniting counterflow diffusion flamelets. The results show that the EDC/DRM19 and PDF/FGM models predict the experimentally observed decreasing trend of lift-off height with increase of the coflow temperature. Although more detailed chemistry is used with EDC, the temperature fluctuations at the coflow inlet (approximately 100K) cannot be included resulting in a significant overprediction of the flame temperature. Only the PDF modeling results with temperature fluctuations predict the correct mean temperature profiles of the biogas case and compare well with the experimental temperature distributions.

Interaction of Acoustic Waves with a Laminar Line-Flame

A systematic study was conducted to elucidate the effects of acoustic perturbations on laminar diffusion line-flames and the conditions required to cause acoustically-driven extinction. Flames were produced from the fuels n-pentane, n-hexane, n-heptane, n-octane, and JP-8, using fuel-laden wicks. The wicks were housed inside of a burner whose geometry produced flames that approximated a two dimensional flame sheet. The acoustics utilized ranged in frequency between 30-50 Hz and acoustic pressures between 5-50 Pa. The unperturbed mass loss rate and flame height of the alkanes were studied, and they were found to scale in a linear manner consistent with Burke-Schumann. The mass loss rate of hexane-fueled flames experiencing acoustic perturbations was then studied. It was found that the strongest influence on the mass loss rate was the magnitude of oscillatory air movement experienced by the flame. Finally, acoustic perturbations were imposed on flames using all fuels to determine acoustic extinction criterion. Using the data collected, a model was developed which characterized the acoustic conditions required to cause flame extinction. The model was based on the ratio of an acoustic Nusselt Number to the Spalding B Number of the fuel, and it was found that at the minimum speaker power required to cause extinction this ratio was a constant. Furthermore, it was found that at conditions where the ratio was below this constant, a flame could still exist; at conditions where the ratio was greater than or equal to this constant, flame extinction always occurred.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.

A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.

Implicit difference approximation for the time fractional diffusion equation

In this paper, we consider a time fractional diffusion equation on a finite domain. The equation is obtained from the standard diffusion equation by replacing the first-order time derivative by a fractional derivative (of order $0<\alpha<1$ ). We propose a computationally effective implicit difference approximation to solve the time fractional diffusion equation. Stability and convergence of the method are discussed. We prove that the implicit difference approximation (IDA) is unconditionally stable, and the IDA is convergent with $O(\tau+h^2)$, where $\tau$ and $h$ are time and space steps, respectively. Some numerical examples are presented to show the application of the present technique.