SEARCH · Search NASA

Results for “Lift coefficient”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Aerodynamic Rotor Design for a 25 MW Offshore Downwind Turbine

Continuously increasing offshore wind turbine scales require rotor designs that maximize power and performance. Downwind rotors offer advantages in lower mass due to reduced potential for tower strike, and is especially true at large scales, e.g., for a 25 MW turbine. In this study, three 25 MW downwind rotors, each with different prescribed lift coefficient distributions were designed (chord, geometry, and twist) and compared to maximize power production at unprecedented scales and Reynolds numbers, including a new approach to optimize rotor tilt and coning based on aeroelastic effects. To achieve this objective the design process was focused on achieving high power coefficients, while maximizing swept area and minimizing blade mass. Maximizing swept area was achieved by prescribing pre-cone and shaft tilt angles to ensure the aeroelastic orientation when the blades point upwards was nearly vertical at nearly rated conditions. Maximizing the power coefficient was achieved by prescribing axial induction factor and lift coefficient distributions which were then used as inputs for an inverse rotor design tool. The resulting rotors were then simulated to compare performance and subsequently optimized for minimum rotor mass. To achieve these goals, a high Reynolds number design space was developed using computational predictions as well as new empirical correlations for flatback airfoil drag and maximum lift. Within this design space, three rotors of small, medium and large chords were considered for clean airfoil conditions (effects of premature transition were also considered but did not significantly modify the design space). The results indicated that the medium chord design provided the best performance, producing the highest power in Region 2 from simulations while resulting in the lowest rotor mass, both of which support minimum LCOE. The methodology developed herein can be used for the design of other extreme-scale (upwind and downwind) turbines.

downwind rotors↗

DOE OSTI2024

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗

DOE OSTI2025

Fuel cells for single-aisle regional aircraft: System configuration, performance and cost

A hydrogen fuel cell propelled electric aircraft can compete with incumbent turbofan technologies for single-aisle regional aircraft by coupling design of stack, air handling, thermal management, propulsion, and airframe to optimize performance. The stack operates at 95°C to facilitate heat rejection during take-off and below 75°C during cruise to extend lifetime and is oversized to satisfy power requirements at end of life. A multi-stage turbocompressor with a compression ratio >10 is selected to reach high stack power density at 11,300-m cruise altitude. The propulsion system is configured to accommodate air handling within the core duct, an inclined heat exchanger in the outer duct to limit the nacelle size, and variable area nozzles to independently control mass flows through the core and bypass ducts. The airframe is modified for maximum lift coefficient and longer balanced field length for dramatically reduced thrust during take-off, and the fuselage is stretched by 20% to store liquid hydrogen (LH 2 ). Modularization of power systems promotes safety in one engine inoperative scenarios and allows reaching specific power metrics for stack, balance-of-plant and fuel cell system (FCS), necessary for acceptable take-off weight. In conclusion, cost parity requires increase in FCS lifetime, LH 2 cost reduction, and improved FCS specific power.

Catalyst durability↗

DOE OSTI2021

Influence of turbulence intensity on vortex pattern for a rigid cylinder with turbulence stimulation in flow induced oscillations

The effect of ambient turbulence intensity on the flow induced oscillations of rigid circular cylinders with symmetric and asymmetric local turbulence stimulation is studied. Cylinders oscillate in one degree of freedom transversely to a steady uniform flow. Two-dimensional unsteady Reynolds-Average Navier-Stokes equations with the Spalart-Allmaras turbulence model are used to solve the problem numerically. Three freestream turbulence intensity values (0.2%, 1%, 5%) are used to study the influence of turbulence intensity on vortex pattern. Simulation results are compared with experimental data measured in the Marine Renewable Energy Laboratory of the University of Michigan in the range of 30,000 ≤ Re ≤ 110,000. The amplitude ratio, lift coefficient, vortex modes, and the interactions between vortices and cylinders are observed and discussed. In conclusion, major conclusions are: (a) Vortex patterns strongly depend on the freestream ambient turbulence intensity. (b) Low turbulence intensity can generate multi-vortex patterns at high Re. (c) High turbulence intensity results in fewer separations of a shed vortex into multiple vortices inducing emergence of elongated vortex patterns.

42 ENGINEERING↗

DOE OSTI2023

TEAMER: Water Tunnel Data from Testing the Pterofin Skimmer Concept

