Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite-volume”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the acceleration and separation regions of the flow. An unstructured finite-volume solver is used along with an equilibrium wall model. The Mach number of the oncoming flow selected for the simulation is 0.176, and two Reynolds numbers (Re) are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Spanwise periodic results for the lower Re case are compared with the available DNS data, while the full three-dimensional simulation results for the higher Re case are compared to available experimental data. Sensitivity of WMLES results are assessed for a number of factors including grid resolution, wall model exchange location, type of wall model, and unstructured grid topology. The WMLES results for the lower Re case agree well with available DNS data in terms of the wall pressure variation, velocity and turbulent stress profile comparisons. The skin-friction predictions show a reduced tendency to separate. This appears to be accentuated for the higher Re case, for the medium resolution grid used in this study, where WMLES does not predict any flow separation in the mean, which is in contrast to the large separation bubble observed in the experiments.

turbulence↗
NASA NTRS2020

Experimental Investigation of Shock-Shock Interactions Over a 2D Wedge at M = 6

The effects of fin-leading-edge radius and sweep angle on peak heating rates due to shock-shock interactions were investigated in the NASA Langley Research Center 20-Inch Mach 6 Air Tunnel. The fin model leading edges, which represent cylindrical leading edges or struts on hypersonic vehicles, were varied from 0.25 inches to 0.75 inches in radius. A 9° wedge generated a planar oblique shock at 16.7° to the flow that intersected the fin bow shock, producing a shock-shock interaction that impinged on the fin leading edge. The fin angle of attack was varied from 0° (normal to the free-stream) to 15° and 25° swept forward. Global temperature data were obtained from the surface of the fused silica fins through phosphor thermography. Metal oil flow models with the same geometries as the fused silica models were used to visualize the streamline patterns for each angle of attack. High-speed zoom-schlieren videos were recorded to show the features and temporal unsteadiness of the shock-shock interactions. The temperature data were analyzed using one-dimensional semi-infinite as well as one- and two-dimensional finite-volume methods to determine the proper heat transfer analysis approach to minimize errors from lateral heat conduction due to the presence of strong surface temperature gradients induced by the shock interactions. The general trends in the leading-edge heat transfer behavior were similar for the three shock-shock interactions, respectively, between the test articles with varying leading-edge radius. The dimensional peak heat transfer coefficient augmentation increased with decreasing leading-edge radius. The dimensional peak heat transfer output from the two-dimensional code was about 20% higher than the value from a standard, semi-infinite one-dimensional method.

shock-shock interaction↗
NASA NTRS2021

Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities such as line generation for general unstructured grids and optimization of the iterative solver. The efficiency of the new developments is assessed in a massively parallel environment. The assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Three benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere cylinder configuration, a transonic flow around the ONERA M6 wing, and a supersonic flow through a square duct are considered.

Mohagna J Pandya↗
NASA NTRSDate not supplied

Multiscale Mesh Adaptation for Transonic Aeroelastic Flutter Problems

This work applies multiscale mesh adaptation with refine to reduce spatial discretization error of aeroelastic computational fluid dynamics (CFD) simulations. Benchmark flutter models, such as the pitch and plunge NACA64A-010 airfoil and the benchmark supercritical wing, are studied with both a linearized frequency-domain solver and time-marching CFD coupled to a modal structural solver in FUN3D. The undeformed NASA Common Research Model (CRM), an aeroelastic jig shape variant of the CRM, is also studied with the linearized frequency-domain approach. For these cases, the adaptation process converges to comparable flutter predictions to hand-generated meshes but with smaller node counts. However the additional disciplines of the linearized frequency-domain analysis, the mesh deformation, and the unsteady finite-volume solver create robustness challenges that need to be addressed before it can be applied as a fully automated process for complex transonic aeroelastic problems. In particular, negative volumes are observed to be an issue for FUN3D’s linear elasticity mesh deformation solver when moving the adapted meshes.

Aeroelasticity↗
NASA NTRSDate not supplied

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗
NASA NTRSDate not supplied

Wall-Modeled Large-Eddy Simulations for High-Lift Configurations using FUN3D

Wall-modeled large-eddy simulation (WMLES) capability has recently been implemented into FUN3D, an unstructured, node-centered, finite-volume solver developed at the NASA Langley Research Center. In this paper, WMLES is assessed for two configurations that are representative for high-lift applications. The first configuration is a nominal two-dimensional multielement airfoil that has been extensively studied in the literature. WMLES solutions are computed for four angles of attack and compared with previously reported solutions. Good agreement of integrated forces, surface pressures, and boundary-layer velocity profiles is shown with available experimental data especially at lower angles of attack. WMLES solutions are also computed for the NASA High-Lift Common Research Model over a large range of angles of attack. Forces, pitching moments, and pressure distributions are favorably compared with the experimental data up to the maximum lift, including the angle of attack where the maximum lift is obtained experimentally. Eddy visualization techniques of q-criterion and density-gradient magnitude illustrate the resolved content.

