The actual impedance of non-reflecting boundary conditions : implications for the computation of resonators

Non-reflecting boundary conditions are essential elements in the computation of many compressible flows: such simulations are very sensitive to the treatment of acoustic waves at boundaries. Non-reflecting conditions allow acoustic waves to propagate through boundaries with zero or small levels of reflection into the domain. However, perfectly non-reflecting conditions must be avoided because they can lead to ill-posed problems for the mean flow. Various methods have been proposed to construct boundary conditions which can be sufficiently non-reflecting for the acoustic field while still making the mean-flow problem well posed. This paper analyses a widely-used technique for non-reflecting outlets (Rudy and Strikwerda, Poinsot and Lele). It shows that the correction introduced by these authors can lead to large reflection levels and non-physical resonant behaviors. A simple scaling is proposed to evaluate the relaxation coefficient used in theses methods for a non-reflecting outlet. The proposed scaling is tested for simple cases (ducts) both theoretically and numerically

Triple flame structure and diffusion flame stabilization

The stabilization of diffusion ñames is studied using asymptotic techniques and numerical tools. The configuration studied corresponda to parallel streams of cold oxidizer and fuel initially separated by a splitter píate. It is shown that stabilization of a diffusion flame may only occur in this situation by two processes. First, the flame may be stabilized behind the flame holder in the wake of the splitter píate. For this case, numerical simulations confirm scalings previously predicted by asymptotic analysis. Second, the flame may be lifted. In this case a triple flame is found at longer distanees downstream of the flame holder. The structure and propagation speed of this flame are studied by using an actively controlled numerical technique in which the triple flame is tracked in its own reference frame. It is then possible to investigate the triple flame structure and velocity. It is shown, as suggested from asymptotic analysis, that heat reléase may induce displacement speeds of the triple flame larger than the laminar flame speed corresponding to the stoichiometric conditions prevailing in the mixture approaching the triple flame. In addition to studying the characteristics of triple flames in a uniform flow, their re-sistance to turbulence is investigated by subjecting triple flames to different vortical configurations

Effects of pressure gradients on turbulent premixed flames

The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them

DNS and modeling of the interaction between turbulent premixed flames and walls

The interaction between turbulent premixed flames and walls is studied using a two-dimensional full Navier-Stokes solver with simple chemistry. The effects of wall distance on the local and global flame structure are investigated. Quenching distances and maximum wall heat fluxes during quenching are computed in laminar cases and are found to be comparable to experimental and analytical results. For turbulent cases, it is shown that quenching distances and maximum heat fluxes remain of the same order as for laminar flames. Based on simulation results, a 'law-of-the-wall' model is derived to describe the interaction between a turbulent premixed flame and a wall. This model is constructed to provide reasonable behavior of flame surface density near a wall under the assumption that flame-wall interaction takes place at scales smaller than the computational mesh. It can be implemented in conjunction with any of several recent flamelet models based on a modeled surface density equation, with no additional constraints on mesh size or time step

Large-eddy simulation and experimental study of heat transfer, nitric oxide emissions and combustion instability in a swirled turbulent high-pressure burner

Nitric oxide formation in gas turbine combustion depends on four key factors: flame stabilization, heat transfer, fuel-air mixing and combustion instability. The design of modern gas turbine burners requires delicate compromises between fuel efficiency, emissions of oxides of nitrogen (NOx) and combustion stability. Burner designs allowing substantial NOx reduction are often prone to combustion oscillations. These oscillations also change the NOx fields. Being able to predict not only the main species field in a burner but also the pollutant and the oscillation levels is now a major challenge for combustion modelling. This must include a realistic treatment of unsteady acoustic phenomena (which create instabilities) and also of heat transfer mechanisms (convection and radiation) which control NOx generation. In this work, large-eddy simulation (LES) is applied to a realistic gas turbine combustion chamber configuration where pure methane is injected through multiple holes in a cone-shaped burner. In addition to a non-reactive simulation, this article presents three reactive simulations and compares them to experimental results. The first reactive simulation neglects effects of cooling air on flame stabilization and heat losses by radiation and convection. The second reactive simulation shows how cooling air and heat transfer affect nitric oxide emissions. Finally, the third reactive simulation shows the effects of combustion instability on nitric oxide emissions. Additionally, the combustion instability is analysed in detail, including the evaluation of the terms in the acoustic energy equation and the identification of the mechanism driving the oscillation. Results confirm that LES of gas turbine combustion requires not only an accurate chemical scheme and realistic heat transfer models but also a proper description of the acoustics in order to predict nitric oxide emissions and pressure oscillation levels simultaneousl

