Search NASA⌕ Search

SEARCH · Search NASA

Results for “deep 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.

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↗

High-Fidelity Simulation Aerodynamics Dataset of NACA 0012 and 0021 Airfoils

The repository contains time series data of pressure, viscous and moment forces for the NACA 0012 and 0021 airfoils at angles of attack (AOA), 5, 17, 30, 45, 60 and 90 degrees. The data was generated during the high-fidelity simulations of the airfoils using the open-source CFD code, Nalu-Wind (https://github.com/Exawind/nalu-wind). The research objective was to study the effects of mesh resolution and turbulence model on the three-dimensional deep stall aerodynamics of the airfoils, which has been published in the Journal of Turbulence (https://doi.org/10.1080/14685248.2023.2225141). Several meshes of the O-grid type were generated using the modified Pointwise Glyph script of Carrigan (https://github.com/pointwise/AirfoilMesh) by varying the wall normal and spanwise resolutions. Note that, the geometry was extruded by four chord lengths in the spanwise direction. In the paper, we showed that Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS/LES turbulence model with a minimum of 24 cells per chord length in the spanwise direction is necessary to correctly predict the loads in the deep stall regime. Under the NACA 0012 airfoil directory, the data for several combinations of wall-normal and spanwise resolutions are presented in each of the AOA subdirectories. The corresponding experimental data for the lift and drag polars can be found in the "exp_data" subdirectory. The directory for the NACA 0021 airfoil consists of CFD force data and airfoil polars for Reynolds numbers 2.7e5 and 2.0e6.

17 WIND ENERGY↗

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗

Three-Dimensional Aerodynamics and Vortex-Shedding Characteristics of Wind Turbine Airfoils over 360-Degree Angles of Attack

In this work, we present the first three-dimensional (3D) computational investigation of wind turbine airfoils over 360° angles of attack to predict unsteady aerodynamic loads and vortex-shedding characteristics. To this end, static–airfoil simulations are performed for the FFA-W3 airfoil family at a Reynolds number of 107 with the Improved Delayed Detached Eddy Simulation turbulence model. Aerodynamic forces reveal that the onset of boundary-layer instabilities and flow separation does not necessarily coincide with the onset of stall. In addition, a comparison with two-dimensional simulation data and flat plate theory extension of airfoil polars, suggest that, in the deep stall regime, 3D effects remain critical for predicting both the unsteady loads and the vortex-shedding dynamics. For all airfoils, the vortex-shedding frequencies are found to be inversely proportional to the wake width. In the case of slender airfoils, the frequencies are nearly independent of the airfoil thickness, and their corresponding Strouhal number St is approximately 0.15. Based on the calculated St, the potential for shedding frequencies to coincide with the natural frequencies of the International Energy Agency 15 MW reference wind turbine blades is investigated. The analysis shows that vortex-induced vibrations occur primarily at angles of attack of around ±90° for all airfoils.

17 WIND ENERGY↗

A DNN surrogate unsteady aerodynamic model for wind turbine loads calculations

A recurrent deep-neural network (DNN) surrogate model capable of modeling the unsteady aerodynamic response and dynamic stall behavior of wind turbine blades has been developed and validated for use in engineering design codes. The model is trained using a subset of the oscillating airfoil experiments conducted at the Ohio State University wind tunnel. The predictions from our DNN model show excellent agreement with the measured data and, in all cases, a marked improvement over the state-of-the-art unsteady aerodynamic models. The DNN-based unsteady aerodynamics model was integrated with OpenFAST to perform full-turbine load computations for the NREL-5MW rotor. The largest differences are observed for the inboard stations, particularly in the pitching moment response, when using the new surrogate model compared to the other models available in OpenFAST.

17 WIND ENERGY↗

Formation of a Large Cold Groundwater Mantle Helium Anomaly and High Temperature Geothermal Resources in Response to Bimodal Magmatism Near Roosevelt Hot Springs and Utah FORGE, Milford Valley, Southwest Utah

Abstract A large mantle helium anomaly and separate domains of high heat flow are the predominant manifestations of bimodal magmatic activity in the Milford valley. The mantle helium anomaly (1.9–2.6 R/Ra) covers 270 km 2 and is subdivided into two separated domains: a cold shallow groundwater regime and high temperature hydrothermal activity. The zone of anomalous heat flow covers >100 km 2 and is also subdivided into two adjacent domains, comprising hydrothermal activity at Roosevelt Hot Springs (RHS) (3–7 W/m 2 ) and conductive heat flow (100–180 mW/m 2 ). While the transfer of heat and mantle helium at RHS are coupled, heat and helium transfer are decoupled in the adjacent cold groundwater regime to the west. Both the mantle helium and geothermal anomalies are attributed to recent mafic‐felsic magmatic intrusions of >400 km 3 , however, the absence of volcanic eruptions <500,000 years indicates magmas stall before rising to shallow crustal level <10 km depth. Deep level magmatism produces a felsic composition melt, which is inferred to be responsible for the widespread and near uniform range of diluted mantle helium values. A thick and impermeable mass of crystalline granitic basement rock at the mid‐crustal level divides the ascent of mantle helium into separate flow paths. It may also impede the rise of buoyant magma trapping thermal energy that facilitates partial melting, slow cooling, and development of a thick thermal aureole. Partitioning of convective and conductive thermal regimes and independent flow paths supplying deeply derived helium characterize the development of a large long‐lived magma‐related geothermal system.

Geochemistry & Geophysics↗

Nonenergy Biomass Carbon Removal and Storage (BiCRS): Assessing Durability of Nongaseous Carbon Products Across Terrestrial Storage Fates

Biomass Carbon Removal and Storage, or BiCRS, pathways use plants or algae that remove carbon dioxide from the atmosphere through photosynthesis and store it underground or in long-lived products. While some BiCRS approaches generate an energy product, all BiCRS approaches generate a carbon product. A new subset of BiCRS approaches focus on the storage of these raw or converted carbon products for generation of carbon credits. However, the durability of these approaches is highly variable as carbon products vary widely in their “form” and the conditions of their “fate.” We organize our thinking about carbon products and their durability around these two primary axes. The durability of carbon product “forms” is mediated by chemical recalcitrance and ranges substantially across agricultural residues, municipal solid waste, woody biomass, and nongaseous products of thermochemical conversion (e.g., biochars and bio-oils). Meanwhile, terrestrial storage “fates” vary in the mechanism employed to stall decay, including surface storage, dry storage, shallow anoxic storage, and deep or geologic anoxic storage (or injection). Each mechanism has different implications for suitability with different feedstock forms as well as long-term risks. We present a framework for assessing durability of solid or liquid raw and conversion carbon products under terrestrial storage fates, highlighting knowns, unknowns, and research priorities moving forward.

09 BIOMASS FUELS↗

Intracycle RPM control for vertical axis wind turbines

The wind energy market is currently dominated by horizontal axis wind turbines (HAWTs); however, vertical axis wind turbines (VAWTs) are emerging as a design alternative, especially for deep-water offshore siting due to their low center of gravity, ease of access to drivetrain components, and overall simplicity. Due to the absence of a pitch mechanism in large-scale Darrieus VAWTs, stall control has often been used to manage power and loads. Introducing a pitching mechanism in H-type VAWTs has been studied, but this diminishes the mechanical simplicity advantage, and the use of a pitching mechanism in a large-scale Darrieus-type VAWT is not practical. This work examines an innovative, alternative method to control the rotor dynamics of a large-scale 5 MW VAWT to maximize power while constraining loads without introducing any new or complex mechanical elements. This control strategy is termed intracycle revolution per minute (RPM) control, where the rotational speed of the turbine is allowed to vary in an optimal fashion with the azimuthal location of blades as opposed to typical constant RPM operation. An optimization framework is formulated for an open-loop optimal control problem and solved to maximize power subject to constraints on aerodynamic design loads. Results are presented to demonstrate the benefits and the performance limits of intracycle RPM control for large-scale 5 MW Darrieus VAWTs, namely, (1) power production (quantified in terms of AEP) that can be increased subject to baseline load limits and (2) opportunities to significantly increase AEP or decrease loads via intracycle RPM control that are examined for both two-bladed and three-bladed VAWTs.

17 WIND ENERGY↗

Trans-crustal structural control of CO 2 -rich extensional magmatic systems revealed at Mount Erebus Antarctica

Erebus volcano, Antarctica, with its persistent phonolite lava lake, is a classic example of an evolved, CO 2 -rich rift volcano. Seismic studies provide limited images of the magmatic system. Here we show using magnetotelluric data that a steep, melt-related conduit of low electrical resistivity originating in the upper mantle undergoes pronounced lateral re-orientation in the deep crust before reaching shallower magmatic storage and the summit lava lake. The lateral turn represents a structural fault-valve controlling episodic flow of magma and CO 2 vapour, which replenish and heat the high level phonolite differentiation zone. This magmatic valve lies within an inferred, east-west structural trend forming part of an accommodation zone across the southern termination of the Terror Rift, providing a dilatant magma pathway. Unlike H 2 O-rich subduction arc volcanoes, CO 2 -dominated Erebus geophysically shows continuous magmatic structure to shallow crustal depths of < 1 km, as the melt does not experience decompression-related volatile supersaturation and viscous stalling.

58 GEOSCIENCES↗

Submersible Multistage Centrifugal Pump for Versatile Test Reactor Cartridge Test Loop

Submersible multistage centrifugal pumps are ideal for pumping in narrow confined spaces and achieving necessary head pressures and flow rates. Once a diameter is determined then manipulation of the number of stages and motor speed are all that are required to meet desired flow conditions. The Versatile Test Reactor (VTR) closed loop cartridge systems will need forced convection cooling independent of the main reactor. A multistage centrifugal pump can meet the necessary flow rates and pumping pressures while minimizing space taken. The pump considered for this work was based off a deep well submersible pump, a variation of a multistage centrifugal pump. We experimented with two pump sizes, 5 cm (2 inch) and 7.5 cm (3 inch) diameters. These diameters were chosen to fit into the inner diameter of standard 5 and 7.5 (2 and 3 inch) Schedule 40 pipe, respectively. This made the design for the test loop both simpler and less expensive as the need for an engineered pump housing was eliminated. Initial test cartridge planning indicated space for only a 5 cm (2 inch) diameter pump, though early testing of this size showed the need for an abnormally high-speed and high-power motor. Fine tuning of the cartridge design allowed a pump size increase to 7.5 cm (3 inches), which was the pump size most extensively tested in this work. The test loop is composed of various sizes of PVC and aluminum piping components in a loop configuration. The pump is driven by a Pittman 250 W (1/3 horsepower) electric motor with maximum speed of 3,450 RPM. Testing consisted of running the pump at a constant motor speed while varying a control valve to restrict flow through the loop, with differential pressure and flow rate recorded. This was done for one and two stage configurations for the 5 cm (2 inch) diameter impeller design and one, two, and three stage configurations for the 7.5 cm (3 inch diameter) impeller design, respectively. Due to pumping power requirements, two and three stage 7.5 cm (3 inch) diameter impeller testing at higher flowrates lowered the motor speed substantially. In regions where motor speed could not be maintained constant, the data were discarded. The test loop was also reconfigured to allow for the pump to be tested for pressure drop in a stalled or inoperable (0 RPM) flow condition. Demonstration of adequate natural convection cooling of the test cartridge fuel type is necessary under accident conditions, and this will depend upon the flow resistance through the impeller assembly when the pump is not operating. Thus, accurate knowledge of the effective impeller assembly loss coefficient is important for safety evaluations. The test loop was modified to provide water inlet and outlets on either side of the pump impeller stack, and a metered flow of lab water was provided in order to measure the pressure drop across the cartridges as a function of flowrate. Data from the pump head curve testing developed as part of this work and supported by analysis using pump head affinity laws indicates that a three stage 7.5 cm (3 inch) pump impeller design will meet target requirements for coolant flow within the VTR cartridge sodium cartridge at full power conditions [1] with margin; this corresponds to a flowrate of 45 l/min (12 gpm) at a pressure drop of 6.1 m (20 feet) of water head. The results of the pressure loss measurements across the impeller assembly when the pump is stationary (i.e., at 0 RPM) indicate that the pressure loss coefficient is 0.921 for a two impeller stack configuration; this value is calculated based on the flow velocity through the minimum available flow area within a single stage of the impeller which corresponds to 1.4 cm2.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