Search NASA⌕ Search

SEARCH · Search NASA

Results for “FUN3D”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 343 records · Page 19

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

Three new supersonic jet cases were built 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-CFD, all using the Menter Shear Stress Transport k- turbulence model with vorticity source term (SST-V). The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, meaning that they solve the posed CFD equations as intended. The three codes generally provided very close agreement with each other for jet plume quantities, with the only exception being static temperature. This discrepancy was determined to be the result of whether or not turbulent kinetic energy was considered in the definition of total energy. Validation of the SST-V turbulence model for this case (its level of agreement/disagreement with experimental data) is also addressed.

turbulance↗

Study of Trailing Edge Mini Flaps Interactions on a Wind Tunnel Scale Aspect Ratio 13.5 Common Research Model

This paper presents a computation study of an aspect ratio 13.5 wind tunnel scale Common Research Model (CRM) with trailing edge mini-flaps. A vortex-lattice model of the CRM is developed with transonic small disturbance and integral boundary-layer corrections coupled with NASTRAN equivalent beam model for the rapid CRM aeroelastic simulation. The simulations are performed with individual trailing edge mini-flap deflected and a series of groups of two adjacent trailing edge mini-flaps deflected at the same time. A surrogate model is developed to estimate the aerodynamic interactions between the flaps. To validate the surrogate model, static aeroelastic simulations are performed using FUN3D CFD solver.

Interaction↗

Error Estimate of the Ares I Vehicle Longitudinal Aerodynamic Characteristics Based on Turbulent Navier-Stokes Analysis

Numerical predictions of the longitudinal aerodynamic characteristics for the Ares I class of vehicles, along with the associated error estimate derived from an iterative convergence grid refinement, are presented. Computational results are based on the unstructured grid, Reynolds-averaged Navier-Stokes flow solver USM3D, with an assumption that the flow is fully turbulent over the entire vehicle. This effort was designed to complement the prior computational activities conducted over the past five years in support of the Ares I Project with the emphasis on the vehicle s last design cycle designated as the A106 configuration. Due to a lack of flight data for this particular design s outer mold line, the initial vehicle s aerodynamic predictions and the associated error estimates were first assessed and validated against the available experimental data at representative wind tunnel flow conditions pertinent to the ascent phase of the trajectory without including any propulsion effects. Subsequently, the established procedures were then applied to obtain the longitudinal aerodynamic predictions at the selected flight flow conditions. Sample computed results and the correlations with the experimental measurements are presented. In addition, the present analysis includes the relevant data to highlight the balance between the prediction accuracy against the grid size and, thus, the corresponding computer resource requirements for the computations at both wind tunnel and flight flow conditions. NOTE: Some details have been removed from selected plots and figures in compliance with the sensitive but unclassified (SBU) restrictions. However, the content still conveys the merits of the technical approach and the relevant results.

Aerodynamic characteristics↗

Ares I Vehicle Computed Turbulent Ascent Aerodynamic Data Development and Analysis

An overview of the computational ascent aerodynamic data development from a Reynolds-averaged Navier-Stokes flow solver for the Ares I vehicle design is presented. The computed results are assessed for grid and turbulence model effects and verified against the results obtained from other flow solvers. The numerical predictions are analyzed for the surface pressure, sectional line loads, longitudinal aerodynamic and rolling moment coefficients, and their trends with respect to angle of attack, Mach number, and the vehicle's roll angle. The results revealed that the solution development with the Spalart-Allmaras turbulence model was the most robust, stable, and efficient. These predictions were generally found to compare well with experimental data. Relative to wind-tunnel flow conditions, the results at flight Reynolds number showed the largest reduction of about 7% in the computed longitudinal aerodynamic coefficients. Protuberance size and relative position were found to have a significant effect on the vehicle's force and moment coefficients. The strake concept was proposed and shown to reduce the vehicle's maximum rolling moment coefficient Results and analyses have demonstrated the feasibility of the present numerical method, as an engineering tool, for predicting the external flow aerodynamic characteristics of this class of vehicle designs.

Aerodynamic characteristics↗

Experimental and Computational Sonic Boom Assessment of Boeing N+2 Low Boom Models

