On the formation of nitrogen oxides in the combustion processes of hydrogen in air

A research programme is being carried out at the Institute Nacional de Tecnica Aeroespacial of Spain, on several aspects of the formation of nitrogen oxides in continuous flow combustion systems, considering hydrogen and hydrocarbons as fuels. The research programme is fundamentally oriented on the basic aspects of the problem, although it also includes the study of the influence on the formation process of several operational and design variables of the combusters, such as type of fuels, fuel/air ratio, degree of mixing in premixed type flames, existence of droplets as compared with homogeneous combustion.This problem of nitrogen oxides formation is receiving lately great attention, specially in connection with automobile reciprocating engines and aircraft gas turbines. This is due to the fact of the increasing frequency and intensity of photochemical hazes or smog, typical of urban areas submitted to strong solar radiation, which are originated by the action on organic compounds of the oxidants resulting from the photochemical decomposition of nitrogen dioxide N02. In the combustion process almost all nitrogen oxides are in form of NO. This nitric oxide reacts with the oxygen of the air and forms N02, this reaction only taking place in or near the exhaust of tne motors, since the N0-02 reaction becomes frozen for the concentration existing in the atmosphere.

Flammability conditions for ultra-lean hydrogen premixed combustion based on flame-ball analyses

It has been reasoned that the structures of strongly cellular flames in very lean mixtures approach an array of flame balls, each burning as if it were isolated, thereby indicating a connection between the critical conditions required for existence of steady flame balls and those necessary for occurrence of self-sustained premixed combustion. This is the starting assumption of the present study, in which structures of near-limit steady sphericosym-metrical flame balls are investigated with the objective of providing analytic expressions for critical combustion conditions in ultra-lean hydrogen-oxygen mixtures diluted with N2 and water vapor. If attention were restricted to planar premixed flames, then the lean-limit mole fraction of H2 would be found to be roughly ten percent, more than twice the observed flammability limits, thereby emphasizing the relevance of the flame-ball phenomena. Numerical integrations using detailed models for chemistry and radiation show that a onestep chemical-kinetic reduced mechanism based on steady-state assumptions for all chemical intermediates, together with a simple, optically thin approximation for water-vapor radiation, can be used to compute near-limit fuel-lean flame balls with excellent accuracy. The previously developed one-step reaction rate includes a crossover temperature that determines in the first approximation a chemical-kinetic lean limit below which combustión cannot occur, with critical conditions achieved when the diffusion-controlled radiation-free peak temperature, computed with account taken of hydrogen Soret diffusion, is equal to the crossover temperature. First-order corrections are found by activation-energy asymptotics in a solution that involves a near-field radiation-free zone surrounding a spherical flame sheet, together with a far-field radiation-conduction balance for the temperature profile. Different scalings are found depending on whether or not the surrounding atmosphere contains wáter vapor, leading to different analytic expressions for the critical conditions for flame-ball existence, which give results in very good agreement with those obtained by detailed numerical computations.

Numerical simulation of group combustion of pulverized coal

A mathematical model for the group combustion of pulverized coal particles was developed in a previous work. It includes the Lagrangian description of the dehumidification, devolatilization and char gasification reactions of the coal particles in the homogenized gaseous environment resulting from the three fuels, CO, H2 and volatiles, supplied by the gasification of the particles and their simultaneous group combustion by the gas phase oxidation reactions, which are considered to be very fast. This model is complemented here with an analysis of the particle dynamics, determined principally by the effects of aerodynamic drag and gravity, and its dispersion based on a stochastic model. It is also extended to include two other simpler models for the gasification of the particles: the first one for particles small enough to extinguish the surrounding diffusion flames, and a second one for particles with small ash content when the porous shell of ashes remaining after gasification of the char, non structurally stable, is disrupted. As an example of the applicability of the models, they are used in the numerical simulation of an experiment of a non-swirling pulverized coal jet with a nearly stagnant air at ambient temperature, with an initial region of interaction with a small annular methane flame. Computational algorithms for solving the different stages undergone by a coal particle during its combustion are proposed. For the partial differential equations modeling the gas phase, a second order finite element method combined with a semi-Lagrangian characteristics method are used. The results obtained with the three versions of the model are compared among them and show how the first of the simpler models fits better the experimental results.

