Flame acceleration and DDT in large-scale obstructed channels filled with methane-air mixtures

We study the evolution of premixed methane-air flames in large-scale obstructed channels using reactive CFD simulations. We vary the channel height d, the blockage ratio br, and the scaled distance between obstacles L/d to study their effects on the distance to DDT, L_DDT, and the distance to the shock-flame complex, L_SF. The results of simulations show two main effects. On one hand, the increase of br and decrease of L/d promote the flame acceleration and reduce L_SF and L_DDT. On the other hand, some configurations with higher br and smaller L/d prevent the detonation development. As a result, the leading shock and the flame never merge, and continue to propagate as a quasi-steady-state shock-flame complex. A collision of this complex with solid structures generates high pressures and strong reflected shocks that can ignite a detonation. This detonation would propagate in a shock-compressed material and result in extremely high pressures exceeding pressures of a regular detonation. Thus, the distance L_SF at which the shock-flame complex forms provides an important measure of a destructive potential in addition to L_DDT. This is particularly relevant for channels with high br where L_SF can be significantly shorter than L_DDT.     

Characterization of the effects of pressure and hydrogen concentration on laminar burning velocities of methane–hydrogen–air mixtures

The aim of the present work was to characterize both the effects of pressure and of hydrogen addition on methane/air premixed laminar flames. The experimental setup consists of a spherical combustion chamber coupled to a classical shadowgraphy system. Flame pictures are recorded by a high speed camera. Global equivalence ratios were varied from 0.7 to 1.2 for the initial pressure range from 0.1 to 0.5 MPa. The mole fraction of hydrogen in the methane + hydrogen mixture was varied from 0 to 0.2. Experimental results were compared to calculations using a detailed chemical kinetic scheme (GRIMECH 3.0). First, the results for atmospheric laminar CH₄/air flames were compared to the literature. Very good agreements were obtained both for laminar burning velocities and for burned gas Markstein length. Then, increasing the hydrogen content in the mixture was found to be responsible for an increase in the laminar burning velocity and for a reduction of the flame dependence on stretch. Transport effects, through the reduction of the fuel Lewis number, play a role in reducing the sensitivity of the fundamental flame velocity to the stretch. Finally, when the pressure was increased, the laminar burning velocity decreased for all mixtures. The pressure domain was limited to 0.5 MPa due to the onset of instabilities at pressures above this value.     

Laminar and unstable burning velocities and Markstein lengths of hydrogen–air mixtures at engine-like conditions

Hydrogen offers an attractive alternative to conventional fuels for use in spark ignition engines. It can be burned over a very wide range of equivalence ratios and with considerable exhaust gas recirculation. These help to minimise pumping losses through throttleless operation and oxides of nitrogen (NO_x) production through reduced temperature. Full understanding of hydrogen-fuelled engine operation requires data on the laminar burning rate of hydrogen–air residuals under a wide range of conditions. However, such data are sparse. The present work addresses this need for experimental data. Spherically expanding H₂–air flames were measured at a range of temperatures, pressures, and equivalence ratios and with varying concentrations of residuals of combustion. Unstretched burning velocities, u_l, and Markstein lengths, L_b, were determined from stable flames. At the higher pressures, hydrodynamic and diffusional-thermal instabilities caused the flames to be cellular from inception and prohibited the derivation of values of u_l and L_b. The effect of pressure on the burning rate was demonstrated to have opposing trends when comparing stoichiometric and lean mixtures. The present measurements were compared with those available in the literature, and discrepancies were attributed to neglect, in some works, the effects of stretch and instabilities. From the present measurements, the effects of pressure, temperature, and residual gas concentration on burning velocity are quantified for use in a first step towards a general correlation.     

Premixed flame oscillations in obstructed channels with both ends open

