Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic stall”

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.

At least 307 records · Page 17

Numerical Experiment with Time and Spatial Accuracy of Navier-Stokes Computation For Helicopter Problems

Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other blades. This interaction, known as blade-vortex interactions (BVI), results in unsteady loading of the blades and can cause a distinctive acoustic signature. Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interactions can help rotor designers and engineers to predict rotor performance and to develop designs for low acoustic signature. Such a predictive method must preserve a blade's shed vortex for several blade revolutions before being dissipated. A number of researchers have explored the requirements for this task. This paper will outline some new capabilities that have been added to the NASA Ames' OVERFLOW code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these improvements, a number of case studies will be presented. These case studies consist of free convection of a 2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a full 3-D unsteady viscous solution of a helicopter rotor in forward flight In this study both central and upwind difference schemes are modified to be more accurate. Central difference scheme is chosen for this simulation because the flowfield is not dominated by strong shocks. The feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interaction. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in fully-implicit manner at each time-step. The spatial accuracy is either second and fourth-order central difference or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and demonstrate the methods for several sample cases and for a helicopter rotor. Preliminary computations on a rotor were performed by using this method and are in the process of documentation.

Ahmad, Jasim↗

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗

An Automated Stall-Speed Warning System

The NASA Dryden Flight Research Facility embarked upon a project with the United States Army Aviation Engineering Flight Activity (USAAEFA) to develop and test a stall-speed warning system. NASA designed and built an automated stall-speed warning system which presents both airspeed and stall speed to the pilot. The airspeed and stall speed are computed in real time by monitoring the basic aerodynamic parameters (dynamic pressure, horizontal and vertical accelerations, and pressure altitude) and other parameters (elevator and flap positions, engine torques, and fuel flow). In addition, an aural warning at predetermined stall margins is presented to the pilot through a voice synthesizer. Once the system was designed and installed in the aircraft, a flight-test program of less than 20 hr was anticipated to determine the stall-speed software coefficients. These coefficients would then be inserted in the system's software and then test flown over a period of about 10 hr for the purposes of evaluation.

Dynamic Pressure↗

Wind-tunnel free-flight investigation of a model of a spin-resistant fighter configuration

An investigation was conducted to provide some insight into the features affecting the high-angle-of-attack characteristics of a high-performance twin-engine fighter airplane which in operation has exhibited excellent stall characteristics with a general resistance to spinning. Various techniques employed in the study included wind-tunnel free-flight tests, flow-visualization tests, static force tests, and dynamic (forced-oscillation) tests. In addition to tests conducted on the basic configuration tests were made with the wing planform and the fuselage nose modified. The results of the study showed that the model exhibited good dynamic stability characteristics at angles of attack well beyond that for wing stall. The directional stability of the model was provided by the vertical tail at low and moderate angles of attack and by the fuselage forebody at high angles of attack. The wing planform was found to have little effect on the stability characteristics at high angles of attack. The tests also showed that although the fuselage forebody produced beneficial contributions to static directional stability at high angles of attack, it also produced unstable values of damping in yaw. Nose strakes located in a position which eliminated the beneficial nose contributions produced a severe directional divergence.

Grafton, S. B.↗

An Automated Stall-Speed Warning System

The NASA Dryden Flight Research Facility embarked upon a project with the United States Army Aviation Engineering Flight Activity (USAAEFA) to develop and test a stall-speed warning system. NASA designed and built an automated stall-speed warning system which presents both airspeed and stall speed to the pilot. The airspeed and stall speed are computed in real time by monitoring the basic aerodynamic parameters {dynamic pressure, horizontal and vertical accelerations, and pressure altitude) and other parameters (elevator and flap positions, engine torques, and fuel flow). In addition, an aural warning at predetermined stall margins is presented to the pilot through a voice synthesizer. Once the system was designed and installed in the aircraft, a flight-test program of less than 20 hr was anticipated to determine the stall-speed software coefficients. These coefficients would then be inserted in the system's software and then test flown over a period of about 10 hr for the purposes of evaluation.

Pilot↗

Fully integrated aerodynamic/dynamic optimization of helicopter rotor blades