Near-field pressure signatures were measured and computational predictions made for several sonic boom models representing Boeing's Quiet Experimental Validation Concept (QEVC) supersonic transport, as well as three axisymmetric calibration models. Boeing developed the QEVC under a NASA Research Announcement (NRA) contract for Experimental Systems Validations for N+2 Supersonic Commercial Transport Aircraft, which was led by the NASA High Speed Project under the Fundamental Aeronautics Program. The concept was designed to address environmental and performance goals given in the NRA, specifically for low sonic boom loudness levels and high cruise efficiency, for an aircraft anticipated to enter service in the 2020 timeframe. Wind tunnel tests were conducted on the aircraft and calibration models during Phases I and II of the NRA contract from 2011 to 2013 in the NASA Ames 9- by 7-Foot and NASA Glenn 8- by 6-Foot Supersonic Wind Tunnels. Sonic boom pressure signatures were acquired primarily at Mach 1.6 and 1.8, and force and moment data were acquired from Mach 0.8 to 1.8. The sonic boom test data were obtained using a 2-in. flat-top pressure rail and a 14-in. round-top tapered "reflection factor 1" (RF1) pressure rail. Both rails capture an entire pressure signature in one data point, and successive signatures at varying positions along or above the rail were used to improve data quality through spatial averaging. The sonic boom data obtained by the rails were validated with high-fidelity numerical simulations of off-body pressures using the CFD codes USM3D, Cart3D, and OVERFLOW. The test results from the RF1 rail showed good agreement between the computational and experimental data when a variety of testing techniques including spatial averaging of a series of pressure signatures were employed, however, reflections off the 2-in. flat-top rail caused distortions in the signatures that did not agree with the CFD predictions. The 9 x 7 and 8 x 6 wind tunnels generally produced comparable data.

pressure signature↗

Kestrel Results at Liftoff Conditions for a Space Launch System Configuration in Proximity to the Launch Tower

Aerodynamic data books for Space Launch System vehicles require databases for the integrated forces and moments and section loads during liftoff and transition to the ascent phase of flight. While the force and moment database can be generated from wind tunnel results, computational analyses are necessary to provide the extensive surface information required to generate proper lineloads. Of the two flight regimes, the liftoff problem is the more costly and complex situation to simulate, as it requires modeling of the vehicle in proximity to the launch tower. The effects of massive separation on the leeward pressure fields of both the tower and vehicle are not well captured with RANS methods, necessitating the use of more advanced methods, such as Delayed Detached Eddy Simulation, in conjunction with computational grids sufficiently refined to resolve the wakes. Details on the computational setup for employing the Kestrel flow solver to address the liftoff problem are presented. The methodology involves the use of independent unstructured near-body grids for the vehicle and the tower, overset by a solution adaptive Cartesian off-body grid. Results from the simulations are compared to experimental results from a test in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel.

Computational fluid dynamics↗

Edge Based Viscous Method for Node-Centered Formulations

This paper presents a novel, efficient, conservative, edge-based method for evaluation of mean flow viscous fluxes and turbulence-model diffusion terms of the Reynolds-averaged Navier-Stokes equations on tetrahedral grids. The new method is implemented in a practical, node-centered, finite-volume computational fluid dynamics solver. The baseline finite-volume scheme that is equivalent to a second-order accurate finite-element Galerkin approximation of viscous stresses is reformulated. The order of operations to compute the cell-based Green-Gauss gradients is changed to combine the operations by edge, which leads to an equivalent formulation on tetrahedral grids, improves efficiency, and preserves the compact discretization stencil based on the nearest neighbors. The computational results presented in this paper verify the implementation of this edge-based method by comparing its accuracy and iterative convergence with those of the well verified and validated baseline formulation. Efficiency gains for residual and Jacobian evaluations result in significant reduction of time to solution. This novel edge-based formulation on tetrahedra can be seamlessly combined with the baseline formulation on cells of other types for computing solutions on mixed-element grids.

Edge Based↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗

Simulations and case study of X-59 low-booms propagated through measured atmospheric profiles

NASA's X-59 Quiet Supersonic Technology aircraft will soon be used to collect data to support the development of a dose-response relationship between low-boom level and human perception. The X-59's low-boom level will depend on aircraft conditions and trajectory, which can be controlled, and on atmospheric conditions, which cannot be controlled. To assess variability in low-boom levels produced by realistic atmospheres, NASA's PCBoom code was used to simulate propagation of an X-59 nearfield pressure condition through atmospheric profiles measured during NASA's Quiet Supersonic Flights 2018 (QSF18) test. Despite QSF18 lasting only 11 days, substantial weather variability occurred including snow and record high temperatures. A PL range of about 8.5 dB was predicted due to the QSF18 atmospheric variability. These results demonstrate the necessity for X-59's flight condition to be adjusted based on atmospheric conditions in order to achieve desired loudness levels during community surveys. Undertrack booms' Perceived Levels (PL) were predicted not to exceed 75 dB, X-59's target level in a standard atmosphere. Attenuation rate, ray tube area, path length, and other quantities are presented throughout propagation for the atmospheres that produce the loudest and quietest booms. Humidity differences below 15kft were a primary driver of the PL differences.