An initially laminar premixed flame front accelerates extremely fast and may even trigger a detonation when propagating in a semi-open obstructed channel (one end of the channel is closed; the flame is ignited at the closed end and moves towards the open one). However, industrial and laboratory conduits oftentimes have both ends open, or vented, with a flame ignited at one of these ends. The latter constitutes the focus of the present work. Specifically, premixed flame propagation through a comb-shaped array of obstacles, in-built in a channel with both ends open, is studied by means of computational simulation of the reacting flow equations with fully-compressible hydrodynamics and an Arrhenius chemical kinetics. The parametric study includes various blockage ratios and spacing as well as the thermal expansion ratios, with oscillations of the burning rate observed in the majority of the cases, which conceptually differs from fast flame acceleration in semi-open channels. Such a difference is devoted to the fact that while the entire flame-generated jet-flow is pushed towards a single exit in a semi-open channel, in a channel with two ends open, this jet-flow is distributed between the upstream and downstream flows, thereby moderating flame propagation. The flame oscillations are nonlinear in all cases where they are observed. The oscillation period grows with the blockage ratio but decreases with the thermal expansion. The present results also support the recent experiments, modeling and theory of flames in obstructed channels with both ends open, which all yielded steady or quasi-steady flame propagation prior to the onset of flame acceleration. Indeed, the present oscillations can be treated as the fluctuations around a quasi-steady solution.     

Numerical simulations of the flow field ahead of an accelerating flame in an obstructed channel

The development of the unburned gas flow field ahead of a flame front in an obstructed channel was investigated using large eddy simulation (LES). The standard Smagorinsky–Lilly and dynamic Smagorinsky–Lilly subgrid models were used in these simulations. The geometry is essentially two-dimensional. The fence-type obstacles were placed on the top and bottom surfaces of a square cross-section channel, equally spaced along the channel length at the channel height. The laminar rollup of a vortex downstream of each obstacle, transition to turbulence, and growth of a recirculation zone between consecutive obstacles were observed in the simulations. By restricting the simulations to the early stages of the flame acceleration and by varying the domain width and domain length, the three-dimensionality of the vortex rollup process was investigated. It was found that initially the rollup process was two-dimensional and unaffected by the domain length and width. As the recirculation zone grew to fill the streamwise gap between obstacles, the length and width of the computational domain started to affect the simulation results. Three-dimensional flow structures formed within the shear layer, which was generated near the obstacle tips, and the core flow was affected by large-scale turbulence. The simulation predictions were compared to experimental schlieren images of the convection of helium tracer. The development of recirculation zones resulted in the formation of contraction and expansion regions near the obstacles, which significantly affected the centerline gas velocity. Oscillations in the centerline unburned gas velocity were found to be the dominate cause for the experimentally observed early flame-tip velocity oscillations. At later simulation times, regular oscillations in the unburned streamwise gas velocity were not observed, which is contrary to the experimental evidence. This suggests that fluctuations in the burning rate might be the source of the late flame-tip velocity oscillations. The effect of the obstacle blockage ratio (BR) on the development of the unburned gas flow field was also investigated by varying the obstacle height. Simulation predictions show favorable agreement with the experimental results and indicate that turbulence production increases with increasing obstacle BR.     

Turbulent flame propagation limits of ammonia/methane/air premixed mixture in a constant volume vessel

Ammonia is one of promising energy carriers that can be directly used as carbon-neutral fuel for combustion applications. However, because of the low-burning velocity of ammonia, it is challenging to introduce ammonia to practical combustors those are designed for general hydrocarbon fuels. One of ways to enhance the combustibility of ammonia is by mixing it with other hydrocarbon fuels, such as methane, with a burning velocity is much higher than the burning velocity of ammonia. In this study, we conducted flame propagation experiments of ammonia/methane/air using a fan-stirred constant volume vessel to clarify the effect of methane addition to ammonia on the turbulent flame propagation limit. From experimental results, we constructed the flame propagation maps and clarified the flame propagation limits. The results show that the flame propagation limits were extended with an increase in mixing a fraction of methane to ammonia. Additionally, ammonia/methane/air mixtures with the equivalence ration of 0.9 can propagate at the highest turbulent intensity, even though the peak of the laminar burning velocity is the fuel-rich side because of the diffusional-thermal instability of the flame surface. Furthermore, the Markstein number of the mixture obtained in this research successfully expressed the strength of the diffusional-thermal instability effect on the flame propagation capability. The turbulence Karlovitz number at the flame propagation limit monotonically increases with the decreasing Markstein number.     

