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.

Vaporization and collision modeling of liquid fuel sprays in a co-axial fuel and air pre-mixer

Droplet collision occurs frequently in regions where the droplet number density is high. Even for Lean Premixed and Pre-vaporized (LPP) liquid sprays, the collision effects can be very high on the droplet size distributions, which will in turn affect the droplet vaporization process. Hence, in conjunction with vaporization modeling, collision modeling for such spray systems is also essential. The standard O'Rourke's collision model, usually implemented in CFD codes, tends to generate unphysical numerical artifact when simulations are performed on Cartesian grid and the results are not grid independent. Thus, a new collision modeling approach based on no-time-counter method (NTC) proposed by Schmidt and Rutland is implemented to replace O'Rourke's collision algorithm to solve a spray injection problem in a cylindrical coflow premixer. The so called ``four-leaf clover'' numerical artifacts are eliminated by the new collision algorithm and results from a diesel spray show very good grid independence. Next, the dispersion and vaporization processes for liquid fuel sprays are simulated in a coflow premixer. Two liquid fuels under investigation are jet-A and Rapeseed Methyl Esters (RME). Results show very good grid independence in terms of SMD distribution, droplet number distribution and fuel vapor mass flow rate. A baseline test is first established with a spray cone angle of 90 degrees and injection velocity of 3 m/s and jet-A achieves much better vaporization performance than RME due to its higher vapor pressure. To improve the vaporization performance for both fuels, a series of simulations have been done at several different combinations of spray cone angle and injection velocity. At relatively low spray cone angle and injection velocity, the collision effect on the average droplet size and the vaporization performance are very high due to relatively high coalescence rate induced by droplet collisions. Thus, at higher spray cone angle and injection velocity, the results expectedly show improvement in fuel vaporization performance since smaller droplet has a higher vaporization rate. The vaporization performance and the level of homogeneity of fuel-air mixture can be significantly improved when the dispersion level is high, which can be achieved by increasing the spray cone angle and injection velocity. (C) 2012 Elsevier Ltd. All rights reserved.

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.

CFD Analysis of two phase flow with a prediction of torque demand of a Bi-lobe rotor of a mixing chamber

An attempt to study the fluid dynamic behavior of two phase flow comprising of solid and liquid with nearly equal density in a geometrical case that has an industrial significance in theareas like processing of polymers, food, pharma ceutical, paints. In this work,crystalline silica is considered as the dispersed medium in glycerin. In the CFD analysis carried out,the two phase components are considered to be premixed homogeneously at the initial state. The flow in a cylinder that has an axially driven bi-lobe rotor, a typical blender used in polymer industry for mixing or kneading to render the multi-component mixture to homogeneous condition is considered. A viscous, incompressible, isothermal flow is considered with an assumption that the components do not undergo any physical change and the solids are rigid and mix in fully wetting conditions. Silica with a particle diameter of 0.4 mm is considered and flow is analyzed for different mixing fractions. An industry standard CFD code is used for solving 3D-RANS equations. As the outcome of the study the torque demand by the bi-lobe rotor for different mixture fractions which are estimated show a behavioral consistency to the expected physical phenomena occurring in the domain considered.

Strategies for high efficiency and stability in biogas-fuelled small engines

This work assesses the performance of small biogas-fuelled engines and explores high-efficiency strategies for power generation in the very low power range of less than 1000 W. Experiments were performed on a small 95-cc, single-cylinder, four-stroke spark-ignition engine operating on biogas. The engine was operated in two modes, i.e., `premixed' and `fuel injection' modes, using both single and dual spark plug configurations. Measurements of in-cylinder pressure, crank angle, brake power, air and fuel flow rates, and exhaust emissions were conducted. Cycle-to-cycle variations in engine in-cylinder pressure and power were also studied and assessed quantitatively for various loading conditions. Results suggest that biogas combustion can be fairly sensitive to the ignition strategies thereby affecting the power output and efficiency. Further, results indicate that continuous fuel injection shows superior performance compared to the premixed case especially at low loads owing to possible charge stratification in the engine cylinder. Overall, this study has demonstrated for the first time that a combination of technologies such as lean burn, fuel injection, and dual spark plug ignition can provide highly efficient and stable operation in a biogas-fuelled small S.I. engine, especially in the low power range of 450-1000W. (C) 2014 Elsevier Inc. All rights reserved.

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.

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.

Optical diagnostics of fuel-air mixing and vortex formation in a cavity combustor

The current work reports optical diagnostic measurements of fuel-air mixing and vortex structure in a single cavity trapped vortex combustor (TVC). Specifically, the mixture fraction using acetone PLIF technique in the non-reacting flow, and PIV measurements in the reacting flow are reported for the first time in trapped vortex combustors. The fuel-air momentum flux ratio, where the air momentum corresponds to that entering the cavity through a specially-incorporated flow guide vane, is used to characterize the mixing. The acetone PLIF experiments show that at high momentum flux ratios, the fuel-air mixing in the cavity is very minimal and is enhanced as the momentum flux ratio reduces, due to a favourable vortex formation in the cavity. Stoichiometric mixture fraction surfaces show that the mixing causes the reaction surfaces to shift from non-premixed to partially-premixed stratified mixtures. PIV measurements conducted in the non-reacting flow in the cavity further reinforce this observation. The scalar dissipation rates of mixture fraction were compared with the contours of RMS of fluctuating velocity and showed very good agreement. The regions of maximum mixing are observed to be along the fuel air interface. Reacting flow Ply measurements which differ substantially from the non-reacting cases primarily because of the heat release from combustion and the resulting gas expansion show that the vortex is displaced from the centre of the cavity towards the guide vane. Overall, the measurements show interesting features of the flow including the presence of the dual cavity structure and lead to a clear understanding of the underlying physics of the cavity flow highlighting the importance of the fuel-air momentum ratio parameter. (C) 2014 Elsevier Inc. All rights reserved.

