Search NASA⌕ Search

SEARCH · Search NASA

Results for “Blade element method”

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.

50 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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

Hybrid Data‐Driven Discovery of High‐Performance Silver Selenide‐Based Thermoelectric Composites

Optimizing material compositions often enhances thermoelectric performances. However, the large selection of possible base elements and dopants results in a vast composition design space that is too large to systematically search using solely domain knowledge. To address this challenge, a hybrid data-driven strategy that integrates Bayesian optimization (BO) and Gaussian process regression (GPR) is proposed to optimize the composition of five elements (Ag, Se, S, Cu, and Te) in AgSe-based thermoelectric materials. Data is collected from the literature to provide prior knowledge for the initial GPR model, which is updated by actively collected experimental data during the iteration between BO and experiments. Within seven iterations, the optimized AgSe-based materials prepared using a simple high-throughput ink mixing and blade coating method deliver a high power factor of 2100 µW m −1 K −2 , which is a 75% improvement from the baseline composite (nominal composition of Ag 2 Se 1 ). In conclusion, the success of this study provides opportunities to generalize the demonstrated active machine learning technique to accelerate the development and optimization of a wide range of material systems with reduced experimental trials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

High-order methods for wind turbine simulations and moving to next-generation platforms

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. This milestone represents an effort to increase the fidelity of Nalu-Wind at a fixed resolution through the implementation of a tensor-product based, matrix-free high order scheme. High order finite element methods have increased local work per datum communicated and have the potential to provide significantly more accurate solutions at a fixed number of degrees of freedom. Previous to this milestone, Nalu-Wind had an arbitrary order Control Volume Finite Element Method discretization as a solver option, but it required too much memory and was too slow to be of practical use. The work in this milestone addresses these issues by first implementing an implicit, high order solver that only partially assembles the global system. This reduces the memory footprint of the high-order scheme by orders of magnitude for higher polynomial orders. Second, a faster, tensor-product based method for evaluating the action of the left-hand side was implemented. This reduces the amount of computational work required by the scheme and dramatically enhanced the time-to-solution on example problems. Finally, this milestone is an evaluation of the value of high order methods in the wind application space. With the enhancements to memory and computational cost, accuracy vs. time-to-solution was evaluated for several resolutions on an under-resolved Taylor Green vortex test case. Results show that the high order scheme is cost-competitive with the production low-order schemes in Nalu-Wind, being moderately more expensive than the production edge-based vertex centered finite volume scheme. The evaluation of accuracy on the test case shows a potential benefit to high order at the highest resolution while not deteriorating accuracy on the lowest tested resolution. More work is needed to show value in the wind application, but positive strides have been made.

17 WIND ENERGY↗
DOE OSTI2024

Enhanced joining strength in additive-manufactured polylactic-acid structures fused by embedded heated metallic meshes

Additively manufactured thermoplastic polymers, such as polylactic acid (PLA), hold significant promise for sustainable engineering structures, including wind turbine blades. Upscaling these structures beyond the limitations of 3D printer build volumes is a challenge; fusion joining presents a potential solution. This paper introduces a displacement-controlled resistance welding process for PLA, as an alternative to the typical force-controlled methods. Here, we investigated the bonding quality of resistance-welded and adhesive-bonded PLA beams through three-point bending and measured the surface deformations using digital image correlation. Different metal meshes (30%/0.11 mm Ni—Cu, 34%/0.07 mm Ni—Cu, and 36%/0.25 mm Co—Ni) served as heating elements. The process parameters were varied for the 34%/0.07 mm Ni—Cu mesh to identify an optimum set of parameters. Results showed that this optimized displacement-controlled welding achieved 94% of the original strength of monolithic samples. This indicates that the new welding process not only ensures high-quality bonding and fine surface finishing but also promotes sustainability, recyclability, and economic efficiency in various polymer and composite structural applications.

Additive manufacturing↗
DOE OSTI2018

Creating Pathways to Success for Supersized Wind Turbine Blades: 2018 Workshop Summary Report

