Nonlinear thermoacoustics of ducted premixed flames: The influence of perturbation convection speed

When a premixed flame is placed within a duct, acoustic waves induce velocity perturbations at the flame's base. These travel down the flame, distorting its surface and modulating its heat release. This can induce self-sustained thermoacoustic oscillations. Although the phase speed of these perturbations is often assumed to equal the mean flow speed, experiments conducted in other studies and Direct Numerical Simulation (DNS) conducted in this study show that it varies with the acoustic frequency. In this paper, we examine how these variations affect the nonlinear thermoacoustic behaviour. We model the heat release with a nonlinear kinematic G-equation, in which the velocity perturbation is modelled on DNS results. The acoustics are governed by linearised momentum and energy equations. We calculate the flame describing function (FDF) using harmonic forcing at several frequencies and amplitudes. Then we calculate thermoacoustic limit cycles and explain their existence and stability by examining the amplitude-dependence of the gain and phase of the FDF. We find that, when the phase speed equals the mean flow speed, the system has only one stable state. When the phase speed does not equal the mean flow speed, however, the system supports multiple limit cycles because the phase of the FDF changes significantly with oscillation amplitude. This shows that the phase speed of velocity perturbations has a strong influence on the nonlinear thermoacoustic behaviour of ducted premixed flames. © 2013 The Combustion Institute.

A novel theoretical approach to passive control of thermo-acoustic oscillations: Application to ducted heat sources

In this paper, we develop a linear technique that predicts how the stability of a thermo-acoustic system changes due to the action of a generic passive feedback device or a generic change in the base state. From this, one can calculate the passive device or base state change that most stabilizes the system. This theoretical framework, based on adjoint equations, is applied to two types of Rijke tube. The first contains an electrically-heated hot wire and the second contains a diffusion flame. Both heat sources are assumed to be compact so that the acoustic and heat release models can be decoupled. We find that the most effective passive control device is an adiabatic mesh placed at the downstream end of the Rijke tube. We also investigate the effects of a second hot wire and a local variation of the cross-sectional area but find that both affect the frequency more than the growth rate. This application of adjoint sensitivity analysis opens up new possibilities for the passive control of thermo-acoustic oscillations. For example, the influence of base state changes can be combined with other constraints, such as that the total heat release rate remains constant, in order to show how an unstable thermo-acoustic system should be changed in order to make it stable. Copyright © 2013 by ASME.

Computations of turbulent lean premixed combustion using conditional moment closure

Conditional Moment Closure (CMC) is a suitable method for predicting scalars such as carbon monoxide with slow chemical time scales in turbulent combustion. Although this method has been successfully applied to non-premixed combustion, its application to lean premixed combustion is rare. In this study the CMC method is used to compute piloted lean premixed combustion in a distributed combustion regime. The conditional scalar dissipation rate of the conditioning scalar, the progress variable, is closed using an algebraic model and turbulence is modelled using the standard k-e{open} model. The conditional mean reaction rate is closed using a first order CMC closure with the GRI-3.0 chemical mechanism to represent the chemical kinetics of methane oxidation. The PDF of the progress variable is obtained using a presumed shape with the Beta function. The computed results are compared with the experimental measurements and earlier computations using the transported PDF approach. The results show reasonable agreement with the experimental measurements and are consistent with the transported PDF computations. When the compounded effects of shear-turbulence and flame are strong, second order closures may be required for the CMC. © 2013 Copyright Taylor and Francis Group, LLC.

Laminar propagation of lean premixed flames ignited in stratified mixture

