The internal structure of igniting turbulent sprays as revealed by complex chemistry DNS

A parametric study of spark ignition in a uniform monodisperse turbulent spray is performed with complex chemistry three-dimensional Direct Numerical Simulations in order to improve the understanding of the structure of the ignition kernel. The heat produced by the kernel increases with the amount of fuel evaporated inside the spark volume. Moreover, the heat sink by evaporation is initially higher than the heat release and can have a negative effect on ignition. With the sprays investigated, heat release occurs over a large range of mixture fractions, being high within the nominal flammability limits and finite but low below the lean flammability limit. The burning of very lean regions is attributed to the diffusion of heat and species from regions of high heat release, and from the spark, to lean regions. Two modes of spray ignition are reported. With a relatively dilute spray, nominally flammable material exists only near the droplets. Reaction zones are created locally near the droplets and have a non-premixed character. They spread from droplet to droplet through a very lean interdroplet spacing. With a dense spray, the hot spark region is rich due to substantial evaporation but the cold region remains lean. In between, a large surface of flammable material is generated by evaporation. Ignition occurs there and a large reaction zone propagates from the rich burned region to the cold lean region. This flame is wrinkled due to the stratified mixture fraction field and evaporative cooling. In the dilute spray, the reaction front curvature pdf contains high values associated with single droplet combustion, while in the dense spray, the curvature is lower and closer to the curvature associated with gaseous fuel ignition kernels. © 2011 The Combustion Institute.

Experimental aspects of partially premixed pulsating combustion

The pulsating combustion process has won interest in current research due to indications that its application in energy generation can offer several advantages, such as: fuel economy, reduced pollutants formation, increased rate of convective heat transfer and reduced investment, when compared with conventional techniques. An experimental study has been conducted with the objective of investigating the effects of combustion driven acoustic oscillations in the emission rates of combustion gases, especially carbon monoxide and nitrogen oxides. The experiments were conducted in a water-jacketed 1-m long by 25-cm internal diameter stainless steel vertical tube. The combustor operated with liquefied petroleum gas (LPG) in both oscillatory and non oscillatory conditions, under the same input conditions. Part of the reactant mixture was excited acoustically, before the burner exit, by a speaker positioned strategically. The burner was aligned with the chamber longitudinal axis and positioned at its bottom. The experiments were conducted for 0.16 g/s of LPG burning in stoichiometric equivalence ratio. The main conclusions were: a) the pulsating combustion process produces more uniform fuel/air profile than the non pulsating process, b) close to stoichiometric equivalence ratio the pulsating combustion process generates higher rates of NO x; c) the frequency has a strong influence in NO x emission, but the pressure amplitude has a weak influence; d) the presence of the acoustic field may change drastically the combustion gas emissions in diffusion flames, but in pre-mixed flames the influence is not as strong.

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.

Prediction of Upward Flame Spread over Polymers

In this work, the existing understanding of flame spread dynamics is enhanced through an extensive study of the heat transfer from flames spreading vertically upwards across 5 cm wide, 20 cm tall samples of extruded Poly (Methyl Methacrylate) (PMMA). These experiments have provided highly spatially resolved measurements of flame to surface heat flux and material burning rate at the critical length scale of interest, with a level of accuracy and detail unmatched by previous empirical or computational studies. Using these measurements, a wall flame model was developed that describes a flame’s heat feedback profile (both in the continuous flame region and the thermal plume above) solely as a function of material burning rate. Additional experiments were conducted to measure flame heat flux and sample mass loss rate as flames spread vertically upwards over the surface of seven other commonly used polymers, two of which are glass reinforced composite materials. Using these measurements, our wall flame model has been generalized such that it can predict heat feedback from flames supported by a wide range of materials. For the seven materials tested here – which present a varied range of burning behaviors including dripping, polymer melt flow, sample burnout, and heavy soot formation – model-predicted flame heat flux has been shown to match experimental measurements (taken across the full length of the flame) with an average accuracy of 3.9 kW m-2 (approximately 10 – 15 % of peak measured flame heat flux). This flame model has since been coupled with a powerful solid phase pyrolysis solver, ThermaKin2D, which computes the transient rate of gaseous fuel production of constituents of a pyrolyzing solid in response to an external heat flux, based on fundamental physical and chemical properties. Together, this unified model captures the two fundamental controlling mechanisms of upward flame spread – gas phase flame heat transfer and solid phase material degradation. This has enabled simulations of flame spread dynamics with a reasonable computational cost and accuracy beyond that of current models. This unified model of material degradation provides the framework to quantitatively study material burning behavior in response to a wide range of common fire scenarios.