Hydrodynamic instabilities in gaseous detonations: comparison of Euler, Navier–Stokes, and large-eddy simulation

A large-eddy simulation is conducted to investigate the transient structure of an unstable detonation wave in two dimensions and the evolution of intrinsic hydrodynamic instabilities. The dependency of the detonation structure on the grid resolution is investigated, and the structures obtained by large-eddy simulation are compared with the predictions from solving the Euler and Navier–Stokes equations directly. The results indicate that to predict irregular detonation structures in agreement with experimental observations the vorticity generation and dissipation in small scale structures should be taken into account. Thus, large-eddy simulation with high grid resolution is required. In a low grid resolution scenario, in which numerical diffusion dominates, the structures obtained by solving the Euler or Navier–Stokes equations and large-eddy simulation are qualitatively similar. When high grid resolution is employed, the detonation structures obtained by solving the Euler or Navier–Stokes equations directly are roughly similar yet equally in disagreement with the experimental results. For high grid resolution, only the large-eddy simulation predicts detonation substructures correctly, a fact that is attributed to the increased dissipation provided by the subgrid scale model. Specific to the investigated configuration, major differences are observed in the occurrence of unreacted gas pockets in the high-resolution Euler and Navier–Stokes computations, which appear to be fully combusted when large-eddy simulation is employed

Numerical simulations of turbulent premixed H2/O2/N2 flames with complex chemistry

Premixed stoichiometric H2/O2/N2 flames propagating in two-dimensional turbulence were studied using direct numerical simulation (simulations in which all fluid and thermochemical scales are fully resolved) including realistic chemical kinetics and molecular transport. Results are compared with earlier zero-chemistry (flame sheet) and one-step chemistry simulations. Consistent with the simpler models, the turbulent flame with realistic chemistry aligns preferentially with extensive strain rates in the tangent plane and flame curvature probability density functions are close to symmetric with near-zero means. By contrast to simple-chemistry results with non-unity Lewis numbers (ratio of thermal to species diffusivity), local flame structure does not correlate with curvature but rather with tangential strain rate. Turbulent straining results in substantial thinning of the flame relative to the steady unstrained laminar case. Heat release and H2O2 contours remain thin and connected ('flamelet-like') while species including H-atom and OH are more diffuse. Peak OH concentration occurs well behind the peak heat-release zone. The feasibility of incorporating realistic chemistry into full turbulence simulations to address issues such as pollutant formation in hydrocarbon-air flames is suggested

Simulation of a turbulent flame in a channel

The interaction between turbulent premixed flames and channel walls is studied. Combustion is represented by a simple irreversible reaction with a large activation temperature. Feedback to the flowfield is suppressed by invoking a constant density assumption. The effect of wall distance on local and global flame structure is investigated. Quenching distances and maximum wall heat fluxes computed in laminar cases are compared to DNS results. It is found that quenching distances decrease and maximum heat fluxes increase relative to laminar flame values. It is shown that these effects are due to large coherent structures which push flame elements towards to wall. The effect of wall strain is studied in flame-wall interaction in a stagnation line flow; this is used to explain the DNS results. It is also shown that 'remarkable' flame events are produced by interaction with a horseshoe vortex: burnt gases are pushed towards the wall at high speed and induce quenching and high wall heat fluxes while fresh gases are expelled from the wall region and form finger-like structures. Effects of the wall on flame surface density are investigated, and a simple model for flame-wall interaction is proposed; its predictions compare well with the DNS results