Search NASA⌕ Search

SEARCH · Search NASA

Results for “RANS simulations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The myth of URANS

Since the 1990s, RANS practitioners have observed spontaneous unsteadiness in RANS simulations. Some have suggested deliberately using this as a method of resolving large turbulent structures. However, to date, no one has produced a theoretical justification for this unsteady RANS (URANS) approach. Here, we extend the dynamical system fixed point analysis to create a theoretical model for URANS dynamics. The results are compared to URANS simulations for homogeneous isotropic decaying turbulence. The model shows that URANS can predict incorrect decay rates and that the solution tends towards steady RANS over time. Similar analysis for forced turbulence shows a fixed modelled energy of about 30% of total energy, regardless of the model parameters. The same analysis can be used to show how hybrid type models can begin to address these issues.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zero-pressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient flat plate case.

Transition Model↗
NASA NTRSDate not supplied

Transonic Lift and Drag Predictions using Wall Modelled Large Eddy Simulations

Wall-modelled Large Eddy Simulations (WMLES) of the NASA Common Research Model (CRM) at transonic conditions and various angles of attacks leading up to and including shock-induced flow separation are performed using the LAVA computational framework. The simulations are shown to accurately predict the lift curve slope and the onset of separation characterized by the break in the pitching moment. Furthermore, careful assessment of skin-friction drag at cruise condition and its subsequent decrease with increasing angles of attack are shown to be in agreement with viscous sublayer resolving Reynolds Averaged Navier Stokes(RANS) simulations. The small differences between WMLES and RANS appear to be of the same order as differences seen between two RANS models at the cruise-point. Some sensitivity is observed to the coefficient used in the subgrid scale model, although this can be reconciled by noting the low chord- and shock-incidence Reynolds numbers (Rex≈106) seen in the outboard regions of the wing, along with the uncertainties associated with tripping and the numerical transition that occurs near the leading edge. Among the primary drawbacks seen in the WMLES predictions is the subdued increase in wave-drag with increasing angle-of-attack and the shock intensity when compared with experimental data. Preliminary assessment suggests that additional span- and stream-wise grid resolution is likely needed in mid- and out-board portions of the wing to better resolve the shock-induced separated flow and to further investigate prediction accuracy of unsteady temporal characteristics of the problem.

ARMD↗
DOE OSTI2026

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↗
NASA NTRS1978

Real time digital propulsion system simulation for manned flight simulators

A real time digital simulation of a STOL propulsion system was developed which generates significant dynamics and internal variables needed to evaluate system performance and aircraft interactions using manned flight simulators. The simulation ran at a real-to-execution time ratio of 8.8. The model was used in a piloted NASA flight simulator program to evaluate the simulation technique and the propulsion system digital control. The simulation is described and results shown. Limited results of the flight simulation program are also presented.

Mihaloew, J. R.↗
NASA NTRSDate not supplied

Computational Analysis of a Boundary-Layer Ingesting Tailcone Thruster Configuration Within the National Transonic Facility Using the LAVA Solver

Boundary-layer ingesting (BLI) propulsion systems are one of the many technologies currently under investigation within the aerospace community to meet NASA’s Advanced Air Transport Technology project goal of a sustainable future in aviation. To that end, an experimental campaign was conducted in the National Transonic Facility (NTF) to study the propulsion-airframe integration effects of the Boundary-Layer Ingesting Tailcone System installed in the aft portion of a 2.7% scale design of the NASA Common Research Model (CRM). This test was accompanied by a numerical simulation effort using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) solver, to validate the applicability of current best-practice Reynolds-averaged Navier Stokes (RANS) models in accurately predicting the relevant flow-physics in free-air. A total of 205 steady RANS simulations were performed using LAVA, covering the range of conditions present in the NTF test matrix. Flow conditions ranged from 5- to 15-million Reynolds number, Mach numbers of 0.75, 0.80 and 0.85, and angles-of-attack between –3° and 4°. Four different mass-flow plugs were also tested to assess propulsor operating condition effects. Sensitivity to these flow conditions are analyzed in detail, especially in terms of the nacelle flow distortion which was the main subject of the study. Aerodynamic load coefficients, aftbody boundary-layer measurements and fuselage pressure tap results are also discussed, providing a wide range of validation results from this experimental campaign. This validation effort shows that current best-practices using RANS models within the LAVA solver flow solver are well suited for predicting the inlet flow distortion characteristics of novel aircraft configurations employing BLI at the aft end of the fuselage.

AATT↗
NASA NTRSDate not supplied

Computational Analysis of a Boundary Layer Ingesting Tailcone Thruster Configuration Using the LAVA Curvilinear Solver

