Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes simulation”

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 253 records · Page 14

USM3D-ME Contributions to the 5th AIAA High Lift Prediction Workshop

This paper presents the results of Reynolds-averaged Navier-Stokes (RANS) simulations conducted by NASA’s flow solver, mixed-element USM3D (USM3D-ME), for the 5th AIAA High-Lift Prediction Workshop. As part of the Fixed-Grid RANS Technology Focus Group (TFG), these simulations were performed to assess the accuracy and efficiency of the USM3D-ME solutions in predicting high-lift flows. The High-Lift Common Research Model (CRM-HL) served as the primary geometry. Several CRM-HL configurations were used for three case studies: a verification study (Case 1), a configuration buildup study (Case 2), and a Reynolds-number variation study (Case 3). Overall, USM3D-ME RANS results aligned with the solutions selected by the Fixed-Grid RANS TFG and available wind tunnel data Simulations for Cases 1 and Configuration 2.1 achieved machine-zero residual convergence, with aerodynamic coefficients converging to steady-state values. However, Configurations 2.2-2.4 and Case 3 encountered iterative- and grid-convergence challenges, particularly at high angles of attack. Compared with the experimental data available for Configurations 2.2-2.4, close agreement was demonstrated at low angles of attack. However, for angles of attack approaching the maximum lift conditions, the predicted lift coefficient and pitching moment deviated from experimental values. The drag-coefficient predictions were in a relatively good agreement, however, slight overpredictions were observed at the highest angle of attack corresponding to the maximum-lift condition. Although iterative convergence for Configurations 2.2-2.4 at high angles of attack remains a persistent challenge, averaging aerodynamic coefficients over the last 5000 iterations yielded satisfactory agreement with the available wind tunnel experimental data. During the workshop, the lack of iterative convergence was attributed to the vortex structures emanating from the slat brackets. To investigate this issue further, post-workshop simulations were conducted on Configuration 2.2. In one study, RANS simulations were performed on a simplified geometry with the slat brackets removed. The second study focused on performing URANS simulations on the original Configuration 2.2 geometry. Preliminary results from both studies are presented and compared with wind tunnel data for Configuration 2.2. Consistent with the findings of other participants in the Fixed-Grid RANS TFG, this study emphasizes the necessity for further exploration and advancement in RANS technology for predicting high-lift flows.

CFD↗

USM3D-ME Contributions to the 5th AIAA High Lift Prediction Workshop

This paper presents the results of Reynolds-averaged Navier-Stokes (RANS) simulations conducted by the NASA flow solver, mixed-element USM3D (USM3D-ME), for the 5th AIAA High-Lift Prediction Workshop. As part of the Fixed-Grid RANS Technology Focus Group (TFG), these simulations were performed to assess the accuracy and efficiency of the USM3D-ME solutions in predicting high-lift flows. The High-Lift Common Research Model (CRM-HL) served as the primary geometry. Several CRM-HL configurations were used for three case studies: a verification study (Case 1), a configuration buildup study (Case 2), and a Reynolds-number variation study (Case 3). Overall, USM3D-ME RANS results aligned with the solutions selected by the Fixed-Grid RANS TFG and available wind tunnel data. Simulations for Case 1 and Configuration 2.1 achieved machine-zero residual convergence, with aerodynamic coefficients converging to steady-state values. However, Configurations 2.2-2.4 and Case 3 encountered iterative- and grid-convergence challenges, particularly at high angles of attack. Compared with the experimental data available for Configurations 2.2-2.4, close agreement was demonstrated at low angles of attack. However, for angles of attack approaching the maximum lift conditions, the predicted lift coefficient and pitching moment deviated from experimental values. The drag-coefficient predictions were in relatively good agreement, however, slight overpredictions were observed at the highest angle of attack corresponding to the maximum-lift condition. Although iterative convergence for Configurations 2.2-2.4 at high angles of attack remains a persistent challenge, averaging aerodynamic coefficients over the last 5000 iterations yielded satisfactory agreement with the available wind tunnel experimental data. During the workshop, the lack of iterative convergence was attributed to the vortex structures emanating from the slat brackets. To investigate this issue further, post-workshop simulations were conducted on Configuration 2.2. In one study, RANS simulations were performed on a simplified geometry with the slat brackets removed. The second study focused on performing unsteady RANS (URANS) simulations on the original Configuration 2.2 geometry. Preliminary results from both studies are presented and compared with wind tunnel data for Configuration 2.2. Consistent with the findings of other participants in the Fixed-Grid RANS TFG, this study emphasizes the necessity for further exploration and advancement in RANS technology for predicting high-lift flows.