Sandwich propellant combustion: Modeling and experimental comparison

This paper reports reacting fluid dynamics calculations for an ammonium percholrate binder sandwich and extracts experimentally observed features including surface profiles and maximum regression rates as a function of pressure and binder thickness. These studies have been carried out by solving the two-dimensional unsteady Navier-Stokes equations with energy and species conservation equations and a kinetic model of three reaction steps (ammonium perchlorate decomposition flame, primary diffusion flame, and final diffusion flame) in the gas phase. The unsteady two-dimensional conduction equation is solved in the condensed phase. The regressing surface is unsteady and two dimensional. Computations have been carried out for a binder thickness range of 25-125 mum and a pressure range of 1.4 to 6.9 MPa. Good comparisons at several levels of detail are used to demonstrate the need for condensed-phase two-dimensional unsteady conduction and three-step gas-phase reactions. The choice of kinetic and thermodynamic parameters is crucial to good comparison with experiments. The choice of activation energy parameters for ammonium percholrate combustion has been made with stability of combustion in addition to experimentally determined values reported in literature. The choice of gas-phase parameters for the diffusion flames are made considering that (a) primary diffusion flame affects the low-pressure behavior and (b) final diffusion flame affects high-pressure behavior. The predictions include the low-pressure deflagration limit of the sandwich apart from others noted above. Finally, this study demonstrates the possibility of making meaningful comparisons with experimental observations on sandwich propellant combustion.

Effect of port gas injection on in-cylinder stratification in a CNG fuelled SI engine

The concept of barrel stratification of air-fuel mixture is evaluated for a port gas injection (PGI) single cylinder spark ignition (SI) internal combustion (IC) engine using a transient three-dimensional computational fluid dynamic (CFD) model. The gaseous fuel used in the study is compressed natural gas (CNG). It is observed that compared to the premixed gas carburettor case, a substantial amount of in-cylinder stratification can be achieved with port gas injection system. A detailed parametric study is reported to understand the effect of the various injection parameters such as injection location, injection orientation, start of injection (SOT) and its duration, and injection rate. Furthermore, the best injection timing is evaluated for various load and speed cases. It is observed that the best stratification pattern can be achieved at 50% engine load. The injection location is observed to have a profound effect on the in-cylinder stratification pattern, and injection towards the side of the spark plug is observed to give a rich fuel-air mixture near the spark plug. It is also shown that there exists an optimal injection pressure.

First and second law thermodynamic analysis of air and oxy-steam biomass gasification

