Elimination of Numerical Damping in the Stability Analysis of Non-Compact Thermoacoustic Systems With Linearized Euler Equations

2021 ◽  
Author(s):  
Thomas Hofmeister ◽  
Tobias Hummel ◽  
Frederik Magnus Berger ◽  
Noah Klarmann ◽  
Thomas Sattelmayer
Author(s):  
Thomas Hofmeister ◽  
Tobias Hummel ◽  
Frederik Berger ◽  
Noah Klarmann ◽  
Thomas Sattelmayer

Abstract The hybrid Computational Fluid Dynamics/Computational AeroAcoustics (CFD/CAA) approach represents an effective method to assess the stability of non-compact thermoacoustic systems. This paper summarizes the state-of-the-art of this method, which is currently applied for the stability prediction of a lab-scale configuration of a perfectly-premixed, swirl-stabilized gas turbine combustion chamber at the Thermodynamics institute of the Technical University of Munich. Specifically, 80 operational points, for which experimentally observed stability information is readily available, are numerically investigated concerning their susceptibility to develop thermoacoustically unstable oscillations at the first transversal eigenmode of the combustor. Three contributions are considered in this work: (1) flame driving due the deformation and displacement of the flame, (2) visco-thermal losses in the acoustic boundary layer and (3) damping due to acoustically induced vortex shedding. The analysis is based on eigenfrequency computations of the Linearized Euler Equations with the stabilized Finite Element Method (sFEM). One main advancement presented in this study is the elimination of the non-physical impact of artificial diffusion schemes, which is necessary to produce numerically stable solutions, but falsifies the computed stability results.


2021 ◽  
Vol 143 (3) ◽  
Author(s):  
Thomas Hofmeister ◽  
Tobias Hummel ◽  
Frederik Berger ◽  
Noah Klarmann ◽  
Thomas Sattelmayer

Abstract The hybrid computational fluid dynamics/computational aeroacoustics (CFD/CAA) approach represents an effective method to assess the stability of noncompact thermoacoustic systems. This paper summarizes the state-of-the-art of this method, which is currently applied for the stability prediction of a lab-scale configuration of a perfectly premixed, swirl-stabilized gas turbine combustion chamber at the Thermodynamics institute of the Technical University of Munich. Specifically, 80 operational points, for which experimentally observed stability information is readily available, are numerically investigated concerning their susceptibility to develop thermoacoustically unstable oscillations at the first transversal eigenmode of the combustor. Three contributions are considered in this work: (1) flame driving due the deformation and displacement of the flame, (2) visco-thermal losses in the acoustic boundary layer and (3) damping due to acoustically induced vortex shedding. The analysis is based on eigenfrequency computations of the Linearized Euler Equations with the stabilized finite element method (sFEM). One main advancement presented in this study is the elimination of the nonphysical impact of artificial diffusion schemes, which is necessary to produce numerically stable solutions, but falsifies the computed stability results.


2014 ◽  
Vol 2014 ◽  
pp. 1-10 ◽  
Author(s):  
Stefan Balint ◽  
Agneta M. Balint

This paper considers the stability of constant solutions to the 1D Euler equation. The idea is to investigate the effect of different function spaces on the well-posedness and stability of the null solution of the 1D linearized Euler equations. It is shown that the mathematical tools and results depend on the meaning of the concepts “perturbation,” “small perturbation,” “solution of the propagation problem,” and “small solution, that is, solution close to zero,” which are specific for each function space.


2016 ◽  
Vol 138 (10) ◽  
Author(s):  
Santosh Konangi ◽  
Nikhil K. Palakurthi ◽  
Urmila Ghia

The goal of this paper is to derive the von Neumann stability conditions for the pressure-based solution scheme, semi-implicit method for pressure-linked equations (SIMPLE). The SIMPLE scheme lies at the heart of a class of computational fluid dynamics (CFD) algorithms built into several commercial and open-source CFD software packages. To the best of the authors' knowledge, no readily usable stability guidelines appear to be available for this popularly employed scheme. The Euler equations are examined, as the inclusion of viscosity in the Navier–Stokes (NS) equation serves to only soften the stability limits. First, the one-dimensional (1D) Euler equations are studied, and their stability properties are delineated. Next, a rigorous stability analysis is carried out for the two-dimensional (2D) Euler equations; the analysis of the 2D equations is considerably more challenging as compared to analysis of the 1D form of equations. The Euler equations are discretized using finite differences on a staggered grid, which is used to achieve equivalence to finite-volume discretization. Error amplification matrices are determined from the stability analysis, stable and unstable regimes are identified, and practical stability limits are predicted in terms of the maximum allowable Courant–Friedrichs–Lewy (CFL) number as a function of Mach number. The predictions are verified using the Riemann problem, and very good agreement is obtained between the analytically predicted and the “experimentally” observed CFL values. The successfully tested stability limits are presented in graphical form, as compared to complicated mathematical expressions often reported in published literature. Since our analysis accounts for the solution scheme along with the full system of flow equations, the conditions reported in this paper offer practical value over the conditions that arise from analysis of simplified 1D model equations.