A fully integrated aerodynamic/dynamic optimization procedure is described for helicopter rotor blades. The procedure combines performance and dynamic analyses with a general purpose optimizer. The procedure minimizes a linear combination of power required (in hover, forward flight, and maneuver) and vibratory hub shear. The design variables include pretwist, taper initiation, taper ratio, root chord, blade stiffnesses, tuning masses, and tuning mass locations. Aerodynamic constraints consist of limits on power required in hover, forward flight and maneuvers; airfoil section stall; drag divergence Mach number; minimum tip chord; and trim. Dynamic constraints are on frequencies, minimum autorotational inertia, and maximum blade weight. The procedure is demonstrated for two cases. In the first case, the objective function involves power required (in hover, forward flight and maneuver) and dynamics. The second case involves only hover power and dynamics. The designs from the integrated procedure are compared with designs from a sequential optimization approach in which the blade is first optimized for performance and then for dynamics. In both cases, the integrated approach is superior.

Walsh, Joanne L.↗

Fully integrated aerodynamic/dynamic optimization of helicopter rotor blades

This paper describes a fully integrated aerodynamic/dynamic optimization procedure for helicopter rotor blades. The procedure combines performance and dynamics analyses with a general purpose optimizer. The procedure minimizes a linear combination of power required (in hover, forward flight, and maneuver) and vibratory hub shear. The design variables include pretwist, taper initiation, taper ratio, root chord, blade stiffnesses, tuning masses, and tuning mass locations. Aerodynamic constraints consist of limits on power required in hover, forward flight and maneuver; airfoil section stall; drag divergence Mach number; minimum tip chord; and trim. Dynamic constraints are on frequencies, minimum autorotational inertia, and maximum blade weight. The procedure is demonstrated for two cases. In the first case the objective function involves power required (in hover, forward flight, and maneuver) and dynamics. The second case involves only hover power and dynamics. The designs from the integrated procedure are compared with designs from a sequential optimization approach in which the blade is first optimized for performance and then for dynamics. In both cases, the integrated approach is superior.

Walsh, Joanne L.↗

Surge and Stall Detection Using Acoustic Analysis FOR GAS Turbine Hybrid Cycles

In this article compressor dynamics were studied in a gas turbine - fuel cell hybrid power system having a larger compressor volume than traditionally found in gas turbine systems. This larger compressor volume adversely affects the surge margin of the gas turbine. Industrial acoustic sensors were placed near the compressor to identify when the equipment was getting close to the surge line. Fast Fourier transform (FFT) mathematical analysis was used to obtain spectra representing the probability density across the frequency range (0 - 5000 Hz). Comparison between FFT spectra for nominal and transient operations revealed that higher amplitude spikes were observed during incipient stall at three different frequencies, 900, 1020, and 1800 Hz. These frequencies were compared to the natural frequencies of the equipment and the frequency for the rotating turbomachinery to identify more general nature of the acoustic signal typical of the onset of compressor surge. The primary goal of this acoustic analysis was to establish an online methodology to monitor compressor stability that can be anticipated and avoided.

42 ENGINEERING↗

Computational Study of a Generic T-tail Transport

This paper presents a computational study on the static and dynamic stability characteristics of a generic transport T-tail configuration under a NASA research program to improve stall models for civil transports. The NASA Tetrahedral Unstructured Software System (TetrUSS) was used to obtain both static and periodic dynamic solutions at low speed conditions for three Reynolds number conditions up to 60 deg angle of attack. The computational results are compared to experimental data. The dominant effects of Reynolds number for the static conditions were found to occur in the stall region. The pitch and roll damping coefficients compared well to experimental results up to up to 40 deg angle of attack whereas yaw damping coefficient agreed only up to 20 deg angle of attack.

McMillin, S. Naomi↗

Flutter Prediction Report in Support of the High Angle Working Group at the Third Aeroelastic Prediction Workshop

This paper presents a summary of the computational flutter results associated with the AIAA third Aeroelastic Prediction Workshop (AePW-3) High Angle Working Group (HAWG). The computational results are compared against the experimental data collected during the Pitch and Plunge Apparatus (PAPA) Benchmark Supercritical Wing (BSCW) test campaign conducted in the Transonic Dynamics Tunnel (TDT) at NASA Langley Research Center thirty years ago. During that test several flutter points were identified at transonic conditions. One of these points, specifically near Mach 0.8 and 5° angle of attack, became a focal point of the computational challenge with in HAWG. Various fidelity time-domain, Reduced Order Model (ROM) and Linearized Frequency Domain (LFD) methods were used by seven participating teams. While there are encouraging trends in the computational results, the range of the predicted flutter dynamic pressure is still quite large due to the stall flutter mechanism. A description of each participating team’s software and methods is included.

Pawel Chwalowski↗

Dynamic surface measurements on a model helicopter rotor during blade slap at high angles of attack