Boundary layer ingesting (BLI) propulsion systems are one of the many technologies currently under investigation within the aerospace community to meet NASA's Advanced Air Transport Technology project goal of a sustainable future in aviation. To that end, an experimental campaign was conducted in the National Transonic Facility (NTF) wind tunnel to study the propulsion-airframe integration effects of the Boundary Layer Ingesting Tailcone System installed in the aft portion of a 2.7% scale design of the NASA Common Research Model (CRM). This test was accompanied by a numerical simulation effort using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) curvilinear solver to validate the applicability of current best-practice Reynolds-Averaged Navier-Stokes (RANS) models in accurately predicting the relevant flow physics in free-air. A total of 205 steady RANS simulations were performed using LAVA, covering the range of conditions present in the NTF test matrix. Flow conditions ranged from 5- to 15-million Reynolds number, Mach numbers of 0.75, 0.80 and 0.85, and angles of attack between -3 and 4 degrees. Four different mass-flow plugs were also tested to assess propulsor operating condition effects. Sensitivity to these flow conditions is analyzed in detail, especially in terms of the nacelle flow distortion which was the main subject of the study. Aerodynamic load coefficients, aft body boundary layer measurements and fuselage pressure tap results are also discussed, providing a wide range of validation results from this experimental campaign. This validation effort shows that current best-practices using RANS models within the LAVA curvilinear solver flow solver are well suited for predicting the inlet flow distortion characteristics of novel aircraft configurations employing BLI at the aft end of the fuselage. This talk will cover the combined experimental and numerical efforts that resulted from this BLI-focused investigation.

ARMD↗
DOE OSTI2024

Effects of boundary hydraulics, dissolved oxygen, and dissolved organic carbon on growth and death dynamics of aerobic microbes in riverbed dune-induced hyporheic zones

Surface and groundwater interact in the hyporheic zone beneath and adjacent to rivers in the presence of a diverse microbial community. Heterotrophic bacteria mediate a range of environmentally important reactions, yet few studies have quantified bacterial growth and death dynamics in the hyporheic zone, and none have systematically analyzed their response to variations in hydraulic or chemical conditions. Here we used MODFLOW and SEAM3D to simulate hydraulics; dissolved oxygen (DO) and dissolved organic carbon (DOC) transport; and aerobic microbial metabolism, growth, and death in hyporheic zones induced by riverbed dunes. We ran simulations both with and without growth/death processes, and varied hydraulic parameters and DO/DOC boundary concentrations. Microbial biomass reached steady state (t = 3 days) in every simulation, at which time there was greater biomass and DOC biodegradation rates in the hyporheic flowcell (300% and 85% higher for the base case, respectively) when accounting for microbial growth dynamics. This occurred as microbial biomass tailored its spatial distribution to the availability of DO and DOC, demonstrating the importance of simulating growth/death processes. Biomass generally increased with hyporheic flow cell area as upwelling groundwater decreased. When varying surface water DO and DOC source concentrations relative to the base case, the greatest effect on biomass occurred when increasing DOC and decreasing DO. We determined minimum DO and DOC steady-state concentrations required for microbial growth, but the minimums were not absolute or related by stoichiometry. Increasing DOC created a smaller area of microbes with higher concentrations relative to the base case. Increasing DO slightly increased the area occupied by microbes while keeping the total biomass nearly constant. Overall, microbial growth and death dynamics depend on DO and DOC availability in the hyporheic zone, which is dependent on DOC/DO boundary concentrations and hyporheic flow paths, and in turn the hydraulic interaction between surface water and groundwater.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Verification of Anisotropic Mesh Adaptation for Turbulent Simulations over ONERA M6 Wing

Unstructured anisotropic mesh adaptation is known to be an efficient way to control discretization errors in Computational Fluid Dynamics (CFD) simulations. Method verification is required to provide the confidence for routine use in production analysis. The current work aims at verification of anisotropic mesh adaptation for RANS simulations over the ONERA M6 wing. The present verification study is performed using four different flow solvers, three different implementations of the metric field, and three mesh mechanics packages. Two of the flow solvers use stabilized finite-element discretizations (FUN3D-SFE and GGNS), one uses finite-volume discretization (FUN3D-FV), and the last one uses mixed finite-volume and finite element discretizations (Wolf). The mesh adaptation is based on an error estimator that aims to control the quadratic error term in the linear interpolation of Mach number. Two sets of adaptations were performed; the first one controls the interpolation error in L2 norm and the second one controls the interpolation error in L4 norm. Convergence studies were performed on the forces and the pitching moment using all four solvers, and the results are compared with previously verified convergence studies on fixed (nonadapted) meshes. Both forces and pitching moment on adapted meshes are found to be converging to the fine mesh values faster than those on fixed meshes. In addition to forces and moments, convergence of surface pressure and skin friction coefficients at various measurement locations on the wing are also presented. Adapted-mesh surface pressure distributions agree with the fine fixed mesh pressure distributions. Adapted-mesh skin friction distributions contain high frequency noise with mean values approaching the fixed mesh pressure skin friction distributions.

