Modelling soot formation in a DISI engine

In this work, the formation of soot in a Direct Injection Spark Ignition (DISI) engine is simulated using the Stochastic Reactor Model (SRM) engine code. Volume change, convective heat transfer, turbulent mixing, direct injection and flame propagation are accounted for. In order to simulate flame propagation, the cylinder is divided into an unburned, entrained and burned zone, with the rate of entrainment being governed by empirical equations but combustion modelled with chemical kinetics. The model contains a detailed chemical mechanism as well as a highly detailed soot formation model, however computation times are relatively short. The soot model provides information on the morphology and chemical composition of soot aggregates along with bulk quantities, including soot mass, number density, volume fraction and surface area. The model is first calibrated by simulating experimental data from a Gasoline Direct Injection (GDI) Spark Ignition (SI) engine. The model is then used to simulate experimental data from the literature, where the numbers, sizes and derived mass particulate emissions from a 1.83 L, 4-cylinder, 4 valve production DISI engine were examined. Experimental results from different injection and spark timings are compared with the model and the qualitative trends in aggregate size distribution and emissions match the exhaust gas measurements well. © 2010 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Effects of flow-geometry and Damkohler number on turbulence-flame interaction

DNS of turbulent hydrogen-air premixed flame is conducted for freely propagating and V-flames, using complex chemical kinetics. The results are analysed to study the influence of flame configuration on the turbulence-scalar interaction, which is critical for the scalar gradient generation process. The result suggests that this interaction process is not influenced by the flame configuration and the flame normal is found to predominantly align with the most extensive strain in the region of intense heat release.

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.

DNS of EGR-type turbulent flame in MILD condition

Three-dimensional direct numerical simulation (DNS) of exhaust gas recirculation (EGR)-type turbulent combustion operated in moderate and intense low-oxygen dilution (MILD) condition has been carried out to study the flame structure and flame interaction. In order to achieve adequate EGR-type initial/inlet mixture fields, partially premixed mixture fields which are correlated with the turbulence are carefully preprocessed. The chemical kinetics is modelled using a skeletal mechanism for methane-air combustion. The results suggest that the flame fronts have thin flame structure and the direct link between the mean reaction rate and scalar dissipation rate remains valid in the EGR-type combustion with MILD condition. However, the commonly used canonical flamelet is not fully representative for MILD combustion. During the flame-flame interactions, the heat release rate increases higher than the maximum laminar flame value, while the gradient of progress variable becomes smaller than laminar value. It is also proposed that the reaction rate and the scalar gradient can be used as a marker for the flame interaction. © 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

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.

Modelling of turbulent lifted jet flames using flamelets: A priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z-c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z-c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences. © 2014 Taylor & Francis.

Kinetics of the chemical looping oxidation of H2 by a co-precipitated mixture of CuO and Al2O3

Chemical-looping combustion (CLC) has the inherent property of separating 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. When used for the combustion of gaseous fuels, such as natural gas, or synthesis gas from the gasification of coal, the technique gives a stream of CO2 which, on an industrial scale, would be sufficiently pure for geological sequestration. An important issue is the form of the metal oxide, since it must retain its reactivity through many cycles of complete reduction and oxidation. Here, we report on the rates of oxidation of one constituent of synthesis gas, H2, by co-precipitated mixtures of CuO+Al2O3 using a laboratory-scale fluidised bed. To minimise the influence of external mass transfer, and also of errors in the measurement of [H2], particles sized to 355-500μm were used at low [H2], with the temperature ranging from 450 to 900°C. Under such conditions, the reaction was slow enough for meaningful measurements of the intrinsic kinetics to be made. The reaction was found to be first order with respect to H2. Above ∼800°C, the reaction of CuO was fast and conformed to the shrinking core mechanism, proceeding via the intermediate, Cu2O, in: 2CuO+H2→Cu2O+H2O, ΔH1073 K0=- 116.8 kJ/mol; Cu2O+H2→2Cu+H2O, ΔH1073 K0-80.9 kJ/mol. After oxidation of the products Cu and Cu2O back to CuO, the kinetics in subsequent cycles of chemical looping oxidation of H2 could be approximated by those in the first. Interestingly, the carrier was found to react at temperatures as low as 300°C. The influence of the number of cycles of reduction and oxidation is explored. Comparisons are drawn with previous work using reduction by CO. Finally, these results indicate that the kinetics of reaction of the oxygen carrier with gasifier synthesis gases is very much faster than rates of gasification of the original fuel. © 2010 The Institution of Chemical Engineers.