Gasification is an energy transformation process in which solid fuel undergoes thermochemical conversion to produce gaseous fuel, and the two most important criteria involved in such process to evaluate the performance, economics and sustainability of the technology are: the total available energy (exergy) and the energy conserved (energy efficiency). Current study focuses on the energy and exergy analysis of the oxy-steam gasification and comparing with air gasification to optimize the H-2 yield, efficiency and syngas energy density. Casuarina wood is used as a fuel, and mixture of oxygen and steam in different proportion and amount is used as a gasifying media. The results are analysed with respect to varying equivalence ratio and steam to biomass ratio (SBR). Elemental mass balance technique is employed to ensure the validity of results. First and second law thermodynamic analysis is used towards time evaluation of energy and exergy analysis. Different component of energy input and output has been studied carefully to understand the influence of varying SBR on the availability of energy and irreversibility in the system to minimize the losses with change in input parameters for optimum performance. The energy and exergy losses (irreversibility) for oxy-steam gasification system are compared with the results of air gasification, and losses are found to be lower in oxy-steam thermal conversion; which has been argued and reasoned due to the presence of N-2 in the air-gasification. The maximum exergy efficiency of 85% with energy efficiency of 82% is achieved at SBR of 0.75 on the molar basis. It has been observed that increase in SBR results in lower exergy and energy efficiency, and it is argued to be due to the high energy input in steam generation and subsequent losses in the form of physical exergy of steam in the product gas, which alone accounts for over 18% in exergy input and 8.5% in exergy of product gas at SBR of 2.7. Carbon boundary point (CBP), is identified at the SBR of 1.5, and water gas shift (WGS) reaction plays a crucial role in H-2 enrichment after carbon boundary point (CBP) is reached. Effects of SBR and CBP on the H-2/CO ratio is analysed and discussed from the perspective of energy as well as the reaction chemistry. Energy density of syngas and energy efficiency is favoured at lower SBR but higher SBR favours H-2 rich gas at the expense of efficiency. Copyright (C) 2014, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Carboxydothermus hydrogenoformans comme catalyseur biologique pour la conversion du monoxyde de carbone en hydrogène simultanément a la minéralisation de calcium et phosphate