A workshop on Pathways to Success for Supersized Wind Turbine Blades was conducted by the U.S. Department of Energy’s (DOE’s) office of Energy Efficiency and Renewable Energy (EERE) at the Kimpton Hotel Palomar in Washington D.C., on March 6-7, 2018. Approximately 40-50 experts and industry stakeholders came together for the event, including manufacturers, transportation specialists, project developers, operators, engineering firms, consultants, and university researchers. Technical experts from the national laboratories and WETO were also present to engage in discussions about solving the challenges faced by supersized wind turbine blades. The workshop attendees participated in evaluating the current status of wind turbine blade design, manufacture, transportation, erection and operation, identifying constraints to cost-effective application of current technologies and methods for blades of increasing size, and discussing needs and opportunities for research, development and deployment of materials, manufacturing, structural configuration, and transportation. The workshop was one step within a larger initiative to identify specific R&D opportunities DOE could pursue to address technical barriers or implementation challenges faced by the U.S. wind energy industry to achieve continued decreases in LCOE. Following a plenary session, the 2-day workshop featured three group discussion sessions, with each session focusing on a specific “pathway” to enabling LCOE reductions for rotor blades of increasing size. The three pathways considered were “on-site manufacture,” “transport,” and “hybrid and alternative,” which included various options involving central manufacturing of sub-elements following on-site assembly. Each of the pathway group discussions was opened with a short presentation from one or more invited speakers, followed by an open discussion with balanced input from stakeholder groups and individuals. Participation among the workshop attendees was considered highly productive. Experts and stakeholders were engaged throughout the sessions and offered significant insights into the challenges and potential enabling technologies for supersized blades. Discussion highlights and take-aways for the three pathways are described.

17 WIND ENERGY↗
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2021

Alternative Damage Tolerant Materials for Wind Turbine Blades: An Overview

Current wind turbine blade materials may not be damage tolerant to the extent necessary to optimize the Levelized Cost of Energy (LCOE) of wind energy plants. Traditionally, wind turbine blades have been designed using a safe-life approach, but advances in inspection techniques and structural health monitoring solutions give rise to the opportunity to design wind turbine blades using a damage tolerant approach. Materials selection is a key element of da mage tolerant design, so the extent of the damage tolerance of alternative materials has been analyzed through a literature review and discussions with industry leaders. Fabrics and resin selection significantly affect the damage tolerance of composites. Changes to fabric architecture may include through-the-thickness (TTT) fibers, stretch-broken carbon fiber (SBCF) composites, and aligned discontinuous fiber reinforced composites (ADFRCs). Previous research has demonstrated that using TTT fibers in creases damage tolerance, but additional research is necessary to demonstrate the effectiveness of SBCFs and ADFRCs in mitigating damage. Several studies have demonstrated increased damage tolerance when toughened resin systems are used. In addition to toughened resin systems, thermoplastics have been shown to be tougher than thermosets. However, thermosets have been traditionally preferred in wind turbine blade manufacturing due to ease of manufacturing. Thermoplastic resin system s have been developed that can be used with conventional manufacturing methods but have yet to be studied for its damage tolerant capabilities. Furthermore, cost and stress analyses on where to effectively implement TTT fibers, SBCF composites, ADFRCs, and toughened resin systems must be executed prior to incorporating new materials into wind turbine blade manufacturing.

17 WIND ENERGY↗
DOE OSTI2022

Multiple Stream Low-Cost Recycling Method

