Effects of the biodiesel fuel physical properties on the swirl stabilised spray combustion characteristics

An increasin g interest in biofuel applications in modern engines requires a better understanding of biodiesel combustion behaviour. Many numerical studies have been carried out on unsteady combustion of biodiesel in situations similar to diesel engines, but very few studies have been done on the steady combustion of biodiesel in situations similar to a gas turbine combustor environment. The study of biodiesel spray combustion in gas turbine applications is of special interest due to the possible use of biodiesel in the power generation and aviation industries. In modelling spray combustion, an accurate representation of the physical properties of the fuel is a first important step, since spray formation is largely influenced by fuel properties such as viscosity, density, surface tension and vapour pressure. In the present work, a calculated biodiesel properties database based on the measured composition of Fatty Acid Methyl Esters (FAME) has been implemented in a multi-dimensional Computational Fluid Dynamics (CFD) spray simulation code. Simulations of non-reacting and reacting atmospheric-pressure sprays of both diesel and biodiesel have been carried out using a spray burner configuration for which experimental data is available. A pre-defined droplet size probability density function (pdf) has been implemented together with droplet dynamics based on phase Doppler anemometry (PDA) measurements in the near-nozzle region. The gas phase boundary condition for the reacting spray cases is similar to that of the experiment which employs a plain air-blast atomiser and a straight-vane axial swirler for flame stabilisation. A reaction mechanism for heptane has been used to represent the chemistry for both diesel and biodiesel. Simulated flame heights, spray characteristics and gas phase velocities have been found to compare well with the experimental results. In the reacting spray cases, biodiesel shows a smaller mean droplet size compared to that of diesel at a constant fuel mass flow rate. A lack of sensitivity towards different fuel properties has been observed based on the non-reacting spray simulations, which indicates a need for improved models of secondary breakup. By comparing the results of the non-reacting and reacting spray simulations, an improvement in the complexity of the physical modelling is achieved which is necessary in the understanding of the complex physical processes involved in spray combustion simulation. Copyright © 2012 SAE International.

Modal dynamics of self-excited azimuthal instabilities in an annular combustion chamber

In this paper we describe the time-varying amplitude and its relation to the global heat release rate of self-excited azimuthal instabilities in a simple annular combustor operating under atmospheric conditions. The combustor was modular in construction consisting of either 12, 15 or 18 equally spaced premixed bluff-body flames around a fixed circumference, enabling the effect of large-scale interactions between adjacent flames to be investigated. High-speed OH* chemiluminescence imaged from above the annulus and pressure measurements obtained at multiple locations around the annulus revealed that the limit cycles of the modes are degenerate in so much as they undergo continuous transitions between standing and spinning modes in both clockwise (CW) and anti-clockwise (ACW) directions but with the same resonant frequency. Similar behaviour has been observed in LES simulations which suggests that degenerate modes may be a characteristic feature of self-excited azimuthal instabilities in annular combustion chambers. By modelling the instabilities as two acoustic waves of time-varying amplitude travelling in opposite directions we demonstrate that there is a statistical prevalence for either standing m=1 or spinning m=±1 modes depending on flame spacing, equivalence ratio, and swirl configuration. Phase-averaged OH* chemiluminescence revealed a possible mechanism that drives the direction of the spinning modes under limit-cycle conditions for configurations with uniform swirl. By dividing the annulus into inner and outer annular regions it was found that the spin direction coincided with changes in the spatial distribution of the peak heat release rate relative to the direction of the bulk swirl induced along the annular walls. For standing wave modes it is shown that the globally integrated fluctuations in heat release rate vary in magnitude along the acoustic mode shape with negligible contributions at the pressure nodes and maximum contributions at the pressure anti-nodes. © 2013.

Oxygen carrier dispersion in inert packed beds to improve performance in chemical looping combustion

Various packed beds of copper-based oxygen carriers (CuO on Al2O3) were tested over 100 cycles of low temperature (673K) Chemical Looping Combustion (CLC) with H2 as the fuel gas. The oxygen carriers were uniformly mixed with alumina (Al2O3) in order to investigate the level of separation necessary to prevent agglomeration. It was found that a mass ratio of 1:6 oxygen carrier to alumina gave the best performance in terms of stable, repeating hydrogen breakthrough curves over 100 cycles. In order to quantify the average separation achieved in the mixed packed beds, two sphere-packing models were developed. The hexagonal close-packing model assumed a uniform spherical packing structure, and based the separation calculations on a hypergeometric probability distribution. The more computationally intensive full-scale model used discrete element modelling to simulate random packing arrangements governed by gravity and contact dynamics. Both models predicted that average 'nearest neighbour' particle separation drops to near zero for oxygen carrier mass fractions of x≥0.25. For the packed bed systems studied, agglomeration was observed when the mass fraction of oxygen carrier was above this threshold. © 2013 Elsevier B.V.

Model-based control of thermoacoustic instabilities in partially premixed lean combustion - A design case study

Self-excited oscillation is becoming a major issue in low-emission, lean partially premixed combustion systems, and active control has been shown to be a feasible method to suppress such instabilities. A number of robust control methods are employed to obtain a feedback controller and it is observed that the robustness to system uncertainty is significantly better for a low complexity controller in spite of the norms being similar. Moreover, we demonstrate that closed-loop stability for such a complex system can be proved via use of the integral quadratic constraint method. Open- and closed-loop nonlinear simulations are provided. © 2013 Copyright Taylor and Francis Group, LLC.

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.

