A two-way coupling for modeling thermoacoustic instabilities in a flat flame Rijke tube

Thermoacoustic instabilities are a serious problem for lean premixed combustion systems. Due to different time and length scales associated with the flow field, combustion, and acoustics, numerical computations of thermoacoustic phenomena are conceptually challenging. This work presents a coupled method for the simulation of thermoacoustic instabilities in low Mach number reacting flows. The acoustics are represented by a reduced order model that can be obtained from network techniques or finite element computations. A detailed chemistry finite-difference zero Mach number solver is used for the small scale flame dynamics. Under the assumption that the pressure is continuous across the flame, the acoustic model can be reduced to a time-domain relation mapping the velocity perturbation downstream of the flame to that upstream. Closure is obtained by the flame code, which delivers the jump in velocity across the combustion zone. The method is applied to an experimental laminar premixed burner-stabilized flat flame Rijke tube, that exhibits strong thermoacoustic oscillations associated with the 5k=4 mode of the geometrical set-up. In addition to the fundamental oscillation, a significant subharmonic response of the flame is observed. Results from the coupled simulation are compared to the experimental data. Good qualitative and quantitative agreement is found

Unsteady Thick Airfoil Aerodynamics: Experiments, Computation, and Theory

An experimental, computational and theoretical investigation was carried out to study the aerodynamic loads acting on a relatively thick NACA 0018 airfoil when subjected to pitching and surging, individually and synchronously. Both pre-stall and post-stall angles of attack were considered. Experiments were carried out in a dedicated unsteady wind tunnel, with large surge amplitudes, and airfoil loads were estimated by means of unsteady surface mounted pressure measurements. Theoretical predictions were based on Theodorsen's and Isaacs' results as well as on the relatively recent generalizations of van der Wall. Both two- and three-dimensional computations were performed on structured grids employing unsteady Reynolds-averaged Navier-Stokes (URANS). For pure surging at pre-stall angles of attack, the correspondence between experiments and theory was satisfactory; this served as a validation of Isaacs theory. Discrepancies were traced to dynamic trailing-edge separation, even at low angles of attack. Excellent correspondence was found between experiments and theory for airfoil pitching as well as combined pitching and surging; the latter appears to be the first clear validation of van der Wall's theoretical results. Although qualitatively similar to experiment at low angles of attack, two-dimensional URANS computations yielded notable errors in the unsteady load effects of pitching, surging and their synchronous combination. The main reason is believed to be that the URANS equations do not resolve wake vorticity (explicitly modeled in the theory) or the resulting rolled-up un- steady flow structures because high values of eddy viscosity tend to \smear" the wake. At post-stall angles, three-dimensional computations illustrated the importance of modeling the tunnel side walls

Is the Blade Element Momentum theory overestimating wind turbine loads? – An aeroelastic comparison between OpenFAST's AeroDyn and QBlade's Lifting-Line Free Vortex Wake method

Abstract. Load calculations play a key role in determining the design loads of different wind turbine components. To obtain the aerodynamic loads for these calculations, the industry relies heavily on the Blade Element Momentum (BEM) theory. BEM methods use several engineering correction models to capture the aerodynamic phenomena present in Design Load Cases (DLCs) with turbulent wind. Because of this, BEM methods can overestimate aerodynamic loads under challenging conditions when compared to higher-order aerodynamic methods – such as the Lifting-Line Free Vortex Wake (LLFVW) method – leading to unnecessarily high design loads and component costs. In this paper, we give a quantitative answer to the question of load overestimation of a particular BEM implementation by comparing the results of aeroelastic load calculations done with the BEM-based OpenFAST code and the QBlade code, which uses a particular implementation of the LLFVW method. We compare extreme and fatigue load predictions from both codes using sixty-six 10 min load simulations of the Danish Technical University (DTU) 10 MW Reference Wind Turbine according to the IEC 61400-1 power production DLC group. Results from both codes show differences in fatigue and extreme load estimations for the considered sensors of the turbine. LLFVW simulations predict 9 % lower lifetime damage equivalent loads (DELs) for the out-of-plane blade root and the tower base fore–aft bending moments compared to BEM simulations. The results also show that lifetime DELs for the yaw-bearing tilt and yaw moments are 3 % and 4 % lower when calculated with the LLFVW code. An ultimate state analysis shows that extreme loads of the blade root out-of-plane bending moment predicted by the LLFVW simulations are 3 % lower than the moments predicted by BEM simulations. For the maximum tower base fore–aft bending moment, the LLFVW simulations predict an increase of 2 %. Further analysis reveals that there are two main contributors to these load differences. The first is the different way both codes treat the effect of the nonuniform wind field on the local blade aerodynamics. The second is the higher average aerodynamic torque in the LLFVW simulations. It influences the transition between operating modes of the controller and changes the aeroelastic behavior of the turbine, thus affecting the loads

