Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational aerodynamics”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Code Artifact for: Clustering Analysis of Commercial Vehicles Using Automatically Extracted Features from Time Series Data [SWR-21-96]

This repository contains data ingestion, feature extraction, and analysis code used in NREL Technical report "Clustering Analysis of Commercial Vehicles Using Automatically Extracted Features from Time Series Data." The code is written in Python. The ETL and feature extraction code must be run in a Spark context. The analysis code can be run without Spark, provided you have pre-computed features in a CSV file. Analysis code related to the NREL Technical Report NREL/TP-2C00-74212. Includes PySpark functions to perform trip segmentation and feature extraction over big time series data in Apache Spark. Includes "domain specific" features such as Aerodynamic Speed (ft/s), Characteristic Acceleration (ft/s2), Percent Below 55 (%), Percent Zero (%), Stops Per Mile, Average Speed (mph), Maximum Speed (mph), and Speed Standard Deviation (mph). Includes Pyspark UDF to compute "domain agnostic" features using the TSFresh library. This software record also includes the analysis notebooks and code to generate the results in the previously mentioned technical report.

Perr-Sauer, Jordan↗
DOE OSTI2025

Wind Turbine Rotor Design Using High-Fidelity Aerostructural Optimization

Large wind turbines yield more energy but demand careful aeroelastic blade design. Coupled multiphysics design strategies can reduce wind energy costs by exploiting fluid-structure interactions. This work presents the first high-fidelity aerostructural optimization study of a large wind turbine rotor. We use blade-resolved fluid dynamics and structural solvers in a monolithic gradient-based optimization framework to explore steady-state torque and blade mass tradeoffs. The coupled-adjoint approach computes gradients efficiently, enabling the optimization of over 100 structural and geometric parameters simultaneously. Our optimization study modifies a DTU 10 MW benchmark with a simplified structure and isotropic material properties. The tightly coupled optimizations increase torque by 14% while reducing rotor mass by 9% or reduce blade mass by 27% while maintaining torque. Blade-resolved models provide greater design freedom, enabling 5% higher mass reductions than conventional parameterizations at equal torque. This framework paves the way for more detailed high-fidelity optimization studies to complement conventional design approaches.

17 WIND ENERGY↗
DOE OSTI2025

Unsteady aerodynamic loads on pitching aerofoils represented by Gaussian body force distributions

The actuator line model (ALM) is an approach commonly used to represent lifting and dragging devices like wings and blades in large-eddy simulations (LES). The crux of the ALM is the projection of the actuator point forces onto the LES grid by means of a Gaussian regularisation kernel. The minimum width of the kernel is constrained by the grid size; however, for most practical applications like LES of wind turbines, this value is an order of magnitude larger than the optimal value that maximises accuracy. This discrepancy motivated the development of corrections for the actuator line, which, however, neglect the effect of unsteady spanwise shed vorticity. In this work we develop a model for the impact of spanwise shed vorticity on the unsteady loading of an aerofoil modelled as a Gaussian body force distribution, where the model is applicable within the regime of unsteady attached flow. The model solution is derived both in the time and frequency domain and features an explicit dependence on the Gaussian kernel width. We verify the model with ALM-LES for both pitch steps and periodic pitching. The model solution is compared with Theodorsen theory and validated with both computational fluid dynamics using body fitted grids and experiment. It is concluded that the optimal kernel width for unsteady aerodynamics is approximately 40 % of the chord. The ALM is able to predict the magnitude of the unsteady loading up to a reduced frequency of 𝑘 ≈ 0.2.

17 WIND ENERGY↗
DOE OSTI2022

Passive Ultrasonic Deterrents to Reduce Bat Mortality in Wind Farms (Final Report)