An investigation of the kinetics of CO2 uptake by a synthetic calcium based sorbent

This study examines the kinetics of carbonation by CO2 at temperatures of ca. 750 °C of a synthetic sorbent composed of 15 wt% mayenite (Ca12Al14O33) and CaO, designated HA-85-850, and draws comparisons with the carbonation of a calcined limestone. In-situ XRD has verified the inertness of mayenite, which neither interacts with the active CaO nor does it significantly alter the CaO carbonation–calcination equilibrium. An overlapping grain model was developed to predict the rate and extent of carbonation of HA-85-850 and limestone. In the model, the initial microstructure of the sorbent was defined by a discretised grain size distribution, assuming spherical grains. The initial input to the model – the size distribution of grains – was a fitted parameter, which was in good agreement with measurements made with mercury porosimetry and by the analysis of SEM images of sectioned particles. It was found that the randomly overlapping spherical grain assumption offered great simplicity to the model, despite its approximation to the actual porous structure within a particle. The model was able to predict the performance of the materials well and, particularly, was able to account for changes in rate and extent of reaction as the structure evolved after various numbers of cycles of calcination and carbonation.

An investigation of the kinetics of CO 2 uptake by a synthetic calcium based sorbent

This study examines the kinetics of carbonation by CO 2 at temperatures of ca. 750°C of a synthetic sorbent composed of 15wt% mayenite (Ca 12Al 14O 33) and CaO, designated HA-85-850, and draws comparisons with the carbonation of a calcined limestone. In-situ XRD has verified the inertness of mayenite, which neither interacts with the active CaO nor does it significantly alter the CaO carbonation-calcination equilibrium. An overlapping grain model was developed to predict the rate and extent of carbonation of HA-85-850 and limestone. In the model, the initial microstructure of the sorbent was defined by a discretised grain size distribution, assuming spherical grains. The initial input to the model - the size distribution of grains - was a fitted parameter, which was in good agreement with measurements made with mercury porosimetry and by the analysis of SEM images of sectioned particles. It was found that the randomly overlapping spherical grain assumption offered great simplicity to the model, despite its approximation to the actual porous structure within a particle. The model was able to predict the performance of the materials well and, particularly, was able to account for changes in rate and extent of reaction as the structure evolved after various numbers of cycles of calcination and carbonation. © 2011 Elsevier Ltd.

The effect of addition of ZrO2 to Fe2O3 for hydrogen production by chemical looping

In this paper, a synthetic mixture of ZrO2 and Fe 2O3 was prepared by coprecipitation for use in chemical looping and hydrogen production. Cycling experiments in a fluidized bed showed that a material composed of 30 mol % ZrO2 and 70 mol % Fe 2O3 was capable of producing hydrogen with a consistent yield of 90 mol % of the stoichiometric amount over 20 cycles of reduction and oxidation at 1123 K. Here, the iron oxide was subjected to cycles consisting of nearly 100% reduction to Fe followed by reoxidation (with steam or CO 2 and then air) to Fe2O3. There was no contamination by CO of the hydrogen produced, at a lower detection limit of 500 ppm, when the conversion of Fe3O4 to Fe was kept below 90 mol %. A preliminary investigation of the reaction kinetics confirmed that the ZrO2 support does not inhibit rates of reaction compared with those observed with iron oxide alone. © 2012 American Chemical Society.