Aerodynamics↗

Using 3D Measurements from MRI to Improve VVUQ

Magnetic resonance velocimetry (MRV) is a diagnostic that can measure 3D, three-component turbulent velocity fields in arbitrarily complex flow configurations. The approach uses magnetic resonance imaging (MRI), which is a system commonly available in radiology departments or medical research centers. MRV acquires signal from hydrogen protons in water channel flows. Despite the obvious utility of the measurements for investigating a variety of flows, comparisons with computational fluid dynamic (CFD) simulations largely retain traditional metrics, comparing velocity profiles and planes at discrete regions within the flow rather than using the 3D nature of the measured field. In this effort, MRV was conducted in a simple water channel with six spanwise-centered periodic obstacles roughly shaped as cubes. At a fully turbulent Reynolds number of 15,000, the channel includes two flows: a main streamwise flow and, in the wake of the second obstacle, an injected flow oriented perpendicular to the streamwise flow. The flow geometry includes partial obstacles on the side walls and is inherently 3D because of the interaction of the two streams and the wake features, making the flow challenging to completely measure with traditional optical techniques. Aside from line profiles and planar comparisons, two 3D metrics are used to compare the data with a steady Reynolds-averaged Navier–Stokes (RANS) simulation result. A brief discussion about the comparison is provided, including comments about uncertainty.

Benson, Mike [ORNL] (ORCID:000000023210116X)↗

Development and assessment of models for turbulent Rayleigh-Taylor mixing using the macroscopic forcing method

Reynolds-Averaged Navier Stokes (RANS) simulations are a popular method for designing ICF experiments, and accurate mixing models are crucial for these simulations to give good predictions. To this end, the present work seeks to demonstrate the Macroscopic Forcing Method (MFM) as a tool for both improving existing RANS models as well as assessing RANS model forms. First, MFM analysis from Lavacot et al. (Phys. Rev. Fluids, 2025) is used to develop the k–L–F model, an extension of the k–L model of Dimonte and Tipton (Phys. Fluids, 2006) that incorporates nonlocality through addition of a turbulent species flux transport equation. MFM is then applied to the k–L–F model along with the k–L and BHR–4 models to assess their forms and compare the model-implied eddy diffusivity moments to those measured from high-fidelity simulations. Furthermore, the analysis reveals that models incorporating nonlocality (k–L–F and BHR–4) match the high-fidelity simulation data better than purely local models (k–L), both in terms of mean fields and eddy diffusivity moments. However, all of the considered RANS models struggle to match temporal moments at high Atwood numbers, highlighting the importance of temporal nonlocality in these regimes and the need for additional improvement even among models incorporating nonlocality.

general physics↗

Three dimensional hypersonic flow simulations with the CSCM implicit upwind Navier-Stokes method

The Bardina and Lombard (1985) three-dimensional CSCM Navier-Stokes method is presently extended to the simulation of complex hypersonic reentry vehicle external flows at angle of attack. The robust stability of the method derives from the combination of conservative implicit upwind flux difference splitting with a three-dimensional diagonally-dominant approximate factorization and relaxation scheme and characteristic-based implicit boundary approximations. The method's efficiency derives from an implicit symmetric Gauss-Seidel 'method of planes' relaxation scheme with alternating directional space marching sweeps along the flow coordinate direction.

Bardina, Jorge↗

Transonic Flows About A Fighter Airplane

Report describes computational study of transonic flows about airplane having shape similar to F-16A fighter. Navier-Stokes calculations simulate complicated flow fields and interactions.

Reznick, Steve↗

Two-equation turbulence modeling for 3-D hypersonic flows

An investigation to verify, incorporate and develop two-equation turbulence models for three-dimensional high speed flows is presented. The current design effort of hypersonic vehicles has led to an intensive study of turbulence models for compressible hypersonic flows. This research complements an extensive review of experimental data and the current development of 2D turbulence models. The review of experimental data on 2D and 3D flows includes complex hypersonic flows with pressure profiles, skin friction, wall heat transfer, and turbulence statistics data. In a parallel effort, turbulence models for high speed flows have been tested against flat plate boundary layers, and are being tested against the 2D database. In the present paper, we present the results of 3D Navier-Stokes numerical simulations with an improved k-omega two-equation turbulence model against experimental data and empirical correlations of an adiabatic flat plate boundary layer, a cold wall flat plate boundary layer, and a 3D database flow, the interaction of an oblique shock wave and a thick turbulent boundary layer with a free stream Mach number = 8.18 and Reynolds number = 5 x 10 to the 6th.

