Results for “flow solver” — Page 53

NASA NTRSDate not supplied

Mixed-Element USM3D Contributions to the 4th AIAA High-Lift Prediction Workshop

This paper discusses results of the mixed-element USM3D (USM3D-ME) simulations performed for the 4th AIAA High-Lift Prediction Workshop. The workshop was separated into six Technical Focus Groups to investigate the impact of geometry modeling, grid, and computational methods for predicting high lift flows. This work was performed under the Fixed Grid RANS Technical Focus Group. The primary geometry selected for the workshop was the High-Lift Common Research Model. The performed simulations included a flap deflection study and both a grid refinement study and pitch sweep for the nominal flap deflection configuration. The results show that USM3D-ME RANS solutions, generally, tends to underpredict the lift coefficient and to predict a less negative pitching moment relative to the experimental data. The predicted drag coefficient values agree better with experiment for smaller angles of attack but were observed to be larger than experiment for the largest angle of attack simulated. The results of the grid refinement study demonstrated a lack of grid convergence for the provided grid family. The results of the grid refinement study are consistent with the submissions to the Fixed Grid and Mesh Adaptation Technical Focus Groups. Grid convergence for the provided grid family remains elusive for the international community. A 2D Multielement Airfoil configuration was included to enable a turbulence model verification study, which illustrated favorable agreement between USM3D-ME and the solutions provided by other flow solvers.

CFD↗

NASA NTRSDate not supplied

High Speed Data Investigation of an LOX/LNG Combustor Assisted by Non-linear Acoustic Simulations

During testing of a regenerative cooled liquid oxygen/liquid natural gas combustor large amplitude combustion instabilities were observed in half of the series conducted. These instabilities produced an unsteady pressure Root Mean Square (RMS) value greater than the chamber steady state pressure. Quick look investigations of the test data presented some standard indications of non-linear tangential instability such as a non-sinusoidal time series, and larger peaks than trough. Also observed in the data quick look were some atypical findings including a first harmonic amplitude larger than the fundamental, and numerous side bands on the fundamental and subsequent harmonics. A simplified model of the combustor was created using a commercially available finite element software. The first tangential acoustic mode was then simulated using a non-linear wave equation in conjunction with a time explicit discontinuous Gelarkin flow solver. Numerous simulations were evaluated providing valuable insight and reproducing the atypical findings in the data. Side bands were created through non-linear signal modulation like an AM radio. This happens due to the non-linear steepening and relaxing that occurs in large amplitude oscillations. The simulation was also able to reproduce a signal with a first harmonic amplitude larger than the fundamental. Investigating the pressure signal at multiple locations across the chamber head end revealed radial dependent frequency content. A literature review of the topic revealed corroboration of these effect in idealized geometries.

Combustion Stability↗

NASA NTRSDate not supplied

High Speed Data Investigation of an LOX/LNG Combustor Assisted by Non-linear Acoustic Simulations

During testing of a regenerative cooled liquid oxygen/liquid natural gas combustor large amplitude combustion instabilities were observed in half of the series conducted. These instabilities produced an unsteady pressure Root Mean Square (RMS) value greater than the chamber steady state pressure. Quick look investigations of the test data presented some standard indications of non-linear tangential instability such as a non-sinusoidal time series, and larger peaks than trough. Also observed in the data quick look were some atypical findings including a first harmonic amplitude larger than the fundamental, and numerous side bands on the fundamental and subsequent harmonics. A simplified model of the combustor was created using a commercially available finite element software. The first tangential acoustic mode was then simulated using a non-linear wave equation in conjunction with a time explicit discontinuous Gelarkin flow solver. Numerous simulations were evaluated providing valuable insight and reproducing the atypical findings in the data. Side bands were created through non-linear signal modulation like an AM radio. This happens due to the non-linear steepening and relaxing that occurs in large amplitude oscillations. The simulation was also able to reproduce a signal with a first harmonic amplitude larger than the fundamental. Investigating the pressure signal at multiple locations across the chamber head end revealed radial dependent frequency content. A literature review of the topic revealed corroboration of these effect in idealized geometries.

Combustion Stability↗

NASA NTRSDate not supplied

Implementation and Verification of the SST-𝜸 and SA-AFT Transition Models in FUN3D

