Results for “FUN3D” — Page 18

NASA NTRSDate not supplied

Simulations of Mach 0.8 Transonic Truss-Braced Wing Aircraft Aerodynamics at High Angles of Attack

This paper presents numerical simulations of the Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft at high angle of attack using CFD solver FUN3D. Three different approaches - steady Reynolds Averaged Navier-Stokes (RANS), Unsteady Reynolds Averaged Navier-Stokes (URANS), and Delayed Detached Eddy Simulation (DDES) are used to simulation the TTBW aircraft at Mach number 0.8 and altitude 40, 000 ft. The pitch break phenomenon is observed at this flight condition. The simulation results show that the pitch break angle of attack occurs at 2.75˚ by DDES, 2.81◦ by URANS, and 2.89˚ by RANS. The DDES unsteady pressure begins to grow at angle of attack of 3.5˚, which might indicate the buffet onset, and rises one order of magnitude larger at angle of attack of 4.0˚. However, the URANS unsteady pressure begins to grow at angle of attack of 4.0˚. Both unsteady simulation results show the pitch break occurs before buffet onset.

TTBW↗

NASA NTRSDate not supplied

Simulation and Regression Modeling of X-59 Low-Boom Carpets Across America

The NASA X-59 aircraft is predicted to produce a significantly quieter cruise sonic boom than traditional N-wave-producing aircraft. A propagation simulation study was undertaken to quantify loudness levels, exposure size, and variability of the X-59 low-boom carpet using realistic atmospheric profiles across the contiguous United States of America (CONUS). Near-field pressure data of the X-59 in supersonic cruise from NASA’s fully unstructured Navier–Stokes three-dimensional (known as FUN3D) computational fluid dynamics code were propagated using NASA’s PCBoom code, which solves an enhanced Burgers equation along acoustic rays. Atmospheric profiles from the National Oceanic and Atmospheric Administration’s Climate Forecast System Version 2 database were used for propagation at 138 locations across the CONUS. Carpets at each location were generated for aircraft headings in the four cardinal directions. Over one million X-59 carpets were generated in total. The effects of the heading, season, geography, and climate zone on boom levels and exposure size are presented. Multiple linear regression models were developed to estimate carpet width and loudness metrics across the CONUS. These results inform regulators and mission planners on expected variations in boom levels and carpet extent from atmospheric variations. Understanding potential carpet variability is important when planning community noise surveys using the X-59.

X-59↗

NASA NTRSDate not supplied

Transition Analysis for the Pazy Wing

The Pazy wing, designed as a benchmark for highly flexible wings that experience large deformations, presents intriguing nonlinear phenomena due to its unique geometry with spanwise varying static deformation and low Reynolds number flow conditions. These phenomena include limit cycle oscillation related to laminar separation bubbles, boundary-layer transition, and dynamic stall. The primary focus of this study is to investigate transition behavior on the sagged Pazy wing under static loading by using the shear-stress-transport (SST)-based Langtry-Menter 𝜸 − 𝑹𝒆𝜽𝒕 transition model in the NASA OVERFLOW and FUN3D solvers. Based on the steady-state results obtained via different variants of the Langtry-Menter model, we investigate the likely transition behavior over a range of incidence angles. The predicted flowfield on the sagged Pazy wing is strongly three dimensional and the application of turbulence index as a viable indicator of the transition front in such flows is also explored. The study demonstrates the limitations of the Langtry-Menter model when applied to low Reynolds number flows. We introduce a spanwise periodic configuration to simulate a narrow-span section of the Pazy wing, specifically centered on one of the 14 nearly equidistant ribs spanning the wing. The computational results demonstrate that the spanwise periodic configuration provides a straightforward environment for studying transition characteristics as well as demonstrating systematic grid convergence and successful code-to-code comparisons. The research offers useful insights into the transition characteristics of flexible wings under various loading conditions, contributing to the field of aeroelastic simulations.

CFD Modeling↗

NASA NTRSDate not supplied

Linearized Frequency-Domain Gust Analysis and Adjoint-Based Sensitivities

