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

Title: Response analysis of a laminar premixed M-flame to flow perturbations using a linearized compressible Navier-Stokes solver

Abstract

The response of a laminar premixed methane-air flame subjected to flow perturbations around a steady state is examined experimentally and using a linearized compressible Navier-Stokes solver with a one-step chemistry mechanism to describe combustion. The unperturbed flame takes an M-shape stabilized both by a central bluff body and by the external rim of a cylindrical nozzle. This base flow is computed by a nonlinear direct simulation of the steady reacting flow, and the flame topology is shown to qualitatively correspond to experiments conducted under comparable conditions. The flame is then subjected to acoustic disturbances produced at different locations in the numerical domain, and its response is examined using the linearized solver. This linear numerical model then allows the componentwise investigation of the effects of flow disturbances on unsteady combustion and the feedback from the flame on the unsteady flow field. It is shown that a wrinkled reaction layer produces hydrodynamic disturbances in the fresh reactant flow field that superimpose on the acoustic field. This phenomenon, observed in several experiments, is fully interpreted here. The additional perturbations convected by the mean flow stem from the feedback of the perturbed flame sheet dynamics onto the flow field by a mechanism similar tomore » that of a perturbed vortex sheet. The different regimes where this mechanism prevails are investigated by examining the phase and group velocities of flow disturbances along an axis oriented along the main direction of the flow in the fresh reactant flow field. It is shown that this mechanism dominates the low-frequency response of the wrinkled shape taken by the flame and, in particular, that it fully determines the dynamics of the flame tip from where the bulk of noise is radiated.« less

Authors:
 [1];  [2];  [3];  [4];  [5]
  1. LadHyX, CNRS and Ecole Polytechnique, 91128 Palaiseau (France)
  2. CNRS, UPR 288, Laboratoire d’Energétique Moléculaire et Macroscopique Combustion (EM2C), Grande Voie des Vignes, 92290 Châtenay-Malabry (France)
  3. (France)
  4. ONERA-DAFE, 8 rue des Vertugadins, 92190 Meudon (France)
  5. Department of Mathematics, Imperial College London, London SW7 2AZ (United Kingdom)
Publication Date:
OSTI Identifier:
22403225
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physics of Fluids (1994); Journal Volume: 27; Journal Issue: 4; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICAL METHODS AND COMPUTING; AIR; CHEMISTRY; COMBUSTION; CYLINDRICAL CONFIGURATION; DISTURBANCES; FLAMES; FLOW RATE; METHANE; NAVIER-STOKES EQUATIONS; NONLINEAR PROBLEMS; NOZZLES; STEADY-STATE CONDITIONS; TOPOLOGY; UNSTEADY FLOW

Citation Formats