Combustion in stratified mixtures is envisaged in practical energy systems such as direct-injection spark-ignited (DISI) car engines, gas turbines, for reducing CO2 and pollutant emissions while protecting their efficiency. The mixture gradients change the fundamental properties of the flame, especially by a difference in temperature and composition between the burnt gases and those of a flame consuming a homogeneous mixture. This paper presents an investigation of the properties of the flame propagating in a lean homogeneous mixture after ignition in a richer mixture according to the magnitude of the stratification. Three magnitudes of stratification are investigated. The local flame burning velocity is determined by an original PIV algorithm developed previously. The local equivalence ratio in the fresh gases is measured from anisole PLIF. From the simultaneous PIV-PLIF measurements, the flame burning velocities conditioned on the local stretch rate and equivalence ratio in fresh gases are measured. The flame propagating through the homogeneous lean mixture has properties depending on the ignition conditions in the stratified layer. The flame propagating in the lean mixture is back-supported longer for ignition under the richer condition. The change of stretch sensitivity and burning velocity of the flame in the lean mixture is measured over time for the three magnitudes of mixture stratification investigated. The ignition in richer mixtures compensates for the nonequidiffusion effect of lean propane flame and sustains its robustness to stretch. The flame propagation in the lean homogeneous mixture is enhanced by ignition in a richer stratified layer, as much by their robustness to stretch as by an increase in the flame speed or the burning velocity. The decay time of this influence of the stratification, called memory effect, is determined. © 2013 The Combustion Institute.

Scalar gradient behaviour in MILD combustion

The results of three-dimensional Direct Numerical Simulation (DNS) of Moderate, Intense Low-oxygen Dilution (MILD) and conventional premixed turbulent combustion conducted using a skeletal mechanism including the effects of non-unity Lewis numbers and temperature dependent transport properties are analysed to investigate combustion characteristics using scalar gradient information. The DNS data is also used to synthesise laser induced fluorescence (LIF) signals of OH, CH2O, and CHO. These signals are analysed to verify if they can be used to study turbulent MILD combustion and it has been observed that at least two (OH and CH2O) LIF signals are required since the OH increase across the reaction zone is smaller in MILD combustion compared to premixed combustion. The scalar gradient PDFs conditioned on the reaction rate obtained from the DNS data and synthesised LIF signals suggests a strong gradient in the direction normal to the MILD reaction zone with moderate reaction rate implying flamelet combustion. However, the PDF of the normal gradient is as broad as for the tangential gradient when the reaction rate is high. This suggests a non-flamelet behaviour, which is due to interaction of reaction zones. The analysis of the conditional PDFs for the premixed case confirms the expected behaviour of scalar gradient in flamelet combustion. It has been shown that the LIF signals synthesised using 2D slices of DNS data also provide very similar insights. These results demonstrate that the so-called flameless combustion is not an idealised homogeneous reactive mixture but has common features of conventional combustion while containing distinctive characteristics. © 2013 The Combustion Institute.

Adaptive combustion instability control with saturation

The control of a class of combustion systems, suceptible to damage from self-excited combustion oscillations, is considered. An adaptive stable controller, called Self-Tuning Regulator (STR), has recently been developed, which meets the apparently contradictory challenge of relying as little as possible on a particular combustion model while providing some guarantee that the controller will cause no harm. The controller injects some fuel unsteadily into the burning region, thereby altering the heat release, in response to an input signal detecting the oscillation. This paper focuses on an extension of the STR design, when, due to stringent emission requirements and to the danger of flame extension, the amount of fuel used for control is limited in amplitude. A Lyapunov stability analysis is used to prove the stability of the modified STR when the saturation constraint is imposed. The practical implementation of the modified STR remains straightforward, and simulation results, based on the nonlinear premixed flame model developed by Dowling, show that in the presence of a saturation constraint, the self-excited oscillations are damped more rapidly with the modified STR than with the original STR. © 2001 by S. Evesque. Published by the American Institute of Aeronautics and Astronautics, Inc.

Comparison of acoustic velocity perturbation measurements using PIV vs. two-microphone technique