Aravind Balan↗
DOE OSTI2024

CFD simulations of Molten Salt Fast Reactor core cavity flows

Computational Fluid Dynamics (CFD) has become increasingly important in the research and development of advanced nuclear reactors. Here, in the current study, extensive CFD simulations were conducted for the coolant flow in Molten Salt Fast Reactor (MSFR) core models using the state-of-the-art spectral element flow solver Nek5000 and multiscale coarse-mesh thermal-hydraulic software Pronghorn. The underlying motivation is to seek an in-depth understanding of how the internal velocity distribution can be influenced by the MSFR core cavity shape, the Reynolds number, turbulence modeling options and the inlet boundary conditions. The CFD techniques involved in this investigation range from coarse-mesh CFD, RANS modeling, to the high-fidelity LES calculations. Specifically, a series of RANS simulations were performed for the 2-D axisymmetric core model and 3-D wedge domains to study the flow distribution inside the MSFR core. It is observed that a proper representation of the MSFR inlet channel duct is important for the prediction of internal flow distribution. It is also showcased here how researchers can leverage the Nek5000 CFD results to calibrate more efficient coarse-mesh CFD tools, like Pronghorn, for the actual MSFR design needs. Moreover, this paper highlights a 3-D LES model for an entire MSFR core using the spectral element method and demonstrates the feasibility of this modeling approach. The readiness and potential limitations of the RANS approach are examined with respect to the high-fidelity LES simulations. The present investigation lays a solid foundation as we are leveraging the high-fidelity CFD capabilities to inform MSFR design efforts.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

A General Relativistic Magnetohydrodynamic Simulation of Jet Formation

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation ofjet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity approx.0.3c) is created, as shown by previous two-dimensional axi- symmetric simulations with mirror symmetry at the equator. The three-dimensional simulation ran over 100 light crossing time units (T(sub s) = r(sub s)/c, where r(sub s = 2GM/c(sup 2), which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted owing in part to magnetic pressure from the twisting of the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r(sub s). At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface ofthe thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outward with a wider angle than the initial jet. The widening of the jet is consistent with the outward-moving torsional Alfven waves. This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk because of the iack of streaming materiai from an accompanying star.

Nishikawa, K.-I.↗
NASA NTRS2023

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

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 airfoil 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 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985).1The 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

Transonic Buffet↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

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 airfoil 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 million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). 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.

Transonic Buffet↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

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 airfoil 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 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). 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 span wise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗
DOE OSTI2022

SST ..kappa..-..omega.. Simulations of the Atmospheric Boundary Layer Including the Coriolis Effect

For large-scale structures in the atmospheric boundary layer (ABL), the Coriolis effect and near-wall behavior can have a meaningful impact. For example, both the Coriolis effect and blade boundary layer impact how muchpower wind farms produce. RANS simulations of the ABL typically use the ..kappa..-epsilon turbulence model, which has been developed to accurately capture the Coriolis effect but typically does not perform well near walls. The SST ..kappa..-..omega.. turbulence model performs well near walls but does not accurately capture the Coriolis effect. This work modifies SST ..kappa..-..omega.. to accurately model the Coriolis effect. We discuss the similarities and differences in how to modify ..kappa..-epsilon and SST ..kappa..-..omega.. for the Coriolis effect. Finally, we compare ..kappa..-epsilon and SST ..kappa..-..omega.. simulations of the ABL, including the Coriolis effect, with large eddy simulations and measurements.

atmospheric boundary layer↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2017

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations (LES) were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver was used to compute the flow from a 2-inch round nozzle. Three different flow conditions of varying jet Mach numbers and temperature ratios were examined. The LES results showed that the temperature field behaves similar to the velocity field, but with a more rapidly spreading mixing layer. Predictions of mean, mu-bar(sub i), and fluctuating, mu'(sub i), velocities were compared to particle image velocimetry data. Predictions of mean, T-bar, and fluctuating, T', temperature were compared to data obtained using Rayleigh scattering and Raman spectroscopy. Very good agreement with experimental data was demonstrated for the mean and fluctuating velocities. The LES correctly predicts the behavior of the turbulent temperature field, but over-predicts the levels of the fluctuations. The turbulent heat flux was examined and compared to Reynolds-averaged Navier-Stokes (RANS) results. The LES and RANS simulations produced very similar results for the radial heat flux. However, the axial heat flux obtained from the LES differed significantly from the RANS result in both structure and magnitude, indicating that the gradient diffusion type model in RANS is inadequate. Finally, the LES data was used to compute the turbulent Prandtl number and verify that a constant value of 0.7 used in the RANS models is a reasonable assumption.

hot jets↗