Pterofin's Skimmer concept relies on a flapping and pitching hydrofoil to extract hydrokinetic energy from water flows. The concept aims to utilize unsteady fluid dynamics phenomena (added mass, shed vorticity, and unsteady boundary layer development) to achieve higher lift coefficients, enabling increased power density of the hydrokinetic device and a fundamental shift in the rpm/torque scaling of the power take off compared with turbines. The Applied Research Laboratory at Penn State, in collaboration with Pterofin, designed and built a proof-of-concept flapping/pitching mechanism which was subsequently tested in ARL's 12-inch water tunnel facility. The mechanical power supplied to or extracted from the mechanism was measured for a range of hydrofoils provided by Pterofin over operating conditions including reduced frequency, Reynolds number, and the ratio between pitching and flapping amplitudes. The power lost to friction in the mechanism was removed from the net power measurement by means of a bare hub tare, with the resultant hydrodynamic power being used to calculate a mechanism-independent and non-dimensional power coefficient. The product of this effort is a dataset describing the power coefficient of a hydrofoil having simultaneous pitching and flapping motions, both of which are approximately sinusoidal. Power coefficients were collected for a range of primary design variables including: - Reduced frequency: 0.01 to 0.95 - Pitching/flapping peak angle ratio: 1.5 to 3.0 - Chord-based Reynolds number: 60,000 to 560,000 Secondary design variables relating to the hydrofoil geometry were explored including: - Aspect ratio - Planform shape - Section thickness distribution - Hydrofoil position relative to the pitching axis - Hydrofoil sweep angle relative to the pitching axis Measured data are provided in mean and time series formats. MATLAB scripts are provided which can be used to generate figures of time-averaged and phase-averaged hydrodynamic power coefficients calculated from the measured data. A complete description of the experiment and data reduction can be found in the Post Access Report for the Pterofin Skimmer test effort which will be available on the TEAMER website. This work was supported by the Pacific Energy Ocean Trust via a TEAMER award.

16 TIDAL AND WAVE POWER↗

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 OSTI2025

Effect of Surface Roughness on Dynamic Stall in Pitching Motion

Dynamic stall plays a critical role in determining the performance and stability of a wide range of fluid-dynamics systems in various engineering applications. This unsteady aerodynamic phenomenon is particularly significant for maneuvering aircraft wings, jet aircraft subjected to gust encounters, helicopter rotor blades, and wind turbine blades [1–3]. The prediction of the dynamic stall vortex (DSV) is challenging due to factors such as unsteady aerodynamics, three-dimensional (3-D) effects, turbulence and flow separation, incoming gust, and surface impact effects [4–6]. Hence, advanced computational techniques and modeling approaches in computational fluid dynamics (CFD) would be required to enhance the accuracy and reliability of DSV predictions in dynamic stall scenarios. In the past, Batther and Lee [7] employed delayed detached eddy simulations (DDES) to understand the flow physics associated with the onset of dynamic stall. Their approach demonstrated that DDES achieves results comparable to those obtained from large-eddy simulations at a reduced computational cost. In another study, Khalifa et al. [8] examined the 3-D aspects of dynamic stall on a NACA 0012 airfoil using DES solvers. In conclusion, the findings underlined the superiority of 3-D simulations over two-dimensional approaches, particularly in predicting the lift coefficient values and capturing dynamic stall stages more precisely.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

High Efficiency 2-Bladed Utility Wind Turbine Enhanced by CoFlow Jet Active Flow Control

