skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Oscillatory response of an idealized two-dimensional diffusion flame: Analytical and numerical study

Abstract

The heat release response of a two-dimensional (2-D) co-flow diffusion flame is theoretically investigated for the mass fraction fluctuations of reactant concentration at the inlet boundary and for time-varying spatially uniform flow velocity field in the domain. The present work is an extension of the Burke-Schumann steady flame model to an oscillatory situation, but in a 2-D framework. The governing equation of the problem is the scalar advection equation for the Schwab-Zel'dovich variable. An exact solution is found in terms of an infinite series for the case of mixture fraction fluctuations at the inlet boundary. The cases of time-varying uniform flow velocity field and a combination of velocity and mixture fraction fluctuations are investigated numerically. The temperature and heat release rate of the flame are thermodynamically calculated utilizing the mixed-is-burnt approach. The main results of the paper are the response functions. The nonlinearity in the calculation of the heat release rate and the convective nonlinearity due to velocity fluctuations result in the generation of higher harmonics in the response function for a given sinusoidal excitation. Therefore, the response function is decomposed and obtained for each of the significant harmonics. The results show that, in general, the response function decreases withmore » increase in the excitation frequency as reported with premixed flames in the literature. However, in the present case, the decrease occurs when the excitation time-scale is less than the diffusion time-scale. Interesting flame shape variations such as flame clip-off, flipping between overexpanded and underexpanded conditions, and flame wrinkling are observed in the case of mixture fraction oscillations. (author)« less

Authors:
 [1]; ;  [2]
  1. Department of Mechanical Engineering, Indian Institute of Technology Madras, Chennai 600036 (India)
  2. Department of Aerospace Engineering, Indian Institute of Technology Madras, Chennai 600036 (India)
Publication Date:
OSTI Identifier:
20880663
Resource Type:
Journal Article
Resource Relation:
Journal Name: Combustion and Flame; Journal Volume: 149; Journal Issue: 3; Other Information: Elsevier Ltd. All rights reserved
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; RESPONSE FUNCTIONS; FLUCTUATIONS; FLAMES; MIXTURES; NUMERICAL SOLUTION; DIFFUSION; VELOCITY; HARMONICS; EXCITATION; HEAT; TWO-DIMENSIONAL CALCULATIONS; NONLINEAR PROBLEMS; EXACT SOLUTIONS; EQUATIONS; ADVECTION; SCALARS; SHAPE; SERIES EXPANSION; MASS; OSCILLATIONS; COMBUSTION KINETICS; FREQUENCY DEPENDENCE

Citation Formats

