Statistics of reaction progress variable and mixture fraction gradients from DNS of turbulent partially premixed flames

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios = 0.7 and 1.0 using three-dimensional, simplified chemistry-based direct numerical simulations (DNS). The simulation parameters are chosen such that the flames represent the thin reaction zones regime combustion. A random bimodal distribution of equivalence ratio is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity. The results suggest that the probability density functions (PDFs) of the mixture fraction gradient magnitude |Δξ| (i.e., P(|Δξ|)) can be reasonably approximated using a log-normal distribution. However, this presumed PDF distribution captures only the qualitative nature of the PDF of the reaction progress variable gradient magnitude |Δc| (i.e., P(|Δc|)). It has been found that a bivariate log-normal distribution does not sufficiently capture the quantitative behavior of the joint PDF of |Δξ| and |Δc| (i.e., P(|Δξ|, |Δc|)), and the agreement with the DNS data has been found to be poor in certain regions of the flame brush, particularly toward the burned gas side of the flame brush. Moreover, the variables |Δξ| and |Δc| show appreciable correlation toward the burned gas side of the flame brush. These findings are corroborated further using a DNS data of a lifted jet flame to study the flame geometry dependence of these statistics. © 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.

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.

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.

Effect of drying methods on the structure and catalytic combustion activity of Mn-substituted hexaaluminate catalysts

Conventional oven drying (COD) and supercritical drying (SCD) methods were applied to the preparation of Mn-substituted hexaaluminate (BaMnA(11)O(19-alpha)) catalysts. The effect of drying methods on phase composition, specific surface area, pore structure and combustion activity of the samples was investigated. The samples obtained by SCD have higher surface area, narrower pore size distribution, and higher combustion activity than those obtained by COD.

Polypropylene as a carbon source for the synthesis of multi-walled carbon nanotubes via catalytic combustion

Multi-walled carbon nanotubes (MWCNTs) were efficiently synthesized by catalytic combustion of polypropylene (PP) using nickel compounds (such as Ni2O3, NiO, Ni(OH)(2) and NiCO3 (.) 2Ni(OH)(2)) as catalysts in the presence of organic-modified montmorillonite (OMMT) at 630-830 degrees C. Morphologies of the sample undergoing different combustion times were observed to investigate actual process producing MWCNTs by this method. The obtained MWCNTs were characterized by X-ray diffraction (XRD), transmission electron microscope and Raman spectroscopy. The yield of MWCNTs was affected by the composition of PP mixtures with OMMT and nickel compounds and the combustion temperature. The proton acidic sites from the degraded OMMT layers due to the Hoffman reaction of the modifiers at high temperature played an important role in the catalytic degradation of PP to supply carbon sources that are easy to be catalyzed by nickel catalyst for the growth of MWCNTs. The XRD measurements demonstrated that the nickel compounds were in situ reduced into the Ni(0) state with the aid of hydrogen gas and/or hydrocarbons in the degradation products of PP, and the Ni(O) was really the active site for the growth of MWCNTs. The combination of nickel compounds with OMMT was a key factor to efficiently synthesize MWCNTs via catalytic combustion of PP.

Nitrate-citrate combustion synthesis and properties of Ce1-xSmxO2-x/2 solid solutions

A nitrate-citrate combustion route to synthesize nanocrystalline samarium-doped ceria powders for solid electrolyte ceramics is presented. This route is based on the gelling of nitrate solutions by the addition of citric acid and ammonium hydroxide, followed by an intense combustion process due to an exothermic redox reaction between nitrate and citrate ions. The influence of ignition temperature on the characteristics of the powders was studied. The change of the crystal structure with the content of doped Sm was investigated. High temperature X-ray, and Raman scattering were used to characterize the sample. The lattice constant and unit volume increase with doping level and increasing temperature. Dense ceramic samples prepared by uniaxial pressing and sintering in air were also studied.

Current status and perspectives of accelerated carbonation processes on municipal waste combustion residues

The increasing volumes of municipal solid waste produced worldwide are encouraging the development of processes to reduce the environmental impact of this waste stream. Combustion technology can facilitate volume reduction of up to 90%, with the inorganic contaminants being captured in furnace bottom ash, and fly ash/APC residues. The disposal or reuse of these residues is however governed by the potential release of constituent contaminants into the environment. Accelerated carbonation has been shown to have a potential for improving the chemical stability and leaching behaviour of both bottom ash and fly ash/APC residues. However, the efficacy of carbonation depends on whether the method of gas application is direct or indirect. Also important are the mineralogy, chemistry and physical properties of the fresh ash, the carbonation reaction conditions such as temperature, contact time, CO2 partial pressure and relative humidity. This paper reviews the main issues pertaining to the application of accelerated carbonation to municipal waste combustion residues to elucidate the potential benefits on the stabilization of such residues and for reducing CO2 emissions. In particular, the modification of ash properties that occur upon carbonation and the CO2 sequestration potential possible under different conditions are discussed. Although accelerated carbonation is a developing technology, it could be introduced in new incinerator facilities as a "finishing step" for both ash treatment and reduction of CO2 emissions.