The global composite industry generates large quantities of waste and which mostly end as landfill due to lack of meaningful end-use applications for the multiple waste streams. In a recent report by the Electric Power Research Institute (EPRI), waste generated by the wind industry could reach 370,000 tons a year of composite wind energy blades being decommissioned and scrapped. Wind energy is just one major industry utilizing composite materials. The waste generated by industry includes End-of-Life (EoL) materials and manufacturing process scrap. GreenTex Solutions has developed a unique and innovative technology to recycle the composite waste streams in a range of forms made from production waste and EoL materials. This includes manufacturing waste materials such as dry chopped fiber tow, loose fibers, shredded fibers from reinforcement fabrics, cured/semi-cured prepregs, and it also includes fully-cured composite structure waste (such as edge trims from cured parts) from manufacturing aircraft, automobiles, wind blades, boats, and composite cylinders (tanks). Current recycling methods involve recovering the structural fiber by removing the matrix resin through methods such as pyrolysis. The resulting fibers are used in injection molding or wet laid nonwoven mats and other usable forms. The GreenTex technology bypasses these intermediate steps to create the lowest possible recycling processing costs and the lowest embodied energy/CO2 emissions. The end-product from the GreenTex technology is a finished industrial composite part/application versus intermediate fibers or fabrics. The GreenTex manufacturing process enables cross-industry reuse of recycled feedstock by taking waste from multiple industries (wind energy, aerospace, marine, etc.) and recycles the waste into a product used in other industries. The initial target market application is structural flooring for intermodal shipping containers and truck bodies. One of the team’s key partners is Wabash National which produced 29,000 truck bodies in 2019. The current flooring system is comprised of solid oak “butcher board” laminated panels. Additionally, Mediterranean Shipping Containers (MSC) transports over 1.8 million twenty-foot equivalent units per year in intermodal shipping containers. The current container floor is laminated hardwood that is harvested from the rain forests of Central and South America. The project is to develop a flooring system made from recycled composites that can be qualified for both companies (Wabash and MSC). Initial prototypes validated that the recycled composites panels are lighter and thinner with much higher mechanical strength. These results suggest a typical truck trailer would have 20% lower tare weight. The GreenTex technology is not limited to flooring and is widely applicable to other transportation elements such as walls, roof elements, cab areas and related structural components. Under this project different composite waste streams were evaluated and then combined to develop a formulation that would meet the targeted performance criteria for a flooring system. Wet compression molding was used to fabricate plaques at different tonnage using various composite waste streams. The plaques were tested for flexure and impact.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Distinguishing Research Reactor Fuels Using Neutron Multiplicity - 25615

Neutron measurements are one of the reliable means of interrogation of fissile material. Neutron measurements are used for total counting (pure metal fissile samples) or counting coincidence/multiplicity (for impure fissile samples). Various flavors of neutron measurement are used routinely by the IAEA to determine the fissile mass of an unknown sample. The objective of this work is to design an experiment to distinguish the two different fuel types using neutron multiplicity method. The design of experiment is performed using Monte Carlo N- Particle (MCNP) transport code. General Electric (GE) designed 1-MWt pool-type research reactor at UMass Lowell (UMLRR) uses Materials Testing Reactor (MTR) plate-type fuel consisting of 19 full and 2 partial fuel assemblies, controlled by four boron control blades and one regulating control rod. Two different fuel types are available at the UMLRR facility. In 2008, UMass Lowell obtained slightly burnt fuel assemblies from the decommissioned WPI research reactor. The simulation results discussed in this paper show that within the limits of statistical error, the single and double moments for the WPI and UMLRR fuels can be measured. The two fuels have very similar composition and enrichments resulting in very similar neutron multiplicity moments. This implies that neutron multiplicity measurements alone would not be sufficient to distinguish between the two fuel types at low burnups. The change in the observed single-double ratio between each fuel type is less than 6 % which falls below the uncertainty in the efficiency of the JCC71 slab detectors. The neutron signature from low burnup and low enriched uranium fuel are not strong, however, these can be measured under low background and long measurement times. This would be an important consideration in designing the experiments and parallel gamma measurements will also be planned. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Roughness Sensitivity Comparisons of Wind Turbine Blade Sections

One explanation for wind turbine power degradation is insect roughness. Historical studies on insect-induced power degradation have used simulation methods which are either un- representative of actual insect roughness or too costly or time-consuming to be applied to wide-scale testing. Furthermore, the role of airfoil geometry in determining the relations between insect impingement locations and roughness sensitivity has not been studied. To link the effects of airfoil geometry, insect impingement locations, and roughness sensitivity, a simulation code was written to determine representative insect collection patterns for different airfoil shapes. Insect collection pattern data was then used to simulate roughness on an NREL S814 airfoil that was tested in a wind tunnel at Reynolds numbers between 1.6 x 10 6 and 4.0 x 10 6 . Results are compared to previous tests of a NACA 63 3 -418 airfoil. Increasing roughness height and density results in decreased maximum lift, lift curve slope, and lift-to-drag ratio. Increasing roughness height, density, or Reynolds number results in earlier bypass transition, with critical roughness Reynolds numbers lying within the historical range. Increased roughness sensitivity on the 25% thick NREL S814 is observed compared to the 18% thick NACA 63 3 -418. Blade-element-momentum analysis was used to calculate annual energy production losses of 4.9% and 6.8% for a NACA 63 3 -418 turbine and an NREL S814 turbine, respectively, operating with 200 μm roughness. These compare well to historical field measurements.