The transition modeling capability in the NASA unstructured FUN3D suite of codes has been augmented by incorporating two additional transport-equation-based transition models, namely, the 𝜸 transition model coupled with Menter’s 2003 Shear-Stress Transport (SST) turbulence model and Coder’s Amplification Factor Transport (AFT) transition model combined with the Spalart-Allmaras (SA) turbulence model. Both of these models, SST-𝜸 and SA-AFT, are used to compute transition characteristics of several test cases from the literature. The selected test cases cover a range of transition scenarios including bypass transition over a flat plate (ERCOFTAC T3A and T3A- flow configurations), transition in the presence of separation bubbles (NLF-0416 airfoil at selected angles of attack), and natural transition due to Tollmien-Schlichting instabilities (flat-plate experiment by Schubauer and Klebanoff). To address the urgent need for model verification as highlighted in recent workshops on transition modeling, the drag, lift, surface pressure, and skin-friction coefficients from the FUN3D solutions for a sequence of meshes are compared with the results obtained with the same models implemented in the NASA structured grid flow solver OVERFLOW. Comparisons between the local and global surface coefficients from the FUN3D and OVERFLOW solvers for all test cases resulted in good agreement for the finest meshes.

Verification↗

NASA NTRSDate not supplied

Validation of a Mid-Fidelity Approach for Aircraft Stability and Control Characterization

This paper describes an efficient computational approach for aircraft aerodynamic model development intended for use in flight dynamics simulations. A commercial surface-vorticity flow solver called FlightStream is employed to predict performance, stability, and control characteristics for a NASA subscale electric vertical takeoff and landing aircraft in its isolated-airframe configuration. A framework is described to rapidly convert OpenVSP geometry into a form compatible with the FlightStream software. FlightStream predictions for the variation of aerodynamic force and moment coefficients with airflow angles and control surface deflection angles are compared to static wind-tunnel data to illustrate the process and assess the accuracy of FlightStream solutions. The results and approach presented in this paper are encouraging for use in rapid aerodynamic modeling and flight dynamics simulation development early in the aircraft design process.

aerodynamic modeling↗

NASA NTRSDate not supplied

Multi-Fidelity Computational Analysis of a Quiet Single-Main Rotor Helicopter for Air Taxi Operations

A multi-fidelity computational study has been conducted on NASA’s quiet single-main rotor helicopter concept for urban air mobility using computational fluid dynamics and rotorcraft comprehensive analysis tools. Two approaches have been followed in order to simulate the rotor flows: (1) a rotor disk model, which models the rotors as source terms, and (2) a high-fidelity approach using high-order accurate schemes and dual-time stepping, which simulates the rotor with its individual rotating blade grids. The delayed detached-eddy simulation model has been employed. A loose-coupling approach between the flow solver and a rotorcraft comprehensive code is utilized to include vehicle trim and motion variables. The quiet single-main rotor helicopter is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

RVLT↗

NASA NTRSDate not supplied

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

CFD↗

NASA NTRSDate not supplied

Simulations of Coupled Shoulder Ablation of A Conceptual Aeroshell

The current study explores coupled numerical simulations of flow dynamics and material response affect ablation of the shoulder region of an entry capsule used in space missions like the Mars Sample Return (MSR) project. A new computational tool named Ares, developed at NASA Ames Research Center, is applied to carry out the coupled simulations by interconnecting well-established US3D (flow solver) and Icarus (material response solver), and the material considered is three-dimensional mid-density carbon phenolic (3MDCP). Particularly, this work extends the uncoupled analysis of Shrestha et al. 2022 by incorporating the coupling mechanism across the ablating surface between fluid and the capsule geometry. For simplicity, an axisymmetric geometry of the entry system is considered. A few mission-relevant parametric studies such as angles of orientation of the material layers with respect to the streamwise direction, backshell boundary conditions, and alternation of convective heat-transfer coefficients are investigated to essentially explore Ares capability with the primary focus on shoulder location.

Coupled Ablation Modeling↗

NASA NTRSDate not supplied

Developing New Tools for Modeling Rocket Plume-Surface Interactions

With NASA’s goal to land the next human on the lunar surface in the next few years, it has become vitally important that we have a better understanding of how future landing spacecraft will interact with the unique properties of regolith¬¬––the layer of loose, unconsolidated dust and rock on the lunar surface¬¬––which can cause hazards like visual obstructions, particulate clouds, and cratering of the landing zone. Researchers from the Fluid Dynamics Branch at NASA’s Marshall Space Flight Center are performing plume-surface interaction (PSI) simulations between lander engine plumes and unprepared regolith surfaces, and have developed new tools to provide predictive PSI environments for various NASA projects and missions, including the Human Lander System (HLS), Commercial Lunar Payload Services (CLPS), and future Mars landers. These tools allow the researchers to determine how to best meet the simulation and time requirements for each project by varying model fidelity. The highest fidelity tool is the Gas Granular Flow Solver (Loci/GGFS) that models gas-particle multi-phase interactions to predict regolith cratering and ejection of particles into the immediate surroundings of the lander. At its highest fidelity, it can model microscopic regolith particle interactions with a particle size/shape distribution that statistically replicates actual regolith, however, to be most effective with today’s computing resources, it is currently run using only one to three equivalent particle sizes/shapes. The team also incorporated engineering models into their software suite to create production-ready hybrid tools with reduced fidelity. At the lowest fidelity, the computational fluid dynamics (CFD) code Loci/CHEM+DIGGEM can predict crater depth over time by relating local CFD-predicted surface shear stresses to a model of erosion mass flux.