Bardina, J. E.↗

Viscous hypersonic shock-shock interaction on a blunt body at high altitude

The shock interaction produced when an incident oblique shock impinges on a blunt body, such as an engine inlet cowl lip of a hypersonic vehicle, is investigated for high altitude flight conditions. A perfect gas, Navier-Stokes numerical simulation of this problem at various altitudes representing continuum through transitional conditions is performed using the modified flux vector splitting method of Steger and Warming (1979). Two series of solutions are produced. First, a number-of-shock-positions are studied at a particular altitude and Mach number. Second, given a fixed shock position and Mach number, the interaction is investigated at several altitudes ranging from continuum to transitional flow conditions. It is shown that the interaction becomes fundamentally different as the density is lowered, and its effect on the overheating problem is progressively diminished. The maximum stagnation point heating at the highest altitude is reached only when the incident shock misses the cowl lip completely, and any interaction with the cowl bow shock that does occur takes place downstream and thus has little effect on the conditions at the stagnation point.

Comeaux, Keith A.↗

Compressor Stability Enhancement Using Discrete Tip Injection

Mass injection upstream of the tip of a high-speed axial compressor rotor is a stability enhancement approach known to be effective in suppressing small in tip-critical rotors. This process is examined in a transonic axial compressor rotor through experiments and time-averaged Navier-Stokes CFD simulations. Measurements and simulations for discrete injection are presented for a range of injection rates and distributions of injectors around the annulus. The simulations indicate that tip injection increases stability by unloading the rotor tip and that increasing injection velocity improves the effectiveness of tip injection. For the tested rotor, experimental results demonstrate that at 70 percent speed the stalling flow coefficient can be reduced by 30 percent using an injected mass- flow equivalent to 1 percent of the annulus flow. At design speed, the stalling flow coefficient was reduced by 6 percent using an injected mass-fiow equivalent to 2 percent of the annulus flow. The experiments show that stability enhancement is related to the mass-averaged axial velocity at the tip. For a given injected mass-flow, the mass-averaged axial velocity at the tip is increased by injecting flow over discrete portions of the circumference as opposed to full-annular injection. The implications of these results on the design of recirculating casing treatments and other methods to enhance stability will be discussed.

Suder, Kenneth L.↗

Development of a Flexible Framework of Common Hypersonic Navier-Strokes Meshes for the Space Shuttle Orbiter

A flexible framework for the development of block structured volume grids for hypersonic Navier-Stokes flow simulations was developed for analysis of the Shuttle Orbiter Columbia. The development of the flexible framework, resulted in an ability to quickly generate meshes to directly correlate solutions contributed by participating groups on a common surface mesh, providing confidence for the extension of the envelope of solutions and damage scenarios. The framework draws on the experience of NASA Langely and NASA Ames Research Centers in structured grid generation, and consists of a grid generation process that is implemented through a division of responsibilities. The nominal division of labor consisted of NASA Johnson Space Center coordinating the damage scenarios to be analyzed by the Aerothermodynamics Columbia Accident Investigation (CAI) team, Ames developing the surface grids that described the computational volume about the orbiter, and Langely improving grid quality of Ames generated data and constructing the final volume grids. Distributing the work among the participants in the Aerothermodynamic CIA team resulted in significantly less time required to construct complete meshes than possible by any individual participant. The approach demonstrated that the One-NASA grid generation team could sustain the demand for new meshes to explore new damage scenarios within a aggressive timeline.

Alter, Stephen J.↗

Studying Turbulence Using Numerical Simulation Databases - X Proceedings of the 2004 Summer Program

This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.

Moin, Parviz↗

Opposed Jet Burner Approach for Characterizing Flameholding Potentials of Hydrocarbon Scramjet Fuels