The overall objective of this project was to develop active and passive ultrasonic deterrent concepts to mitigate bat mortality at wind turbines. The research was supported by the US Department of Energy under the contract DE-EE0008731. Novel ultrasonic bat deterrents are investigated in this study. The deterrents are based on the idea of aerodynamic whistles wherein flow-acoustic resonance is used to produce high-amplitude tonal sound at desired ultrasonic frequencies. The deterrents can be grouped into active and passive. The active deterrents/whistles are driven by pressurized air supplied by an external source (e.g., an air compressor). The passive deterrents/whistles are designed to be mounted on wind turbine blades and are “powered” by the energy in the air moving past the rotor blades. Multiple active and passive deterrent ideas are proposed and investigated numerically and experimentally. The computational analysis involves solving the compressible unsteady Reynolds-averaged Navier-Stokes (URANS) equations and coupling the near-field aerodynamic solution with an integral method based on the Ffowcs Williams-Hawkings acoustic analogy to predict ultrasound in the farfield. Experiments are performed in the anechoic chamber at ISU (for the active whistles) and in the Stability wind tunnel at Virginia Tech (for the passive whistles).

17 WIND ENERGY↗
DOE OSTI2020

Preliminary Introduction of a Free Vortex Wake Method Into OpenFAST

Increasing wind turbine size can have substantial economic benefits. However, increasing wind turbine turbine size requires the use of models that accurately capture the aerodynamic response to large deflections and nonstraight blade geometries. Additionally, the installation of offshore floating turbines requires turbine models that can handle large rotor motion. Free vortex wake (FVW) methods are capable of modeling this type of complex physics while remaining computationally tractable to perform the many simulations necessary for turbine design. In this work, a FVW model is added to the National Renewable Energy Laboratory engineering tool OpenFAST to allow for the aerodynamic modeling of highly flexible turbines along with the aero-hydro-servo-elastic response capabilities of OpenFAST. Code implementation is not complete and thus no results are presented. However, areas of potential modeling improvement are discussed, such as modifications to the circulation calculation and the need for different FVW input model specifications when high-thrust cases are being modeled.

17 WIND ENERGY↗
DOE OSTI2024

Improvements to the Blade Element Momentum Formulation of OpenFAST for Skewed Inflows

In this work, we modify the blade element momentum algorithm of OpenFAST to improve its predictions under large skewed inflow conditions. We use the well-known Glauert's skew correction and introduce continuous extension of the model for high-thrust conditions. We present the rationale behind Glauert's empirical model and discuss the different conventions possible for the axial induction factor. We verify the model against the higher-fidelity lifting-line vortex method and blade-resolved computational fluid dynamics, and we observe that the new implementation enhances the accuracy and reliability of OpenFAST's aerodynamic modeling capabilities in conditions involving large skew angles. For the parametric studies run using the different codes, we find that the power changes with the skew angle as cos 1.7 (θ skew ) and the thrust as cos 0.65 (θ skew ). An analysis of the azimuthal variation of the induced velocities in the rotor plane reveals that current redistribution models used in blade element momentum codes may need to be refined.

17 WIND ENERGY↗
DOE OSTI2022

Enabling Innovation in Wind Turbine Design Using Artificial Intelligence

The Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements (INTEGRATE) project is developing a new inverse-design capability for wind turbine rotors using invertible neural networks. This artificial intelligence (AI)-based technology can capture complex nonlinear aerodynamic effects 100 times faster than alternative design approaches.

aerodynamics↗
DOE OSTI2022

RAFT (Response Amplitudes of Floating Turbines) [SWR-22-40]

RAFT is a Python code for frequency-domain analysis of floating wind turbines. RAFT works in the frequency domain to provide efficient computations of a floating wind system's response spectra accounting for platform hydrodynamics, quasi-static mooring reactions, rotor aerodynamics, and linearized control. It can be used to calculate response amplitude operators (RAO's), power spectral densities, mean/static properties, and response metrics based on mean and root-mean-square values of each output. RAFT serves as the Level 1 model in the WEIS (Wind Energy with Integrated Servo-control) toolkit, and features an OpenMDAO wrapper.