Understanding combustion instabilities requires accurate measurements of the acoustic velocity perturbation into injectors. This is often accomplished via the use of the two microphone technique, as this only requires two pressure transducers. However, measurements of the actual velocities emerging from the injectors are not often taken, leaving questions regarding the assumptions about the acoustic velocity. A comparison of velocity measured at downstream of the injector with that of two-microphone technique can show the accuracy and limitations of two-microphone technique. In this paper, velocity measurements are taken using both particle image velocimetry (PIV) and the two-microphone technique in a high pressure facility designed for aeroengine injector measurements. The flow is excited using an area modulation device installed on the choked end of the combustion chamber, with PIV measurements enabled by optical access downstream of the injector through a quartz tube and windows. Acoustic velocity perturbations at the injector are determined by considering the Fourier transformed pressure fluctuations for two microphones installed at a known distance upstream of the injector. PIV measurements are realized by seeding the air flow with micrometric water particles under 2.5 bar pressure at ambient temperature. Phase locked velocity fields are realized by synchronizing the acquisition of PIV images with the revolution of the acoustic modulator using the pressure signal measured at the face of injector. The mean velocity fluctuation is calculated as the difference between maximum and minimum velocities, normalized by the mean velocity of the unforced case. Those values are compared with the peak-to-peak velocity fluctuation amplitude calculated by the two-microphone technique. Although the ranges of velocity fluctuations for both techniques are similar, the variation of fluctuation with forcing frequencies diverges significantly with frequency. The differences can be attributed to several limitations associated with of both techniques, such as the quality of the signal, the signal/noise ratio, the accuracy of PIV measurements and the assumption of isentropic flow of the particle velocity from the plenum through the injector. We conclude that two-microphone methods can be used as a reference value for the velocity fluctuation in low order applications such as flame transfer functions, but not for drawing conclusions regarding the absolute velocity fluctuations in the injector. Copyright © 2013 by ASME.

Structure of cryogenic flames at elevated pressures

This paper presents new experimental results on cryogenic jet flames formed by a coaxial injector at a pressure of 70 bar, which approaches the pressures found in rocket engines. This element, fed with liquid oxygen and gaseous hydrogen, is placed in a square combustion chamber equipped with quartz windows. The flame is examined via spectroscopy, OH* emission, and backlighting, the aim being to provide basic information on the flame structure. It is found that some of the OH* emission is absorbed by the OH radicals present in the flame. A detailed examination of this effect is presented, in which it is shown that, for this turbulent flame, the Abel transform gives the position of the intense reaction region, whether or not absorption is signficant. The flame is attached to the oxygen injector, as at low pressure. At high pressure, flame expansion is reduced compared with low pressure and is also less dependent on the momentum flux ratio between the hydrogen and the oxygen streams. An analysis of the relevant Damköhler numbers suggests that this is because the rate of combustion is mainly controlled by large-scale turbulent mixing at high pressure, and it is dominated by jet break-up, atomization, and vaporization at low pressures. Jet break-up is particularly dependent on the momentum flux ratio. Finally, the mean volumetric heat release rates and flame surface density in the experimental facility are estimated.

Hermetically-coated superparamagnetic Fe2O3 particles with SiO2 nanofilms

Magnetic nanoparticles are frequently coated with SiO2to improve their functionality and bio-compatibility in a range of biomedical and polymer nanocomposile applications. In this paper, a scalable flame aerosol technology is used to produce highly dispersible, superparamagnetic iron oxide nanoparticles hermetically coaled with silica to retain full magnetization performance. Iron oxide particles were produced by flame spray pyrolysis (FSP) of iron acelylacetonale in xylene/acetonitrile solutions, and the resulting aerosol was in situ coaled with SiO2 by oxidation of swirling hexamethlydisiloxane vapor. The process allows independent control of the core Fe2O3, particle properties and the thickness of their silica coaling film. This ensures that the non-magnetic SiO2 layer can be closely controlled and minimized. The optimal SiO2 content for complete (hermetic) encapsulation of the magnetic core particles was determined by isopropanol chemisorption. The magnetization of Fe2O3 coated with about 2 nm thin SiO2 layers was nearly identical lo that of uncoated, pure Fe2O3 nanoparlicles.