Gust analysis is added to a linearized frequency-domain method in FUN3D, a NASA computational fluid dynamics solver. The method linearizes about a nonlinear static equilibrium condition and is therefore appropriate for problems with small perturbations such as transonic stochastic gust analysis. In addition to the gust analysis, adjoint-based sensitivities of stochastic gust constraints are implemented for multidisciplinary design optimization. The linearized frequency-domain gust model and adjoint-based sensitivities are described and verified. The method is applied to an optimization for mass minimization of the AGARD 445.6 wing subject to a stochastic gust constraint limiting the displacement of the wing tip.

Aeroelasticity↗

NASA NTRSDate not supplied

Assessment of UMUSCL Scheme for DNS of Turbulent Flows

Direct numerical simulation (DNS) is performed using an unstructured node-based finite-volume methodology for canonical turbulent flow configurations to assess the effect of the upwinding parameter (κ) in the unstructured monotonic upstream scheme for conservation laws (UMUSCL) reconstruction scheme. The NASA FUN3D solver is used, and four turbulent validation cases are considered: periodic Taylor-Green vortex, (2) flow through a plane channel, (3) flow through a channel with a constriction, and (4) flow over a flat plate. The results are presented for κ = 0.5, 0.9 and 1 on hexahedral grids, with the last value corresponding to a zero-numerical dissipation scheme. Detailed comparisons with available reference data indicate that for the smallest length-scales simulated, κ = 0.5 is too dissipative while κ = 1 could lead to an energy build-up. Overall, κ =0.9 yields satisfactory results in terms of accuracy and robustness for all the cases considered in this study.

Direct Numerical Simulation↗

NASA NTRSDate not supplied

Deep Neural Network Based Convergence Classification for Computational Fluid Dynamics

A supervised deep learning approach is coupled with heuristic convergence criteria to construct a classification model for detecting the completion (convergence) of computational fluid dynamics (CFD) simulations. Heuristic convergence criteria alone are not always sufficient and more complex decisions are often left to a human analyst. The proposed approach leverages heuristic convergence criteria as well as two deep neural network (DNN) models, one binary and one multi-class, to improve the efficiency and consistency of convergence classification across a wide range of flight regimes. The DNN models presented are each trained on a subset of ascent aerodynamic CFD simulations for NASA’s Space Launch System and were produced using NASA’s unstructured Navier-Stokes solver FUN3D. Individual solutions are analyzed intermittently and are classified as sufficiently converged, further iterations required, or switch from steady Reynolds Averaged Navier-Stokes (RANS) to unsteady RANS CFD based on the iterative histories of four aerodynamic coefficients. The implemented classification model is shown to produce solutions that closely correlate to solutions produced by a human analyst. This work lays groundwork for expanding the capabilities of DNNs for automating and improving more of the CFD process.

SLS↗

NASA NTRSDate not supplied

Aeroelastic Simulation of Transonic Truss-Braced Wing Aircraft for Flight Dynamic Stability Analysis

This paper presents aeroelastic simulation of the Mach 0.8 Transonic Truss-Braced Wing (TTBW) aeroelastic model using CFD solver FUN3D for the flight dynamic and control stability analysis. A jig twist optimization is performed to obtain a flight optimized jig twist for the flexible Mach 0.8 TTBW model. The developed aeroelastic model steady-state stability derivatives with respect to the angle of attack, angle of sideslip, and airspeed are calculated with perturbations in the angle of attack, angle of sideslip, and Mach number, respectively. A series of unsteady simulations is conducted for the developed aeroelastic model in pitch oscillation at various reduced frequencies. The dynamic stability derivatives are estimated using a frequency domain estimation method. The control derivatives of the elevator, rudder, and ailerons of the developed aeroelastic model are estimated for control stability analysis.

TTBW↗

NASA NTRSDate not supplied

Aeroelastic Simulation of Wind Tunnel Scale Aspect Ratio 13.5 Common Research Model with Trailing Edge Mini-Flaps

This paper presents a computation study of an aspect ratio 13.5 wind tunnel scale Common Research Model (CRM) with trailing edge mini-flaps using FUN3D. The aeroelastic simulations are conducted under free-air and inside wind tunnel conditions. The NASA Langley Research Center Transonic Dynamic Tunnel (TDT) is a closed circuit and slotted-test-section with a 16 ft ×16 ft test section. The structure mode shapes are extracted from the latest updated NASTRAN Finite-Element Model (FEM) model. A dynamic pressure sweep is performed to determine the control reversal boundary of the trailing edge mini-flaps. The wind tunnel wall interference effects on the mini-flap control derivatives are analyzed based on the simulations.

