Analysis of the onset of dynamic stall
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The Boeing-Vertol VR-7 airfoil was experimentally studied with steady and pulsed upper-surface blowing for sinusoidal pitching oscillations described by alpha = alpha(sub m) + 10 deg sin(omega t). The tests were conducted in the U.S. Army Aeroflightdynamics Directorate's Water Tunnel at NASA Ames Research Center. The experiment was performed at a Reynolds number of 100,000. Pitch oscillations with alpha(sub m) = 10 deg and 15 deg and with reduced frequencies ranging from k = 0.005 to 0.15 were examined. Blowing conditions ranged from C(sub mu) = 0.03 to 0.66 and F(+) = 0 to 3. Unsteady lift, drag, and pitching-moment loads were measured, and fluorescent-dye flow visualizations were obtained. Steady, upper-surface blowing was found to be capable of trapping a separation bubble near the leading edge during a portion of the airfoil's upward rotation. When this occurred, the lift was increased significantly and stall was averted. In all cases, steady blowing reduced the hysteresis amplitudes present in the loads and produced a large thrust force. The benefits of steady blowing diminished as the reduced frequency and mean angle of oscillation increased. Pulsed blowing showed only marginal benefits for the conditions tested. The greatest gains from pulsed blowing were achieved at F(+) = 0.9.
The first major area of this study was to develop a vorticity-velocity formulation and numerical solution algorithms suitable for the analyses of incompressible as well as low-to- moderate-speed compressible flows. Research performed towards contributing to the determination of the appropriate vorticity and dilation creation boundary conditions suggested to temporarily set aside this approach and use a primitive-variable approach other than the pseudo-compressibility approach used. The second major area of study was initiated to comprehensively examine the INS-2D and INS-3D programs from the point of view of boundary conditions. The research carried out was documented in the form of two technical papers which are included in Appendices A and B; the boundary-condition related issues for INS-3D are briefly mentioned.
The importance of transition modeling in the computation of compressible, unsteady separated flows is discussed. The study showed that it is critical to predict the experimentally attained transition point properly in order to obtain good agreement with data it the same Mach number and Reynolds number.
An experiment was performed in a water tunnel on a Boeing-Vertol VR-7 airfoil to study the effects of tangential blowing over the upper surface. Blowing was applied at the quarter-chord location during sinusoidal pitching oscillations described by alpha = alpha(sub m) + 10 deg sin omega t. Results were obtained for a Reynolds number of 1 x 10(exp 5), mean angles of 10 and 15 deg, reduced frequencies ranging from 0.005 to 0.15, and blowing rates from C(sub mu) = 0.16 to 0.66. Unsteady lift, drag, and pitching moment loads are reported, along with fluorescent-dye flow visualizations. Strong steady blowing was found to prevent the bursting of the leading-edge separation bubble at several test points. When this occurred, the lift was increased significantly, stall was averted, and the shape of the moment response showed a positive damping in pitch. In almost all cases, steady blowing reduced the hysteresis amplitudes present in the loads, but the benefits diminished as the reduced frequency and mean angle of oscillation increased. A limited number of pulsed blowing cases indicated that for low blowing rates, the greatest gains were achieved at F(sup +) = 0.9.
Experimental and computational investigations of the dynamic stall phenomenon continue to attract the attention of various research groups in the major aeronautical research laboratories. There are two reasons for this continued research interest. First, the occurrence of dynamic stall on the retreating blade of helicopters imposes a severe performance limitation and thus suggests to search for ways to delay the onset of dynamic stall. Second, the lift enhancement prior to dynamic stall presents an opportunity to achieve enhanced maneuverability of fighter aircraft. A description of the major parameters affecting dynamic stall and lift and an evaluation of research efforts prior to 1988 has been given by Carr. In this paper the authors' recent progress in the development of experimental and computational methods to analyze the dynamic stall phenomena occurring on NACA 0112 airfoils is reviewed. First, the major experimental and computational approaches and results are summarized. This is followed by an assessment of our results and an outlook toward the future.
Abstract Horizontal axis wind turbines (HAWTs) experience three‐dimensional rotational and unsteady aerodynamic phenomena at the rotor blades sections. These highly unsteady three‐dimensional effects have a dramatic impact on the aerodynamic load distributions on the blades, in particular, when they occur at high angles of attack due to stall delay and dynamic stall. Unfortunately, there is no complete understanding of the flow physics yet at these unsteady 3D flow conditions, and hence, the existing published theoretical models are often incapable of modelling the impact on the turbine response realistically. The purpose of this paper is to provide an insight on the combined influence of the stall delay and dynamic stall on the blade load history of wind turbines in controlled and uncontrolled conditions. New dynamic stall vortex and nonlinear tangential force coefficient modules, which integrally take into account the three dimensional rotational effect, are also proposed in this paper. This module along with the unsteady influence of turbulent wind speed and tower shadow is implemented in a blade element momentum (BEM) model to estimate the aerodynamic loads on a rotating blade more accurately. This work presents an important step to help modelling the combined influence of the stall delay and dynamic stall on the load history of the rotating wind turbine blades which is vital to have lighter turbine blades and improved wind turbine design systems.
We evaluate two leading-edge-based dynamic stall-onset criteria (namely, the maximum magnitudes of the leading-edge suction parameter and the boundary enstrophy flux) for mixed and trailing-edge stall. These criteria have been shown to successfully predict the onset of leading-edge stall at Reynolds numbers of [Formula: see text], where the leading-edge suction drops abruptly. However, for mixed/trailing-edge stall, leading-edge suction tends to persist even when there is significant trailing-edge reversed flow and stall is underway, necessitating further investigation into the effectiveness of these criteria. Using wall-resolved large-eddy simulations and the unsteady Reynolds-averaged Navier–Stokes method, we simulate one leading-edge stall and three mixed/trailing-edge stall cases at Reynolds numbers of [Formula: see text] and [Formula: see text]. We contrast the progression of flow features such as trailing-edge separation and vortex formation across different stall types and evaluate the stall-onset criteria relative to critical points in the flow. We find that the criteria nearly coincide with the instance of leading-edge suction collapse and are reached in advance of dynamic stall vortex formation and lift stall for all four cases. We conclude that the two criteria effectively signal dynamic stall onset in cases where the dynamic stall vortex plays a prominent role.
Dynamic and static stall data in relation to airfoil stall at subsonic speeds
The predictions of regressive lag-mode damping levels are correlated with the database of an isolated, soft-inplane, three-blade rotor operated untrimmed. The database was generated at the Army Aeroflightdynamics Directorate at Ames. The correlation covers a broad range of data, from near-zero thrust conditions in hover to high-thrust and highly stalled conditions in forward flight with advance ratio as high as 0.55 and shaft angle as high as 20 degrees. In the experimental rotor, the airfoil or blade portion has essentially uniform mass and stiffness distributions, but the root flexure has highly nonuniform mass and stiffness distributions. Accordingly, the structural approximations refer to four models of root-flexure-blade assembly. They range from a rigid flap-lag model to three elastic flap-lag-torsion models, which differ in modeling the root flexure. The three models of root-flexure are: three root springs in which the bending-torsion couplings are fully accounted for; a finite-length beam element with some average mass and stiffness distributions such that the fundamental frequencies match those of the experimental model; and accurate modal representation in which the actual mass and stiffness distributions of the experimental root-flexure-blade assembly are used in calculating the nonrotating mode shapes. The four models of root-flexure-blade assembly are referred to as the rigid flap-lag model, spring model, modified model and modal model. For each of these four models of the root-flexure-blade assembly, the predictions are based on the following five aerodynamic theories: ear theory, which accounts for large angle-of-attack and reverse-flow effects on lift, and has constant drag and pitching moment; quasisteady stall theory, which includes quasisteady stall lift, drag and pitching moment characteristics of the airfoil section, dynamics stall theory, which uses the ONERA dynamic stall models of lift, drag and pitching moment; dynamic wake theory, which is based on a finite-state three-dimensional wake model and includes all wake effects including both shed and trailing vorticity; and dynamics and wake theory, which combines both dynamic stall theory and dynamic wake theory and is a relatively complete aerodynamic representation.
Phenomena that control the flow during the stall portion of a dynamic stall cycle are analyzed, and their effect on blade motion is outlined. Four mechanisms by which dynamic stall may be initiated are identified: (1) bursting of the separation bubble, (2) flow reversal in the turbulent boundary layer on the airfoil upper surface, (3) shock wave-boundary layer interaction behind the airfoil crest, and (4) acoustic wave propagation below the airfoil. The fluid mechanics that contribute to the identified flow phenomena are summarized, and the usefulness of a model that incorporates the required fluid mechanics mechanisms is discussed.
The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.
A comparison is made of the results of three methods for calculating the effects of dynamic stall on the performance, airloads, and blade stresses of a helicopter rotor at high loading. The three dynamic-stall methods considered predict essentially the same performance and trim for the rotor. They give roughly the same mean bending moments, but the peak-to-peak torsion and bending moments differ by 25 to 40 percent, and there are differences in the details of the predicted blade motion and stresses. The latter are due to significant differences in the dynamic stall aerodynamic loads, particularly the aerodynamic pitch moment, predicted by the three methods.
Unsteady aerodynamics will be an important part of the floating wind turbines of the future operating under high shear across the rotor disk coupled with platform motion and atmospheric turbulence. We develop an unsteady aerodynamics and dynamic stall model using a long short-term memory variant of recurrent neural networks. The neural network model is trained using the oscillating airfoil data set from Ohio State University. The predictions from our machine learning (ML)-based model show good agreement with the experimental data and other state-of-the-art dynamic stall models for a wide range of airfoils, Reynolds numbers and reduced frequencies. In some cases the predictions are better than the Beddoes-Leishman model implementation in OpenFAST, when using the default coefficients. The ML-based model is also able to capture the key physics associated with dynamic stall, such as the precedence of moment stall before lift stall and cycle-to-cycle variations in the aerodynamic response. The new unsteady aerodynamics model is expected to improve prediction of fatigue loads for yaw-based wake-steering control scenarios in actuator-line and actuator-disc simulations of wind farms. Our methodology for training the ML-model provides a pathway for improving design level tools using high-fidelity computational fluid dynamics (CFD) simulations in the future.
This work presents two aerodynamic models that are used to perform modal and stability analyses of wind turbines. Including aerodynamic effects in stability analyses is important to capture aerodynamic inertia, stiffness, and damping. The two models consist of a dynamic inflow model based on Øye’s model and an unsteady aerodynamics model based on a simplified Beddoes-Leishman dynamic stall model. Both models are written in a continuous-time state-space form, which is the form needed to perform linearizations within OpenFAST. In this article, we present the new models, perform some simple verifications against HAWC2 results, and then perform stability analyses. We investigate the contribution of the aerodynamic effects on the turbine modes, particularly the flapwise modes, which are the most affected. Because of the development presented in this work, Campbell diagrams accounting for dynamic stall and dynamic wake effects were generated for the first time using OpenFAST.
Time-dependent Navier-Stokes simulations have been carried out for a flexible UH-60A rotor in forward flight, where the rotor wake interacts with the rotor blades. These flow conditions involved blade vortex interaction and dynamic stall, two common conditions that occur as modern helicopter designs strive to achieve greater flight speeds and payload capacity. These numerical simulations utilized high-order spatial accuracy and delayed detached eddy simulation. Emphasis was placed on understanding how improved rotor wake resolution affects the prediction of the normal force, pitching moment, and chord force of the rotor. Adaptive mesh refinement was used to highly resolve the turbulent rotor wake in a computationally efficient manner. Moreover, blade vortex interaction was found to trigger dynamic stall. Time-dependent flow visualization was utilized to provide an improved understanding of the numerical and physical mechanisms involved with three-dimensional dynamic stall.
Rotor behavior in stalled conditions is investigated using wind tunnel test data of a l/l0-scale CH-47B/C type rotor, which provides a set of test conditions extending from unstalled to light stall to some deep stall conditions over a wide range of advance ratios. The rotor performance measured in the wind tunnel is similar to the main rotor performance measured during the NASA/Army UH-60A Airloads Program, although the two rotors are quite different. The analysis CAMRAD II has been used to predict the rotor performance and loads. Full-scale airfoil test data are corrected for Reynolds number effects for comparison with the model-scale rotor test. The calculated power coefficient shows good correlation with the measurements below stall with the Reynolds number-corrected airfoil table. Various dynamic stall models are used in the calculations. The Boeing model shows the lift augmentation at low advance ratios and the Leishman-Beddoes model shows better correlation of torsion moment than the other models at mu = 0.2. However, the dynamic stall models, in general, show only a small influence on the rotor power and torsion moment predictions especially at higher advance ratios.