Hall, Matthew↗
DOE OSTI2023

Comparison of free vortex wake and blade element momentum results against large-eddy simulation results for highly flexible turbines under challenging inflow conditions

Abstract. Throughout wind energy development, there has been a push to increase wind turbine size due to the substantial economic benefits. However, increasing turbine size presents several challenges, both physically and computationally. Modeling large, highly flexible wind turbines requires highly accurate models to capture the complicated aeroelastic response due to large deflections and nonstraight blade geometries. Additionally, the development of floating offshore wind turbines requires modeling techniques that can predict large rotor and tower motion. Free vortex wake methods model such complex physics while remaining computationally tractable to perform key simulations necessary during the turbine design process. Recently, a free vortex wake model – cOnvecting LAgrangian Filaments (OLAF) – was added to the National Renewable Energy Laboratory's engineering tool OpenFAST to allow for the aerodynamic modeling of highly flexible turbines along with the aero-hydro-servo-elastic response capabilities of OpenFAST. In this work, free vortex wake and low-fidelity blade element momentum (BEM) results are compared to high-fidelity actuator-line computational fluid dynamics simulation results via the Simulator fOr Wind Farm Applications (SOWFA) method for a highly flexible downwind turbine for varying yaw misalignment, shear exponent, and turbulence intensity conditions. Through these comparisons, it was found that for all considered quantities of interest, SOWFA, OLAF, and BEM results compare well for steady inflow conditions with no yaw misalignment. For OLAF results, this strong agreement with the SOWFA results was consistent for all yaw misalignment values. The BEM results, however, deviated significantly more from the SOWFA results with increasing absolute yaw misalignment. Differences between OLAF and BEM results were dominated by the yaw misalignment angle, with varying shear exponent and turbulence intensity leading to more subtle differences. Overall, OLAF results were more consistent than BEM results when compared to SOWFA results under challenging inflow conditions.

17 WIND ENERGY↗
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Encoding nonlinear and unsteady aerodynamics of limit cycle oscillations using nonlinear sparse Bayesian learning

This article investigates the applicability of a recently proposed, nonlinear sparse Bayesian learning (NSBL) algorithm to identify and estimate the complex aerodynamics of limit cycle oscillations. NSBL provides a semi-analytical framework for determining the data-optimal sparse model nested within a (potentially) over-parameterized model. This is particularly relevant to nonlinear dynamical systems where modelling approaches involve the use of physics-based and data-driven components. In such cases, the data-driven components, where analytical descriptions of the physical processes are not readily available, are often prone to overfitting, meaning that the empirical aspects of these models will often involve the calibration of an unnecessarily large number of parameters. While an overparameterized model may fit the observed data well, such models may be inadequate for making predictions in regimes that are different from those wherein the data were recorded. In view of this, it is desirable to not only calibrate the model parameters, but also identify the optimal compromise between data fit and model complexity. In this article, we exhibit the optimal model discovery for an aeroelastic system wherein the structural dynamics are well-known and described by a differential equation model, coupled with a semi-empirical aerodynamic model for laminar separation flutter, resulting in low-amplitude limit cycle oscillations (LCO). To illustrate the performance of the algorithm, in this article, we use synthetic data and demonstrate the ability of the algorithm to correctly rediscover the optimal model and model parameters, given a known data-generating model. The synthetic data are generated from a forward simulation of a known differential equation model with parameters selected so as to mimic the dynamics observed in wind-tunnel experiments. Subsequently, we demonstrate the performance of the algorithm for model selection using noisy LCO data from wind tunnel experiments. As there is no ground truth available for the experimental data case, we provide a comparison between NSBL and Bayesian model selection to validate the results, and demonstrate the use of NSBL as an efficient alternative to traditional methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enhancement of Unsteady and 3D Aerodynamics Models using Machine Learning

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.

17 WIND ENERGY↗
DOE OSTI2022