Tyagi, Manav, Jamadar, Nachiket, and Chakravarthy, S.R. Oscillatory response of an idealized two-dimensional diffusion flame: Analytical and numerical study. United States: N. p., 2007. Web. doi:10.1016/J.COMBUSTFLAME.2006.12.020.
Tyagi, Manav, Jamadar, Nachiket, & Chakravarthy, S.R. Oscillatory response of an idealized two-dimensional diffusion flame: Analytical and numerical study. United States. doi:10.1016/J.COMBUSTFLAME.2006.12.020.
Tyagi, Manav, Jamadar, Nachiket, and Chakravarthy, S.R. Tue . "Oscillatory response of an idealized two-dimensional diffusion flame: Analytical and numerical study". United States. doi:10.1016/J.COMBUSTFLAME.2006.12.020.
@article{osti_20880663,
title = {Oscillatory response of an idealized two-dimensional diffusion flame: Analytical and numerical study},
author = {Tyagi, Manav and Jamadar, Nachiket and Chakravarthy, S.R.},
abstractNote = {The heat release response of a two-dimensional (2-D) co-flow diffusion flame is theoretically investigated for the mass fraction fluctuations of reactant concentration at the inlet boundary and for time-varying spatially uniform flow velocity field in the domain. The present work is an extension of the Burke-Schumann steady flame model to an oscillatory situation, but in a 2-D framework. The governing equation of the problem is the scalar advection equation for the Schwab-Zel'dovich variable. An exact solution is found in terms of an infinite series for the case of mixture fraction fluctuations at the inlet boundary. The cases of time-varying uniform flow velocity field and a combination of velocity and mixture fraction fluctuations are investigated numerically. The temperature and heat release rate of the flame are thermodynamically calculated utilizing the mixed-is-burnt approach. The main results of the paper are the response functions. The nonlinearity in the calculation of the heat release rate and the convective nonlinearity due to velocity fluctuations result in the generation of higher harmonics in the response function for a given sinusoidal excitation. Therefore, the response function is decomposed and obtained for each of the significant harmonics. The results show that, in general, the response function decreases with increase in the excitation frequency as reported with premixed flames in the literature. However, in the present case, the decrease occurs when the excitation time-scale is less than the diffusion time-scale. Interesting flame shape variations such as flame clip-off, flipping between overexpanded and underexpanded conditions, and flame wrinkling are observed in the case of mixture fraction oscillations. (author)},
doi = {10.1016/J.COMBUSTFLAME.2006.12.020},
journal = {Combustion and Flame},
number = 3,
volume = 149,
place = {United States},
year = {Tue May 15 00:00:00 EDT 2007},
month = {Tue May 15 00:00:00 EDT 2007}
}
  • The authors consider an idealized two-dimensional model of the sintering of a close packed array of cylindrical grains with equivalent planar grain boundaries (GB) midway between each pair of cylinders. The scale and temperature are assumed to be such that material transport is completely dominated by diffusion along the GB's and on the pore surfaces; lattice diffusion is therefore neglected. This requires the condition D[sub gb]a[sub 0]/D[sub L]L >> 1, where D[sub gb] and D[sub L] are respectively the grain boundary and the lattice diffusion coefficients, a[sub 0] is an atomic dimension and L is the order of the cylindermore » radius. For particles sizes of 1 micrometer or less, of interest in the sintering of very fine powders, this condition is satisfied except possibly very near the melting point. The treatment takes into account an applied hydrostatic pressure during sintering but assumes that it is insufficient to cause plastic flow. Further, they ignore complications such as grain growth and the concomitant dragging and/or abandonment of pores by GB's, as well as the coarsening of pores. Possible extensions of the model to estimate the effect of these complications are discussed in the last section. As material moves from the GB's into the pores, thereby permitting the movement of the cylinders together, it will tend to alter the shape of the pores. They shall show, however, that it is reasonable to assume that surface diffusion is sufficiently rapid to redistribute the material and maintain the pore in the shape it would have in static equilibrium.« less
  • The authors investigate the sintering of two touching circular particles by surface and grain boundary diffusion. Typical examples for the evolution of the shape of the particles, their surface curvatures, and their surface fluxes are given. The sintering kinetics are evaluated as a function of the dihedral angle at the grain boundary-surface junctions and the grain boundary to surface diffusivity ratio. In particular, the growth rates of the neck between the two particles, the growth rate exponents, and the changes in the lengths of the particle pairs are monitored. The times needed to reach certain fractions of the final equilibriummore » neck sizes are tabulated for typical experimental dihedral angles and diffusivity ratios. The authors simulation is based on a rigorous mathematical system modeling the sintering of the two particles, and a rigorous numerical method for solving this system is adopted.« less
  • We study the ground-state properties of trapped inhomogeneous systems of hardcore bosons in two- and three-dimensional lattices. We obtain results both numerically, using quantum Monte Carlo techniques, and via several analytical approximation schemes, such as the Gutzwiller mean-field approach, a cluster mean-field method, and a spin-wave analysis which takes quantum fluctuations into account. We first study the homogeneous case, for which simple analytical expressions are obtained for all observables of interest, and compare the results with the numerical ones. We obtain the equation of state of the system along with other thermodynamic properties such as the free energy, kinetic energy,more » superfluid density, condensate density, and compressibility. In the presence of a trap, there is in general a spatial coexistence of superfluid and insulating domains. We show that the spin-wave-based method reproduces the quantum Monte Carlo results for global as well as for local quantities with a high degree of accuracy. We also discuss the validity of the local density approximation. Our analysis can be used to describe bosons in optical lattices where the onsite interaction U is much larger than the hopping amplitude t.« less
  • An approximate evaluation of the pair distribution and the structure factor is performed analytically for the two-dimensional, one-component plasma at any value of the coupling constant Gamma. The approximate distribution remains positive and satisfies three sum rules, including the compressibility one. When Gamma ..-->.. 0 or Gamma ..-->.. 2, exact results are found. At Gamma = 2 the transition from monotonic (Gamma < 2) to oscillatory (Gamma > 2) decay of correlations takes place. Comparison with the Monte Carlo simulations shows good agreement for 0 < Gamma < 4.
  • The burning characteristics of a premixed, H{sub 2}/air Bunsen-type flame are investigated using a time-dependent, axisymmetric numerical model with variable transport properties and a detailed-chemical-kinetics mechanism. The temperature, species concentration, and velocity fields are investigated under fuel-lean, stoichiometric, and fuel-rich conditions. The calculations show that under fuel-lean conditions the flame exhibits the ``tip-opening`` phenomenon, while under fuel-rich condition the tip of the flame burns intensely. These results are in agreement with the experimental findings of Mizomoto et al. who have suggested that the tip-opening phenomenon results from the nonunity Lewis number. To further investigate the impact of local Lewis numbermore » on the premixed-flame structure, numerical experiments are performed by modifying the local Lewis numbers of the individual species. Results for the fuel-lean condition confirm that the local Lewis number is responsible for the tip-opening phenomenon. Indeed, when the local Lewis number is set equal to 2, the burning pattern of the fuel-lean premixed flame resembles that of a fuel-rich flame with a closed tip. The spatial distributions of NO in the fuel-lean, stoichiometric, and fuel-rich flames are also examined. Under the fuel-lean and stoichiometric conditions, the NO is formed along the high-temperature cone of the flame, as expected. In the fuel-rich case, a dual flame structure is observed. The NO production occurs primarily in the secondary ``diffusion`` flame which is established at the interface of the excess fuel and ambient oxygen. Buoyancy-induced toroidal vortices are found to form in the vertically mounted premixed flames. However, the dynamic characteristics of a premixed flame, in contrast to those of a jet diffusion flame, are observed to be dependent on the inlet velocity profile of the fuel jet.« less