Density Functional Study on Reaction Activity of Steam Decomposition Products Adsorbed on Carbonaceous Surface

A systematic theoretical study on the adsorption of steam and its thermal decomposition products on carbon both zigzag and armchair surface was performed to provide molecular-level understanding of the reaction activity of all these reactants in biomass steam gasification process. All the calculations were carried out using density functional theory (DFT) at the B3LYP/6-31+g(d,p) level. The structures of carbonaceous surfaces, all reactants and surface complexes were optimized and characterized. Based on the value of adsorption heat been obtained from the calculation, the activity of all reactants can be ordered as: O > O2 >H2 >H >OH >H2O for both zigzag and armchair surface, and the adsorption style is physisorption to water molecule and chemisorption to the other dissociated components.

Modèles de flammelette en combustion turbulente avec extinction et réallumage : étude asymptotique et numérique, estimation d’erreur a posteriori et modélisation adaptative

On s’intéresse ici aux erreurs de modélisation liées à l’usage de modèles de flammelette sous-maille en combustion turbulente non prémélangée. Le but de cette thèse est de développer une stratégie d’estimation d’erreur a posteriori pour déterminer le meilleur modèle parmi une hiérarchie, à un coût numérique similaire à l’utilisation de ces mêmes modèles. Dans un premier temps, une stratégie faisant appel à un estimateur basé sur les résidus pondérés est développée et testée sur un système d’équations d’advection-diffusion-réaction. Dans un deuxième temps, on teste la méthodologie d’estimation d’erreur sur un autre système d’équations, où des effets d’extinction et de réallumage sont ajoutés. Lorsqu’il n’y a pas d’advection, une analyse asymptotique rigoureuse montre l’existence de plusieurs régimes de combustion déjà observés dans les simulations numériques. Nous obtenons une approximation des paramètres de réallumage et d’extinction avec la courbe en «S», un graphe de la température maximale de la flamme en fonction du nombre de Damköhler, composée de trois branches et d’une double courbure. En ajoutant des effets advectifs, on obtient également une courbe en «S» correspondant aux régimes de combustion déjà identifiés. Nous comparons les erreurs de modélisation liées aux approximations asymptotiques dans les deux régimes stables et établissons une nouvelle hiérarchie des modèles en fonction du régime de combustion. Ces erreurs sont comparées aux estimations données par la stratégie d’estimation d’erreur. Si un seul régime stable de combustion existe, l’estimateur d’erreur l’identifie correctement ; si plus d’un régime est possible, on obtient une fac˛on systématique de choisir un régime. Pour les régimes où plus d’un modèle est approprié, la hiérarchie prédite par l’estimateur est correcte.

La(2-x)Ce(x)Cu(1-y)Zn(y)O(4) perovskites for high temperature water-gas shift reaction

The performance of La(2-x)Ce(x)Cu(1-y)Zn(y)O(4) perovskites as catalysts for the high temperature water-gas shift reaction (H T-W G S R) was investigated. The catalysts were characterized by EDS, XRD, BET surface area, TPR, and XANES. The results showed that all the perovskites exhibited the La(2)CuO(4) orthorhombic structure, so the Pechini method is suitable for the preparation of pure perovskite. However, the La(1.90)Ce(0.10)CuO(4) perovskite alone, when calcined at 350/700 degrees C, also showed a (La(0.935)Ce(0.065))(2)CuO(4) perovskite with tetragonal structure, which produced a surface area higher than the other perovskites. The perovskites that exhibited the best catalytic performance were those calcined at 350/700 degrees C and, among these, La(1.90)Ce(0.10)CuO(4) was outstanding, probably because of the high surface area associated with the presence of the (La(0.935)Ce(0.065))(2)CuO(4) perovskite with tetragonal structure and orthorhombic La(2)CuO(4) phase.

Effects of the partial replacement of La by M (M=Ce, Ca and Sr) in La(2-x)M(x)CuO(4) perovskites on catalysis of the water-gas shift reaction

The performance of La(2-x)M(x)CuO(4) perovskites (where M = Ce, Ca or Sr) as catalysts for the water-gas shift reaction was investigated at 290 degrees C and 360 degrees C. The catalysts were characterized by EDS, XRD, N(2) adsorption-desorption, XPS and XANES. The XRD results showed that all the perovskites exhibited a single phase (the presence of perovskite structure), suggesting the incorporation of metals in the perovskite structure. The XPS and XANES results showed the presence of Cu(2+) on the surface. The perovskites that exhibited the best catalytic performance were La(2-x)Ce(x)CuO(4) perovslcites, with CO conversions of 85%-90%. Moreover, these perovskites have higher surface areas and larger amounts of Cu on the surface. And Ce has a higher filled energy level than the other metals, increasing the energy of the valence band of Ce and providing more electrons for the reaction. Besides, the La(1.80)Ca(0.20)CuO(4) perovskite showed a good catalytic performance.