Numerical and experimental evaluation of shock dividers

Mitigation of pressure pulsations in the exhaust of a pulse detonation combustor is crucial for operation with a downstream turbine. For this purpose, a device termed the shock divider is designed and investigated. The intention of the divider is to split the leading shock wave into two weaker waves that propagate along separated ducts with different cross sections, allowing the shock waves to travel with different velocities along different paths. The separated shock waves redistribute the energy of the incident shock wave. The shock dynamics inside the divider are investigated using numerical simulations. A second-order dimensional split finite volume MUSCL-scheme is used to solve the compressible Euler equations. Furthermore, low-cost simulations are performed using geometrical shock dynamics to predict the shock wave propagation inside the divider. The numerical simulations are compared to high-speed schlieren images and time-resolved total pressure recording. For the latter, a high-frequency pressure probe is placed at the divider outlet, which is shown to resolve the transient total pressure during the shock passage. Moreover, the separation of the shock waves is investigated and found to grow as the divider duct width ratio increases. The numerical and experimental results allow for a better understanding of the dynamic evolution of the flow inside the divider and inform its capability to reduce the pressure pulsations at the exhaust of the pulse detonation combustor

Determination of the angle of attack on a research wind turbine rotor blade using surface pressure measurements

In this paper, a method to determine the angle of attack on a wind turbine rotor blade using a chordwise pressure distribution measurement was applied. The approach used a reduced number of pressure tap data located close to the blade leading edge. The results were compared with the measurements from three external probes mounted on the blade at different radial positions and with analytical calculations. Both experimental approaches used in this study are based on the 2-D flow assumption; the pressure tap method is an application of the thin airfoil theory, while the probe method applies geometrical and induction corrections to the measurement data. The experiments were conducted in the wind tunnel at the Hermann Föttinger Institut of the Technische Universität Berlin. The research turbine is a three-bladed upwind horizontal axis wind turbine model with a rotor diameter of 3 m. The measurements were carried out at rated conditions with a tip speed ratio of 4.35 and different yaw and pitch angles were tested in order to compare the approaches over a wide range of conditions. Results show that the pressure tap method is suitable and provides a similar angle of attack to the external probe measurements as well as the analytical calculations. This is a significant step for the experimental determination of the local angle of attack, as it eliminates the need for external probes, which affect the flow over the blade and require additional calibration

About the suitability of different numerical methods to reproduce model wind turbine measurements in a wind tunnel with a high blockage ratio

In the present paper, numerical and experimental investigations of a model wind turbine with a diameter of 3.0 m are described. The study has three objectives. The first one is the provision of validation data. The second one is to estimate the influence of the wind tunnel walls by comparing measurements to simulated results with and without wind tunnel walls. The last objective is the comparison and evaluation of methods of high fidelity, namely computational fluid dynamics, and medium fidelity, namely lifting-line free vortex wake. The experiments were carried out in the large wind tunnel of the TU Berlin where a blockage ratio of 40 % occurs. With the lifting-line free vortex wake code QBlade, the turbine was simulated under far field conditions at the TU Berlin. Unsteady Reynolds-averaged Navier–Stokes simulations of the wind turbine, including wind tunnel walls and under far field conditions, were performed at the University of Stuttgart with the computational fluid dynamics code FLOWer. Comparisons among the experiment, the lifting-line free vortex wake code and the computational fluid dynamics code include on-blade velocity and angle of attack. Comparisons of flow fields are drawn between the experiment and the computational fluid dynamics code. Bending moments are compared among the simulations. A good accordance was achieved for the on-blade velocity and the angle of attack, whereas deviations occur for the flow fields and the bending moments

High-fidelity simulations of the mixing and combustion of a technically premixed hydrogen flame

Numerical simulations are used here to obtain further understanding on the flashback mechanism of a technically premixed hydrogen flame operated in lean conditions. Recent work from the authors (Mira et al., 2020) showed that the hydrogen momentum strongly influences the flame dynamics and plays a fundamental role on the stability limits of the combustor. The axial injection influences the vortex breakdown position and therefore, the propensity of the burner to produce flashback. This work is an extension of our previous work where we include a detailed description of the mixing process and the influence of equivalence ratio fluctuations and heat loss on the flame dynamics