Enhancements to Overset Methods for Achieving Higher Accuracy and Solution Convergence: Preprint

Overset methods, when utilized within a multi-solver framework, provide a very flexible and powerful toolset for analyzing moving body problem. However, many overset framework suffer from several deficiencies such as (1) convergence degradation of the non-linear iterative schemes (2) difficulty in achieving higher-order accuracy and (3) high computational costs over long simulation times. In this work, we explore approaches that mitigate these deficiencies and provide a pathway for improved accuracy and fidelity of overset methods. A high-order, discontinuous Galerkin solver is integrated into an existing overset framework and applied to an unsteady ROBIN fuselage case. For convergence acceleration, an Overset-GMRES linear solver with a timestep controller is applied to a DLR-F6 aircraft steady simulation and an unsteady, notional rotor-fuselage. Finally the same rotor-fuselage case as well as the JVX rotor/wing case is simulated using a reduced-order aerodynamic model.

aerospace↗
DOE OSTI2022

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↗
DOE OSTI2022

Grassmannian Shape Representations for Aerodynamic Applications

Airfoil shape design is a classical problem in engineering, science, and manufacturing. Our motivation is to combine principled physics-based considerations for the shape design problem with modern computational techniques informed by a data-driven approach. Traditional analyses of airfoil shapes emphasize a flow-based sensitivity to deformations which can be represented generally by affine transformations (rotation, scaling, shearing, shifting). We present a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data, (ii) improved low-dimensional parameter domain for inferential statistics, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes.

blade representation↗
DOE OSTI2020

Wake Statistics of Different-Scale Wind Turbines under Turbulent Boundary Layer Inflow

Subscale wind turbines can be installed in the field for the development of wind technologies, for which the blade aerodynamics can be designed in a way similar to that of a full-scale wind turbine. However, it is not clear whether the wake of a subscale turbine, which is located closer to the ground and faces different incoming turbulence, is also similar to that of a full-scale wind turbine. In this work we investigate the wakes from a full-scale wind turbine of rotor diameter 80 m and a subscale wind turbine of rotor diameter of 27 m using large-eddy simulation with the turbine blades and nacelle modeled using actuator surface models. The blade aerodynamics of the two turbines are the same. In the simulations, the two turbines also face the same turbulent boundary inflows. The computed results show differences between the two turbines for both velocity deficits and turbine-added turbulence kinetic energy. Such differences are further analyzed by examining the mean kinetic energy equation.

17 WIND ENERGY↗
DOE OSTI2022

Attributes of Bi-Directional Turbomachinery for Pumped Thermal Energy Storage

Abstract In this paper, we (i) present a methodology for determining the aerodynamic performance of bi-directional turbomachines for pumped thermal energy storage, i.e., turbomachines designed to operate as a compressor in one direction, and then as a turbine in the opposite direction, (ii) carry out performance computations for such turbomachines, and (iii) propose principles for conceptual design of these devices. Focus is placed on using the energy storage cycle not only to identify the novel requirements placed on bi-directional turbomachines but also to estimate the effect of these requirements on the efficiency of the energy storage process. In particular, the difference between aerodynamic loading in forward and in backward operation causes the blading to work at incidences leading to the performance below maximum efficiency, resulting in a lower round-trip efficiency. The description of the design principles includes determination of the number of stages, definition of nondimensional parameters for blading selection, and optimization of two-dimensional blading for bi-directional operation. The assessment of stage count shows the relationship between relative Mach number, pressure ratio, and round-trip efficiency. The nondimensional parameters are assessed through a bi-directional analogue to existing “Smith charts,” for the efficiency of single-direction turbomachines, as a function of camber and stagger. The blade shape evaluation and optimization show how the blade profile can be modified to address the requirements of a bi-directional turbomachine, enabling an increase in round-trip efficiency of two percentage points compared to a baseline double circular arc configuration.

Engineering↗