The modern helicopter offers a unique operational capability to both the public and private sectors. However, the use of the helicopter may become severely limited due to the radiated noise generated by the rotor system. A description is presented of some of the experimental results obtained with a model helicopter rotor in an anechoic wind tunnel with regard to blade stall as a source mechanism of blade slap. Attention is given to dynamic rotor blade surface phenomena and the resulting far field impulsive noise from the model helicopter rotor at high angles of attack and low tip speed. The results of the investigation strongly implicates the boundary layer as playing an important role in blade slap due to blade/vortex interaction (BVI) in a highly loaded rotor. Intermittent stall cannot be ruled out as a possible source mechanism for blade slap. This implies that blade surface characteristics, airfoil shape and local Reynolds number may now be used as tools to reduce the resultant far-field sound pressure levels in helicopters.

Hubbard, J. E., Jr.↗

Supersonic stall flutter of high-speed fans

An analytical model is proposed for predicting the onset of supersonic stall bending flutter in high-speed rotors. The analysis is based on a modified two-dimensional, compressible, unsteady actuator disk theory. The stability boundary predicted by the analysis is shown to be in good agreement with the measured boundary of a high speed fan. The prediction that the flutter mode would be a forward traveling wave sensitive to wheel speed and aerodynamic loading is confirmed by experimental measurements. In addition, the analysis shows that reduced frequency and dynamic head also play a significant role in establishing the supersonic stall bending flutter boundary of an unshrouded fan.

Adamczyk, J. J.↗

Mitigation of Motor Stalling and FIDVR via Energy Storage Systems with Temporal Logic Specifications

The fault-induced delayed voltage recovery (FIDVR) phenomenon has been very common from the distribution system through the transmission system. It causes a delay on recovering significantly depressed local voltage after the fault is cleared, and it can also lead to more widespread cascading system failures. Mitigating this event with current control approaches is challenging and becoming a crucial issue. Here, a model predictive control-based strategy employing signal temporal logic specifications is proposed to help mitigate FIDVR. To this end, it investigates and extends a dynamic performance model allowing analytic insights into the system-wide impact of motor stalling and FIDVR. The proposed controller provides richer descriptions of voltage specifications addressing both magnitude and time simultaneously. We consider different control specifications with reactive power support from energy storage systems to prevent the voltage during/after the fault from dropping too low, and reduce the delay time of voltage recovery. The simulation results conducted with the IEEE 57 bus test network validate the proposed method and demonstrate the effectiveness of the mitigation strategy on FIDVR.

25 ENERGY STORAGE↗

Local correlation-based transition models for high-Reynolds-number wind-turbine airfoils

Abstract. Modern wind-turbine airfoil design requires robust performance predictions for varying thicknesses, shapes, and appropriate Reynolds numbers. The airfoils of current large offshore wind turbines operate with chord-based Reynolds numbers in the range of 3–15 million. Turbulence transition in the airfoil boundary layer is known to play an important role in the aerodynamics of these airfoils near the design operating point. While the lack of prediction of lift stall through Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) is well known, airfoil design using CFD requires the accurate prediction of the glide ratio (L/D) in the linear portion of the lift polar. The prediction of the drag bucket and the glide ratio is greatly affected by the choice of the transition model in RANS CFD of airfoils. We present the performance of two existing local correlation-based transition models – one-equation model (γ− SA) and two-equation model (γ-Reθt‾- SA) coupled with the Spalart–Allmaras (SA) RANS turbulence model – for offshore wind-turbine airfoils operating at a high Reynolds number. We compare the predictions of the two transition models with available experimental and CFD data in the literature in the Reynolds number range of 3–15 million including the AVATAR project measurements of the DU00-W-212 airfoil. Both transition models predict a larger L/D compared to fully turbulent results at all Reynolds numbers. The two models exhibit similar behavior at Reynolds numbers around 3 million. However, at higher Reynolds numbers, the one-equation model fails to predict the natural transition behavior due to early transition onset. The two-equation transition model predicts the aerodynamic coefficients for airfoils of various thickness at higher Reynolds numbers up to 15 million more accurately compared to the one-equation model. As a result, the two-equation model predictions are more comparable to the predictions from eN transition model. However, a limitation of this model is observed at very high Reynolds numbers of around 12–15 million where the predictions are very sensitive to the inflow turbulent intensity. The combination of the two-equation transition model coupled with the Spalart–Allmaras (SA) RANS turbulence model is a good method for performance prediction of modern wind-turbine airfoils using CFD.