Author(s):  
Moritz Schulze ◽  
Tobias Hummel ◽  
Noah Klarmann ◽  
Frederik Berger ◽  
Bruno Schuermans ◽  
...  

A novel methodology for linear stability analysis of high-frequency thermoacoustic oscillations in gas turbine combustors is presented. The methodology is based on the linearized Euler equations, which yield a high-fidelity description of acoustic wave propagation and damping in complex, non-uniform, reactive mean flow environments, such as encountered in gas turbine combustion chambers. Specifically, this work introduces three novelties to the community: (1) Linear stability analysis on the basis of linearized Euler equations. (2) Explicit consideration of three-dimensional, acoustic oscillations at screech level frequencies, particularly the first transversal mode. (3) Handling of non-compact flame coupling with LEE, that is, the spatially varying coupling dynamics between perturbation and unsteady flame response due to small acoustic wavelengths. Two different configurations of an experimental model combustor in terms of thermal power and mass flow rates are subject of the analysis. Linear flame driving is modeled by prescribing the unsteady heat release source term of the linearized Euler equations by local flame transfer functions, which are retrieved from first principles. The required steady state flow field is numerically obtained via CFD, which is based on an extended Flamelet-Generated Manifold combustion model, taking into account heat transfer to the environment. The model is therefore highly suitable for such types of combustors. The configurations are simulated, and thermoacoustically characterized in terms of eigenfrequencies and growth rates associated with the first transversal mode. The findings are validated against experimentally observed thermoacoustic stability characteristics. On the basis of the results, new insights into the acoustic field are discussed.


Author(s):  
Moritz Schulze ◽  
Tobias Hummel ◽  
Noah Klarmann ◽  
Frederik Berger ◽  
Bruno Schuermans ◽  
...  

A novel methodology for linear stability analysis of high-frequency thermoacoustic oscillations in gas turbine combustors is presented. The methodology is based on the linearized Euler equations (LEEs), which yield a high-fidelity description of acoustic wave propagation and damping in complex, nonuniform, reactive mean flow environments, such as encountered in gas turbine combustion chambers. Specifically, this work introduces three novelties to the community: (1) linear stability analysis on the basis of linearized Euler equations. (2) Explicit consideration of three-dimensional, acoustic oscillations at screech level frequencies, particularly the first-transversal mode. (3) Handling of noncompact flame coupling with LEE, that is, the spatially varying coupling dynamics between perturbation and unsteady flame response due to small acoustic wavelengths. Two different configurations of an experimental model combustor in terms of thermal power and mass flow rates are subject of the analysis. Linear flame driving is modeled by prescribing the unsteady heat release source term of the linearized Euler equations by local flame transfer functions, which are retrieved from first principles. The required steady-state flow field is numerically obtained via computational fluid dynamics (CFD), which is based on an extended flamelet-generated manifold (FGM) combustion model, taking into account heat transfer to the environment. The model is therefore highly suitable for such types of combustors. The configurations are simulated, and thermoacoustically characterized in terms of eigenfrequencies and growth rates associated with the first-transversal mode. The findings are validated against experimentally observed thermoacoustic stability characteristics. On the basis of the results, new insights into the acoustic field are discussed.


2019 ◽  
Vol 1 (1) ◽  
pp. 49-60
Author(s):  
Simon Heru Prassetyo ◽  
Ganda Marihot Simangunsong ◽  
Ridho Kresna Wattimena ◽  
Made Astawa Rai ◽  
Irwandy Arif ◽  
...  

This paper focuses on the stability analysis of the Nanjung Water Diversion Twin Tunnels using convergence measurement. The Nanjung Tunnel is horseshoe-shaped in cross-section, 10.2 m x 9.2 m in dimension, and 230 m in length. The location of the tunnel is in Curug Jompong, Margaasih Subdistrict, Bandung. Convergence monitoring was done for 144 days between February 18 and July 11, 2019. The results of the convergence measurement were recorded and plotted into the curves of convergence vs. day and convergence vs. distance from tunnel face. From these plots, the continuity of the convergence and the convergence rate in the tunnel roof and wall were then analyzed. The convergence rates from each tunnel were also compared to empirical values to determine the level of tunnel stability. In general, the trend of convergence rate shows that the Nanjung Tunnel is stable without any indication of instability. Although there was a spike in the convergence rate at several STA in the measured span, that spike was not replicated by the convergence rate in the other measured spans and it was not continuous. The stability of the Nanjung Tunnel is also confirmed from the critical strain analysis, in which most of the STA measured have strain magnitudes located below the critical strain line and are less than 1%.


Sign in / Sign up

Export Citation Format

Share Document