La gazéification est aujourd'hui l'une des stratégies les plus prometteuses pour valoriser les déchets en énergie. Cette technologie thermo-chimique permet une réduction de 95 % de la masse des intrants et génère des cendres inertes ainsi que du gaz de synthèse (syngaz). Le syngaz est un combustible gazeux composé principalement de monoxyde de carbone (CO), d'hydrogène (H2) et de dioxyde de carbone (CO2). Le syngaz peut être utilisé pour produire de la chaleur et de l'électricité. Il est également la pierre angulaire d'un grand nombre de produits à haute valeur ajoutée, allant de l'éthanol à l'ammoniac et l'hydrogène pur. Les applications en aval de la production de syngaz sont dictées par son pouvoir calorifique, lui-même dépendant de la teneur du gaz en H2. L’augmentation du contenu du syngaz en H2 est rendu possible par la conversion catalytique à la vapeur d’eau, largement répandu dans le cadre du reformage du méthane pour la production d'hydrogène. Au cours de cette réaction, le CO est converti en H2 et CO2 selon : CO + H2O → CO2 + H2. Ce processus est possible grâce à des catalyseurs métalliques mis en contact avec le CO et de la vapeur. La conversion catalytique à la vapeur d’eau a jusqu'ici été réservé pour de grandes installations industrielles car elle nécessite un capital et des charges d’exploitations très importantes. Par conséquent, les installations de plus petite échelle et traitant des intrants de faible qualité (biomasse, déchets, boues ...), n'ont pas accès à cette technologie. Ainsi, la seule utilisation de leur syngaz à faible pouvoir calorifique, est limitée à la génération de chaleur ou, tout au plus, d'électricité. Afin de permettre à ces installations une gamme d’application plus vaste de leurs syngaz, une alternative économique à base de catalyseur biologique est proposée par l’utilisation de bactéries hyperthermophiles hydrogénogènes. L'objectif de cette thèse est d'utiliser Carboxydothermus hydrogenoformans, une bactérie thermophile carboxydotrophe hydrogénogène comme catalyseur biologique pour la conversion du monoxyde de carbone en hydrogène. Pour cela, l’impact d'un phénomène de biominéralisation sur la production d’H2 a été étudié. Ensuite, la faisabilité et les limites de l’utilisation de la souche dans un bioréacteur ont été évaluées. Tout d'abord, la caractérisation de la phase inorganique prédominante lorsque C. hydrogenoformans est inoculé dans le milieu DSMZ, a révélé une biominéralisation de phosphate de calcium (CaP) cristallin en deux phases. L’analyse par diffraction des rayons X et spectrométrie infrarouge à transformée de Fourier de ce matériau biphasique indique une signature caractéristique de la Mg-whitlockite, alors que les images obtenues par microscopie électronique à transmission ont montré l'existence de nanotiges cristallines s’apparentant à de l’hydroxyapatite. Dans les deux cas, le mode de biominéralisation semble être biologiquement induit plutôt que contrôlé. L'impact du précipité de CaP endogène sur le transfert de masse du CO et la production d’H2 a ensuite été étudié. Les résultats ont été comparés aux valeurs obtenues dans un milieu où aucune précipitation n'est observée. Dans le milieu DSMZ, le KLa apparent (0.22 ± 0.005 min-1) et le rendement de production d’H2 (89.11 ± 6.69 %) étaient plus élevés que ceux obtenus avec le milieu modifié (0.19 ± 0.015 min-1 et 82.60 ± 3.62% respectivement). La présence du précipité n'a eu aucune incidence sur l'activité microbienne. En somme, le précipité de CaP offre une nouvelle stratégie pour améliorer les performances de transfert de masse du CO en utilisant les propriétés hydrophobes de gaz. En second lieu, la conversion du CO en H2 par la souche Carboxydothermus hydrogenoformans fut étudiée et optimisée dans un réacteur gazosiphon de 35 L. Parmi toutes les conditions opérationnelles, le paramètre majeur fut le ratio du débit de recirculation du gaz sur le débit d'alimentation en CO (QR:Qin). Ce ratio impacte à la fois l'activité biologique et le taux de transfert de masse gaz-liquide. En effet, au dessus d’un ratio de 40, les performances de conversion du CO en H2 sont limitées par l’activité biologique alors qu’en dessous, elles sont limitées par le transfert de masse. Cela se concrétise par une efficacité de conversion maximale de 90.4 ± 0.3 % et une activité spécifique de 2.7 ± 0.4 molCO·g–1VSS·d–1. Malgré des résultats prometteurs, les performances du bioréacteur ont été limitées par une faible densité cellulaire, typique de la croissance planctonique de C. hydrogenoformans. Cette limite est le facteur le plus contraignant pour des taux de charge de CO plus élevés. Ces performances ont été comparées à celles obtenues dans un réacteur à fibres creuses (BRFC) inoculé par la souche. En dépit d’une densité cellulaire et d’une activité volumétrique plus élevées, les performances du BRFC à tout le moins cinétiquement limitées quand elles n’étaient pas impactées par le transfert de masse, l'encrassement et le vieillissement de la membrane. Afin de parer à la dégénérescence de C. hydrogenoformans en cas de pénurie de CO, la croissance de la bactérie sur pyruvate en tant que seule source de carbone a été également caractérisée. Fait intéressant, en présence simultanée de pyruvate et de CO, C. hydrogenoformans n’a amorcé la consommation de pyruvate qu’une fois le CO épuisé. Cela a été attribué à un mécanisme d'inhibition du métabolisme du pyruvate par le CO, faisant ainsi du pyruvate le candidat idéal pour un système in situ de secours.

Balanço energético de um forno túnel de cerâmica estrutural convertido de lenha para gás natural

In energy systems, the balance of entrances, exits and losses are fundamental to rationalize the energy consumption, independently of the source (sun, natural gas, wind, water, firewood or oil). This estimate is important so much in the phase of project of the facilities, as in the exploration or operation. In the project phase it indicates the energy needs of the process and the contribution of the energy in the cost of the product and the capacity of storage of the fuel and in the operation phase it allows to evaluate the use of the energy in the process of it burns, showing the weak points that should suffer intervention to improve the efficiency. With this tool, it can be implemented routines of calculation of thermal balances in ovens of it burns of structural ceramic, in way to generate an optimized mathematical model for application in the current and promising structural ceramic brazilian industry. The ceramic oven in study is located in the metropolitan area of Natal (Rio Grande do Norte) and it is a continuous oven of the type wagons tunnel, converted of firewood for natural gas and it produces blocks of red ceramic. The energy balance was applied in the oven tunnel before and after the conversion and made the comparisons of the energy efficiencies (it burns to the firewood and it burns to natural gas), what showed that the gaseous fuel is more efficient when we burn structural ceramic in ovens tunnels. When we burn natural gas, the requested energy is smaller and better used. Tests were accomplished in the burned product that showed the best quality of the burned brick with natural gas. That quality improvement makes possible to accomplish new interventions for the most rational use of the energy in the oven tunnel of the Ceramic in study and in the industries of structural ceramic of the whole Brazil, that need control tools of burning and of quality