High-Lift flows↗
NASA NTRSDate not supplied

Numerical Investigation of Flow Over Periodic Hill

Hybrid turbulence modeling strategies that blend between large-eddy simulation and Reynolds-averaged Navier-Stokes promise to make optimal use of the available grid resolution. A delayed detached eddy simulation hybrid turbulence modeling strategy, based on the shear-stress transport turbulence model, was implemented into a finite-volume compressible Navier-Stokes code. Then the flow through a two-dimensional periodic hill channel was simulated. This geometry is a widely accepted benchmark test case for validating hybrid turbulence modeling strategies. The computational grid for the simulation was downloaded from a hybrid turbulence model validation website. The present hybrid turbulence model simulation predicts the size of the separated flow region as well as velocity and Reynolds stress profiles in good agreement with large-eddy simulation reference data.

Computational Fluid Dynamics↗
NASA NTRS2022

Monte-Carlo Analysis of Minimum Thermocouple Depths using Icarus

Icarus is a three-dimensional, unstructured, finite-volume material response solver developed at NASA Ames Research Center and has been verified against other NASA material response tools like FIAT, which have a long history of successfully designing thermal protection system (TPS). Icarus solves a set of conservation equations for mass and energy and uses Darcy’s Law in place of momentum conservation. An ecosystem of material response tools has been built around a general-purposed Icarus library that in addition to the typical material response analysis also supports TPS sizing (1-D and multi-dimensional), uncertainty quantification, and has been successfully integrated into a multi-physics architecture built around US3D. In this paper, a brief overview of Icarus and its capabilities will be presented using an illustrative Monte Carlo analysis of the one-dimensional, in-depth material response of a representative Dragonfly trajectory.

Material Response↗
NASA NTRSDate not supplied

Monte-Carlo Analysis of Minimal Thermocouple Depths using Icarus

Icarus is a three-dimensional, unstructured, finite-volume material response solver developed at NASA Ames Research Center \cite{Schulz_2017} and has been verified against other NASA material response tools like FIAT, which have a long history of successfully designing thermal protection system (TPS). Icarus solves a set of conservation equations for mass and energy and uses Darcy’s Law in place of momentum conservation. An ecosystem of material response tools has been built around a general-purposed Icarus library that in addition to the typical material response analysis also supports TPS sizing (1-D and multi-dimensional), uncertainty quantification, and has been successfully integrated into a multi-physics architecture built around US3D \cite{Schroeder_2021}. In this paper, a brief overview of Icarus and its capabilities will be presented using an illustrative Monte Carlo analysis of the one-dimensional, in-depth material response of a representative Dragonfly trajectory.

Material Response↗
NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗
NASA NTRSDate not supplied

Wall-modeled LES of the Three-dimensional Speed Bump Experiment

We evaluate the performance of wall-modeled large eddy simulation (WMLES) in predicting smooth-body turbulent flow separation for a three-dimensional Gaussian-shaped speed bump geometry. The Reynolds number based on the bump length is 2 million, and an unstructured compressible finite-volume solver is used with an equilibrium wall model and dynamic subgrid scale model. Spanwise periodic simulations of the centerline two-dimensional bump were used to assess grid resolution requirements, and resolving the thin internal layer in the accelerating region was found to be necessary to correctly capture the downstream separated flow region. Based on these insights, an optimized grid that was smaller by a factor of two provided comparable accuracy to a finer grid that has been used by us and other researchers in past studies. Using insights gained from the spanwise periodic simulation, an unstructured polyhderal grid with about 250 million cells was used for the three-dimensional (3D) configuration with inviscid tunnel side and top walls. Detailed comparisons of wall skin-friction coefficient, wall pressure, velocity and turbulent stresses with available experimental data indicated excellent agreement. While the low-Reynolds Number Spalart-Allmaras RANS model with rotation/curvature correction gave qualitatively good agreement with experiments, WMLES showed significantly more accurate quantitative predictions of the flowfield. The 3D WMLES showed good agreement with experiments in the separated region, and the centerplane results indicated a different separation topology compared to the spanwise periodic simulation.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRSDate not supplied

On Computationally Efficient Radiative Transfer Calculations for Three-dimensional Entry Problems

