Influence of a combustion-driven oscillation on global mixing in the flame from a refinery flare

An assessment of the influence of strong combustion-driven oscillations on mixing rates and visible radiation in the flame from a full-scale refinery flare is reported. Importantly, the oscillations were generated naturally, with no external forcing, and at a high Reynolds number of 4 × 10⁶. These conditions differentiate this study from those of previous investigations, which all involved some external forcing and were at a Re too low to ensure fully turbulent flow within the flame. A frame-by-frame analysis of video footage, providing good resolution of the instantaneous edge of each flame, was used to assess flame dimensions, and so to determine a global residence time. Since the flames are in the fast-chemistry regime, the visual imagers can be used to determine a global mixing rate. The analysis reveals a consistent picture that the combustion-driven oscillations do not result in a significant change to the global mixing rate, but do increase the visible radiation. This is in contrast to previous investigations, using externally forced jets, where forcing at the preferred mode has been found to increase mixing rates and reduce radiation.     

Three-dimensional direct simulations and structure of expanding turbulent methane flames

Direct numerical simulations (DNS) are ideally suited to investigate in detail turbulent reacting flows in simple geometries. For an increasing number of applications, detailed models must be employed to describe the chemical processes with sufficient accuracy. Despite the huge cost of such simulations, recent progress has allowed the direct numerical simulation of turbulent premixed flames while employing complete reaction schemes. We briefly describe our own developments in this field and use the resulting DNS code to investigate more extensively the structure of premixed methane flames expanding in a three-dimensional turbulent velocity field, initially homogeneous and isotropic. This situation typifies, for example, the initial flame development after spark ignition in a gas turbine or an internal combustion engine. First investigation steps have been carried out at low turbulence levels on this same configuration in the past Symposium, and we build on top of these former results. Here, a considerably higher Reynolds number is considered, the simulation has been repeated twice in to limit the possibility of spurious, very specific results, and several complementary post-processing steps are carried out. Characteristic features concerning the observed combustion regime are presented. We then investigate in a quantitative manner the evolution of flame surface area, global stretch-rate, flame front curvature, flame thickness, and correlation between thickness and curvature. The possibility of obtaining reliable information on flame front curvature from two-dimensional slices is checked by comparison with the exact procedure.     

Interaction of a premixed flame with a liquid fuel film on a wall

In piston engines and in gas turbines, the injection of liquid fuel often leads to the formation of a liquid film on the combustor wall. If a flame reaches this zone, undesired phenomena such as coking may occur and diminish the lifetime of the engine. Moreover, the effect of such an interaction on maximum wall heat fluxes, flame quenching, and pollutant formation is largely unknown. This paper presents a numerical study of the interaction of a premixed flame with a cold wall covered with a film of liquid fuel. Simulations show that the presence of the film leads to a very rich zone at the wall in which the flame cannot propagate. As a result, the flame wall distance remains larger with liquid fuel than it is for a dry wall, and maximum heat fluxes are smaller. The nature of the interaction of flame wall interaction with a liquid fuel is also different from the classical flame/dry wall interaction: it is controlled mainly by chemical mechanisms and not by the thermal quenching effect observed for flames interacting with dry walls: the existence of a very rich zone created above the liquid film is the main mechanism controlling quenching.     

High‐order upwind compact scheme and simulation of turbulent premixed V‐flame

A high‐order accurate upwind compact difference scheme with an optimal control coefficient is developed to track the flame front of a premixed V‐flame. In multi‐dimensional problems, dispersion effect appears in the form of anisotropy. By means of Fourier analysis of the operators, anisotropic effects of the upwind compact difference schemes are analysed. Based on a level set algorithm with the effect of exothermicity and baroclinicity, the flame front is tracked. The high‐order accurate upwind compact scheme is employed to approximate the level set equation. In order to suppress numerical oscillations, the group velocity control technique is used and the upwind compact difference scheme is combined with the random vortex method to simulate the turbulent premixed V‐flame. Distributions of velocities and flame brush thickness are obtained by this technique and found to be comparable with experimental measurement. Copyright © 2005 John Wiley & Sons, Ltd.     

Cellular instabilities of expanding hydrogen/propane spherical flames at elevated pressures: theory and experiment