A Generalized Methodology to Characterize Composite Materials for Pyrolysis Models

The predictive capabilities of computational fire models have improved in recent years such that models have become an integral part of many research efforts. Models improve the understanding of the fire risk of materials and may decrease the number of expensive experiments required to assess the fire hazard of a specific material or designed space. A critical component of a predictive fire model is the pyrolysis sub-model that provides a mathematical representation of the rate of gaseous fuel production from condensed phase fuels given a heat flux incident to the material surface. The modern, comprehensive pyrolysis sub-models that are common today require the definition of many model parameters to accurately represent the physical description of materials that are ubiquitous in the built environment. Coupled with the increase in the number of parameters required to accurately represent the pyrolysis of materials is the increasing prevalence in the built environment of engineered composite materials that have never been measured or modeled. The motivation behind this project is to develop a systematic, generalized methodology to determine the requisite parameters to generate pyrolysis models with predictive capabilities for layered composite materials that are common in industrial and commercial applications. This methodology has been applied to four common composites in this work that exhibit a range of material structures and component materials. The methodology utilizes a multi-scale experimental approach in which each test is designed to isolate and determine a specific subset of the parameters required to define a material in the model. Data collected in simultaneous thermogravimetry and differential scanning calorimetry experiments were analyzed to determine the reaction kinetics, thermodynamic properties, and energetics of decomposition for each component of the composite. Data collected in microscale combustion calorimetry experiments were analyzed to determine the heats of complete combustion of the volatiles produced in each reaction. Inverse analyses were conducted on sample temperature data collected in bench-scale tests to determine the thermal transport parameters of each component through degradation. Simulations of quasi-one-dimensional bench-scale gasification tests generated from the resultant models using the ThermaKin modeling environment were compared to experimental data to independently validate the models.

Effect of diffusion-layer porosity on the performance of polymer electrolyte fuel cells

The gas-diffusion layer (GDL) influences the performance of electrodes employed with polymer electrolyte fuel cells (PEFCs). A simple and effective method for incorporating a porous structure in the electrode GDL using sucrose as the pore former is reported. Optimal (50 w/o) incorporation of a pore former in the electrode GDL facilitates the access of the gaseous reactants to the catalyst sites and improves the fuel cell performance. Data obtained from permeability and porosity measurements, single-cell performance, and impedance spectroscopy suggest that an optimal porosity helps mitigating mass-polarization losses in the fuel cell resulting in a substantially enhanced performance.

Prediction of flame liftoff height of diffusion/partially premixed jet flames and modeling of mild combustion burners

In this article, a new flame extinction model based on the k/epsilon turbulence time scale concept is proposed to predict the flame liftoff heights over a wide range of coflow temperature and O-2 mass fraction of the coflow. The flame is assumed to be quenched, when the fluid time scale is less than the chemical time scale ( Da < 1). The chemical time scale is derived as a function of temperature, oxidizer mass fraction, fuel dilution, velocity of the jet and fuel type. The present extinction model has been tested for a variety of conditions: ( a) ambient coflow conditions ( 1 atm and 300 K) for propane, methane and hydrogen jet flames, ( b) highly preheated coflow, and ( c) high temperature and low oxidizer concentration coflow. Predicted flame liftoff heights of jet diffusion and partially premixed flames are in excellent agreement with the experimental data for all the simulated conditions and fuels. It is observed that flame stabilization occurs at a point near the stoichiometric mixture fraction surface, where the local flow velocity is equal to the local flame propagation speed. The present method is used to determine the chemical time scale for the conditions existing in the mild/ flameless combustion burners investigated by the authors earlier. This model has successfully predicted the initial premixing of the fuel with combustion products before the combustion reaction initiates. It has been inferred from these numerical simulations that fuel injection is followed by intense premixing with hot combustion products in the primary zone and combustion reaction follows further downstream. Reaction rate contours suggest that reaction takes place over a large volume and the magnitude of the combustion reaction is lower compared to the conventional combustion mode. The appearance of attached flames in the mild combustion burners at low thermal inputs is also predicted, which is due to lower average jet velocity and larger residence times in the near injection zone.

