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

Active Management of Flap-Edge Trailing Vortices

The vortex hazard produced by large airliners and increasingly larger airliners entering service, combined with projected rapid increases in the demand for air transportation, is expected to act as a major impediment to increased air traffic capacity. Significant reduction in the vortex hazard is possible, however, by employing active vortex alleviation techniques that reduce the wake severity by dynamically modifying its vortex characteristics, providing that the techniques do not degrade performance or compromise safety and ride quality. With this as background, a series of experiments were performed, initially at NASA Langley Research Center and subsequently at the Berlin University of Technology in collaboration with the German Aerospace Center. The investigations demonstrated the basic mechanism for managing trailing vortices using retrofitted devices that are decoupled from conventional control surfaces. The basic premise for managing vortices advanced here is rooted in the erstwhile forgotten hypothesis of Albert Betz, as extended and verified ingeniously by Coleman duPont Donaldson and his collaborators. Using these devices, vortices may be perturbed at arbitrarily long wavelengths down to wavelengths less than a typical airliner wingspan and the oscillatory loads on the wings, and hence the vehicle, are small. Significant flexibility in the specific device has been demonstrated using local passive and active separation control as well as local circulation control via Gurney flaps. The method is now in a position to be tested in a wind tunnel with a longer test section on a scaled airliner configuration. Alternatively, the method can be tested directly in a towing tank, on a model aircraft, a light aircraft or a full-scale airliner. The authors believed that this method will have significant appeal from an industry perspective due to its retrofit potential with little to no impact on cruise (devices tucked away in the cove or retracted); low operating power requirements; small lift oscillations when deployed in a time-dependent manner; and significant flexibility with respect to the specific devices selected

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

Wind Turbine Tip Vortices under the influence of Wind Tunnel Blockage Effects

The current paper describes the characteristics of the tip vortex in the near wake of a three-bladed upwind horizontal axis wind turbine with a rotor diameter of 3 m. Phase-locked stereo particle image velocimetry measurements were carried out under the influence of the wind tunnel walls that create a high blockage ratio. The location of the vortex, convection velocity, core radius, and strength were investigated and compared with similar investigations, including different blockages cases. Additionally, the same performance of the wind turbine model was simulated in the open source wind turbine tool QBlade, using the lifting line free vortex wake module in the absence of the walls. The results showed that the location of the tip vortices was more inboard the tip and more downstream the tunnel compared to the simulations and similar experiments. The convection velocity remained similar in the axial direction and changed in the lateral direction, contributing to the delay of the movement of the tip vortex outboard the tip. The strength, based on the circulation, was found with a difference of 4% between simulation and experiment

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