Absolute/Convective Instability Transition in a Backward Facing Step Combustor: Fundamental Mechanism and Influence of Density Gradient

Hydrodynamic instabilities of the flow field in lean premixed gas turbine combustors can generate velocity perturbations that wrinkle and distort the flame sheet over length scales that are smaller than the flame length. The resultant heat release oscillations can then potentially result in combustion instability. Thus, it is essential to understand the hydrodynamic instability characteristics of the combustor flow field in order to understand its overall influence on combustion instability characteristics. To this end, this paper elucidates the role of fluctuating vorticity production from a linear hydrodynamic stability analysis as the key mechanism promoting absolute/convective instability transitions in shear layers occurring in the flow behind a backward facing step. These results are obtained within the framework of an inviscid, incompressible, local temporal and spatio-temporal stability analysis. Vorticity fluctuations in this limit result from interaction between two competing mechanisms-(1) production from interaction between velocity perturbations and the base flow vorticity gradient and (2) baroclinic torque in the presence of base flow density gradients. This interaction has a significant effect on hydrodynamic instability characteristics when the base flow density and velocity gradients are colocated. Regions in the space of parameters characterizing the base flow velocity profile, i.e., shear layer thickness and ratio of forward to reverse flow velocity, corresponding to convective and absolute instability are identified. The implications of the present results on understanding prior experimental studies of combustion instability in backward facing step combustors and hydrodynamic instability in other flows such as heated jets and bluff body stabilized flames is discussed.

On flame-turbulence interaction in constant-pressure expanding flames

In this paper we present one of the first high-speed particle image velocimetry measurements to quantify flame-turbulence interaction in centrally-ignited constant-pressure premixed flames expanding in nearisotropic turbulence. Measurements of mean flow velocity and rms of fluctuating flow velocity are provided over a range of conditions both in the presence and absence of the flame. The distributions of stretch rate contributions from different terms such as tangential straining, normal straining and curvature are also provided. It is found that the normal straining displays non-Gaussian pdf tails whereas the tangential straining shows near Gaussian behavior. We have further tracked the motion of the edge points that reside and co-move with the edge of the flame kernel during its evolution in time, and found that within the measurement conditions, on average the persistence time scales of stretch due to pure curvature exceed that due to tangential straining by at least a factor of two. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

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.

ABSOLUTE/CONVECTIVE INSTABILITY TRANSITION IN A BACKWARD FACING STEP COMBUSTOR: FUNDAMENTAL MECHANISM AND INFLUENCE OF DENSITY GRADIENT

Hydrodynamic instabilities of the flow field in lean premixed gas turbine combustors can generate velocity perturbations that wrinkle and distort the flame sheet over length scales that are smaller than the flame length. The resultant heat release oscillations can then potentially result in combustion instability. Thus, it is essential to understand the hydrodynamic instability characteristics of the combustor flow field in order to understand its overall influence on combustion instability characteristics. To this end, this paper elucidates the role of fluctuating vorticity production from a linear hydrodynamic stability analysis as the key mechanism promoting absolute/convective instability transitions in shear layers occurring in the flow behind a backward facing step. These results are obtained within the framework of an inviscid, incompressible, local temporal and spatio-temporal stability analysis. Vorticity fluctuations in this limit result from interaction between two competing mechanisms - (1) production from interaction between velocity perturbations and the base flow vorticity gradient and (2) baroclinic torque in the presence of base flow density gradients. This interaction has a significant effect on hydrodynamic instability characteristics when the base flow density and velocity gradients are co-located. Regions in the space of parameters characterizing the base flow velocity profile, i.e. shear layer thickness and ratio of forward to reverse flow velocity, corresponding to convective and absolute instability are identified. The implications of the present results on prior observations of flow instability in other flows such as heated jets and bluff-body stabilized flames is discussed.

Mitigating instability by actuating the swirler in a combustor

Thermoacoustic instability in a lean premixed combustor is a major impediment towards reliable operation of gas turbine engines for both aerospace and land based applications. In this communication, we investigate the following concept: in a laboratory combustor, could the otherwise static swirler be actuated to a rotary motion, such that the higher intensity turbulence and higher swirl number generated in the flame stabilization region might alter the flame position, structure and thereby assist in mitigating thermoacoustic instabilities? Results obtained using microphone and high speed imaging, show prominent reductions in the amplitudes of the first mode of the thermoacoustically unstable flame, with increased rotation rate of the swirler. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Temperature measurement of reflected shock wave by using chemical indicator

This report describes a new method for measuring the temperature of the gas behind the reflected shock wave in shock tube, corresponding to the reservoir temperature of a shock tunnel, based on the chemical reaction of small amount of CF4 premixed in the test gas. The final product C2F4 is used as the temperature indicator, which is sampled and detected by a gas chromatography in the experiment. The detected concentration of C2F4 is correlated to the temperature of the reflected shock wave with the initial pressure P-1 and test time tau as parameters in the temperature range 3 300 K < T < 5 600 K, pressure range 5 kPa < P1 <12 kPa and tau similar or equal to 0.4 ms.

障碍物扰动对预混气点火的影响

在一根半开口管道中,以乙炔为燃料,改变当量比和进气速度,实验研究了障碍物扰动对预混气点火的影响.实验表明,在小当量比情况下,在没有障碍物的光管中,随着进气速度的增加,点火所需要的当量比逐渐减小,回火区域逐渐增大.而添加障碍物后,回火区域的位置和点火所需要的当量比相对来说比较稳定,基本不随进气速度的增加而大幅度变化.并且由于障碍物引起的扰动使得点火所需要的当量比要比无障碍物时大得多.