Simulation of an Integrated Gasification Combined Cycle with Chemical-Looping Combustion and CO2 sequestration

Chemical-looping combustion allows an integration of CO2 capture in a thermal power plant without energy penalty; secondly, a less exergy destruction in the combustion chemical transformation is achieved, leading to a greater overall thermal efficiency. This paper focus on the study of the energetic performance of this concept of combustion in an integrated gasification combined cycle power plant when synthesis gas is used as fuel for the gas turbines. After thermodynamic modelling and optimization of some cycle parameters, the power plant performance is evaluated under diverse working conditions and compared to a conventional integrated gasification combined cycle with precombustion capture. Energy savings in CO2 capture and storage has been quantified. The overall efficiency increase is found to be significant and even notable, reaching values of around 7%. In order to analyze the influence of syngas composition on the results, different H2-content fuels are considered.

A model for turbulent combustion simulation of large scale hydrogen explosions

El agotamiento, la ausencia o, simplemente, la incertidumbre sobre la cantidad de las reservas de combustibles fósiles se añaden a la variabilidad de los precios y a la creciente inestabilidad en la cadena de aprovisionamiento para crear fuertes incentivos para el desarrollo de fuentes y vectores energéticos alternativos. El atractivo de hidrógeno como vector energético es muy alto en un contexto que abarca, además, fuertes inquietudes por parte de la población sobre la contaminación y las emisiones de gases de efecto invernadero. Debido a su excelente impacto ambiental, la aceptación pública del nuevo vector energético dependería, a priori, del control de los riesgos asociados su manipulación y almacenamiento. Entre estos, la existencia de un innegable riesgo de explosión aparece como el principal inconveniente de este combustible alternativo. Esta tesis investiga la modelización numérica de explosiones en grandes volúmenes, centrándose en la simulación de la combustión turbulenta en grandes dominios de cálculo en los que la resolución que es alcanzable está fuertemente limitada. En la introducción, se aborda una descripción general de los procesos de explosión. Se concluye que las restricciones en la resolución de los cálculos hacen necesario el modelado de los procesos de turbulencia y de combustión. Posteriormente, se realiza una revisión crítica de las metodologías disponibles tanto para turbulencia como para combustión, que se lleva a cabo señalando las fortalezas, deficiencias e idoneidad de cada una de las metodologías. Como conclusión de esta investigación, se obtiene que la única estrategia viable para el modelado de la combustión, teniendo en cuenta las limitaciones existentes, es la utilización de una expresión que describa la velocidad de combustión turbulenta en función de distintos parámetros. Este tipo de modelos se denominan Modelos de velocidad de llama turbulenta y permiten cerrar una ecuación de balance para la variable de progreso de combustión. Como conclusión también se ha obtenido, que la solución más adecuada para la simulación de la turbulencia es la utilización de diferentes metodologías para la simulación de la turbulencia, LES o RANS, en función de la geometría y de las restricciones en la resolución de cada problema particular. Sobre la base de estos hallazgos, el crea de un modelo de combustión en el marco de los modelos de velocidad de la llama turbulenta. La metodología propuesta es capaz de superar las deficiencias existentes en los modelos disponibles para aquellos problemas en los que se precisa realizar cálculos con una resolución moderada o baja. Particularmente, el modelo utiliza un algoritmo heurístico para impedir el crecimiento del espesor de la llama, una deficiencia que lastraba el célebre modelo de Zimont. Bajo este enfoque, el énfasis del análisis se centra en la determinación de la velocidad de combustión, tanto laminar como turbulenta. La velocidad de combustión laminar se determina a través de una nueva formulación capaz de tener en cuenta la influencia simultánea en la velocidad de combustión laminar de la relación de equivalencia, la temperatura, la presión y la dilución con vapor de agua. La formulación obtenida es válida para un dominio de temperaturas, presiones y dilución con vapor de agua más extenso de cualquiera de las formulaciones previamente disponibles. Por otra parte, el cálculo de la velocidad de combustión turbulenta puede ser abordado mediante el uso de correlaciones que permiten el la determinación de esta magnitud en función de distintos parámetros. Con el objetivo de seleccionar la formulación más adecuada, se ha realizado una comparación entre los resultados obtenidos con diversas expresiones y los resultados obtenidos en los experimentos. Se concluye que la ecuación debida a Schmidt es la más adecuada teniendo en cuenta las condiciones del estudio. A continuación, se analiza la importancia de las inestabilidades de la llama en la propagación de los frentes de combustión. Su relevancia resulta significativa para mezclas pobres en combustible en las que la intensidad de la turbulencia permanece moderada. Estas condiciones son importantes dado que son habituales en los accidentes que ocurren en las centrales nucleares. Por ello, se lleva a cabo la creación de un modelo que permita estimar el efecto de las inestabilidades, y en concreto de la inestabilidad acústica-paramétrica, en la velocidad de propagación de llama. El modelado incluye la derivación matemática de la formulación heurística de Bauwebs et al. para el cálculo de la incremento de la velocidad de combustión debido a las inestabilidades de la llama, así como el análisis de la estabilidad de las llamas con respecto a una perturbación cíclica. Por último, los resultados se combinan para concluir el modelado de la inestabilidad acústica-paramétrica. Tras finalizar esta fase, la investigación se centro en la aplicación del modelo desarrollado en varios problemas de importancia para la seguridad industrial y el posterior análisis de los resultados y la comparación de los mismos con los datos experimentales correspondientes. Concretamente, se abordo la simulación de explosiones en túneles y en contenedores, con y sin gradiente de concentración y ventilación. Como resultados generales, se logra validar el modelo confirmando su idoneidad para estos problemas. Como última tarea, se ha realizado un analisis en profundidad de la catástrofe de Fukushima-Daiichi. El objetivo del análisis es determinar la cantidad de hidrógeno que explotó en el reactor número uno, en contraste con los otros estudios sobre el tema que se han centrado en la determinación de la cantidad de hidrógeno generado durante el accidente. Como resultado de la investigación, se determinó que la cantidad más probable de hidrogeno que fue consumida durante la explosión fue de 130 kg. Es un hecho notable el que la combustión de una relativamente pequeña cantidad de hidrogeno pueda causar un daño tan significativo. Esta es una muestra de la importancia de este tipo de investigaciones. Las ramas de la industria para las que el modelo desarrollado será de interés abarca la totalidad de la futura economía de hidrógeno (pilas de combustible, vehículos, almacenamiento energético, etc) con un impacto especial en los sectores del transporte y la energía nuclear, tanto para las tecnologías de fisión y fusión. ABSTRACT The exhaustion, absolute absence or simply the uncertainty on the amount of the reserves of fossil fuels sources added to the variability of their prices and the increasing instability and difficulties on the supply chain are strong incentives for the development of alternative energy sources and carriers. The attractiveness of hydrogen in a context that additionally comprehends concerns on pollution and emissions is very high. Due to its excellent environmental impact, the public acceptance of the new energetic vector will depend on the risk associated to its handling and storage. Fromthese, the danger of a severe explosion appears as the major drawback of this alternative fuel. This thesis investigates the numerical modeling of large scale explosions, focusing on the simulation of turbulent combustion in large domains where the resolution achievable is forcefully limited. In the introduction, a general description of explosion process is undertaken. It is concluded that the restrictions of resolution makes necessary the modeling of the turbulence and combustion processes. Subsequently, a critical review of the available methodologies for both turbulence and combustion is carried out pointing out their strengths and deficiencies. As a conclusion of this investigation, it appears clear that the only viable methodology for combustion modeling is the utilization of an expression for the turbulent burning velocity to close a balance equation for the combustion progress variable, a model of the Turbulent flame velocity kind. Also, that depending on the particular resolution restriction of each problem and on its geometry the utilization of different simulation methodologies, LES or RANS, is the most adequate solution for modeling the turbulence. Based on these findings, the candidate undertakes the creation of a combustion model in the framework of turbulent flame speed methodology which is able to overcome the deficiencies of the available ones for low resolution problems. Particularly, the model utilizes a heuristic algorithm to maintain the thickness of the flame brush under control, a serious deficiency of the Zimont model. Under the approach utilized by the candidate, the emphasis of the analysis lays on the accurate determination of the burning velocity, both laminar and turbulent. On one side, the laminar burning velocity is determined through a newly developed correlation which is able to describe the simultaneous influence of the equivalence ratio, temperature, steam dilution and pressure on the laminar burning velocity. The formulation obtained is valid for a larger domain of temperature, steam dilution and pressure than any of the previously available formulations. On the other side, a certain number of turbulent burning velocity correlations are available in the literature. For the selection of the most suitable, they have been compared with experiments and ranked, with the outcome that the formulation due to Schmidt was the most adequate for the conditions studied. Subsequently, the role of the flame instabilities on the development of explosions is assessed. Their significance appears to be of importance for lean mixtures in which the turbulence intensity remains moderate. These are important conditions which are typical for accidents on Nuclear Power Plants. Therefore, the creation of a model to account for the instabilities, and concretely, the acoustic parametric instability is undertaken. This encloses the mathematical derivation of the heuristic formulation of Bauwebs et al. for the calculation of the burning velocity enhancement due to flame instabilities as well as the analysis of the stability of flames with respect to a cyclic velocity perturbation. The results are combined to build a model of the acoustic-parametric instability. The following task in this research has been to apply the model developed to several problems significant for the industrial safety and the subsequent analysis of the results and comparison with the corresponding experimental data was performed. As a part of such task simulations of explosions in a tunnel and explosions in large containers, with and without gradient of concentration and venting have been carried out. As a general outcome, the validation of the model is achieved, confirming its suitability for the problems addressed. As a last and final undertaking, a thorough study of the Fukushima-Daiichi catastrophe has been carried out. The analysis performed aims at the determination of the amount of hydrogen participating on the explosion that happened in the reactor one, in contrast with other analysis centered on the amount of hydrogen generated during the accident. As an outcome of the research, it was determined that the most probable amount of hydrogen exploding during the catastrophe was 130 kg. It is remarkable that the combustion of such a small quantity of material can cause tremendous damage. This is an indication of the importance of these types of investigations. The industrial branches that can benefit from the applications of the model developed in this thesis include the whole future hydrogen economy, as well as nuclear safety both in fusion and fission technology.

