Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scale resolving 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 55 records · Page 3

Wall Modeled Large Eddy Simulations for NASA’s Jet Noise Consensus Database of Single-Flow, Round, Convergent Jets

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics(LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction(PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology(CST) project. The goal of PUR is to quantify and reduce uncertainties from scale-resolving simulations to assess noise characteristics of next generation quiet supersonic commercial jets during takeoff and landing conditions where the noise from the exhaust jet dominates. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several conditions both for static (no ambient co-flow), which is the focus of this article, and in-flight (flight stream co-flow) conditions, which are beyond the current scope. Nine different flow conditions ranging in jet exit Mach number from 0.38 to 1.0 with nozzle temperature ratios (NTR) from 0.84 to 2.7 have been conducted. Details of the structured overset grids, numerical discretization and wall-model are provided. Near-field comparisons to PIV data show great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the lip line shear-layer. A permeable Ffowcs Williams Hawkings (FWH) surface, enclosing the jet, is used to predict far-field noise from the simulated flow-field. Comparison of CFD predictions to microphone array measurements demonstrate excellent agreement within the resolved frequency range. A systematic under-prediction of far-aft observer angles larger than 150 degrees has been observed across all simulations. We achieved a cost reduction of an order of magnitude for these WMLES compared to an earlier study of this configuration due to algorithmic and software improvements. The accuracy of the results and short turnaround time demonstrate that WMLES within the LAVA framework is a cost-effective approach for jet noise predictions that could soon be incorporated into the design cycle of jet noise reduction technologies.

CST↗

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

Exploring the Accuracy of RANS Simulations for Mars Entry Vehicles

Accurate yet inexpensive predictions of aerodynamic coefficients for Mars entry vehicles have remained a consistent challenge over the past five decades. Below Mach 6, drag on the backshell becomes significant and must be accurately predicted. Steady Reynolds-averaged Navier-Stokes (RANS) models are commonly used, despite their poor predictions of backshell pressure. While scale-resolving simulations have shown promise in the past decade, there is still a need for cheap, accurate RANS predictions for large aerodynamic databases. The Mars Science Laboratory (MSL) is used as a case study to examine predictive accuracy and known shortcomings for RANS predictions of Mars entry vehicles. Several different grid generation techniques are compared, including a comparison between prismatic boundary layer grids and fully unstructured, tetrahedral grids. Comparisons are made to experimental data for Mach 2.5, 3.5, and 4.5. The accuracy of predicted aerodynamic coefficients is examined. Overpredictions in axial force and drag are explained by a closer examination of the surface pressure. These findings document sensitivities and best practices for future RANS database development of Mars entry vehicles.

RANS↗

Exploring the Accuracy of RANS Simulations for Mars Entry Vehicles

Accurate yet inexpensive predictions of aerodynamic coefficients for Mars entry vehicles have remained a consistent challenge over the past five decades. Below Mach 6, drag on the backshell becomes significant and must be accurately predicted. Steady Reynolds-averaged Navier-Stokes (RANS) models are commonly used, despite their poor predictions of backshell pressure. While scale-resolving simulations have shown promise in the past decade, there is still a need for cheap, accurate RANS predictions for large aerodynamic databases. The Mars Science Laboratory (MSL) is used as a case study to examine predictive accuracy and known shortcomings for RANS predictions of Mars entry vehicles. Several different grid generation techniques are compared, including a comparison between prismatic boundary layer grids and fully unstructured, tetrahedral grids. Comparisons are made to experimental data for Mach 2.5, 3.5, and 4.5. The accuracy of predicted aerodynamic coefficients is examined. Overpredictions in axial force and drag are explained by a closer examination of the surface pressure. These findings document sensitivities and best practices for future RANS database development of Mars entry vehicles.

CFD↗

Implicit Thermochemical Nonequilibrium Flow Simulations on Unstructured Grids using GPUs

Thermochemical nonequilibrium flow simulation capabilities have been previously implemented, verified, and validated for central processing unit (CPU) systems in NASA’s unstructured-grid computational fluid dynamics solver FUN3D. Many exascale-class high-performance computing systems will rely on graphics processing unit (GPU) architectures for high throughput and energy efficiency; thus, CPU-based scientific computing software unable to effectively utilize these systems must be updated. In this work, we present a CUDA C++ implementation of FUN3D’s thermochemical nonequilibrium flow simulation capabilities targeting NVIDIA Tesla GPUs. An overview of the porting and optimization strategy is described and performance comparisons with other recent architectures are presented. Scaling to thousands of GPUs is demonstrated, yielding computational performance equivalent to that of several million CPU cores. The implementation enables efficient, high-fidelity, scale-resolving simulations of thermochemical nonequilibrium flows for many applications including atmospheric entry, hypersonics, and combustion.

GPU↗

Implicit Thermochemical Nonequilibrium Compressible Flow Simulations on Unstructured Grids Using GPUs

As next-generation exascale-class systems arrive, existing software must be updated accordingly to effectively utilize these systems. For high concurrency and energy efficiency, many of these systems utilize GPU architectures. In this work, we present a CUDA C++ implementation of FUN3D's thermochemical nonequilibrium capability for turbulent flows. Efficiency is demonstrated at scale using the Summit system at the Oak Ridge Leadership Computing Facility which is representative of future exascale systems. This work enables faster, higher fidelity, and scale-resolving simulations of thermochemical nonequilibrium flows including reentry, hypersonics, and combustion.

CFD, GPU, HPC, Hypersonics, Chemistry↗

The First Stars: A Low-Mass Formation Mode

We perform numerical simulations of the growth of a Population III stellar system under photodissociating feedback. We start from cosmological initial conditions at z = 100, self-consistently following the formation of a minihalo at z = 15 and the subsequent collapse of its central gas to high densities. The simulations resolve scales as small as approx. 1 AU, corresponding to gas densities of 10(exp 16)/cu cm. Using sink particles to represent the growing protostars, we evolve the stellar system for the next 5000 yr. We find that this emerging stellar group accretes at an unusually low rate compared with minihalos which form at earlier times (z = 20-30), or with lower baryonic angular momentum. The stars in this unusual system will likely reach masses ranging from <1Stellar Mass to approx. 5 Stellar Mass by the end of their main-sequence lifetimes, placing them in the mass range for which stars will undergo an asymptotic giant branch (AGB) phase. Based upon the simulation, we predict the rare existence of Population III stars that have survived to the present day and have been enriched by mass overflow from a previous AGB companion.

Low-Mass↗

Aeroacoustic Computations of a Transonic Truss-Braced Wing Aircraft: Part 1 – Aero-dynamic and Airframe Noise Simulations

Computational results are presented for an 8%-scale, full-span, transonic truss-braced wing (TTBW) model simulated as installed in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel and in free-air conditions. The simulations were conducted with the lat-tice Boltzmann solver PowerFLOW® to capture the time-accurate characteristics of the flow. The aerodynamic behavior of the aircraft was investigated in the landing configuration, with high-lift devices and landing gear deployed, as well as in the clean (cruise) configuration, with these components stowed. Analyses were performed on local flow quantities, global forces, and time-averaged surface pressures. Aerodynamic quantities were shown to be sensitive to mesh resolution levels, driven by small geometric features inherent to the TTBW model. Flow features of the TTBW model were examined, with the wing/strut configuration of this model presenting unique behaviors generally not found in conventional transport aircraft. Near-field, time-dependent flow quantities obtained from the scale-resolving simulations were used in conjunction with a Ffowcs-Williams and Hawkings integral approach to predict the far-field airframe noise signature of this advanced concept. The effects of permeable data sur-face end caps on the far-field noise spectrum in the flyover direction were determined to be negligible.

Transonic Truss-Braced Wing↗

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗

A New Recycling Method to Generate Turbulent Inflow Profiles

The accuracy of the scale-resolving simulations for practical geometries strongly depends on the inflow boundary conditions. Imposing experimentally observed turbulent inflow profiles for the numerical simulations is a major challenge. Existing methods available in the literature assume self-similar behavior, which is not true for most of the experiments. In the present work, we formulate the turbulent inflow profile generation technique as an optimization problem. An adjoint technique is exploited to evaluate the sensitivities of multiple input parameters for the present problem. The present formulation is then tested to generate a laminar boundary layer profile, turbulent boundary layer profile, and turbulent jet profile.

Garai, Anirban↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

Direct Numerical Simulations↗

Structure and Dynamics of Bleed-Controlled Impinging Shock/Turbulent-Boundary-Layer Interactions

Porous boundary-layer bleed is often employed to alleviate the detrimental effects of shock/turbulent-boundary-layer interactions (STBLI) in high-speed vehicles. However, outstanding questions remain on the underlying mechanisms through which bleed influences STBLI dynamics. This work addresses this knowledge gap by performing a scale-resolved Large-Eddy Simulation of an impinging STBLI controlled by discrete bleed holes. The simulation has a freestream Mach number of 2.5, an oblique shock generator angle of 8 degrees, and a Reynolds number based on the incoming momentum thickness of 2800. Bleed substantially reduces the separation size and upstream influence of the STBLI compared to the canonical, uncontrolled case. Consequently, the characteristic, low-frequency motion of STBLI-induced separation is mitigated, as revealed by wall-pressure spectra. Energy is shifted from the low frequency separation behavior to a higher-frequency signature caused by the bleed shocks and expansions. The scale-resolved simulations also reveal that these small, localized structures severely limit the turbulence production traditionally seen in the uncontrolled interaction.

computational fluid dynamics↗

Direct simulation of compressible wall-bounded turbulence

When analyzing many turbulent flows, the effects of compressibility can be neglected. Even some relatively high-speed flows, such as boundary layers generated by a supersonic aircraft, produce turbulent statistics that are similar to those found for the incompressible case. There are other situations, however, in which the non-zero divergence of the turbulence leads to behavior that is fundamentally different from that found at constant density. Examples include flows created by internal combustion engines, hypersonic flight, and supersonic combustion. It is with instances such as these that this project is concerned. In particular, we are interested in the effects of compressibility on turbulence near a smooth solid constant temperature surface; our primary objective is an increased physical understanding that can be used to improve turbulence models of wall-bounded compressible flows. With this in mind, we have begun a direct numerical simulation (DNS) study of turbulence in a plane channel. Because all of the relevant spatial and temporal scales are to be resolved, the simulations require no subgrid scale parameterization. The DNS code developed by Buell to study compressible plane Couette flow was modified to solve the compressible Navier-Stokes equations in the plane channel. The channel was chosen over the Couette flow for two reasons: (1) to avoid using the very large streamwise domains needed to adequately capture the large Couette vortical structures; and (2) to make use of previous experience by considering the compressible version of a well established case and isolate finite Mach number effects by comparing to the incompressible channel.

Coleman, Gary N.↗