Computational engineering design for micro-scale combustion devices: A thermally improved configuration

A multi-objective design optimisation study has been carried out with the objectives to improve the overall efficiency of the device and to reduce the fuel consumption for the proposed micro-scale combustor design configuration. In a previous study we identified the topology of the combustion chamber that produced improved behaviour of the device in terms of the above design criteria. We now extend our design approach, and we propose a new configuration by the addition of a micro-cooling channel that will improve the thermal behaviour of the design as previously suggested in literature. Our initial numerical results revealed an improvement of 2.6% in the combustion efficiency when we applied the micro-cooling channel to an optimum design configuration we identified from our earlier multi-objective optimisation study, and under the same operating conditions. The computational modelling of the combustion process is implemented in the commercial computational fluid dynamics package ANSYS-CFX using Finite Rate Chemistry and a single step hydrogen-air reaction. With this model we try to balance good accuracy of the combustion solution and at the same time practicality within the context of an optimisation process. The whole design system comprises also the ANSYS-ICEM CFD package for the automatic geometry and mesh generation and the Multi-Objective Tabu Search algorithm for the design space exploration. We model the design problem with 5 geometrical parameters and 3 operational parameters subject to 5 design constraints that secure practicality and feasibility of the new optimum design configurations. The final results demonstrate the reliability and efficiency of the developed computational design system and most importantly we assess the practicality and manufacturability of the revealed optimum design configurations of micro-combustor devices. Copyright © 2013 by ASME.

Prediction of combustion noise for an aeroengine combustor

Combustion noise may become an important noise source for lean-burn gas turbine engines, and this noise is usually associated with highly unsteady flames. This work aims to compute the broadband combustion noise spectrum for a realistic aeroengine combustor and to compare with available measured noise data on a demonstrator aeroengine. A low-order linear network model is applied to a demonstrator engine combustor to obtain the transfer function that relates to unsteadiness in the rate of heat release, acoustic, entropic, and vortical fluctuations. A spectral model is used for the heat release rate fluctuation, which is the source of the noise. The mean flow of the aeroengine combustor required as input data to this spectral model is obtained from Reynolds-averaged Navier-Stokes simulations. The computed acoustic field for a low-to-medium power setting indicates that the models used in this study capture the main characteristics of the broadband spectral shape of combustion noise. Reasonable agreement with the measured spectral level is achieved. © 2012 AIAA.

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.

A flame surface density approach to large-eddy simulation of premixed turbulent combustion

The flame surface density approach to the modeling of premixed turbulent combustion is well established in the context of Reynolds-averaged simulations. For the future, it is necessary to consider large-eddy simulation (LES), which is likely to offer major advantages in terms of physical accuracy, particularly for unsteady combustion problems. LES relies on spatial filtering for the removal of unresolved phenomena whose characteristic length scales are smaller than the computational grid scale. Thus, there is a need for soundly based physical modeling at the subgrid scales. The aim of this paper is to explore the usefulness of the flame surface density concept as a basis for LES modeling of premixed turbulent combustion. A transport equation for the filtered flame surface density is presented, and models are proposed for unclosed terms. Comparison with Reynolds-averaged modeling is shown to reveal some interesting similarities and differences. These were exploited together with known physics and statistical results from experiment and from direct numerical stimulation in order to gain insight and refine the modeling. The model has been implemented in a combustion LES code together with standard models for scalar and momentum transport. Computational results were obtained for a simple three-dimensional flame propagation test problem, and the relative importance of contributing terms in the modeled equation for flame surface density was assessed. Straining and curvature are shown to have a major influence at both the resolved and subgrid levels.

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.

CO2-gasification of a lignite coal in the presence of an iron-based oxygen carrier for chemical-looping combustion

Chemical-looping combustion (CLC) has the inherent property of separating the product CO2 from flue gases. Instead of air, it uses an oxygen carrier, usually in the form of a metal oxide, to provide oxygen for combustion. All techniques so far proposed for chemical looping with solid fuels involve initially the gasification of the solid fuel in order for the gaseous products to react with the oxygen carrier. Here, the rates of gasification of coal were compared when gasification was undertaken in a fluidised bed of either (i) an active Fe-based oxygen carrier used for chemical looping or (ii) inert sand. This enabled an examination of the ability of chemical looping materials to enhance the rate of gasification of solid fuels. Batch gasification and chemical-looping combustion experiments with a German lignite and its char are reported, using an electrically-heated fluidised bed reactor at temperatures from 1073 to 1223 K. The fluidising gas was CO2 in nitrogen. The kinetics of the gasification were found to be significantly faster in the presence of the oxygen carrier, especially at temperatures above 1123 K. A numerical model was developed to account for external and internal mass transfer and for the effect of the looping agent. The model also included the effects of the evolution of the pore structure at different conversions. The presence of Fe2O3 led to an increase in the rate of gasification because of the rapid oxidation of CO by the oxygen carrier to CO2. This resulted in the removal of CO and maintained a higher mole fraction of CO2 in the mixture of gas around the particle of char, i.e. within the mass transfer boundary layer surrounding the particle. This effect was most prominent at about 20% conversion when (i) the surface area for reaction was at its maximum and (ii) because of the accompanying increase in porosity and pore size, intraparticle resistance to gas mass transfer within the particle of char had fallen, compared with that in the initial particle. Excellent agreement was observed between the rates predicted by the numerical model and those observed experimentally. ©2013 Elsevier Ltd. All rights reserved.