Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Effect of H-2/CO Composition on Extinction Strain Rates of Counterflow Syngas Flames

This work presents a detailed experimental and numerical investigation of the effect of H-2/CO composition on extinction characteristics of premixed and nonpremixed syngas flames. Experimental measurements of local and global extinction strain rates in counterflow diffusion flames have been reported at atmospheric pressure for six different compositions of syngas fuel. The concentration of H-2 was varied from 5 to 20% with a 3% increment, and correspondingly, CO was decreased from 35 to 20% in steps of 3%. Particle imaging velocimetry has been used to determine the local extinction strain rates. Local extinction strain rates increased with an increase in the H-2/CO ratio in both nonpremixed and premixed flames. The predicted extinction strain rates for both nonpremixed and premixed counterflow flames using five different mechanisms available in the literature were compared with measurements. The Davis H-2/CO and Ranzi H-2/CO mechanisms predicted extinction strain rates within 10% of experimental values irrespective of the H-2/CO ratio. In the nonpremixed case, the Cl mechanism by Li et al., GRI 3.0, and the Ranzi H-2/CO mechanism predicted extinction strain rates well for low H-2/CO ratios (from 5:35 to 14:26) but deviated from experiments for higher H-2/CO values (17:23 and 20:20). In addition to kinetics, preferential diffusion effects were found to affect the reaction zone significantly and create distinct localized reaction zone structures in nonpremixed flames, which could contribute to discrepancies in extinction predictions.

Studies on a new high-intensity low-emission burner

This paper presents computational and experimental results on a new burner configuration with a mild combustion concept with heat release rates up to 10 MW/m(3). The burner configuration is shown to achieve mild combustion by using air at ambient temperature at high recirculation rates (similar to250%-290%) both experimentally and computationally. The principal features of the configuration are: (1) a burner with forward exit for exhaust gases; (2) injection of gaseous fuel and air as multiple, alternate, peripheral highspeed jets at the bottom at ambient temperature, thus creating high enough recirculation rates of the hot combustion products into fresh incoming reactants; and (3) use of a suitable geometric artifice-a frustum of a cone to help recirculation. The computational studies have been used to reveal the details of the flow and to optimize the combustor geometry based on recirculation rates. Measures, involving root mean square temperature fluctuations, distribution of temperature and oxidizer concentration inside the proposed burner, and a classical turbulent diffusion jet flame, are used to distinguish between them quantitatively. The system, operated at heat release rates of 2 to 10 MW/m(3) (compared to 0.02 to 0.32 MW/m(3) in the earlier studies), shows a 10-15 dB reduction in noise in the mild combustion mode compared to a simple open-top burner and exhaust NOx emission below 10 ppm for a 3 kW burner with 10% excess air. The peak temperature is measured around 1750 K, approximately 300 K lower than the peak temperature in a conventional burner.

A detailed numerical study of NOx kinetics in low calorific value H-2/CO syngas flames

The present study provides an extensive and detailed numerical analysis of NO chemical kinetics in low calorific value H-2/CO syngas flames utilizing predictions by five chemical kinetic mechanisms available out of which four deal with H-2/CO while the fifth mechanism (GRI 3.0) additionally accounts for hydrocarbon chemistry. Comparison of predicted axial NO profiles in premixed flat flames with measurements at 1 bar, 3.05 bar and 9.15 bar shows considerably large quantitative differences among the various mechanisms. However, at each pressure, the quantitative reaction path diagrams show similar NO formation pathways for most of the mechanisms. Interestingly, in counterflow diffusion flames, the quantitative reaction path diagrams and sensitivity analyses using the various mechanisms reveal major differences in the NO formation pathways and reaction rates of important reactions. The NNH and N2O intermediate pathways are found to be the major contributors for NO formation in all the reaction mechanisms except GRI 3.0 in syngas diffusion flames. The GRI 3.0 mechanism is observed to predict prompt NO pathway as the major contributing pathway to NO formation. This is attributed to prediction of a large concentration of CH radical by the GRI 3.0 as opposed to a relatively negligible value predicted by all other mechanisms. Also, the back-conversion of NNH into N2O at lower pressures (2-4 bar) was uniquely observed for one of the five mechanisms. The net reaction rates and peak flame temperatures are used to correlate and explain the differences observed in the peak NO] at different pressures. This study identifies key reactions needing assessment and also highlights the need for experimental data in syngas diffusion flames in order to assess and optimize H-2/CO and nitrogen chemistry. Copyright (C) 2014, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Influence of fuel hydrogen fraction on syngas fueled SI engine: Fuel thermo-physical property analysis and in-cylinder experimental investigations