Preparation, characterisation and testing of CuO/Ce0.8Zr0.2O2 catalysts for NO oxidation to NO2 and mild temperature diesel soot combustion

CuO/ceria-zirconia catalysts have been prepared, deeply characterised (N2 adsorption–desorption isotherms at −196 °C, XRD, Raman spectroscopy, XPS, TEM and H2-TPR) and tested for NO oxidation to NO2 in TPR conditions, and for soot combustion at mild temperature (400 °C) in a NOx/O2 stream. The behaviour has been compared to that of a reference Pt/alumina commercial catalyst. The ceria-zirconia support was prepared by the co-precipitation method, and different amounts of copper (0.5, 1, 2, 4 and 6 wt%) were loaded by incipient wetness impregnation. The results revealed that copper is well-dispersed onto the ceria-zirconia support for the catalysts with low copper loading and CuO particles were only identified by XRD in samples with 4 and 6% of copper. A very low loading of copper increases significantly the activity for the NO oxidation to NO2 with regard to the ceria-zirconia support and an optimum was found for a 4% CuO/ceria-zirconia composition, showing a very high activity (54% at 348 °C). The soot combustion rate at 400 °C obtained with the 2% CuO/ceria-zirconia catalyst is slightly lower to that of 1% Pt/alumina in terms of mass of catalyst but higher in terms of price of catalyst.