Aeroelastic↗

NASA NTRSDate not supplied

Transition Analysis for the Pazy Wing

The Pazy wing, designed as a benchmark for highly flexible wings that experience large deformations, presents intriguing nonlinear phenomena due to its unique geometry with spanwise varying static deformation and low Reynolds number flow conditions. These phenomena include limit cycle oscillation related to laminar separation bubbles, boundary-layer transition, and dynamic stall. The primary focus of this study is to investigate transition behavior on the sagged Pazy wing under static loading by using the shear-stress-transport (SST)-based Langtry-Menter 𝜸 − 𝑹𝒆𝜽𝒕 transition model in the NASA OVERFLOW and FUN3D solvers. Based on the steady-state results obtained via different variants of the Langtry-Menter model, we investigate the likely transition behavior over a range of incidence angles. The predicted flowfield on the sagged Pazy wing is strongly three dimensional and the application of turbulence index as a viable indicator of the transition front in such flows is also explored. The study demonstrates the limitations of the Langtry-Menter model when applied to low Reynolds number flows. We introduce a spanwise periodic configuration to simulate a narrow-span section of the Pazy wing, specifically centered on one of the 14 nearly equidistant ribs spanning the wing. The computational results demonstrate that the spanwise periodic configuration provides a straightforward environment for studying transition characteristics as well as demonstrating systematic grid convergence and successful code-to-code comparisons. The research offers useful insights into the transition characteristics of flexible wings under various loading conditions, contributing to the field of aeroelastic simulations.

CFD Modeling↗

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

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

This paper will describe three new supersonic jet cases for the NASA Turbulence Modeling Resource (TMR). These cases were taken from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test problem, which utilized data taken at NASA GRC consisting of non-intrusive measurements of jet plume velocities and temperatures, including both mean values and turbulence statistics. The three jet cases all used a Mach 1.63 axisymmetric nozzle: A heated jet for each at on-design conditions, a temperature-matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle. Computational fluid dynamics (CFD) solutions are obtained using three established Reynolds-averaged Navier-Stokes (RANS) codes: Wind-US, FUN3D, and Vulcan, all using the Menter Shear Stress Transport (SST) k-ω turbulence model. The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, that is - providing solutions that may not necessarily provide close comparisons with experimental data, but solve the posed CFD equations as intended.

nozzle↗

NASA NTRSDate not supplied

GlennICE Simulation of 24, 28 and 36 Inch Diameter eVTOL Propellers in Forward Flight

Flow solutions and water collection results are presented for three different propeller sizes. Flow solutions were generated using the NASA FUN3D software, and water collection was simulated using NASA GlennICE software. The geometries considered represent experimental articles that have been tested in the NASA Icing Research Tunnel. Each configuration has four blades. For the experiment, the same spinner centerbody is used for each propeller size. The different sized blades are similar for radii larger than about 40% of total span. Operating conditions were chosen to match operating points where experimental data was collected. Results are presented for both single bin simulations, as well as 7-bin Langmuir D distributions. For 80 microns MVD the Langmuir D distribution showed small differences compared to the single bin. For 15 microns MVD the Langmuir D distribution had a larger impact, especially on the spinner. In the experiment, the centerbody was unchanged as propeller size was changed. In this paper, additional simulations were completed with the 36 inch case perfectly scaled down to 24 and 28 inch. The perfect geometric scaling allows for an assessment of the effect of constant centerbody size on the results. While slight differences in pressure distribution on the propeller are visible between scaled and non-scaled centerbody, the effect seems quite minimal regarding water collection. All simulations were run very cold to ensure rime ice, since handling runback water in the non-inertial frame is beyond the scope of the present work. Ice shapes were generated using a single time step growth. The total mass on the blade and spinner was generally underpredicted for 15 microns MVD, and overpredicted for 80 microns MVD. When comparing ice shapes at operating conditions the numerical results produce very similar ice shapes.

Aircraft Icing↗

NASA NTRSDate not supplied

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

NASA NTRSDate not supplied