Hermetically coated superparamagnetic Fe2O3 particles with SiO2 nanofilms

Magnetic nanoparticles are frequently coated with SiO2 to improve their functionality and biocom-patibility in a range of biomedical and polymer nanocomposite applications. In this paper, a scalable flame aerosol technology is used to produce highly dispersible, superparamagnetic iron oxide nanoparticles hermetically coated with silica to retain full magnetization performance. Iron oxide particles were produced by flame spray pyrolysis of iron acetylacetonate in xylene/acetonitrile solutions and the resulting aerosol was in situ coated with silicon dioxide by oxidation of swirling hexamethlydisiloxane vapor. The process allows independent control of the core Fe2O3 (maghemite) particle properties and the thickness of their silica coating film. This ensures that the nonmagnetic SiO2 layer can be closely controlled and minimized. The optimal SiO2 content for complete (hermetic) encapsulation of the magnetic core particles was determined by isopropanol chemisorption. The magnetization of Fe 2O3 coated with about 2 nm thin SiO2 layers was nearly identical to that of uncoated, pure Fe2O3 nanoparticles. © 2009 American Chemical Society.

Effect of size, composition, and morphology on magnetic performance: First-order reversal curves evaluation of iron oxide nanoparticles

© 2014 AIP Publishing LLC. Superparamagnetic nanoparticles are employed in a broad range of applications that demand detailed magnetic characterization for superior performance, e.g., in drug delivery or cancer treatment. Magnetic hysteresis measurements provide information on saturation magnetization and coercive force for bulk material but can be equivocal for particles having a broad size distribution. Here, first-order reversal curves (FORCs) are used to evaluate the effective magnetic particle size and interaction between equally sized magnetic iron oxide (Fe2O3) nanoparticles with three different morphologies: (i) pure Fe2O3, (ii) Janus-like, and (iii) core/shell Fe2O3/SiO2synthesized using flame technology. By characterizing the distribution in coercive force and interaction field from the FORC diagrams, we find that the presence of SiO2in the core/shell structures significantly reduces the average coercive force in comparison to the Janus-like Fe2O3/SiO2and pure Fe2O3particles. This is attributed to the reduction in the dipolar interaction between particles, which in turn reduces the effective magnetic particle size. Hence, FORC analysis allows for a finer distinction between equally sized Fe2O3particles with similar magnetic hysteresis curves that can significantly influence the final nanoparticle performance.

Discrete multicomponent model for biodiesel spray combustion simulation

An important first step in spray combustion simulation is an accurate determination of the fuel properties which affects the modelling of spray formation and reaction. In a practical combustion simulation, the implementation of a multicomponent model is important in capturing the relative volatility of different fuel components. A Discrete Multicomponent (DM) model is deemed to be an appropriate candidate to model a composite fuel like biodiesel which consists of four components of fatty acid methyl esters (FAME). In this paper, the DM model is compared with the traditional Continuous Thermodynamics (CTM) model for both diesel and biodiesel. The CTM model is formulated based on mixing rules that incorporate the physical and thermophysical properties of pure components into a single continuous surrogate for the composite fuel. The models are implemented within the open-source CFD code OpenFOAM, and a semi-quantitative comparison is made between the predicted spray-combustion characteristics and optical measurements of a swirl-stabilised flame of diesel and biodiesel. The DM model performs better than the CTM model in predicting a higher magnitude of heat release rate in the top flame brush region of the biodiesel flame compared to that of the diesel flame. Using both the DM and CTM models, the simulation successfully reproduces the droplet size, volume flux, and droplet density profiles of diesel and biodiesel. The DM model predicts a longer spray penetration length for biodiesel compared to that of diesel, as seen in the experimental data. Also, the DM model reproduces a segregated biodiesel fuel vapour field and spray in which the most abundant FAME component has the longest vapour penetration. In the biodiesel flame, the relative abundance of each fuel component is found to dominate over the relative volatility in terms of the vapour species distribution and vice versa in the liquid species distribution. © 2014 Elsevier Ltd. All rights reserved.

