In premixed turbulent combustion, the modelling of the turbulent flux of the mean reaction progress variable remains somewhat controversial. Classical gradient transport assumptions based on the eddy viscosity concept are often used while both experimental data and theoretical analysis have pointed out the existence of countergradient turbulent diffusion. Direct numerical simulation (DNS) is used in this paper to provide basic information on the turbulent flux of and study the occurrence of counter-gradient transport. The numerical configuration corresponds to twoor three-dimensional premixed flames in isotropic turbulent flow. The simulations correspond to various flame and flow conditions that are representative of flamelet combustion. They reveal that different flames will feature different turbulent transport properties and that these differences can be related to basic dynamical differences in the flame-flow interactions: counter-gradient diffusion occurs when the flow field near the flame is dominated by thermal dilatation due to chemical reaction, whereas gradient diffusion occurs when the flow field near the flame is dominated by the turbulent motions. The DNS-based analysis leads to a simple expression to describe the turbulent flux of , which in turn leads to a simple criterion to delineate between the gradient and counter-gradient turbulent diffusion regimes. This criterion suggests that the occurrence of one regime or the other is determined primarily by the ratio of turbulence intensity divided by the laminar flame speed, and by the flame heat release factor, τ ≡ (Tb — Tu)/Tu, where Tu and Tb are respectively the temperature within unburnt and burnt gas. Consistent with the Bray-Moss-Libby theory, counter-gradient (gradient) diffusion is promoted by low (high) values and high (low) values of τ. DNS also shows that these results are not restricted to the turbulent transport of . Similar results are found for the turbulent transport of flame surface density, Σ. The turbulent fluxes of and Σ are strongly correlated in the simulated flames and counter-gradient (gradient) diffusion of always coincides with counter-gradient (gradient) diffusion of Σ.
One basic effect of turbulence in turbulent premixed combustion is for the fluctuating velocity field to wrinkle the flame and greatly increase its surface area. In the flamelet theory, this effect is described by the flame surface density. An exact evolution equation for the flame surface density, called the Σ-equation, may be written, where basic physical mechanisms like production by hydrodynamic straining and destruction by propagation effects are described explicitly. Direct numerical simulation (DNS) is used in this paper to estimate the different terms appearing in the Σ-equation. The numerical configuration corresponds to three-dimensional premixed flames in isotropic turbulent flow. The simulations are performed for various mixture Lewis numbers in order to modify the strength and nature of the flame-flow coupling. The DNS-based analysis provides much information relevant to flamelet models. In particular, the flame surface density, and the source and sink terms for the flame surface density, are resolved spatially across the turbulent flame brush. The geometry as well as the dynamics of the flame differ quite significantly from one end of the reaction zone to the other. For instance, contrary to the intuitive idea that flame propagation effects merely counteract the wrinkling due to the turbulence, the role of flame propagation is not constant across the turbulent brush and switches from flame surface production at the front to flame surface dissipation at the back. Direct comparisons with flamelet models are also performed. The Bray-Moss-Libby assumption that the flame surface density is proportional to the flamelet crossing frequency, a quantity that can be measured in experiments, is found to be valid. Major uncertainties remain, however, over an appropriate description of the flamelet crossing frequency. In comparison, the coherent flame model of Marble & Broadwell achieves closure at the level of the Σ-equation and provides a more promising physically based description of the flame surface dynamics. Some areas where the model needs improvement are identified.
Combustion instability is investigated in the case of a multiple inlet combustor with dump. It is shown that low-frequency instabilities are acoustically coupled and occur at the eigenfrequencies of the system. Using spark-schlieren and a special phase-average imaging of the C2-radical emission, the fluid-mechanical processes involved in a vortex-driven mode of instability are investigated. The phase-average images provide maps of the local non-steady heat release. From the data collected on the combustor the processes of vortex shedding, growth, interactions and burning are described. The phases between the pressure, velocity and heat-release fluctuations are determined. The implications of the global Rayleigh criterion are verified and a mechanism for low-frequency vortex-driven instabilities is proposed.
Low-frequency combustion instabilities are studied in a model ramjet combustor facility. The facility is two-dimensional, and is comprised of a long inlet duct, a dump combustor cavity with variable size capability, and an exhaust nozzle. The flame is observed to be unstable over a wide range of operating conditions. Acoustic pressure and velocity measurements are made at various locations in the system. They show that the inlet duct acts as a long-wavelength acoustic resonator. However, the instability frequency does not lock to any particular value. This result suggests that the instability mechanism is not purely acoustic in nature. Schlieren imaging reveals that the instability is associated with large-scale flame-front motions which are driven by periodic vortex shedding at the instability frequency. Vortices are generated at the dump in phase with the acoustic velocity fluctuations in the inlet duct. The unsteady heat addition process closely follows the vortex history: the vortices form, grow in size, convect through the combustor cavity, impinge on the exhaust nozzle, break down to small scales and burn. C2 and CH radical spectroscopy is used to determine the phase relation between heat release and pressure in the reaction zone. Rayleigh's criterion is thereby shown to be satisfied. Next, the crucial question of how the oscillation frequency is determined is addressed. Inlet velocity and combustor length are systematically varied to assess the role of vortices by modification of their characteristic lifetime. The influence of the acoustic feedback time is also studied by shortening the inlet duct. The results show that the instability frequency is controlled by both vortex kinetics in the combustor and acoustic response of the inlet section. Therefore, the instability may be considered as a mixed acoustic-convective mode. Finally, combining Rayleigh's criterion with a global feedback loop equation, it is found that the resonant frequencies are selected according to the restriction \[ \frac{1}{4N-1} < \frac{\tau_{\rm v}}{\tau_{\rm f}} < \frac{3}{4N-3}, \] where N is the mode of oscillation and τv is the time for vortices to be convected from inlet to exhaust with τf being the feedback time taken for a pressure disturbance to travel up the inlet system and back.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.