17 WIND ENERGY↗

Thermally Stable Ohmic Contacts on Silicon Carbide Developed for High- Temperature Sensors and Electronics

The NASA aerospace program, in particular, requires breakthrough instrumentation inside the combustion chambers of engines for the purpose of, among other things, improving computational fluid dynamics code validation and active engine behavioral control (combustion, flow, stall, and noise). This environment can be as high as 600 degrees Celsius, which is beyond the capability of silicon and gallium arsenide devices. Silicon-carbide- (SiC-) based devices appear to be the most technologically mature among wide-bandgap semiconductors with the proven capability to function at temperatures above 500 degrees Celsius. However, the contact metalization of SiC degrades severely beyond this temperature because of factors such as the interdiffusion between layers, oxidation of the contact, and compositional and microstructural changes at the metal/semiconductor interface. These mechanisms have been proven to be device killers. Very costly and weight-adding packaging schemes that include vacuum sealing are sometimes adopted as a solution.

Okojie, Robert S.↗

A Wind-Tunnel Investigation of the Stall-Flutter Characteristics of a Supersonic-Type Propeller at Positive and Negative Thrust

An investigation has been made utilizing a three-blade, 10-foot- diameter, supersonic-ty-pe propeller to determine propeller flutter characteristics. The particular flutter characteristics of interest were (1) the effect of stall flutter on a propeller operating in positive and negative thrust, (2) the effect of stall flutter on a propeller operating with the thrust axis inclined, and (3) the variation of vibratory blade shear stresses as the stall flutter boundary is penetrated and exceeded. Thrust and power measurements were made for all test conditions. Wake and inflow surveys were made when appropriate, to define the thrust and torque distributions and the magnitude of the inflow velocity. Stress measurements were made simultaneously to obtain the propeller flutter and bending response. It was found when operating both in the positive and negative thrust regions that, for most cases after the onset of flutter, the magnitude of the flutter stresses at first increased rapidly with section blade angle P, after which further increases in 0 resulted in only a moderate increase or a reduction in stress. Thrust-axis inclination up to the limit of the tests (angle of attack of 15 deg and dynamic pressure of 40 psf) appeared to have no effect on stall flutter. The stall flutter stresses were found to be directly associated with the section thrust characteristics of the blades. The onset of flutter was found to occur simultaneously with the divergence of the section thrust variation with blade angle from linearity for stations outboard of the blade 0.8-radius station. The maximum flutter stresses appeared to be a function of the maximum section thrust obtained at or in the vicinity of the blade 0.8-radius station. In an attempt to correlate two-dimensional airfoil data with three-dimensional data to predict the stall angle of attack (divergence of the section thrust) of the blade sections, it was found that no consistent correlation could be obtained. Also, a knowledge of the inflow conditions appeared to be insufficient to account for differences in airfoil characteristics between the two-dimensional and the three-dimensional cases.

Rogallo, Vernon L.↗

Symbolic generation of elastic rotor blade equations using a FORTRAN processor and numerical study on dynamic inflow effects on the stability of helicopter rotors

The process of performing an automated stability analysis for an elastic-bladed helicopter rotor is discussed. A symbolic manipulation program, written in FORTRAN, is used to aid in the derivation of the governing equations of motion for the rotor. The blades undergo coupled bending and torsional deformations. Two-dimensional quasi-steady aerodynamics below stall are used. Although reversed flow effects are neglected, unsteady effects, modeled as dynamic inflow are included. Using a Lagrangian approach, the governing equations are derived in generalized coordinates using the symbolic program. The program generates the steady and perturbed equations and writes into subroutines to be called by numerical routines. The symbolic program can operate on both expressions and matrices. For the case of hovering flight, the blade and dynamic inflow equations are converted to equations in a multiblade coordinate system by rearranging the coefficients of the equations. For the case of forward flight, the multiblade equations are obtained through the symbolic program. The final multiblade equations are capable of accommodating any number of elastic blade modes. The computer implementation of this procedure consists of three stages: (1) the symbolic derivation of equations; (2) the coding of the equations into subroutines; and (3) the numerical study after identifying mass, damping, and stiffness coefficients. Damping results are presented in hover and in forward flight with and without dynamic inflow effects for various rotor blade models, including rigid blade lag-flap, elastic flap-lag, flap-lag-torsion, and quasi-static torsion. Results from dynamic inflow effects which are obtained from a lift deficiency function for a quasi-static inflow model in hover are also presented.

Reddy, T. S. R.↗