Formation of brominated pollutants during the pyrolysis and combustion of tetrabromobisphenol A at different temperatures

Tetrabromobisphenol A (TBBPA) is the most widely used brominated flame retardant worldwide. A detailed examination of the degradation products emitted during thermal decomposition of TBBPA is presented in the study. Runs were performed in a laboratory furnace at different temperatures (650 and 800 °C) and in different atmospheres (nitrogen and air). More than one hundred semivolatile compounds have been identified by GC/MS, with special interest in brominated ones. Presence of HBr and brominated light hydrocarbons increased with temperature and in the presence of oxygen. Maximum formation of PAHs is observed at pyrolytic condition at the higher temperature. High levels of 2,4-, 2,6- and 2,4,6- bromophenols were found. The levels of polybrominated dibenzo-p-dioxins and furans have been detected in the ppm range. The most abundant isomers are 2,4,6,8-TeBDF in pyrolysis and 1,2,3,7,8-PeBDF in combustion. These results should be considered in the assessment of thermal treatment of materials containing brominated flame retardants.

Behavior of different soot combustion catalysts under NOx/O2. Importance of the catalyst–soot contact

Four different catalysts (Pt/Al2O3, Ce0.8Zr0.2O2, PrO2−x and SrTiCuO3) have been investigated on a laboratory scale to evaluate their potential as diesel soot combustion catalysts under different experimental conditions, which simulate the situation found in a continuous regeneration technology trap (dual-bed configuration of catalyst and soot) or a catalyst-coated filter system (single-bed configuration, both catalyst and soot particles mixed under loose-contact mode). Under dual-bed configuration, the behavior of the catalysts towards soot combustion are very similar, despite the differences observed in the NO2 production profiles. However, under single-bed configuration, there are important differences in the soot combustion activities and in the NO2 slip profiles. The configurations chosen have an enormous impact on CO/(CO + CO2) ratios of combustion products as well. The most active catalyst under NOx + O2 is PrO2−x combining a high contribution of active oxygen-assisted soot combustion as well as high NO2 production activity along the catalytic bed.