plume surface interaction↗

NASA NTRSDate not supplied

GPU Implementation of the OVERFLOW CFD Code

The high-performance computing (HPC) landscape is quickly changing to systems where most of the performance comes from specialized chips, specifically graphics processing units (GPUs). Such GPU systems are throughput machines, where efficient use of the GPU often requires code refactoring to expose a few orders of magnitude more fine grain parallelism than was previously used on the CPU. Recent modifications to OVERFLOW, an overset, structured grid, computational fluid dynamics flow solver, written in Fortran will be presented. These modifications include both code modernization efforts and algorithmic changes to enable OVERFLOW to efficiently utilize GPUs. Many of these algorithmic changes would likely also be applicable for other structured grid, stencil-based codes wanting to utilize GPUs. The capabilities that have been ported to run on the GPUs are presented, along with the performance gains of the GPU version relative the CPU version of OVERFLOW.

GPU Programming↗

NASA NTRSDate not supplied

Shock-Buffet Prediction Report in Support of the High Angle Working Group at the Third Aeroelastic Prediction Workshop

This paper presents a summary of the computational shock-buffet results associated with the third Aeroelastic Prediction Workshop (AePW-3) High-Angle Working Group. The computational results are compared against experimental data collected during the Oscillating Turn Table (OTT) Benchmark Supercritical Wing (BSCW) test campaign conducted in NASA Langley Research Center Transonic Dynamics Tunnel (TDT) in early 2000. In addition to experimental forced oscillation data, unsteady pressure data was collected at a single spanwise station at transonic Mach numbers, several values of dynamic pressure, and fixed angles of attack. It is acknowledged that the current experimental data are limited, but future experiments are planned to obtain pressure data at two wingspan stations in addition to unsteady Pressure Sensitive Paint (uPSP). With that, the computational results presented here are considered to be semi-blind. The results are presented in the form of time-varying pressure coefficients, running-mean values of the pressure coefficients, Power Spectral Density (PSD) and Strouhal number computations. The results are split into two parts. In the first part, the computational results obtained by each team using their flow solvers and grids are compared. In the second part, a single computational tool was used to compute unsteady pressure on participant-generated grids. In addition, a description of each participating team’s software and methods is included.

Pawel Chwalowski↗

NASA NTRSDate not supplied

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

Computational Fluid Dynamics↗

NASA NTRS2024

Computational Techniques to Generate Space Launch System Aerodynamic Databases

This document describes the reasoning and trade studies used to evaluate tools for constructing the aerodynamic lineload databases for the liftoff and transition phases of flight for the Space Launch System. Three computational fluid dynamics codes (USM3D, FUN3D, Kestrel) were investigated with various turbulence models, as well as detached eddy simulation variants for the launch vehicle in free air and in proximity to the tower. Decisions were made mostly based on results from brief developmental studies performed in response to specific, unforeseen challenges that were encountered in the analysis of a given configuration. The need to develop databases in a timely manner, as well as accurately capture the expected leeward-wake flowfield characteristics, led to the selection of the Kestrel flow solver with its delayed detached eddy simulation method, the Spalart-Allmaras turbulence model, and the adaptive mesh refinement capability in the off-body Cartesian grid region.

Karen A. Deere↗

NASA NTRSDate not supplied

GPU Implementation of the OVERFLOW CFD Code

The high-performance computing (HPC) landscape is quickly changing to systems where most of the performance comes from specialized chips, specifically graphics processing units (GPUs). Such GPU systems are throughput machines, where efficient use of the GPU often requires code refactoring to expose a few orders of magnitude more fine grain parallelism than was previously used on the CPU. Recent modifications to OVERFLOW, an overset, structured grid, computational fluid dynamics flow solver, written in Fortran will be presented. These modifications include both code modernization efforts and algorithmic changes to enable OVERFLOW to efficiently utilize GPUs. Many of these algorithmic changes would likely also be applicable for other structured grid, stencil-based codes wanting to utilize GPUs. The capabilities that have been ported to run on the GPUs are presented, along with the performance gains of the GPU version relative the CPU version of OVERFLOW.