An experimental and theoretical investigation of the onset of cellular instabilities on spherically expanding flames in mixtures of hydrogen and propane in air at elevated pressures was conducted. Critical conditions for the onset of instability were measured and mapped out over a range of pressures and mixture compositions. An asymptotic theory of hydrodynamic and diffusional-thermal cell development on flames in mixtures comprised of two scarce fuels burning in air was also formulated. Predicted values of Peclet number, defined as the flame radius at the onset of instability normalized by the flame thickness, were shown to compare favorably with the experimentally measured values.     

蜡烛火焰三维温度场重建实验研究

本文在实验室中以蜡烛火焰为研究对象,研究采用辐射图像处理技术重建蜡烛三维温度场。并采用细丝热电偶实测火焰局部温度,对蜡烛火焰三维温度场可视化的结果进行了比较。结果表明,火焰高温区的重建温度是合理的,此方法能够很好地重建单峰、双峰三维温度分布。进一步证实了基于火焰图像处理技术实现三维温度场可视化的可行性。     

The effect of hydrogen addition on flammability limit and NOx emission in ultra-lean counterflow CH4/air premixed flames

The effect of hydrogen addition to ultra lean counterflow CH₄/air premixed flames on the extinction limits and the characteristics of NO_x emission was investigated by numerical simulation. Detailed chemistry and complex thermal and transport properties were employed. The results show that the addition of hydrogen can significantly enlarge the flammable region and extend the flammability limit to lower equivalence ratios. If the equivalence ratio is kept constant, the addition of hydrogen increases the emission of NO in a flame due to the enhancement in the rate of the NNH or N₂O intermediate NO formation routes. The addition of hydrogen causes a monotonic decrease in the formation of NO₂ and N₂O, except flames near the extinction limits, where the emission of NO₂ and N₂O first increases, and then decreases with the increase in the fraction of hydrogen. Overall, hydrogen enrichment technology allows stable combustion under ultra lean conditions, resulting in significant CO₂ and NO emission reduction.     

A segregated implicit solution algorithm for 2D and 3D laminar incompressible flows

A segregated algorithm for the solution of laminar incompressible, two- and three-dimensional flow problems is presented. This algorithm employs the successive solution of the momentum and continuity equations by means of a decoupled implicit solution method. The inversion of the coefficient matrix which is common for all momentum equations is carried out through an approximate factorization in upper and lower triangular matrices. The divergence-free velocity constraint is satisfied by formulating and solving a pressure correction equation. For the latter a combined application of a preconditioning technique and a Krylov subspace method is employed and proved more effecient than the approximate factorization method. The method exhibits a monotonic convergence, it is not costly in CPU time per iteration and provides accurate solutions which are independent of the underrelaxation parameter used in the momentum equations. Results are presented in two- and three-dimensional flow problems.     

Flame inhibition by phosphorus-containing compounds in lean and rich propane flames

Chemical inhibition of laminar propane flames by organophosphorus compounds has been studied experimentally and computationally using a detailed chemical kinetic reaction mechanism. Both fuel-lean and fuel-rich propane flames were studied to examine the role of equivalence ratio in flame inhibition. The experiments examined a wide variety of organophosphorus compounds. We report on experimental species flame profiles for tri-methyl phosphate (TMP) and compare them with modeled species flame profile results of TMP and di-methyl methyl phosphonate (DMMP). Both experiments and kinetic modeling indicate that inhibition efficiency is effectively the same for all of the organophosphorus compounds examined, independent of the molecular structure of the initial inhibitor molecule. Chemical inhibition is due to reactions involving small P-bearing species HOPO₂ and HOPO produced by the organophosphorus compounds (OPCs). Ratios of HOPO₂ and HOPO concentrations differ between lean and rich flames, with HOPO₂ dominant in lean flames while HOPO dominates in rich flames. Resulting HOPO₂ and HOPO species profiles do not significantly depend on the initial source of the HOPO₂ and HOPO, and thus are relatively insensitive to the initial OPC inhibitor. A more generalized form of the Twarowski mechanism is developed to account for the results observed, and new theoretical values are determined for heats of formation of the important P-containing species, using the BAC-G2 method.