Kinetics of the combustion of olive oil. A semi-global model

Support for this work was provided by PROMETEO/2009/043/FEDER of Generalitat Valenciana (Spain) and CTQ2008-05520 (Spanish MCI/research).

Evaluation of hectorites, synthesized in different conditions, as soot combustion catalysts after impregnation with copper

Two microporous hectorites were prepared by conventional and microwave heating, and a delaminated mesoporous hectorite by an ultrasound-assisted synthesis. These three hectorites were impregnated with copper. The characterization techniques used were XRD, N2 adsorption, TEM and H2 reduction after selective surface copper oxidation by N2O (to determine copper dispersion). The catalytic activity for soot combustion of the copper-free and the copper-containing hectorites was tested under a gas mixture of 500 ppm NOx/5% O2/N2 (and 5% O2/N2 in some cases), evaluating their stability through three consecutive soot combustion experiments. The delaminated hectorite showed the highest surface area (353 m2/g) allowing the highest dispersion of copper. This copper-containing catalyst was the most active for soot combustion among those prepared and tested in this study. We have also concluded that the Cu/hectorite-catalyzed soot combustion mechanism is based on the activation of the O2 molecule and not on the NO2-assisted soot combustion.

Pollutant emissions during pyrolysis and combustion of waste printed circuit boards, before and after metal removal

The constant increase in the production of electronic devices implies the need for an appropriate management of a growing number of waste electrical and electronic equipment. Thermal treatments represent an interesting alternative to recycle this kind of waste, but particular attention has to be paid to the potential emissions of toxic by-products. In this study, the emissions from thermal degradation of printed circuit boards (with and without metals) have been studied using a laboratory scale reactor, under oxidizing and inert atmosphere at 600 and 850 °C. Apart from carbon oxides, HBr was the main decomposition product, followed by high amounts of methane, ethylene, propylene, phenol and benzene. The maximum formation of PAHs was found in pyrolysis at 850 °C, naphthalene being the most abundant. High levels of 2-, 4-, 2,4-, 2,6- and 2,4,6-bromophenols were found, especially at 600 °C. Emissions of PCDD/Fs and dioxin-like PCBs were quite low and much lower than that of PBDD/Fs, due to the higher bromine content of the samples. Combustion at 600 °C was the run with the highest PBDD/F formation: the total content of eleven 2,3,7,8-substituted congeners (tetra- through heptaBDD/Fs) was 7240 and 3250 ng WHO2005-TEQ/kg sample, corresponding to the sample with and without metals, respectively.