The aerodynamic performance and power output of a high efficiency Co-Flow Jet (CFJ) 2-bladed Wind Turbine is studied in this paper. CFJ is a Zero-Net-Mass-Flux (ZNMF) active flow control method that dramatically increases airfoil lift coefficient and suppresses flow separation at a low energy expenditure. The 3D Reynolds Averaged Navier-Stokes (RANS) equations with one-equation Spalart-Allmaras (SA) turbulence model are solved to simulate the 3D flows of the wind turbines. The CFJ-Wind Turbine in this paper is modified from the NREL 3-bladed 3.4 MW reference wind turbine with the same RPM, blade length, distribution of chord, twist, and bending. The CFJ 2-bladed wind turbine utilizes the CFJ-NACA 64 series airfoils at most of the blade span from 31.88\% to the tip. The CFJ injection and suction slots are implemented along the full blade span to enhance the effectiveness and power output. The study shows that the 2-bladed CFJ turbine significantly increases power output compared with the 3-bladed baseline turbine at all wind speeds by 15.8\% at the wind speed of 9.8 m/s and by 25.9\% at 4 m/s respectively. The efficiency enhancement is substantially more at the low wind speed than at high wind speed. This would be very beneficial to increase the annual power output or capacity factor as wind turbines work at the speeds lower than the rated speed for most of the time. The results show that the flow field around the blade surface presents a significant enhanced suction peak near the blade leading edge, which results in a very high tangential force and thus increased wind turbine net power coefficient. Consistent with the previous study for the small NREL Phase VI research wind turbine, this paper demonstrates that CFJ 2-bladed wind turbine is an effective technology to extract more power from the wind at all speeds. This study is just an initial effort to demonstrate the advantage of the CFJ blade with no design optimization. It is believed that there is more potential to further enhance the efficiency improvement.

active flow control↗

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 OSTI2020

Lift Equivalence and Cancellation for Airfoil Surge–Pitch–Plunge Oscillations

A NACA 0018 airfoil in freestream velocity is oscillated in longitudinal, transverse, and angle-of-attack directions with respect to the freestream velocity, known as surge, plunge, and pitch. The lift-based equivalence method introduces phase shifts between these three motions to construct in-phase sinusoidal components for maximum lift, waveform construction. Lift cancellation is also determined with the exact negative pitch and plunge motion amplitudes found from the equivalence method to achieve out-of-phase wave destruction. Lift cancellation occurs when a combination of these motions is sought to obtain a constant lift magnitude throughout the oscillation cycle. To achieve both equivalence and cancellation of lift, a prescribed pure pitch amplitude through the Theodorsen theory equates the corresponding equivalent plunge amplitude and pitch–plunge phase shift. These Theodorsen, linear superposition findings of pitch–plunge are leveraged toward the Greenberg theory to determine a closed-form, surge–pitch–plunge solution through the addition of a surge–plunge phase shift and optimal surge amplitude for lift cancellation. The lift cancellation surge–pitch–plunge amplitudes define the equivalence amplitude investigated here and theoretically limit the experiment to combinations of the first lift harmonic of the Greenberg theory. The analytical results are then compared with experimental lift force measurements and dye visualization. Finally, the normalized lift differences due to unsteady wake and boundary-layer behavior are examined to explore the extents of the Greenberg theory for these cases of lift-based equivalence and cancellation.

42 ENGINEERING↗

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↗

DOE OSTI2020

Water Energy Resource Innovation on the Cavitation Characteristics

The main purpose of this study is to numerically correlate the amount of generated vapor over a hydrofoil to the lift and drag coefficients acting on it. As such, cavitation characteristics were investigated of a hydrofoil in the cavitating, sub-cavitating, and non-cavitating flows for different angles of attacks (AoA) with the high upstream flow velocity. The hydrofoil was tested in a square water tunnel with water entering the tunnel at various velocities for each AoA ranges from 9.1 m/s to 12.2 m/s. It was found that lift and drag forces acting on the hydrofoil follow the trend of the experimental data quite closely. While the cavitation can be identified by a unique number (averaged vapor volume fraction), the work done created an inverse correlation between this number and the cavitation number at the same angle of attack. The lift force declines with the increase in the vapor content on the hydrofoil surface, meanwhile the drag force peaks at certain vapor volume fraction, and then, a huge reduction occurs with the considerable decrease in the corresponding cavitation number. A fourth-order correlation generated between the lift to drag (L/D) and the cavitation number (σ). It was found the lift-to-drag ratio decreases by the formation of the cavitation over the hydrofoil, thus causing a drop in the efficiency of the turbomachines.

42 ENGINEERING↗

DOE OSTI2024

Note on two formulations of Crank-Nicolson method for Navier-Stokes equations

Here, we consider two formulations of the Crank-Nicolson (CN) method for the Navier-Stokes equations (NSE). The “natural” way of implementing CN for NSE is formally second order accurate in time for both velocity and pressure, whereas another formulation approximates pressure with only first order accuracy in time. Both versions of the method are applied to the benchmark problem of computing drag and lift in the flow around a cylinder. We show that the presumably more accurate version of the CN can create a solution with nonphysical oscillations and give incorrect predictions for the maximal drag coefficient, whereas the other formulation of the method predicts the drag and lift coefficients more accurately and does not introduce nonphysical oscillations. We locate the source of the issue and suggest several remedies.

