Diffusion and chemical reaction in the porous structures of solid oxide fuel cells

The Rolls-Royce Integrated-Planar Solid Oxide Fuel Cell (IP-SOFC) consists of ceramic modules which have electrochemical cells printed on the outer surfaces. The cathodes are the outermost layer of each cell and are supplied with oxygen from air flowing over the outside of the module. The anodes are in direct contact with the ceramic structure and are supplied with fuel from internal gas channels. Natural gas is reformed into hydrogen for use by the fuel cells in a separate reformer module of similar design except that the fuel cells are replaced by a reforming catalyst layer. The performance of the modules is intrinsically linked to the behaviour of the gas flows within their porous structures. Because the porous layers are very thin, a one-dimensional flow model provides a good representation of the flow property variations between fuel channel and fuel cell or reforming catalyst. The multi-component convective-diffusive flows are simulated using a new theory of flow in porous material, the Cylindrical Pore Interpolation Model. The effects of the catalysed methane reforming and water-gas shift chemical reactions are also considered using appropriate kinetic models. It is found that the shift reaction, which is catalysed by the anode material, has certain beneficial effects on the fuel cell module performance. In the reformer module it was found that the flow resistance of the porous support structure makes it difficult to sustain a high methane conversion rate. Although the analysis is based on IP-SOFC geometry, the modelling approach and general conclusions are applicable to other types of SOFC.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Effects of Lewis number on scalar variance transport in premixed flames

The influences of differential diffusion rates of heat and mass on the transport of the variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been studied using three-dimensional simplified chemistry based Direct Numerical Simulation (DNS) data of statistically planar turbulent premixed flames with global Lewis number ranging from Le = 0.34 to 1.2. The Lewis number effects on the statistical behaviours of the various terms of the transport equations of variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been analysed in the context of Reynolds Averaged Navier Stokes (RANS) simulations. It has been found that the turbulent fluxes of the progress variable and temperature variances exhibit counter-gradient transport for the flames with Lewis number significantly smaller than unity whereas the extent of this counter-gradient transport is found to decrease with increasing Lewis number. The Lewis number is also shown to have significant influences on the magnitudes of the chemical reaction and scalar dissipation rate contributions to the scalar variance transport. The modelling of the unclosed terms in the scalar variance equations for the non-unity Lewis number flames have been discussed in detail. The performances of the existing models for the unclosed terms are assessed based on a-priori analysis of DNS data. Based on the present analysis, new models for the unclosed terms of the active scalar variance transport equations are proposed, whenever necessary, which are shown to satisfactorily capture the behaviours of unclosed terms for all the flames considered in this study. © 2010 Springer Science+Business Media B.V.

Turbulent Premixed Flames

An improved understanding of lean fuel turbulent premixed flames must play a central role in the fundamental science of these new concepts.

The nonlinear heat release response of stratified lean-premixed flames to acoustic velocity oscillations

