Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “airfoil”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

G2Aero Database of Airfoils - Curated Airfoils

This dataset contains a curated set of 19,164 airfoil shapes from various applications and the data-driven design space of separable shape tensors (PGA space), which can be used as a parameter space for machine-learning applications focused on airfoil shapes. We constructed the airfoil dataset in two main stages. First, we identified 13 baseline airfoils from the NREL 5MW and IEA 15MW reference wind turbines. We reparameterized these shapes using least-squares fits of 8-order CST parametrizations, which involve 18 coefficients. By uniformly perturbing all 18 CST coefficients by +/-20% around each baseline airfoil, we generated 1,000 unique airfoils. Each airfoil was sampled with 1,001 shape landmarks whose x-coordinates followed a cosine distribution along the chord. This process resulted in a total of 13,000 airfoil shapes, each with 1,001 landmarks. In the second phase, we gathered additional airfoils from the extensive BigFoil database, which consolidates data from sources such as the University of Illinois Urbana-Champaign (UIUC) airfoil database, the JavaFoil database, the NACA-TR-824 database, and others. We undertook a thorough pre-processing step to filter out shapes with sparse, noisy, or incomplete data. We also removed airfoils with sharp leading edge and those exceeding our threshold for trailing edge thickness. Additionally, we thinned out the collection of NACA airfoils-- parametric sweeps of NACA airfoils with increasing thickness and camber present in BigFoil database-- by selecting every fourth step in the parameter sweeps. Finally, we regularized the airfoils by reparametrizing them with an 8-order CST parametrization (with 1,001 shape landmarks with x coordinated following cosine distribution along the chord) and removing airfoils with high reconstruction errors. This data pre-processing resulted in a set of 6,164 airfoils. In total, our curated airfoil dataset comprises 19,164 airfoils, each with 1,001 landmarks, and is stored in the curated_airfoils.npz file. Using this curated airfoil dataset, we utilized the separable shape tensors framework to develop a data-driven parameterization of airfoils based on principal geodesic analysis (PGA) of separable shape tensors. This PGA space is provided in PGAspace.npz file.

airfoils↗
DOE OSTI2022

Effect of the J-Shaped Wind Turbine Airfoil Opening Ratio and Thickness on the Performance of Symmetrical Airfoils

Abstract In this work, the effect of the inner opening ratio on the J-shaped airfoils aerodynamic performance was studied and documented for symmetrical airfoils. Three different airfoil thicknesses were investigated: small (NACA0008), medium (NACA0015), and large (NACA0024). For each airfoil thickness, effects of three inner opening ratios were analyzed: one-third, one-half, and two-thirds. The performance of each opening ratio was compared with the performance of the solid airfoil “zero opening ratio” for different angles of attack between 5 deg and 20 deg. All designs were simulated using the computational fluid dynamics (CFD) technology against experimental results for solid NACA4412 airfoil in the University of Wisconsin-Milwaukee (UWM) wind tunnel facility and other published experimental data. It was found that large eddy simulation yields accurate solutions with a smaller number of mesh cells compared to the k–ω turbulence model but with much longer computational time. The lift-to-drag ratio for all studied airfoils has a maximum value for solid airfoils compared to those equipped with openings. For airfoils equipped with 00.00% opening ratio “solid,” NACA0015 airfoil has the maximum lift-to-drag ratio. Furthermore, it was found that NACA0008 equipped with a 33.33% opening ratio has the best performance of all studied J-shaped airfoils.

Energy & Fuels↗
DOE OSTI2025

Open-Source Offshore Airfoil Summary (V.1.0)

The Open-Source Offshore (OSO) airfoils have been developed for research purposes for offshore wind turbines, offering a set of airfoils that align with modern turbine design requirements and industry design practices without proprietary constraints on research use. The eventual airfoil family will target the IEA 22 MW reference wind turbine, which was originally developed with the FFA airfoils. The two airfoils summarized in Table 1 (OSO-21-WT1 and OSO-30-WT1) started development as part of a family of airfoils being designed to target the IEA 22 MW wind turbine. The criteria used to design these airfoils are summarized in Table 1, which aim to encapsulate requirements of modern airfoils for offshore wind turbine applications, and were developed with feedback from industry and research experts. The airfoils were designed using XFOIL and candidate airfoils were then analyzed in RFOIL, which is considered more accurate than XFoil for high lift predictions of thicker airfoils. The design process for a preliminary family of airfoils is available, including a more detailed explanation of the design requirements and metrics similar to those used for these airfoils. Most of the design criteria are met for these two airfoils, with two exceptions. For both airfoils, the L/D Roughness Loss metric is exceeded (42% > 40% goal) and the desired lift coefficient margin over the design value (“CL_Margin”) was moderately exceeded (0.43 > 0.3) while smooth-stall characteristics (computed) were achieved. Note that all of the metrics were computed using RFOIL, and like other new airfoils, these will need to be experimentally validated at a range of Reynolds numbers. The airfoil coordinates will be shared publicly on Sadia National Laboratories’ public Github repository:

17 WIND ENERGY↗
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 OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗
DOE OSTI2025

Design of a Preliminary Family of Airfoils for High Reynolds Number Wind Turbine Applications

For the past 30 years, offshore wind turbines exhibited a continual pattern of growth that is expected to continue as the industry pushes for higher efficiency. Current designs for the next generation of wind turbines are so large that the chordwise Reynolds number of the blades is well beyond the design range of existing open-source airfoil families. This paper presents a preliminary family of new airfoils designed specifically for the needs of these next-generation offshore turbines, ranging from 21% thick to 30% thick with operating Reynolds numbers between 12 million and 18 million. These airfoils are intended to be alternative to the FFA airfoils that are commonly used on reference turbines such as the IEA 15MW and 22MW designs. In this work, airfoil performance metrics and design targets are developed, the design process is outlined, an optimization scheme is presented, and finally the airfoils and their simulated performance are compared to existing baselines. Lift to drag ratios in a clean condition were improved by up to 49.3% from the baseline FFA airfoil, and rough condition lift to drag ratio was improved by up to 9.3%. It is estimated that the cumulative improvements provided by this airfoil family would result in an approximately 1% increase of Annual Expected Power (AEP) for the 22 MW turbine compared to the current baseline.

airfoils↗
DOE OSTI2024

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

Invertible Neural Networks for Airfoil Design

We report the airfoil design problem, in which an engineer seeks a shape with desired performance characteristics, is fundamental to aerodynamics. Design workflows traditionally rely on iterative optimization methods using low-fidelity integral boundary-layer methods as higher-fidelity adjoint-based computational fluid dynamics methods are computationally expensive. Surrogate-based approaches can accelerate the design process but still rely on some iterative inverse design procedure. In this work, we leverage emerging invertible neural network (INN) tools to enable the rapid inverse design of airfoil shapes for wind turbines. INNs are specialized deep-learning models with well-defined inverse mappings. When trained appropriately, INN surrogate models are capable of forward prediction of aerodynamic and structural quantities for a given airfoil shape as well as inverse recovery of airfoil shapes with specified aerodynamic and structural characteristics. The INN approach offers a roughly 100 times speed-up compared to adjoint-based methods for inverse design. We demonstrate the INN tool for inverse design on three test cases of 100 airfoils each that satisfy the performance characteristics close to those of airfoils used in wind-turbine blades. All generated shapes satisfy the desired aerodynamic characteristics, demonstrating the success of the INN approach for inverse design of airfoils.

17 WIND ENERGY↗
DOE OSTI2022

Blade and airfoil damping configurations

The present embodiments set forth a blade including an airfoil, the airfoil including a tip cap, a pressure sidewall and a suction sidewall extending axially between corresponding leading and trailing edges and radially between the base and the tip cap. The blade, including the airfoil and base, being formed in at least two airfoil parts, each of the two airfoil parts including contacting edges engaging each other respective contacting edges, the contacting edges defining a joint for preloading each of the at least two parts with each other and with the base. The at least two airfoil parts forming the airfoil being retained to each other by an interference fit at the joint. The interference fit providing frictional damping of vibrations in the blade during blade operation.

Hart, Andrew Clifford↗
DOE OSTI2026

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

Regularizing Invertible Neural Networks for Airfoil Design Through Dimension Reduction

This work examines how dimension reduction can improve the performance of invertible neural networks (INN) for airfoil design. Design workflows are typically expensive, relying on many evaluations of high fidelity computational fluid dynamics (CFD) models. Furthermore, the inverse design problem is typically ill-posed. That is, multiple valid solutions exist that satisfy the design criteria. Regularization can reduce this inverse design space and simplify the problem. We study the use of subspace-based input dimension reduction to act as a regularizer for the INN model and improve the recovery of new airfoil shapes with desired performance characteristics. We find that the dimension reduction identifies two dominant modes, relating to airfoil thickness and camber, that optimally determine the airfoil's aerodynamics. We demonstrate the capability of the proposed INN model to generate 100 airfoils that satisfy the specific aerodynamic and structural characteristics.