17 WIND ENERGY↗
DOE OSTI2022

Computational studies of hydraulic stressors for biological performance assessment in a hydropower plant with Kaplan turbine

We report Hydropower is currently one of the preeminent sources of renewable energy in the United States and globally. Hydropower plants also have detrimental impacts on the environment and ecology, including direct impacts to anadromous fish populations. The computational fluid dynamics (CFD) – based Biological Performance Assessment (BioPA) toolset is used for biological evaluations of fish passage through hydropower plants. The hydraulic environment of a hydropower plant was evaluated using CFD coupled with discrete element method (DEM) simulations. The predicted flow field and particle collision rate were validated against the experimental data in a water flume that has an idealized hydroturbine distributor geometry. Flow simulations were conducted to evaluate the hydraulic stressors, such as nadir pressure, fluid shear, runner collision, in a physical scale in a hydropower plant with Kaplan turbine which are responsible for injury and mortality of fish in a downstream migration. The cumulative exposure probability for the nadir pressure and collision with turbine runner was found to decrease with increased discharge rate. The lowest discharge rate shows the higher value of cumulative shear exposure probability. The maximum value of collision velocity increases with increased discharge rate. We offer the conclusions that will help in understanding various hydraulic stressors for biological assessment for environmentally sustainable hydroturbine passage.

13 HYDRO ENERGY↗
DOE OSTI2016

Reactivity Coefficient Calculation for a Low-Power LEU Research Reactor

Core design studies on a 20-MW thermal power low enriched uranium (LEU) fueled research reactor are underway at the NIST Center for Neutron Research (NCNR) with the primary purpose of providing quality neutron beams for scientific research. The 'tank-in-pool' design pattern is adopted for the new reactor: A compact core composed of two horizontally split halves are surrounded by a cylindrical D{sub 2}O reflector tank, and the tank is itself surrounded by a large H{sub 2}O pool, which functions as thermal and biological shielding. The split core consists of 18 MTR-type fuel elements which are evenly distributed into two regions. The LEU fuel used in the study is the U{sub 3}Si{sub 2}-Al dispersion fuel with U-235 enrichment of 19.75%, the only NRC certified LEU fuel in the U.S. The core is cooled and moderated by H{sub 2}O, and separated from D{sub 2}O by zircaloy core boxes. Hafnium control blades are utilized as both criticality and safety control elements for the reactor. They are controlled by a mechanical driver located at the bottom of the core but with the fully withdrawn positions at the top of the core. In a standard reactor calculation, it is necessary to evaluate the reactivity coefficient with respect to some physical parameters such as fuel or moderator temperature. For nearly all reactor designs, negative reactivity coefficients are required for safety. As a light-water cooled thermal reactor, the new reactor is designed to work in an under-moderated mode with the most achievable reactivity. The core is required to have negative reactivity coefficients under any normal or abnormal circumstances such as over power transient, loss of coolant or reflector, etc. Therefore the reactivity coefficient addresses reactor safety concerns as well as provides valuable quantitative parameters for reactor safety analyses. In this summary, three typical reactivity coefficients for the low-power LEU core are calculated using the Monte-Carlo code MCNP-6: the moderator temperature coefficient (MTC), the fuel temperature coefficient (FTC), and the void coefficient (VC). A direct perturbation methodology is applied to calculate all these coefficients. The perturbed case is produced by manually perturbing a single physical property with respect to the reference case (the one with normal operational conditions), as the multiplication factor (i.e. k{sub eff}) of the perturbed case is calculated by MCNP. The corresponding reactivity coefficient is thereby determined by dividing the reactivity change over the associated parameter changes. Core design studies are underway at NIST to develop a low-power LEU fueled beam reactor to advance the neutron source capability at NCNR for the next century. The reactivity effect caused by moderator temperature, fuel temperature and void is investigated though a direct perturbation approach using the MCNP code. The MTC, FTC, and VC for a representative startup core is obtained and the negative effect characteristics of these coefficients are observed. The reactivity effects under several postulated abnormal scenarios are also investigated. The overall negative feedback to the reactivity at each situation is confirmed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
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 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↗