Combined experimental and computational study of laminar, axisymmetric hydrogen–air diffusion flames

We investigate the structure of two-dimensional, axisymmetric, laminar hydrogen–air flames in which a cylindrical fuel stream is surrounded by coflowing air, using laser-diagnostic and computational methods. Spontaneous Raman scattering and coherent anti-Stokes Raman scattering (CARS) are used to measure the distributions of major species and temperature. Computationally, we solve the governing conservation equations for mass, momentum, energy, and species, using detailed chemistry and transport. The fuel is diluted with nitrogen (1:1) to reduce heat transfer to the burner, to match the zero temperature gradient at the fuel exit. Three average fuel exit velocities are studied: 18, 27, and 50 cm/s. Comparisons of the measured and computed results are performed for radial profiles at a number of axial positions, and along the axial centerline. Peak major species mole fractions and temperatures are quantitatively predicted by the computations, and the axial species profiles are predicted to within the experimental uncertainty. In the radial profiles studied, base-case computations excluding thermal diffusion of light species were in excellent agreement with the measurements. While the addition of thermal diffusion led to some discrepancy with the measured results, the magnitude of the differences was no more than 25%. The computations predicted the axial centerline profiles from the burner exit to the maximum temperature well, though the experimental temperatures in the downstream mixing region decreased somewhat faster than the computed profiles. Radiative losses are seen to be negligible in these flames, and changes in transport properties and variations in initial flow velocities generally led to only modest changes in the axial profiles. The results also show that the detailed axial profiles of major species and temperature at different fuel jet velocities scale quantitatively with the jet velocity.     

WIDECARS spectra fitting in a premixed ethylene–air flame

WIDECARS measures temperature and mole fractions of most of the major species in ethylene–air flames. One of the issues in implementing this technique is fitting the experimental spectra to theory to obtain flame conditions (temperature, species mole fractions). Individual spectra contain many species resonances, and theory is slow to compute. Libraries of precalculated spectra can be used, but a library of sufficient density for accurate interpolation is large given the many variables. A new fitting algorithm is presented which utilizes a less‐dense library, and additional spectra are calculated during fitting to maintain accuracy. The iterative convergence method converts the problem of minimizing fit error, which converges slowly, to a zero finding problem, which converges reliably, rapidly, and accurately to best fit. Various practical fitting issues, such as the effects of dye laser mode noise and variability, phase‐matching efficiency, and shifts of the spectrum on the spectrometer are addressed. The technique is demonstrated in the analysis of experimental measurements in an equivalence ratio 2.1 ethylene–air flame above the surface of a McKenna burner. Precision errors because of experimental and fitting effects are discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

A computational analysis of strained laminar flame propagation in a stratified CH4/H2/air mixture

Propagation of a H₂-added strained laminar CH₄/air flame in a rich-to-lean stratified mixture is numerically studied. The back-support effect, which is known to enhance the consumption speed of a flame propagating into a leaner mixture compared to that into a homogeneous mixture, is evaluated. A new method is devised to characterize unsteady reactant-to-reactant counterflow flames under transiently decreasing equivalence ratio, in order to elucidate the influence of flow strain on the back-support effect. In contrast to the conventional reactant-to-product configurations, the current configuration is more relevant to unsteady stratified flames back-supported by their own combustion products. Moreover, since H₂ distribution downstream of the flame is known to play a crucial role in back-supported CH₄/air flames, the influence of H₂ addition in the upstream mixture is examined. The results suggest that a larger strain rate leads to a larger equivalence ratio gradient at the reaction zone through increased flow divergence, which amplifies the back-support. Meanwhile, since H₂ addition in the upstream mixture does not affect the downstream H₂ content, the relative increase in the consumption speed, i.e. the back-support, is suppressed with larger H₂ addition. Especially, when the upstream H₂ content decreases with the equivalence ratio, the H₂ preferentially diffuses toward the unburned gas, which mitigates H₂ accumulation in the preheat zone and further weakens the back-support.     

