Prediction and control of combustion instabilities in real engines

This paper presents recent progress in the field of thermoacoustic combustion instabilities in propulsion engines such as rockets or gas turbines. Combustion instabilities have been studied for more than a century in simple laminar configurations as well as in laboratory-scale turbulent flames. These instabilities are also encountered in real engines but new mechanisms appear in these systems because of obvious differences with academic burners: larger Reynolds numbers, higher pressures and power densities, multiple inlet systems, complex fuels. Other differences are more subtle: real engines often feature specific unstable modes such as azimuthal instabilities in gas turbines or transverse modes in rocket chambers. Hydrodynamic instability modes can also differ as well as the combustion regimes, which can require very different simulation models. The integration of chambers in real engines implies that compressor and turbine impedances control instabilities directly so that the determination of the impedances of turbomachinery elements becomes a key issue. Gathering experimental data on combustion instabilities is difficult in real engines and Large Eddy Simulation (LES) has become a major tool in this field. Recent examples, however, show that LES is not sufficient and that theory, even in these complex systems, plays a major role to understand both experimental and LES results and to identify mitigation techniques

Curvature and confinement effects for flame speed measurements in laminar spherical and cylindrical flames.

This paper discusses methods used to obtain laminar flame speeds in spherical laminar premixed flames. Most recent studies express the laminar flame consumption speed as ρb/ρudR/dt, where R is the flame radius and ρb/ρu is the ratio of the burnt to the fresh gas density (ρb is evaluated at chemical equilibrium and supposed to be constant). This paper investigates the validity of this assumption by reconsidering it in a more general framework. Other formulae are derived and tested on a DNS of cylindrical flames (methane/air and octane/air). Results show that curvature and confinement effects lead to variations of ρb and ρu and to significant errors on the flame speed. Another expression (first proposed by Bradley and Mitcheson in 1976) is derived where no density evaluation is required and only pressure and flame radius evolution are used. It is shown to provide more precise results for the consumption speed than ρb/ρudR/dt because it takes into account curvature and confinement of the flame in the closed bomb

Large-Eddy Simulation of combustion instabilities in a variable-length combustor.

This article presents a simulation of a model rocket combustor with continuously variable acoustic properties thanks to a variable-length injector tube. Fully compressible Large-Eddy Simulations are conducted using the AVBP code. An original flame stabilization mechanism is uncovered where the recirculation of hot gases in the corner recirculation zone creates a triple flame structure. An unstable operating point is then chosen to investigate the mech- anism of the instability. The simulations are compared to experimental results in terms of frequency and mode structure. Two-dimensional axi-symmetric computations are com- pared to full 3D simulations in order to assess the validity of the axi-symmetry assumption for the prediction of mean and unsteady features of this flow. Despite the inaccuracies in- herent to the 2D description of a turbulent flow, for this configuration and the particular operating point investigated, the axi-symmetric simulation qualitatively reproduces some features of the instability

On the relative impact of subgrid-scale modelling and conjugate heat transfer in LES of hot jets in cross-flow over cold plates

This work describes numerical simulations of a hot jet in cross-flow with applications to anti-ice systems of aircraft engine nacelles. Numerical results are compared with experimental measurements obtained at ONERA to evaluate the performances of LES in this industrial context. The combination of complex geometries requiring unstructured meshes and high Reynolds number does not allow the resolution of boundary layers so that wall models must be employed. In this framework, the relative influence of subgrid-scale modelling and conjugate heat transfer in LESs of aerothermal flows is evaluated. After a general overview of the transverse jet simulation results, a LES coupled with a heat transfer solver in the walls is used to show that thermal boundary conditions at the wall have more influence on the results than subgrid scale models. Coupling fluid flow and heat transfer in solids simulations is the only method to specify their respective thermal boundary conditions

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

DNS and modeling of the turbulent boundary layer over an evaporating liquid film

In this paper, Direct Numercial Simulation (DNS) is used to study the boundary layer above an evaporating liquid film in the fully developed turbulent channel flow configuration. First, the numerical results of an isothermal flow are compared to the reference results of Kim et al. [1] to check the validity and the accuracy of the present approach. Next, the calculations are repeated for an anisothermal case where hot gas is flowing between cold walls. The numerical results corroborate those of Nicoud [2] and Huang & Coleman [3], introducing modified dimensionless variables. Finally, an evaporating liquid film is added at the walls. The complexity of the interaction between the evaporation process and the boundary layer structure, as well as its strong dependence on the transient thermophysical properties of the mixture are highlighted. As in the anisothermal case, the classical wall units are no longer adapted to build wall functions and new dimensionless variables are proposed. In addition a wall function must be developed for the evaporating species mass fraction, using a new dimensionless wall variable. It is shown that using these new variables allows to derive new wall functions for momentum, temperature and mass that lead to a correct description of the boundary layer when compared to DNS. These new wall functions may be directly implemented in CFD codes to take into account the impact of an evaporating liquid film

The combustion program at CTR

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct numerical simulation (DNS) appears to be one of the most attractive tools to use in addressing this problem. The general objective of DNS of reacting flows is to improve our knowledge of turbulent combustion but also to use this information for turbulent combustion models. For the foreseeable future, numerical simulation of the full three-dimensional governing partial differential equations with variable density and transport properties as well as complex chemistry will remain intractable; thus, various levels of simplification will remain necessary. On one hand, the requirement to simplify is not necessarily a handicap: numerical simulations allow the researcher a degree of control in isolating specific physical phenomena that is inaccessible in experiments. CTR has pursued an intensive research program in the field of DNS for turbulent reacting flows since 1987. DNS of reacting flows is quite different from DNS of non-reacting flows: without reaction, the equations to solve are clearly the five conservation equations of the Navier Stokes system for compressible situations (four for incompressible cases), and the limitation of the approach is the Reynolds number (or in other words the number of points in the computation). For reacting flows, the choice of the equations, the species (each species will require one additional conservation equation), the chemical scheme, and the configuration itself is more complex