X-59↗

Summary of the First AIAA Stability and Control Prediction Workshop

Results from the First AIAA Stability and Control Prediction Workshop are summarized in this paper. The workshop series was developed in support of three primary objectives: (1) to establish best practices for the prediction of stability and control derivatives using industry-standard Computational Fluid Dynamics (CFD) solvers, (2) to provide an impartial forum for evaluating the effectiveness of Reynolds-averaged-Navier-Stokes- and Detached-Eddy-Simulation-based modeling techniques, and (3) to identify areas in need of additional research and development. To address these objectives, the inaugural workshop focused on generating computational aerodynamic predictions for the ONERA version of the NASA/Boeing Common Research Model. The configuration includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. While longitudinal wind tunnel test data for the model had been previously documented, unpublished lateral test data at small sideslip angles provided a unique opportunity for participants to generate blind computational predictions for wind tunnel data comparisons. Additional test cases included assessments of the Mach number effect on static lateral/directional stability and the impact of the wind tunnel sting on longitudinal stability. Participants were invited to generate solutions using two workshop-provided series of structured overset and unstructured grids, in addition to any participant custom grids of interest. Total and component-level breakdowns for the force and moment coefficients for each test case are presented, as well as sectional pressure distribution data on the wing and tail components, to assess the agreement between several different Reynolds-averaged Navier-Stokes CFD solvers.

CFD↗

Edge-Based Viscous Method for Mixed-Element Node-Centered Finite-Volume Solvers

A novel, efficient, edge-based viscous (EBV) discretization method has been recently developed, implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver, and applied to viscous-kernel computations that include evaluations of meanflow viscous fluxes, turbulence-model and chemistry-model diffusion terms, and the corresponding Jacobian contributions. Initially, the EBV method had been implemented for tetrahedral grids and demonstrated multifold acceleration of all viscous-kernel computations. This paper presents an extension of the EBV method for mixed-element grids. In addition to the primal edges of a given mixed-element grid, virtual edges are introduced to connect cell nodes that are not connected by a primal edge. The EBV method uses an efficient loop over all (primal and virtual) edges and features a compact discretization stencil based on the nearest neighbors. This study verifies the EBV method and assesses its efficiency on mixed-element grids by comparing the EBV solution accuracy and iterative convergence with those of well-established solutions obtained using a cell-based viscous (CBV) discretization method. The EBV solver’s memory footprint is optimized and often smaller than the memory footprint of the CBV solver. A multifold speedup is demonstrated for all viscous-kernel computations resulting in significant reduction of the time to solutions for several benchmark mixed-element-grid computations, including simulations of a flow around NASA’s juncture-flow model and a hypersonic, chemically reacting flow around a blunt body.

Edge-based viscous method↗

Advanced Parametric Lobe Mixer Concepts for Internally Mixed Nozzles

A continuation of efforts to optimize a dual-stream lobe mixer for supersonic commercial flight are detailed, adding scarfing and scalloping features to a baseline lobe mixer previously optimized within a baseline parametric space. These "advanced features" are parameterized by four quantities: scarfing angle 𝜁, scallop height ℎ 𝑠 , scallop depth 𝑑 𝑠 , and scallop conic focal parameter 𝜌 𝑠 . Two baseline mixers, one parallel and one vortical, form the basis for an eighteen-case study to investigate the effects of each parameter in single-axis fashion. Aerodynamic performance is evaluated through velocity coefficient 𝐶 𝑉 . Moderate performance increases are seen with increasing scallop height and depth on the parallel mixer, and with increasing scarfing angle on the vortical mixer. Flowfield analysis is presented through visualization of axial and cross-stream vorticity development through the duct, as well as density variation. A correlation between axial vorticity and density variation decay was found, with earlier introduction of axial vorticity found to result in improved aerodynamic performance and potential noise suppression.

Lobe Mixer↗