High dynamic range thermometry at 5 kHz in hydrogen–air diffusion flame using chirped‐probe‐pulse femtosecond coherent anti‐stokes Raman scattering

Chirped probe pulse femtosecond coherent anti‐Stokes Raman scattering (CPP fs‐CARS) thermometry was performed at 5 kHz in a hydrogen jet diffusion flame with an air co‐flow. Measurements were performed at different heights and radial locations within the jet diffusion flame, up to 16 nozzle exit diameters downstream (x/d = 16). The near‐nozzle measurements were characterized by large, organized, buoyancy‐driven instabilities that become more chaotic at the downstream locations x/d ≥ 4. The diffusion flame results highlight temperature fluctuations characteristic of the buoyancy‐driven Kelvin–Helmholtz‐type instability and provide new insights into the transient structure of these flames. At some measurement locations, the time‐varying temperatures ranged from 300 K to nearly 2400 K. The CPP fs‐CARS signal intensity is a factor of approximately 1000 times lower at 2400 K compared with 300 K. A dual‐channel detection system was used to increase the dynamic range of the CARS measurements. The determination of temperature from the single shot spectra is discussed in detail. Laser and detection system parameters were determined from CPP fs‐CARS spectra obtained from a near‐adiabatic laminar calibration flame apparatus. The temperature precision of the system was determined from these calibration measurements and was found to be better than 2.0% at 2200 K. The influence of an instrument response function on spectral fitting parameters is systematically assessed. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.     

Temperature and carbon source effects on methane–air flame synthesis of CNTs

We have conducted experimental and numerical studies on flame synthesis of carbon nanotubes (CNTs) to investigate the effects of three key parameters – selective catalyst, temperature and available carbon sources – on CNT growth. Two different substrates were used to synthesize CNTs: Ni-alloy wire substrates to obtain curved and entangled CNTs and Si-substrates with porous anodic aluminum oxide (AAO) nanotemplates to grow well-aligned, self-assembled and size-controllable CNTs, each using two different types of laminar flames, co-flow and counter-flow methane–air diffusion flames. An appropriate temperature range in the synthesis region is essential for CNTs to grow on the substrates. Possible carbon sources for CNT growth were found to be the major species CO and those intermediate species C₂H₂, C₂H₄, C₂H₆, and methyl radical CH₃. The major species H₂, CO₂ and H₂O in the synthesis region are expected to activate the catalyst and help to promote catalyst reaction.     

Ignition of a sequential combustor: Evidence of flame propagation in the autoignitable mixture

Ignition of the second stage in a lab-scale sequential combustor is investigated experimentally. A fuel mixing section between jet-in-cross-flow injection and the second stage chamber allows the fuel and vitiated, hot cross-flow to partially mix upstream of the main heat release zone. The focus of the present work is on the transient ignition process leading to a stable flame in the second stage. High-speed OH-PLIF as well as OH chemiluminescence imaging is applied to obtain complementary planar and line-of-sight integrated information on the ignition. We find experimental evidence for the co-existence of two regimes dominating the chamber ignition, i.e. autoignition and flame propagation. As the mass flow of the dilution air injected downstream of the first stage is increased (i.e. mixing temperatures in the fuel mixing section are decreased), we transition from an autoignition to a flame propagation dominated regime. Hysteresis in the ignition behavior is observed indicating that the first stage in a sequential combustor may be operated at leaner conditions than required for ignition of the second stage. The time traces of integral heat release obtained simultaneously with a photomultiplier tube show distinct features depending on the dominating regime, which is important for high-pressure testing with limited optical access.     