airfoil design↗
DOE OSTI2021

Airfoil having internally cooled wall with liner and shell

An airfoil includes an airfoil section that has an airfoil wall that surrounds a cavity. The airfoil wall includes a first wall section comprised of a monolithic ceramic shell, a second wall section comprised of a CMC liner that is bonded along an interface to an interior side of the monolithic ceramic shell, and a ceramic barrier coating that is bonded to an exterior side of the monolithic ceramic shell. The airfoil wall includes at least one of: inter-section cooling channels that extend along the interface and are bounded by the monolithic ceramic shell and the CMC liner, or intra-section cooling channels that extend within the CMC liner.

Croteau, Paul F.↗
DOE OSTI2022

Wind turbine blade design with airfoil shape control using invertible neural networks

Wind turbine blade design is a highly multidisciplinary process that involves aerodynamics, structures, controls, manufacturing, costs, and other considerations. More efficient blade designs can be found by controlling the airfoil cross-sectional shapes simultaneously with the bulk blade twist and chord distributions. Prior work has focused on incorporating panel-based aerodynamic solvers with a blade design framework to allow for airfoil shape control within the design loop in a tractable manner. Including higher fidelity aerodynamic solvers, such as computational fluid dynamics, makes the design problem computationally intractable. In this work, we couple an invertible neural network trained on high-fidelity airfoil aerodynamic data to a turbine design framework to enable the design of airfoil cross sections within a larger blade design problem. We detail the methodology of this coupled framework and showcase its efficacy by aerostructurally redesigning the IEA 15-MW reference wind turbine blade. The coupled approach reduces the cost of energy by 0.9% compared to a more conventional design approach. This work enables the inclusion of high-fidelity aerodynamic data earlier in the design process, reducing cycle time and increasing certainty in the performance of the optimal design.

17 WIND ENERGY↗
DOE OSTI2022

Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach

This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Aerodynamic Sensitivity of a Novel Data-Driven Airfoil Shape Representation Framework

We explore the aerodynamic implications of a novel data-driven separable shape tensor framework used to represent discrete airfoil shapes. In this study, we construct a data-driven parameter space defined by separable shape tensors and informed by tens of thousands of distinct airfoils. We use this design space to generate new airfoil designs to study parametric sensitivities with respect to various aerodynamic responses. We use a HAM2D RANS solver to approximate the lift, drag, and moment coefficients for the generated airfoils at two different angles-of-attack. We analyze the robustness and sensitivities of using the separable shape tensor design space by examining the coverage of the aerodynamic response space, uncovering low-dimensional polynomial ridge approximations, and computing various sensitivity metrics. The results show that the data-driven design space produce significant variation in target aerodynamic quantities and facilitate highly accurate approximations (R^2 > 0.96) of one- and two-dimensional structures in each aerodynamic response. This further reduces the effective dimension to enable simplified design and optimization tasks.

aerodynamics↗
DOE OSTI2025

Stokes-dependent droplet collection efficiency on a NACA 0012 airfoil from droplet-informed simulations with statistical overloading

Accurate modelling of ice accretion on aircraft wings requires analysing droplet impingement on the surface to optimize the design of ice-protection systems. We perform Euler–Lagrange simulations of a droplet-laden flow impinging on a NACA 0012 airfoil. Our study includes water droplets with eight discrete sizes ranging from 1 to 160 microns. We vary the free-stream velocity of the incoming airflow in the range 60 ≤ U ≤ 240 m s −1 and the chord length of the airfoil in the range 0.5 ≤ c ≤ 2 m. Due to the dilute nature of supercooled clouds, one-way coupling is used in the simulations. The effects of droplet breakup and collision are also neglected. To reduce the computational cost, we employ statistical overloading of droplets, allowing us to simulate millions of impinging droplets in a time span on the order of milliseconds. Our results show that the droplet collection efficiency, which measures the likelihood of droplet impingement on the airfoil surface, increases with droplet size and free-stream velocity but decreases with airfoil size. We demonstrate that collection efficiency, impingement velocity and impingement angle are primarily dictated by a single non-dimensional parameter, the droplet Stokes number. We also identify a critical stagnation-streamline Stokes number below which impingements do not occur and use it to estimate the minimum droplet size for impingement. In addition, we observe droplet behaviour to become Stokes number independent at large values of the Stokes number. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