The current work presents an efficient simulation framework for rigorously modeling radiative fields emanating from non-equilibrium planetary entry flows in complex three-dimensional domains. Key to this endeavor is adoption of finite-volume discretization in lieu of brute-force ray tracing. This change in conjunction with mesh sweeping and Lebedev-type quadrature for angular integration allows spatial-angular resolution of radiative transfer to be performed in a computationally tractable manner. Additionally, a new methodology has been established to build standalone reduced-order spectral databases for non-equilibrium radiative properties that can be applied to a broad range of hypersonic planetary entry problems with minimal loss in accuracy. The efficacy of the new framework has been demonstrated on the atomic nitrogen radiative system. The resulting reduced-order model requires three orders-of-magnitude fewer spectral evaluations and results in a 95% decrease in memory footprint. A comparative study for representative forebody and afterbody lines-of-sight from Stardust, FIRE II, and meteor entries into the Earth atmosphere indicates that both total intensity variation and detailed spectra can be retrieved with as few as 625 reduced-order groups (contrasting with the 100,000 frequencies in the original full set model). Similarly, three-dimensional predictions of radiative heating experienced by the Orion forebody are in excellent agreement with legacy radiation solvers while requiring only a sliver (roughly 0.5%) of computing wall time.

3D radiation↗
NASA NTRS2020

Implicit Large-Eddy Simulations of Hot and Cold Supersonic Jets in Loci-CHEM

Turbulence is a critical factor in most aero- and propulsion-related flows. For many years, the industry-standard approach to incorporating the effect of turbulence in computational fluid dynamics (CFD)calculations has been the Reynolds Averaged Navier Stokes (RANS) modeling approach. RANS implicitly time-averages turbulent motion, and models its effect on the mixing of species, momentum, and energy in a flow field. The approach is economical and often allows a CFD calculation to proceed to steady-state. However, it makes several simplifying assumptions, and the results will only be as good as the RANS turbulence model. Meanwhile, the last two decades have witnessed rapid and explosive growth in computational power, the development of CFD codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. CFD engineering calculations of 108– 109 finite-volume cells have become routine. Thus, large eddy simulation (LES), which attempts to directly resolve the unsteady motion of the largest scales of turbulence, is increasingly of interest for many fluids engineering problems. Additionally, hybrid RANS-LES approaches such as Spalart’s Detached Eddy Simulation (DES), which attempt to marry the best strengths of both methods, are in active use and development. Because it seeks to accurately resolve the unsteady motion of the largest scales of turbulence, LES places certain requirements on the numerical method used to discretize and advance the Navier Stokes equations. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide range of length scales.

C I Morris↗
NASA NTRSDate not supplied

Martian B Storm Evolution: Modeling Dust Activity over the Receding South Polar CO 2 Ice Cap at Southern Hemisphere Summer Solstice

Observations of annually recurring, elevated (25 km) atmospheric warming and associated dust concentrations over the south polar cap near Southern Hemisphere (SH) summer solstice are the defining features of the regional dust storm known as the B storm [1]. Our observational analysis of MCS and TES temperature and dust retrievals include all Mars Years (MY) lacking a global dust storm (GDS) between MY 26 and MY 33 (Figure 1). The data indicate the B storm is initiated around perihelion (L s =252°) in the southern midlatitudes and reaches maximum height and intensity over the south pole near SH summer solstice at L s ~267°. At peak intensity, mid-level atmospheric temperatures are elevated by as much as 30 K and mid-level atmospheric dust mixing ratios are more than triple pre-B storm levels (1.5-2 ppm to 7-8 ppm). Cap edge lifting may contribute significantly to initial dust mobilization in the B storm since surface wind stress along the receding cap edge is likely strengthened by the combined effects of the sea breeze circulation and katabatic winds at this time of year [2, 3, 4]. We are investigating dust lifting and lofting near the south pole of Mars during SH summer using the NASA Ames Mars Global Climate Model (MGCM). The NASA Ames MGCM is a finite-volume numerical model with a cubed sphere grid ideal for modeling processes at high latitudes. Initial simulations at low resolution (4x4 degrees) capture some elevated dust and associated warming at the south pole around L s =265°, but the model fails to capture the intensity of the B storm. We will examine and present the model-predicted surface stress and vertical velocities near the receding seasonal CO 2 cap edge and at the south pole during solstice. Finally, we will present higher resolution simulations to better capture small-scale dynamics and features at the pole.

Courtney M. L. Batterson↗
NASA NTRSDate not supplied

Aerothermal Shape Optimization of Actively-Cooled Battery Packs using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

thermal management↗
NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗
DOE OSTI2025

SCM overview & the EBR-II shutdown heat removal tests validation.

Pronghorn is an engineering-scale, coarse-mesh, thermal-hydraulics tool for supporting reactor-core simulations of advanced nuclear reactors. Most of the current efforts in Pronghorn have been devoted in developing porous finite-volume capabilities and adapting closure correlations for coarse-mesh thermal-hydraulics modeling. However, for liquid-metal reactors (LMRs) with wire-wrapped fuel pin assemblies, a pin-level thermal-hydraulic resolution is required for most safety case studies (pin rupture, channel blockage, etc.). For this purpose, a new Subchannel application is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This presentation presents the main features of the SCM code and demonstrates a validation case based on the EBR-II SHRT tests.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