基于OpenMP的飞秒强激光在空气中传输并行计算

轴对称超短强激光在空气中的传输可以用2D+1维非线性薛定谔方程来描述,该方程一般可用FCN方法进行求解,即在时间上应用快速傅里叶变换方法、在空间横截面上采用Crank-Nicholson差分法。但由于计算非常耗时,传统串行的FCN方法只能计算强激光在百m量级距离上传输。基于OpenMP设计了求解2D+1维非线性薛定谔方程的并行方法及其程序。数值模拟结果证明了程序的正确性和较高的并行效率。当线程数为15时并行加速比为12。此并行方法可应用于模拟长距离的超短超强激光的传输。     

Wave propagation in a waveguide with continuous right-angled corners: Numerical simulations and experiment measurements

This paper studied the acoustic wave propagation in a waveguide with continuous right-angled corners, with emphasis on the effect brought by the distance between the corners. The numerical analyses showed that at middle to high frequencies, the transmission loss (TL) of a multi-cornered waveguide was 2–5 dB higher than that of single-cornered and varied with frequency. To explain the performance at peaks and dips in the TL curve, analyses on eigenmodes and sound intensity distribution were conducted. The performance of multi-cornered waveguides was experimentally investigated, which fit well with the numerical results. The present study indicates that, for a waveguide with continuous corners, its acoustic performance is not simply a “summation” of two individual single-cornered ones. Both the standing wave modes and the evanescent modes between the corners lead to its complicated frequency performance.     

Numerical simulations of the nonlinear propagation of femtosecond optical pulses in gases

A two-dimensional axisymmetric model of the propagation of intense femtosecond laser pulses through dispersion-free transparent media is described. The effects of diffraction, nonlinear Kerr effect (instantaneous and retarded) and multiphoton ionisation are included. Numerical results concerning air and other gases are discussed. In particular, time self-compression of femtosecond pulses is predicted. Stable self-guided pulses are simulated, in agreement with recent experimental observations. Received: 19 June 1998 / Received in final form: 14 January 1999     

Comments concerning “Study of solute–solvent and solvent–solvent interactions in pure and mixed binary solvents”

A number of mathematical equations representing the solute solubility in monosolvent and binary solvent mixtures are discussed. This work is a commentary to the article of Maitra and Baghchi that appeared in this journal.     

The effect of air exposure on palladium–copper composite membranes

It was found that when electrolessly deposited thin Pd and Pd–Cu membranes were exposed to air at temperatures above 350 °C, their H₂ flux increased substantially immediately after the air exposure, then decreased to a new steady-state value. While this was a quasi-reversible change for the H₂ flux, the flux of insoluble species, such as N₂, irreversibly increased with every air exposure but by a much smaller extent. The extent of these changes was found to be dependent on the exposure time and the temperature of the tests. Thus, we decided to investigate the effect of gas exposures on the properties of these materials.

Palladium and palladium–copper films, prepared by electroless deposition on ceramic supports, and commercial foils were exposed to air, hydrogen and helium at 500 and 900 °C for times varying from 1 h to 1 week with the objective of determining the effect of the different exposure conditions on the surface morphology, the flux of different penetrants and the crystalline structure of the materials. Atomic force microscopy (AFM) and X-ray diffraction (XRD) were used to study the changes occurring in the films under those conditions.

It was observed that the exposure of both the electroless films and the foils to hydrogen and air markedly modified their surface morphology. The hydrogen exposure tended to smooth the surface features whereas the oxygen exposure created new surface features such holes and large peaks. Additionally it was found that the air exposure produced some oxidation of the film to create PdO.

These results suggested that a common hypothesis stating that air oxidation just cleans the surface of the membrane might not be sufficient to explain all of those changes. A contributing effect of air exposure may be the increase in surface area due to the formation of palladium oxide. However, the extent of the surface area increase was insufficient to explain the increase in steady-state H₂ flux.