Hydrogen, either in pure form or as a gaseous fuel mixture specie enhances the fuel conversion efficiency and reduce emissions in an internal combustion engine. This is due to the reduction in combustion duration attributed to higher laminar flame speeds. Hydrogen is also expected to increase the engine convective heat flux, attributed (directly or indirectly) to parameters like higher adiabatic flame temperature, laminar flame speed, thermal conductivity and diffusivity and lower flame quenching distance. These factors (adversely) affect the thermo-kinematic response and offset some of the benefits. The current work addresses the influence of mixture hydrogen fraction in syngas on the engine energy balance and the thermo-kinematic response for close to stoichiometric operating conditions. Four different bio-derived syngas compositions with fuel calorific value varying from 3.14 MJ/kg to 7.55 MJ/kg and air fuel mixture hydrogen fraction varying from 7.1% to 14.2% by volume are used. The analysis comprises of (a) use of chemical kinetics simulation package CHEMKIN for quantifying the thermo-physical properties (b) 0-D model for engine in-cylinder analysis and (c) in-cylinder investigations on a two-cylinder engine in open loop cooling mode for quantifying the thermo-kinematic response and engine energy balance. With lower adiabatic flame temperature for Syngas, the in-cylinder heat transfer analysis suggests that temperature has little effect in terms of increasing the heat flux. For typical engine like conditions (700 K and 25 bar at CR of 10), the laminar flame speed for syngas exceeds that of methane (55.5 cm/s) beyond mixture hydrogen fraction of 11% and is attributed to the increase in H based radicals. This leads to a reduction in the effective Lewis number and laminar flame thickness, potentially inducing flame instability and cellularity. Use of a thermodynamic model to assess the isolated influence of thermal conductivity and diffusivity on heat flux suggests an increase in the peak heat flux between 2% and 15% for the lowest (0.420 MW/m(2)) and highest (0.480 MW/m(2)) hydrogen containing syngas over methane (0.415 MW/m(2)) fueled operation. Experimental investigations indicate the engine cooling load for syngas fueled engine is higher by about 7% and 12% as compared to methane fueled operation; the losses are seen to increase with increasing mixture hydrogen fraction. Increase in the gas to electricity efficiency is observed from 18% to 24% as the mixture hydrogen fraction increases from 7.1% to 9.5%. Further increase in mixture hydrogen fraction to 14.2% results in the reduction of efficiency to 23%; argued due to the changes in the initial and terminal stages of combustion. On doubling of mixture hydrogen fraction, the flame kernel development and fast burn phase duration decrease by about 7% and 10% respectively and the terminal combustion duration, corresponding to 90%-98% mass burn, increases by about 23%. This increase in combustion duration arises from the cooling of the near wall mixture in the boundary layer attributed to the presence of hydrogen. The enhancement in engine cooling load and subsequent reduction in the brake thermal efficiency with increasing hydrogen fraction is evident from the engine energy balance along with the cumulative heat release profiles. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

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.

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.

Scalar transport in diffusion flame wrapped up by an air and fuel side vortex

The present work involves a computational study of soot (chosen as a scalar which is a primary pollutant source) formation and transport in a laminar acetylene diffusion flame perturbed by a convecting line vortex. The topology of soot contours resulting from flame vortex interactions has been investigated. More soot was produced when vortex was introduced from the air side in comparison to the fuel side. Also, the soot topography was spatially more diffuse in the case of air side vortex. The computational model was found to be in good agreement with the experimental work previously reported in the literature. The computational simulation enabled a study of various parameters like temperature, equivalence ratio and temperature gradient affecting the soot production and transport. Temperatures were found to be higher in the case of air side vortex in contrast to the fuel side one. In case of fuel side vortex, abundance of fuel in the vortex core resulted in fuel-rich combustion zone in the core and a more discrete soot topography. Besides, the overall soot production was observed to be low in the fuel side vortex. However, for the air side vortex, air abundance in the core resulted in higher temperatures and greater soot production. Probability density functions (PDFs) have been introduced to investigate the spatiotemporal variation of soot yield and transport and their dependence on temperature and acetylene concentration from statistical view point. In addition, the effect of flame curvature on soot production is also studied. The regions convex to fuel stream side witnessed thicker soot layer. All numerical simulations have been carried out on Fluent 6.3.26. (C) 2013 Elsevier Ltd. All rights reserved.

Experimental Analysis of Water Vapor Diffusion Through Porous Membranes in a Proton Exchange Membrane Fuel Cell

We report the diffusion characteristics of water vapor through two different porous media, viz., membrane electrode assembly (MEA) and gas diffusion layer (GDL) in a nonoperational fuel cell. Tunable diode laser absorption spectroscopy (TDLAS) was employed for measuring water vapor concentration in the test channel. Effects of the membrane pore size and the inlet humidity on the water vapor transport are quantified through mass flux and diffusion coefficient. Water vapor transport rate is found to be higher for GDL than for MEA. The flexibility and wide range of application of TDLAS in a fuel cell setup is demonstrated through experiments with a stagnant flow field on the dry side.

A Novel Operator-Splitting Technique for One-Dimensional Laminar Premixed Flames

This paper is concerned the calculation of flame structure of one-dimensional laminar premixed flames using the technique of operator-splitting. The technique utilizes an explicit method of solution with one step Euler for chemistry and a novel probabilistic scheme for diffusion. The relationship between diffusion phenomenon and Gauss-Markoff process is exploited to obtain an unconditionally stable explicit difference scheme for diffusion. The method has been applied to (a) a model problem, (b) hydrazine decomposition, (c) a hydrogen-oxygen system with 28 reactions with constant Dρ 2 approximation, and (d) a hydrogen-oxygen system (28 reactions) with trace diffusion approximation. Certain interesting aspects of behaviour of the solution with non-unity Lewis number are brought out in the case of hydrazine flame. The results of computation in the most complex case are shown to compare very favourably with those of Warnatz, both in terms of accuracy of results as well as computational time, thus showing that explicit methods can be effective in flame computations. Also computations using the Gear-Hindmarsh for chemistry and the present approach for diffusion have been carried out and comparison of the two methods is presented.