Opposed Jet Burner (OJB) tools have been used extensively by the authors to measure Flame Strength (FS) extinction limits of laminar H2/N2 air and (recently) hydrocarbon (HC) air Counterflow Diffusion Flames (CFDFs) at one atm. This paper details normalization of FSs of N2- diluted H2 and HC systems to account for effects of fuel composition, temperature, pressure, jet diameter, inflow Reynolds number, and inflow velocity profile (plug, contoured nozzle; and parabolic, straight tube). Normalized results exemplify a sensitive accurate means of validating, globally, reduced chemical kinetic models at approx. 1 atm and the relatively low temperatures approximating the loss of non-premixed idealized flameholding, e.g., in scramjet combustors. Laminar FS is defined locally as maximum air input velocity, U(sub air), that sustains combustion of a counter-jet of g-fuel at extinction. It uniquely characterizes a fuel. And global axial strain rate at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or (sub t)) can be compared directly with computed extinction limits, determined using either a 1-D Navier Stokes stream-function solution, using detailed transport and finite rate chemistry, or (better yet) a detailed 2-D Navier Stokes numerical simulation. The experimental results define an idealized flameholding reactivity scale that shows wide ranging (50 x) normalized FS s for various vaporized-liquid and gaseous HCs, including, in ascending order: JP-10, methane, JP-7, n-heptane, n-butane, propane, ethane, and ethylene. Results from H2 air produce a unique and exceptionally strong flame that agree within approx. 1% of a recent 2-D numerically simulated FS for a 3 mm tube-OJB. Thus we suggest that experimental FS s and/or FS ratios, for various neat and blended HCs w/ and w/o additives, offer accurate global tests of chemical kinetic models at the Ts and Ps of extinction. In conclusion, we argue the FS approach is more direct and fundamental, for assessing, e.g., idealized scramjet flameholding potentials, than measurements of laminar burning velocity or blowout in a Perfectly Stirred Reactor, because the latter characterize premixed combustion in the absence of aerodynamic strain. And FS directly measures a chemical kinetic characteristic of non-premixed combustion at typical flameholding temperatures. It mimics conditions where gfuels are typically injected into a subsonic flameholding recirculation zone that captures air, where the effects of aerodynamic strain and associated multi-component diffusion become important.

Pellett, Gerald L.↗

Aerothermodynamic Environments Definition for the Mars Science Laboratory Entry Capsule

An overview of the aerothermodynamic environments definition status is presented for the Mars Science Laboratory entry vehicle. The environments are based on Navier-Stokes flowfield simulations on a candidate aeroshell geometry and worst-case entry heating trajectories. Uncertainties for the flowfield predictions are based primarily on available ground data since Mars flight data are scarce. The forebody aerothermodynamics analysis focuses on boundary layer transition and turbulent heating augmentation. Turbulent transition is expected prior to peak heating, a first for Mars entry, resulting in augmented heat flux and shear stress at the same heatshield location. Afterbody computations are also shown with and without interference effects of reaction control system thruster plumes. Including uncertainties, analysis predicts that the heatshield may experience peaks of 225 W/sq cm for turbulent heat flux, 0.32 atm for stagnation pressure, and 400 Pa for turbulent shear stress. The afterbody heat flux without thruster plume interference is predicted to be 7 W/sq cm on the backshell and 10 W/sq cm on the parachute cover. If the reaction control jets are fired near peak dynamic pressure, the heat flux at localized areas could reach as high as 76 W/sq cm on the backshell and 38 W/sq cm on the parachute cover, including uncertainties. The final flight environments used for hardware design will be updated for any changes in the aeroshell configuration, heating design trajectories, or uncertainties.

Edquist, Karl T.↗

Navier-Stokes Analysis of a High Wing Transport High-Lift Configuration with Externally Blown Flaps

Insights and lessons learned from the aerodynamic analysis of the High Wing Transport (HWT) high-lift configuration are presented. Three-dimensional Navier-Stokes CFD simulations using the OVERFLOW flow solver are compared with high Reynolds test data obtained in the NASA Ames 12 Foot Pressure Wind Tunnel (PWT) facility. Computational analysis of the baseline HWT high-lift configuration with and without Externally Blown Flap (EBF) jet effects is highlighted. Several additional aerodynamic investigations, such as nacelle strake effectiveness and wake vortex studies, are presented. Technical capabilities and shortcomings of the computational method are discussed and summarized.

Slotnick, Jeffrey P.↗

Aerothermodynamic Design of the Mars Science Laboratory Heatshield