Effects of Aeroelasticity on Buffet Onset of Mach 0.8 Transonic Truss-Braced Wing

This paper presents a buffet simulation of the Mach 0.8 Transonic Truss-Braced Wing (TTBW) aeroelastic model using CFD solver FUN3D. A jig twist optimization is performed to obtain a flight optimized jig twist for the flexible Mach 0.8 TTBW model. A series of Delayed Detached Eddy Simulation (DDES) aeroelastic simulations has been conducted at Mach number 0.8 and altitude 40,000 ft using a structural damping of 2% to compare the pitch break prediction and the buffet behavior for the TTBW aeroelastic model. The unsteady pressure begins to grow at angle of attack of 4.0 ◦ , which might indicate the buffet onset, and rises one order of magnitude larger at angle of attack of 4.5 ◦ . In contrast, the buffet onset for the TTBW 1g model appears to be at the angle of attack of about 3.5 ◦ . The delay in the buffet onset of the aeroelastic model is due to the load relief from the aeroelastic effect.

TTBW↗

NASA NTRSDate not supplied

GlennICE Simulation of 24, 28 and 36 Inch Diameter eVTOL Propellers in Forward Flight

Flow solutions and water collection results are presented for three different propeller sizes. Flow solutions were generated using the NASA FUN3D software, and water collection was simulated using NASA GlennICE software. The geometries considered represent experimental articles that have been tested in the NASA Icing Research Tunnel. Each configuration has four blades. For the experiment, the same spinner centerbody is used for each propeller size. The different sized blades are similar for radii larger than about 40% of total span. Operating conditions were chosen to match operating points where experimental data was collected. Results are presented for both single bin simulations, as well as 7-bin Langmuir D distributions. For 80 microns MVD the Langmuir D distribution showed small differences compared to the single bin. For 15 microns MVD the Langmuir D distribution had a larger impact, especially on the spinner. In the experiment, the centerbody was unchanged as propeller size was changed. In this paper, additional simulations were completed with the 36 inch case perfectly scaled down to 24 and 28 inch. The perfect geometric scaling allows for an assessment of the effect of constant centerbody size on the results. While slight differences in pressure distribution on the propeller are visible between scaled and non-scaled centerbody, the effect seems quite minimal regarding water collection. All simulations were run very cold to ensure rime ice, since handling runback water in the non-inertial frame is beyond the scope of the present work. Ice shapes were generated using a single time step growth. The total mass on the blade and spinner was generally underpredicted for 15 microns MVD, and overpredicted for 80 microns MVD. When comparing ice shapes at operating conditions the numerical results produce very similar ice shapes.

Aircraft Icing↗

NASA NTRSDate not supplied

Computationally Efficient Frequency Domain Method for Dynamic Control Derivative Estimation with Application to Transonic Truss-Braced Wing

This paper presents a computationally efficient method for dynamic control derivative estimation technique via control surface oscillation numerical experiments. Unsteady RANS CFD simulations in FUN3D are performed to simulate the control surface oscillations with a prescribed truncated square wave containing sufficient frequencies of interest. The truncated square wave oscillation offers the computational efficiency which reduces the computational cost by almost an order of magnitude compared to a sine wave oscillation. The nonlinear effect of large control surface oscillation amplitudes creates a spillover effect whereby the frequency response at the same input frequencies contains not only the linear aerodynamic response but also nonlinear aerodynamic response. A correction procedure is developed to remove the spillover effect from the linear aerodynamic response. A frequency-domain regression is performed to estimate the dynamic control derivatives after the correction. The results generally agree with the previous results obtained from the sine wave oscillation.

Stability and Control↗

NASA NTRSDate not supplied

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

NASA NTRS2024

LAURA Users Manual: 5.7

This users manual provides in-depth information concerning installation and execution of Laura, version 5. Laura is a structured, multi-block, compu- tational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 Laura code toward a modular structure afforded by Fortran 2003. The refactoring improved usability and maintain- ability by eliminating the requirement for problem-dependent re-compilations, providing more intuitive distribution of functionality, and simplifying inter- faces required for multi-physics coupling. As a result, Laura now shares gas-physics modules, MPI modules, and other low-level modules with the Fun3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard instal- lation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

CFD hypersonics reentry↗