GPU Programming↗

NASA NTRSDate not supplied

Coaxial Rotor CFD Validation and ML Surrogate Model Generation

A fixed-pitch speed-controlled coaxial rotor system was tested in the NASA Langley Transonic Dynamics Tunnel in September 2022. The rotors have a diameter of 1.35 meters and an inter-rotor spacing of 25% of the diameter. Though this test was focused on the NASA New Frontiers Dragonfly mission, the resulting dataset is relevant to a wide array of applications of multirotor vehicles, especially those using fixed-pitch variable-speed rotors. Most notably beyond Dragonfly, perhaps, is the application of coaxial rotor pair configurations for eVTOL and Urban Air Mobility (UAM) aircraft. This work provides a thorough CFD validation study quantifying coaxial rotor performance estimation with accuracy on order 5-10% using an efficient hybrid BEMT-URANS flow solver over a wide range of operating conditions. This accuracy was achieved using novel approaches for the construction of both the C81 airfoil performance lookup tables and the BEMT rotor model. These novel approaches are combined with advanced scripting to further accelerate the commercial off-the-shelf CFD solver on GPU accelerated machines. Finally, the development of machine learning surrogate models is presented to produce highly efficient and accurate rotor performance predictions for fixed-pitch variable-speed multirotor aircraft over the complete flight regime including scout, cruise, climb, descent, as well as limiting cases of vortex-ring and windmill-brake states.

Coaxial↗

NASA NTRSDate not supplied

Airfoil Table Generation and Validation for the VR-12 and SSC-A09 Airfoils and Quadrotor Performance Prediction

As part of NASA’s Urban Air Mobility (UAM) mission, the Revolutionary Vertical Lift Technology (RVLT) project is creating rotorcraft reference designs that can be used by the rotorcraft community to develop quiet, efficient, and safer air vehicles. At NASA, this effort is made possible by the use of the RVLT toolchain software. Airfoil tables are often a critical part of the conceptual design process as they inform both comprehensive analysis and Computational Fluid Dynamics (CFD) codes of 2D aerodynamic coefficients, such as lift, drag, and moment coefficients. These airfoil tables, typically generated during experimental testing, are often proprietary and thus not widely distributable; further, there is not always clear or available reference documentation that provides information on the test conditions. The airfoil tables from these experiments are referred to in this work as Legacy airfoil tables. The scope of this study is to generate airfoil tables with the NASA OVERset grid CFD FLOW solver (OVERFLOW) coupled with the AirFoil Table Generator (AFTGen) software and validate these airfoil tables based on performance predictions calculated in Comprehensive Hierarchical Aeromechanics Rotorcraft Model (CHARM). The overall goal is to provide CFD generated airfoil tables, validate them through comprehensive analysis, and publish the airfoil tables for public distribution. The RVLT quadrotor reference model was selected as the basis for all analyses. The quadrotor uses the VR-12 and SSC-A09 airfoils for its rotor blades. Simulations were run in OVERFLOW for the VR-12 and SSC-A09 airfoils for ±20° angles of attack, Mach numbers from 0.3 to 1.0, and a Reynolds-Mach proportionality constant of 1.023x107. The C81 format airfoil tables from AFTGen were then blended with Legacy NACA 0012 airfoil table data for ±20° to ±180° angle of attack range. CHARM was used to compute thrust coefficient, power coefficient, and figure of merit for an isolated rotor in hover using proprietary legacy airfoil tables for the VR-12 and SSC-A09, as well as with the OVERFLOW generated airfoil tables for the VR-12 and SSC-A09. Overall, the OVERFLOW and Legacy performance predictions agreed well, with the best correlation between simulated and Legacy data observed with the VR-12 airfoil tables.

Airfoil↗

NASA NTRSDate not supplied

Simulation of a Mach 2 Jet in a Supersonic Retropropulsion Configuration

High-fidelity simulations of nozzle flows relevant to supersonic retropropulsion environments have been conducted using CRoCCo-AMR, a validated flow solver with curvilinear adaptive mesh refinement. The computational approach for simulating retropropulsion utilizing generalized curvilinear adaptive mesh refinement are discussed in detail. Initial results are presented for a single nozzle condition with comparison to available experimental data. Strong agreement is demonstrated in key regions of the flowfield. Mean and statistical analyses are performed to characterize unsteady, turbulence quantities present in the canonical retropropulsion environment. The findings herein suggest the presence of a dominant low-frequency mode that appears across the flowfield originating from the primary vortex cores where the freestream and retropropulsion jet interact.

retropropulsion↗

NASA NTRSDate not supplied

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