Catalytic performance of CuO/Ce0.8Zr0.2O2 loaded onto SiC-DPF in NOx-assisted combustion of diesel soot

This work presents a comparative study between the catalytic performance of the 2% CuO/ceria-zirconia powder catalyst and the same catalyst supported on silicon carbide DPF (Diesel Particulate Filter) towards NO oxidation reaction and soot combustion reaction. The ceria-zirconia catalyst was prepared by the co-precipitation method and 2 wt% copper was incorporated by the incipient wetness impregnation method. The catalyst was incorporated onto the ceramic support using a simple and organic solvent-free procedure by a simply dipping the DPF into an aqueous solution of the catalyst. The powder catalyst has been characterized using N2 adsorption at −196 °C, XRD and Raman Spectroscopy; whereas the catalytic coating morphology has been evaluated by SEM and the mechanical stability by an adherence test. Both catalyst configurations were tested for NO oxidation to NO2 and for soot combustion under NOx/O2. The results revealed that incorporation of the very active copper/ceria-zirconia catalyst onto SiC-DPF has been successfully achieved by a simple coating procedure. Furthermore, the catalytic coating has shown suitable mechanical, chemical and thermal stability. A satisfactory catalytic performance of the catalytic-coated filter was reached towards the NO oxidation reaction. Moreover, it was proved that the catalytic coating is stable and the corresponding coated DPF can be reused for several cycles of NO oxidation without a significant decrease in its activity. Finally, it was verified that the loose-contact mode is a good choice to simulate the catalytic performance of this active phase in a real diesel particulate filter.

Emissions from the Pyrolysis and Combustion of Different Wastes

Paper submitted to the 7th International Symposium on Feedstock Recycling of Polymeric Materials (7th ISFR 2013), New Delhi, India, 23-26 October 2013.

Pyrolysis and combustion study of flexible polyurethane foam

The thermal degradation of flexible polyurethane foam has been studied under different conditions by thermogravimetric analysis (TG), thermogravimetric analysis-infrared spectrometry (TG-IR) and thermogravimetric analysis-mass spectrometry (TG-MS). For the kinetic study, dynamic and dynamic+isothermal runs were performed at different heating rates (5, 10 and 20 °C min−1) in three different atmospheres (N2, N2:O2 4:1 and N2:O2 9:1). Two reaction models were obtained, one for the pyrolysis and another for the combustion degradation (N2:O2 4:1 and N2:O2 9:1), simultaneously correlating the experimental data from the dynamic and dynamic+isothermal runs at different heating rates. The pyrolytic model considered consisted of two consecutive reactions with activation energies of 142 and 217.5 kJ mol−1 and reaction orders of 0.805 and 1.246. Nevertheless, to simulate the experimental data from the combustion runs, three consecutive reactions were employed with activation energies of 237.9, 103.5 and 120.1 kJ mol−1, and reaction orders of 2.003, 0.778 and 1.025. From the characterization of the sample employing TG-IR and TG-MS, the results obtained showed that the FPUF, under an inert atmosphere, started the decomposition breaking the urethane bond to produce long chains of ethers which were degraded immediately in the next step. However, under an oxidative atmosphere, at the first step not only the urethane bonds were broken but also some ether polyols started their degradation which finished at the second step producing a char that was degraded at the last stage.

Analysis of pollutants evolved from combustion and pyrolysis of flexible polyurethane foam in a laboratory furnace

Resumen del póster presentado en PIC2015 – the 14th International Congress on Combustion By-Products and Their Health Effects, Umeå, Sweden, 14-17 June 2015.