Local strain rate and curvature dependences of scalar dissipation rate transport in turbulent premixed flames: A direct numerical simulation analysis

The statistical behaviours of the instantaneous scalar dissipation rate Nc of reaction progress variable c in turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds number Ret. The statistical behaviours of N c and different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation term T1 is positive, whereas the molecular dissipation term (-D2) acts as a leading order sink. The mean contribution of the strain rate term T2 is predominantly negative for the cases considered here. The mean reaction rate contribution T3 is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term (D) assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours of Nc, T1, T2, T 3, (-D2), and f(D) have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rate aT and curvature km. Detailed physical explanations have been provided for the observed behaviour. © 2014 Y. Gao et al.

Simulations of autoignition and laminar premixed flames in methane/air mixtures diluted with hot products

This article considers constant-pressure autoignition and freely propagating premixed flames of cold methane/air mixtures mixed with equilibrium hot products at high enough dilution levels to burn within the moderate to intense low oxygen dilution (MILD) combustion regime. The analysis is meant to provide further insight on MILD regime boundaries and to identify the effect of hot products speciation. As the mass fraction of hot products in the reactants mixture increases, autoignition occurs earlier. Species profiles show that the products/reactants mixture approximately equilibrates to a new state over a quick transient well before the main autoignition event, but as dilution becomes very high, this equilibration transient becomes more prominent and eventually merges with the primary ignition event. The dilution level at which these two reactive zones merge corresponds well with that marking the transition into the MILD regime, as defined according to conventional criteria. Similarly, premixed flame simulations at high dilutions show evidence of significant reactions involving intermediate species prior to the flame front. Since the premixed flame governing equations system demands that the species and temperature gradients be zero at the "cold" boundary, flame speed cannot be calculated above a certain dilution level. Up to this point, which again agrees reasonably well with the transition into the MILD regime according to convention, the laminar burning velocity was found to increase with hot product dilution while flame thickness remained largely unchanged. Some comments on the MILD combustion regime boundary definition for gas turbine applications are included. Copyright © Taylor & Francis Group, LLC.

Influence of turbulence-chemistry interaction for n-heptane spray combustion under diesel engine conditions with emphasis on soot formation and oxidation

The influence of the turbulence-chemistry interaction (TCI) for n-heptane sprays under diesel engine conditions has been investigated by means of computational fluid dynamics (CFD) simulations. The conditional moment closure approach, which has been previously validated thoroughly for such flows, and the homogeneous reactor (i.e. no turbulent combustion model) approach have been compared, in view of the recent resurgence of the latter approaches for diesel engine CFD. Experimental data available from a constant-volume combustion chamber have been used for model validation purposes for a broad range of conditions including variations in ambient oxygen (8-21% by vol.), ambient temperature (900 and 1000 K) and ambient density (14.8 and 30 kg/m3). The results from both numerical approaches have been compared to the experimental values of ignition delay (ID), flame lift-off length (LOL), and soot volume fraction distributions. TCI was found to have a weak influence on ignition delay for the conditions simulated, attributed to the low values of the scalar dissipation relative to the critical value above which auto-ignition does not occur. In contrast, the flame LOL was considerably affected, in particular at low oxygen concentrations. Quasi-steady soot formation was similar; however, pronounced differences in soot oxidation behaviour are reported. The differences were further emphasised for a case with short injection duration: in such conditions, TCI was found to play a major role concerning the soot oxidation behaviour because of the importance of soot-oxidiser structure in mixture fraction space. Neglecting TCI leads to a strong over-estimation of soot oxidation after the end of injection. The results suggest that for some engines, and for some phenomena, the neglect of turbulent fluctuations may lead to predictions of acceptable engineering accuracy, but that a proper turbulent combustion model is needed for more reliable results. © 2014 Taylor & Francis.