Aerothermodynamic design environments are presented for the Mars Science Laboratory entry capsule heatshield. The design conditions are based on Navier-Stokes flowfield simulations on shallow (maximum total heat load) and steep (maximum heat flux, shear stress, and pressure) entry trajectories from a 2009 launch. Boundary layer transition is expected prior to peak heat flux, a first for Mars entry, and the heatshield environments were defined for a fully-turbulent heat pulse. The effects of distributed surface roughness on turbulent heat flux and shear stress peaks are included using empirical correlations. Additional biases and uncertainties are based on computational model comparisons with experimental data and sensitivity studies. The peak design conditions are 197 W/sq cm for heat flux, 471 Pa for shear stress, 0.371 Earth atm for pressure, and 5477 J/sq cm for total heat load. Time-varying conditions at fixed heatshield locations were generated for thermal protection system analysis and flight instrumentation development. Finally, the aerothermodynamic effects of delaying launch until 2011 are previewed.

Edquist, Karl T.↗

Aerothermodynamic Design of the Mars Science Laboratory Backshell and Parachute Cone

Aerothermodynamic design environments are presented for the Mars Science Laboratory entry capsule backshell and parachute cone. The design conditions are based on Navier-Stokes flowfield simulations on shallow (maximum total heat load) and steep (maximum heat flux) design entry trajectories from a 2009 launch. Transient interference effects from reaction control system thruster plumes were included in the design environments when necessary. The limiting backshell design heating conditions of 6.3 W/sq cm for heat flux and 377 J/sq cm for total heat load are not influenced by thruster firings. Similarly, the thrusters do not affect the parachute cover lid design environments (13 W/sq cm and 499 J/sq cm). If thruster jet firings occur near peak dynamic pressure, they will augment the design environments at the interface between the backshell and parachute cone (7 W/sq cm and 174 J/sq cm). Localized heat fluxes are higher near the thruster fairing during jet firings, but these areas did not require additional thermal protection material. Finally, heating bump factors were developed for antenna radomes on the parachute cone

Edquist, Karl T.↗

Inlet Flow Control and Prediction Technologies for Embedded Propulsion Systems

Fail-safe, hybrid, flow control (HFC) is a promising technology for meeting high-speed cruise efficiency, low-noise signature, and reduced fuel-burn goals for future, Hybrid-Wing-Body (HWB) aircraft with embedded engines. This report details the development of HFC technology that enables improved inlet performance in HWB vehicles with highly integrated inlets and embedded engines without adversely affecting vehicle performance. In addition, new test techniques for evaluating Boundary-Layer-Ingesting (BLI)-inlet flow-control technologies developed and demonstrated through this program are documented, including the ability to generate a BLI-like inlet-entrance flow in a direct-connect, wind-tunnel facility, as well as, the use of D-optimal, statistically designed experiments to optimize test efficiency and enable interpretation of results. Validated improvements in numerical analysis tools and methods accomplished through this program are also documented, including Reynolds-Averaged Navier-Stokes CFD simulations of steady-state flow physics for baseline, BLI-inlet diffuser flow, as well as, that created by flow-control devices. Finally, numerical methods were employed in a ground-breaking attempt to directly simulate dynamic distortion. The advances in inlet technologies and prediction tools will help to meet and exceed "N+2" project goals for future HWB aircraft.

McMillan, Michelle L.↗

Computational Analyses of the LIMX TBCC Inlet High-Speed Flowpath

Reynolds-Averaged Navier-Stokes (RANS) simulations were performed for the high-speed flowpath and isolator of a dual-flowpath Turbine-Based Combined-Cycle (TBCC) inlet using the Wind-US code. The RANS simulations were performed in preparation for the Large-scale Inlet for Mode Transition (LIMX) model tests in the NASA Glenn Research Center (GRC) 10- by 10-ft Supersonic Wind Tunnel. The LIMX inlet has a low-speed flowpath that is coupled to a turbine engine and a high-speed flowpath designed to be coupled to a Dual-Mode Scramjet (DMSJ) combustor. These RANS simulations were conducted at a simulated freestream Mach number of 4.0, which is the nominal Mach number for the planned wind tunnel testing with the LIMX model. For the simulation results presented in this paper, the back pressure, cowl angles, and freestream Mach number were each varied to assess the performance and robustness of the high-speed inlet and isolator. Under simulated wind tunnel conditions at maximum inlet mass flow rates, the high-speed flowpath pressure rise was found to be greater than a factor of four. Furthermore, at a simulated freestream Mach number of 4.0, the high-speed flowpath and isolator showed stability for freestream Mach number that drops 0.1 Mach below the design point. The RANS simulations indicate the yet-untested highspeed inlet and isolator flowpath should operate as designed. The RANS simulation results also provided important insight to researchers as they developed test plans for the LIMX experiment in GRC s 10- by 10-ft Supersonic Wind Tunnel.

Dippold, Vance F., III↗