Crank-Nicolson↗

DOE OSTI2022

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A rotorcraft in-flight ice detection framework using computational aeroacoustics and Bayesian neural networks

Abstract This work develops a novel ice detection framework specifically suitable for rotorcraft using computational aeroacoustics and Bayesian neural networks. In an offline phase of the work, the acoustic signature of glaze and rime ice shapes on an oscillating wing are computed. In addition, the aerodynamic performance indicators corresponding to the ice shapes are also monitored. These performance indicators include the lift, drag, and moment coefficients. A Bayesian neural network is subsequently trained using projected Stein variational gradient descent to create a mapping from the acoustic signature generated by the iced wings to predict their performance indicators along with quantified uncertainty that is highly important for time- and safety-critical decision-making scenarios. While the training is carried out fully offline, usage of the Bayesian neural network to make predictions can be conducted rapidly online allowing for an ice detection system that can be used in real time and in-flight.

42 ENGINEERING↗

DOE OSTI2025

Aerodynamic Sensitivities over Separable Shape Tensors

Here, we present a comprehensive aerodynamic sensitivity analysis of airfoil parameterization informed by separable shape tensors. This parameterization approach uniquely benefits the design process by isolating various well-studied shape characteristics, such as airfoil thickness, and providing a well-regulated low-dimensional parameter domain for aerodynamic designs. Exploring the aerodynamic sensitivities of this novel parameterization can provide valuable insights for more robust designs and future manufacturing efforts. We construct a data-driven parameter space of airfoils using principal geodesic analysis of separable shape tensors informed by a curated database containing almost 20,000 suitable engineering airfoils. Analyzing the shape reconstruction error and the maximum mean discrepancy between joint distributions of aerodynamic quantities, we study the dimensionality of the learned parameter space. This simple numerical experiment demonstrates a dramatic dimension reduction that retains design effectiveness and promotes regularity of the shape representations. Finally, we generate new airfoils and use the HAM2D Reynolds-averaged Navier–Stokes solver to predict lift, drag, and moment coefficients. We compute multiple sensitivity metrics to quantify and assert the consistency of parameter influence on the aerodynamic quantities. We also explore low-dimensional polynomial ridge approximations to motivate physical intuitions and offer explanations of the approximated sensitivities.

17 WIND ENERGY↗

DOE OSTI2025

The Fluid Dynamics Uncertainty Quantification Challenge Problem: XFOIL vs. MFOIL

Uncertainty quantification (UQ) has become more critical in aerospace engineering due to the growing dependence on computational tools for design optimization and performance analyses of aerospace vehicles. Even though the significance of UQ in assessing the credibility of computational analyses is well recognized, its costs and complexity impede its integration into standard practices, particularly in computational fluid dynamics (CFD) and other fluid analyses. This paper presents a UQ study for low-fidelity computational aerodynamics analyses with XFOIL and mfoil (i.e., the MATLAB version of XFOIL with several implementation modifications); these tools are utilized widely in both research and education. The main contributions of this paper are as follows: 1) improved precision in quantifying the uncertainty of the baseline Monte Carlo results used to benchmark surrogate modeling techniques for UQ, 2) quantification of the effect of the implementation differences between XFOIL and mfoil on solution quantities of interest (QoIs), such as lift and pitching moment coefficients, and 3) development of an open-source UQ library for use with XFOIL and mfoil, which has educational values and helps promote UQ for fluid analyses with aerospace applications. Results and discussions revolve around cases 1-4 of the challenge problem posed by the AIAA Fluid Dynamics Technical Committee’s Uncertainty Quantification Discussion Group (UQDG). In case 3, this work employs CFDverify, an open-source solution verification software, to quantify the discretization error and evaluate the extrapolated QoIs based on the grid convergence index (GCI). This UQ study differentiates itself from previous studies in the rigor of handling baseline Monte Carlo uncertainty and in including mfoil, which is a more accessible alternative to XFOIL. Finally, despite the growing computing power, low-fidelity computational tools remain valuable, such as for aerodynamic shape optimization at Mach numbers below 0.65 and low-to-mid Reynolds numbers.

Lay, Aidan S [University of Tennessee, Knoxville (↗