Scaling of turbulent flame speed for expanding flames with Markstein diffusion considerations

In this paper we clarify the role of Markstein diffusivity, which is the product of the planar laminar flame speed and the Markstein length, on the turbulent flame speed and its scaling, based on experimental measurements on constant-pressure expanding turbulent flames. Turbulent flame propagation data are presented for premixed flames of mixtures of hydrogen, methane, ethylene, n-butane, and dimethyl ether with air, in near-isotropic turbulence in a dual-chamber, fan-stirred vessel. For each individual fuel-air mixture presented in this work and the recently published iso-octane data from Leeds, normalized turbulent flame speed data of individual fuel-air mixtures approximately follow a Re-T,f(0.5) scaling, for which the average radius is the length scale and thermal diffusivity is the transport property of the turbulence Reynolds number. At a given Re-T,Re-f, it is experimentally observed that the normalized turbulent flame speed decreases with increasing Markstein number, which could be explained by considering Markstein diffusivity as the leading dissipation mechanism for the large wave number flame surface fluctuations. Consequently, by replacing thermal diffusivity with the Markstein diffusivity in the turbulence Reynolds number definition above, it is found that normalized turbulent flame speeds could be scaled by Re-T,M(0.5) irrespective of the fuel, equivalence ratio, pressure, and turbulence intensity for positive Markstein number flames.

Dynamics of thermal ignition of spray flames in mixing layers

Conditions are identified under which analyses of laminar mixing layers can shed light on aspects of turbulent spray combustion. With this in mind, laminar spray-combustion models are formulated for both non-premixed and partially premixed systems. The laminar mixing layer separating a hot-air stream from a monodisperse spray carried by either an inert gas or air is investigated numerically and analytically in an effort to increase understanding of the ignition process leading to stabilization of high-speed spray combustion. The problem is formulated in an Eulerian framework, with the conservation equations written in the boundary-layer approximation and with a one-step Arrhenius model adopted for the chemistry description. The numerical integrations unveil two different types of ignition behaviour depending on the fuel availability in the reaction kernel, which in turn depends on the rates of droplet vaporization and fuel-vapour diffusion. When sufficient fuel is available near the hot boundary, as occurs when the thermochemical properties of heptane are employed for the fuel in the integrations, combustion is established through a precipitous temperature increase at a well-defined thermal-runaway location, a phenomenon that is amenable to a theoretical analysis based on activation-energy asymptotics, presented here, following earlier ideas developed in describing unsteady gaseous ignition in mixing layers. By way of contrast, when the amount of fuel vapour reaching the hot boundary is small, as is observed in the computations employing the thermochemical properties of methanol, the incipient chemical reaction gives rise to a slowly developing lean deflagration that consumes the available fuel as it propagates across the mixing layer towards the spray. The flame structure that develops downstream from the ignition point depends on the fuel considered and also on the spray carrier gas, with fuel sprays carried by air displaying either a lean deflagration bounding a region of distributed reaction or a distinct double-flame structure with a rich premixed flame on the spray side and a diffusion flame on the air side. Results are calculated for the distributions of mixture fraction and scalar dissipation rate across the mixing layer that reveal complexities that serve to identify differences between spray-flamelet and gaseous-flamelet problems.