This paper analyzes the forced response of swirl-stabilized lean-premixed flames to high-amplitude acoustic forcing in a laboratory-scale stratified burner operated with CH4 and air at atmospheric pressure. The double-swirler, double-channel annular burner was specially designed to generate high-amplitude acoustic velocity oscillations and a radial equivalence ratio gradient at the inlet of the combustion chamber. Temporal oscillations of equivalence ratio along the axial direction are dissipated over a long distance, and therefore the effects of time-varying fuel/air ratio on the response are not considered in the present investigation. Simultaneous measurements of inlet velocity and heat release rate oscillations were made using a constant temperature anemometer and photomultiplier tubes with narrow-band OH*/CH* interference filters. Time-averaged and phase-synchronized CH* chemiluminescence intensities were measured using an intensified CCD camera. The measurements show that flame stabilization mechanisms vary depending on equivalence ratio gradients for a constant global equivalence ratio (φg=0.60). Under uniformly premixed conditions, an enveloped M-shaped flame is observed. In contrast, under stratified conditions, a dihedral V-flame and a toroidal detached flame develop in the outer stream and inner stream fuel enrichment cases, respectively. The modification of the stabilization mechanism has a significant impact on the nonlinear response of stratified flames to high-amplitude acoustic forcing (u'/U∼0.45 and f=60, 160Hz). Outer stream enrichment tends to improve the flame's stiffness with respect to incident acoustic/vortical disturbances, whereas inner stream stratification tends to enhance the nonlinear flame dynamics, as manifested by the complex interaction between the swirl flame and large-scale coherent vortices with different length scales and shedding points. It was found that the behavior of the measured flame describing functions (FDF), which depend on radial fuel stratification, are well correlated with previous measurements of the intensity of self-excited combustion instabilities in the stratified swirl burner. The results presented in this paper provide insight into the impact of nonuniform reactant stoichiometry on combustion instabilities, its effect on flame location and the interaction with unsteady flow structures. © 2011 The Combustion Institute.

Effect of flow-geometry on turbulence-scalar interaction in premixed flames

Turbulent combustion of stoichiometric hydrogen-air mixture is simulated using direct numerical simulation methodology, employing complex chemical kinetics. Two flame configurations, freely propagating and V-flames stabilized behind a hot rod, are simulated. The results are analyzed to study the influence of flame configuration on the turbulence-scalar interaction, which is critical for the scalar gradient generation processes. The result suggests that this interaction process is not influenced by the flame configuration and the flame normal is found to align with the most extensive strain in the region of intense heat release. The combustion in the rod stabilized flame is found to be flamelet like in an average sense and the growth of flame-brush thickness with the downstream distance is represented well by Taylor theory of turbulent diffusion, when the flame-brushes are non-interacting. The thickness is observed to saturate when the flame-brushes interact, which is found to occur in the simulated rod stabilized flame with Taylor micro-scale Reynolds number of 97. © 2011 American Institute of Physics.

Effects of preferential transport in turbulent bluff-body-stabilized lean premixed CH 4/air flames

Preferential species diffusion is known to have important effects on local flame structure in turbulent premixed flames, and differential diffusion of heat and mass can have significant effects on both local flame structure and global flame parameters, such as turbulent flame speed. However, models for turbulent premixed combustion normally assume that atomic mass fractions are conserved from reactants to fully burnt products. Experiments reported here indicate that this basic assumption may be incorrect for an important class of turbulent flames. Measurements of major species and temperature in the near field of turbulent, bluff-body stabilized, lean premixed methane-air flames (Le=0.98) reveal significant departures from expected conditional mean compositional structure in the combustion products as well as within the flame. Net increases exceeding 10% in the equivalence ratio and the carbon-to-hydrogen atom ratio are observed across the turbulent flame brush. Corresponding measurements across an unstrained laminar flame at similar equivalence ratio are in close agreement with calculations performed using Chemkin with the GRI 3.0 mechanism and multi-component transport, confirming accuracy of experimental techniques. Results suggest that the large effects observed in the turbulent bluff-body burner are cause by preferential transport of H 2 and H 2O through the preheat zone ahead of CO 2 and CO, followed by convective transport downstream and away from the local flame brush. This preferential transport effect increases with increasing velocity of reactants past the bluff body and is apparently amplified by the presence of a strong recirculation zone where excess CO 2 is accumulated. © 2011 The Combustion Institute.

Low-order modelling of ducted flames with temporally varying equivalence ratio in realistic geometries

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the socalled G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies. Copyright © 2011 by ASME.

The behaviour of laminar stratified methane/air flames in counterflow

Flames propagating through a mixture with a gradient of equivalence ratio have been previously demonstrated to travel faster or slower than their equivalent premixed flames. The present study aims to numerically investigate the response of strained laminar methane-air flames to such gradients. The flames are simulated in a counterflow configuration where a premixed reactant stream at equivalence ratio φR opposes a hot equilibrium stream at equivalence ratio φP. Premixed and stratified flames are compared with respect to the equivalence ratio φ* and the corresponding gradient ∇φ* at the point of peak heat release rate, for three strain rates, a=50, 300 and 500s-1 and a range of φ*. The effect of different stratification levels is also investigated by varying the ratio of φP to φR, Θ. Results indicate that, as long as flames stabilize within the diffusion layer and Θ>1, increased heat release rate Q is seen throughout the progress variable space in comparison to the premixed state. In contrast, an attenuation of heat release rate is seen for Θ<1. The enhancement (or attenuation) of heat release varies monotonically with Θ. The effect of stratification on flame behavior becomes more pronounced as the strain rate increases. The present study reveals the mechanisms for the propagation of quasi-steady stratified flames under lean and rich conditions: stratified flames are primarily dominated by the diffusion of heat under lean conditions, and diffusion of H2 under rich conditions. Thanks to species and thermal support, stratified flames continue to burn beyond the premixed lean and rich flammability limits. Further investigation on the unsteady response of flames to the fluctuating equivalence ratio implies that the steady results represent the unsteady response well, as long as φ* and ∇φ* are similar in both steady and unsteady cases. © 2013 The Combustion Institute.

Experimental investigation of dynamic response of acoustically forced turbulent premixed CH4/CO2/air flames

Increasing demand for energy and continuing increase in environmental as well as financial cost of use of fossil fuels drive the need for utilization of fuels from sustainable sources for power generation. Development of fuel-flexible combustion systems is vital in enabling the use of sustainable fuels. It is also important that these sustainable combustion systems meet the strict governmental emission legislations. Biogas is considered as one of the viable sustainable fuels that can be used to power modern gas turbines: However, the change in chemical, thermal and transport properties as well as change in Wobbe index due to the variation of the fuel constituents can have a significant effect on the performance of the combustor. It is known that the fuel properties have strong influence on the dynamic flame response; however there is a lack of detailed information regarding the effect of fuel compositions on the sensitivity of the flames subjected to flow perturbations. In this study, we describe an experimental effort investigating the response of premixed biogas-air turbulent flames with varying proportions of CH4 and CO2 to velocity perturbations. The flame was stabilized using a centrally placed conical bluff body. Acoustic perturbations were imposed to the flow using loud speakers. The flame dynamics and the local heat release rate of these acoustically excited biogas flames were studied using simultaneous measurements of OH and H2CO planar laser induced fluorescence. OH* chemiluminescence along with acoustic pressure measurements were also recorded to estimate the total flame heat release modulation and the velocity fluctuations. The measurements were carried out by keeping the theoretical laminar flame speed constant while varying the bulk velocity and the fuel composition. The results indicate that the flame sensitivity to perturbations increased with increased dilution of CH4 by CO2 at low amplitude forcing, while at high amplitude forcing conditions the magnitude of the flame response was independent of dilution.