Blanchard, M., E-mail: mathieu.blanchard@ladhyx.polytechnique.fr, Schuller, T., Centrale-Supélec, Grande Voie des Vignes, 92290 Châtenay-Malabry, Sipp, D., and Schmid, P. J. Response analysis of a laminar premixed M-flame to flow perturbations using a linearized compressible Navier-Stokes solver. United States: N. p., 2015. Web. doi:10.1063/1.4918672.
Blanchard, M., E-mail: mathieu.blanchard@ladhyx.polytechnique.fr, Schuller, T., Centrale-Supélec, Grande Voie des Vignes, 92290 Châtenay-Malabry, Sipp, D., & Schmid, P. J. Response analysis of a laminar premixed M-flame to flow perturbations using a linearized compressible Navier-Stokes solver. United States. doi:10.1063/1.4918672.
Blanchard, M., E-mail: mathieu.blanchard@ladhyx.polytechnique.fr, Schuller, T., Centrale-Supélec, Grande Voie des Vignes, 92290 Châtenay-Malabry, Sipp, D., and Schmid, P. J. Wed . "Response analysis of a laminar premixed M-flame to flow perturbations using a linearized compressible Navier-Stokes solver". United States. doi:10.1063/1.4918672.
@article{osti_22403225,
title = {Response analysis of a laminar premixed M-flame to flow perturbations using a linearized compressible Navier-Stokes solver},
author = {Blanchard, M., E-mail: mathieu.blanchard@ladhyx.polytechnique.fr and Schuller, T. and Centrale-Supélec, Grande Voie des Vignes, 92290 Châtenay-Malabry and Sipp, D. and Schmid, P. J.},
abstractNote = {The response of a laminar premixed methane-air flame subjected to flow perturbations around a steady state is examined experimentally and using a linearized compressible Navier-Stokes solver with a one-step chemistry mechanism to describe combustion. The unperturbed flame takes an M-shape stabilized both by a central bluff body and by the external rim of a cylindrical nozzle. This base flow is computed by a nonlinear direct simulation of the steady reacting flow, and the flame topology is shown to qualitatively correspond to experiments conducted under comparable conditions. The flame is then subjected to acoustic disturbances produced at different locations in the numerical domain, and its response is examined using the linearized solver. This linear numerical model then allows the componentwise investigation of the effects of flow disturbances on unsteady combustion and the feedback from the flame on the unsteady flow field. It is shown that a wrinkled reaction layer produces hydrodynamic disturbances in the fresh reactant flow field that superimpose on the acoustic field. This phenomenon, observed in several experiments, is fully interpreted here. The additional perturbations convected by the mean flow stem from the feedback of the perturbed flame sheet dynamics onto the flow field by a mechanism similar to that of a perturbed vortex sheet. The different regimes where this mechanism prevails are investigated by examining the phase and group velocities of flow disturbances along an axis oriented along the main direction of the flow in the fresh reactant flow field. It is shown that this mechanism dominates the low-frequency response of the wrinkled shape taken by the flame and, in particular, that it fully determines the dynamics of the flame tip from where the bulk of noise is radiated.},
doi = {10.1063/1.4918672},
journal = {Physics of Fluids (1994)},
number = 4,
volume = 27,
place = {United States},
year = {Wed Apr 15 00:00:00 EDT 2015},
month = {Wed Apr 15 00:00:00 EDT 2015}
}
  • Combustion instability is a resonance phenomenon that arises due to the coupling between the system acoustics and the unsteady heat release. The constructive feedback between the two processes, which is known to occur as a certain phase relationship between the pressure and the unsteady heat release rate is satisfied, depends on many parameters among which is the acoustic mode, the flame holder characteristics, and the dominant burning pattern. In this paper, the authors construct an analytical model to describe the dynamic response of a laminar premixed flame stabilized on the rim of a tube to velocity oscillation. They consider uniformmore » and nonuniform velocity perturbations superimposed on a pipe flow velocity profile. The model results show that the magnitude of heat release perturbation and its phase with respect to the dynamic perturbation dependent primarily on the flame Strohal number, representing the ratio of the dominant frequency times the tube radius to the laminar burning velocity. In terms of this number, high-frequency perturbations pass through the flame while low frequencies lead to a strong response. The phase with respect to the velocity perturbation behaves in the opposite way. Results of this model are shown to agree with experimental observations and to be useful in determining how the combustion excited model is selected among all the acoustic unstable modes. The model is then used to obtain a time-domain differential equation describing the relationship between the velocity perturbation and the heat release response over the entire frequency range.« less
  • Two versions of a three-dimensional multistage Navier-Stokes code were used to optimize the design of an eleven-stage high-pressure compressor. The first version of the code utilized a mixing plane approach to compute the flow through multistage machines. The effects due to tip clearances and flowpath cavities were not modeled. This code was used to minimize the regions of separation on airfoil and endwall surfaces for the compressor. The resulting compressor contained bowed stators and rotor airfoils with contoured endwalls. Experimental data acquired for the HPC showed that it achieved 2% higher efficiency than a baseline machine, but it had 14%more » lower stall margin. Increased stall margin of the HPC was achieved by modifying the stator airfoils without compromising the gain in efficiency as demonstrated in subsequent rig and engine tests. The modifications to the stators were defined by using the second version of the multistage Navier-Stokes code, which models the effects of tip clearance and endwall flowpath cavities, as well as the effects of adjacent airfoil rows through the use of bodyforces and deterministic stresses. The application of the Navier-Stokes code was assessed to yield up to 50% reduction in the compressor development time and cost.« less
  • The foundation for the development of modern compressible flow solver is based on the Riemann solution of the inviscid Euler equations. The high-order schemes are basically related to high-order spatial interpolation or reconstruction. In order to overcome the low-order wave interaction mechanism due to the Riemann solution, the temporal accuracy of the scheme can be improved through the Runge-Kutta method, where the dynamic deficiencies in the first-order Riemann solution is alleviated through the sub-step spatial reconstruction in the Runge-Kutta process. The close coupling between the spatial and temporal evolution in the original nonlinear governing equations seems weakened due to itsmore » spatial and temporal decoupling. Many recently developed high-order methods require a Navier-Stokes flux function under piece-wise discontinuous high-order initial reconstruction. However, the piece-wise discontinuous initial data and the hyperbolic-parabolic nature of the Navier-Stokes equations seem inconsistent mathematically, such as the divergence of the viscous and heat conducting terms due to initial discontinuity. In this paper, based on the Boltzmann equation, we are going to present a time-dependent flux function from a high-order discontinuous reconstruction. The theoretical basis for such an approach is due to the fact that the Boltzmann equation has no specific requirement on the smoothness of the initial data and the kinetic equation has the mechanism to construct a dissipative wave structure starting from an initially discontinuous flow condition on a time scale being larger than the particle collision time. The current high-order flux evaluation method is an extension of the second-order gas-kinetic BGK scheme for the Navier-Stokes equations (BGK-NS). The novelty for the easy extension from a second-order to a higher order is due to the simple particle transport and collision mechanism on the microscopic level. This paper will present a hierarchy to construct such a high-order method. The necessity to couple spatial and temporal evolution nonlinearly in the flux evaluation can be clearly observed through the numerical performance of the scheme for the viscous flow computations.« less
  • Convergence is proved and error bounds are derived for a finite-difference approximation to discontinuous solutions of the Navier-Stokes equations for nonisentropic, compressible flow in one space dimension. The scheme is fully implicit and can be implemented under reasonable mesh conditions. It is shown that the approximations converge at the rate O([Delta]x[sup 1/2]) when the initial data is in H[sup 1], and O([Delta]x[sup a]) (a < 1/12) when the initial velocity and energy are in L[sup 2], and the initial density is piece wise H[sup 1]. The errors are measured in a norm that dominates the sup-norm of the error inmore » the density, which in general is discontinuous. This choice is indicated by the known continuous-dependence theory and accounts for the low rate of convergence in the second case.« less
  • A multistage compressor performance analysis method based on the three-dimensional Reynolds-averaged Navier-Stokes equations is presented in this paper. This method is an average passage approach where deterministic stresses are used to ensure continuous physical properties across interface planes. The average unsteady effects due to neighboring blades and/or vanes are approximated using deterministic stresses along with the application of bodyforces. Bodyforces are used to account for the potential interaction between closely coupled (staged) rows. Deterministic stresses account for the average wake blockage and mixing effects both axially and radially. The attempt here is to implement an approximate technique for incorporating periodicmore » unsteady flow physics that provides for a robust multistage design procedure incorporating reasonable computational efficiency. The present paper gives the theoretical development of the stress/bodyforce models incorporated in the code, and demonstrates the usefulness of these models in practical compressor applications. Compressor performance prediction capability is then established through a rigorous code/model validation effort using the power of networked workstations. The numerical results are compared with experimental data in terms of one-dimensional performance parameters such as total pressure ratio and circumferentially averaged radial profiles deemed critical to compressor design. This methodology allows the designer to design from hub to tip with a high level of confidence in the procedure.« less