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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Development Status of a 3-D Electron Fluid Model for Hall Thruster Plume Simulations

A 3-D electron fluid model has been developed as a steping stone to fully describe the electron current flow across magnetic fields inside a vacuum chamber and to provide electron flux to solar arrays for spacecraft surface charging model. A detailed description of the numerical treatment of the electric potential solver, including finite-volume formulation, implementation, and the treatment of boundary conditions, are presented in this paper. Verification tests of the model are presented.

Choi, Maria↗
NASA NTRS2018

Unstructured Grid Adaptation and Solver Technology for Turbulent Flows

Unstructured grid adaptation is a tool to control Computational Fluid Dynamics (CFD) discretization error. However, adaptive grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Issues that prevent the use of adaptive grid methods are identified by applying unstructured grid adaptation methods to a series of benchmark cases. Once identified, these challenges to existing adaptive workflows can be addressed. Unstructured grid adaptation is evaluated for test cases described on the Turbulence Modeling Resource (TMR) web site, which documents uniform grid refinement of multiple schemes. The cases are turbulent flow over a Hemisphere Cylinder and an ONERA M6Wing. Adaptive grid force and moment trajectories are shown for three integrated grid adaptation processes with Mach interpolation control and output error based metrics. The integrated grid adaptation process with a finite element (FE) discretization produced results consistent with uniform grid refinement of fixed grids. The integrated grid adaptation processes with finite volume schemes were slower to converge to the reference solution than the FE method. Metric conformity is documented on grid/metric snapshots for five grid adaptation mechanics implementations. These tools produce anisotropic boundary conforming grids requested by the adaptation process.

Park, Michael A.↗
NASA NTRS2018

Development Status of a Three-Dimensional Electron Fluid Model for Hall Thruster Plume Simulations

A 3-D electron fluid model has been developed as a stepping stone to fully describe the electron current flow across magnetic fields inside a vacuum chamber and to provide electron flux to solar arrays for spacecraft surface charging model. A detailed description of the numerical treatment of the electric potential solver, including finite-volume formulation, implementation, and the treatment of boundary conditions, are presented in this paper. Verification tests of the model are presented.

Choi, Maria↗
NASA NTRS2018

Development Status of a 3-D Electron Fluid Model for Hall Thruster Plume Simulations

A 3-D electron fluid model has been developed as a steping stone to fully describe the electron current flow across magnetic fields inside a vacuum chamber and to provide electron flux to solar arrays for spacecraft surface charging model. A detailed description of the numerical treatment of the electric potential solver, including finite-volume formulation, implementation, and the treatment of boundary conditions, are presented in this paper. Verification tests of the model are presented.

Choi, Maria↗
NASA NTRS2018

Material Response Analysis of a Titan Entry Heatshield

Accurate calculation of thermal protection material response is critical to the vehicle design for missions to the Saturn moon Titan. In this study, Icarus, a three-dimensional, unstructured, finite-volume material response solver under active development at NASA Ames Research Center, is used to compute the in-depth material response of the Huygens spacecraft along its November 11 entry trajectory. The heatshield analyzed in this study consists of a five-layer stack-up of Phenolic Impregnated Carbon Ablator (PICA), aluminum honeycomb, adhesive, and face sheetmaterials. During planetary entry, the PICA outer layer is expected to undergo pyrolysis. A surface energy balance boundary condition that captures both time- and spatial-variance of surface properties during entry is used in the simulation.

material response↗
NASA NTRS2018

High-Resolution Modeling of the Dust and Water Cycles with the NASA Ames Mars Global Climate Model

NASA’s Mars Climate Modeling Center at Ames Research Center is currently undergoing an exciting period of growth in personnel, modeling capabilities, and science productivity. We are transitioning from our legacy Arakawa C-grid finite-difference dynamical core to the NOAA/GFDL cubed-sphere finite-volume dynamical core for simulating the climate of Mars in a global framework. This highly parallelized core is scalable and flexible, which allows for significant improvements in the horizontal and vertical resolutions of our simulations. We have implemented the Ames water ice cloud microphysics package described in Haberle et al. (2018) into this new dynamical core. We will present high-resolution simulations of the dust and water cycles that show that sub-degree horizontal resolution improves the agreement between the vertical distribution of dust and water ice and observations. In particular, both water ice clouds and dust are transported to higher altitudes due to stronger topographic circulations at high resolution. Preliminary results suggest that high-resolution global modeling is needed to properly capture critical features of the dust and water cycles, and thus the current Mars climate.

Kahre, Melinda A.↗
NASA NTRS2018

P43K-3884 Characterizing Martian Crater Circulations with the NASA Ames Mars GCM

Observations made in Gale Crater by instruments on the MSL Curiosity Rover show that the diurnal amplitude of the surface pressure is increased and the depth of the Convective Boundary Layer (CBL) is decreased relative to other lander locations on flatter regions of Mars (Haberle et al., 2014; Moores et al., 2015). Mesoscale modeling studies of Gale Crater suggest that crater circulations produce these effects. Tyler & Barnes (2013) show that local upslope/downslope flows along the crater rim and Mt. Sharp amplify the diurnal pressure cycle. These same flows are thought to be at least partly responsible for the suppression of the CBL because upward air flow at the rim and in the center (due to Mt. Sharp) forces subsidence over the lowest regions of the crater during the day. Regional flows, largely due to the location of Gale near the dichotomy boundary, may also play a role in shaping the circulation internal to the crater. Whether the behavior of the CBL and the amplified diurnal pressure cycle are phenomena observed in craters morphologically different from Gale (i.e. bowl-shaped, irregular, degraded) is not yet understood. We will explore these questions by characterizing the behavior of these processes as they are shaped by the morphology of craters greater than 100 km in diameter. We use the NASA Ames Mars Global Circulation Model (GCM) that now utilizes the NOAA/GFDL cubed-sphere finite-volume dynamical core to examine ~100 craters of varying size and shape from a database of known Martian craters (Robbins & Hynek, 2014). Run at 7.5 km resolution, the GCM is capable of resolving surface winds, temperature, and pressure inside craters of this size allowing for the analysis of dozens of craters simulated at various seasons and within the context of synoptic and global-scale phenomena.

Batterson, Courtney MaryLou↗
NASA NTRS2019

Numerical Modeling of Boiling in a Heated Tube

This paper presents numerical models of boiling in a heated tube using the Generalized Fluid System Simulation Program (GFSSP), a finite-volume-based general-purpose flow network code developed at NASA/Marshall Space Flight Center. The heated tube is discretized into a one-dimensional array of nodes and branches to represent the flow of liquid and vapor in a tube with a prescribed pressure differential. The solid wall is also discretized into solid nodes and conductors to allow for heat transfer between the wall and the fluid. The conservation equations of mass, momentum, and energy of the fluid are solved simultaneously with the energy conservation equation for the solid wall. Two experimental configurations of fluid flowing in a vertical tube have been simulated, one with water and the other with liquid hydrogen. This paper compares experimental data with numerical predictions based on four different published correlations for boiling heat transfer coefficients. Three of these correlations are applicable to the saturated vertical flow conditions of the experiments. One of them is applicable to film boiling and has been used for the liquid hydrogen experiment, which was in film boiling regime. For the case of boiling water, the predictions of wall temperatures using the boiling heat transfer correlations agreed well with the experimental results. However, in the case of boiling hydrogen larger discrepancies were observed between the experimental data and numerical predictions.

Majumdar, Alok↗
NASA NTRS2019

Numerical Predictions of Dust-Induced Heat Flux Augmentation in Hypersonic Blunt-Body Flows Using a Discontinuous Galerkin Multiphase Flow Solver

Recent interest in human-scale missions to Mars has motivated the need for high-fidelity simulations of reentry flows. During a dust storm, there can be high levels of suspended dust in the Martian atmosphere, which cannot only enhance erosion of thermal protection systems but also transfer energy and momentum to the shock layer, thereby significantly augmenting the surface heat flux. Second-order finite-volume schemes are typically employed for hypersonic flow simulations, but such schemes suffer from a number of disadvantages. An attractive alternative is discontinuous Galerkin methods, which benefit from arbitrarily high spatial order of accuracy, geometric flexibility, and other properties. To enable accurate computations of high-speed particle-laden flows, an Euler-Lagrange methodology was developed in which the Eulerian field of the carrier gas is calculated using a discontinuous Galerkin scheme while the disperse phase is treated with Lagrangian particle tracking. We discuss challenges associated with coupling these two formulations and how to handle them. Momentum and energy transfer between the carrier gas and the particle phase is considered, and the importance of accounting for interparticle collisions is assessed. In addition, we describe the physical model of the particle phase and examine effects of its uncertainties on the numerical solution. We demonstrate the performance of the Euler-Lagrange method in representative testcases, with focus on the accurate prediction of particle trajectories and heating augmentation. Quantitative comparisons with experiments are provided.

Ching, Eric J.↗
NASA NTRS2019

Free-Flight CFD Simulations of Ballistic Range and Flight Tests

Free-Flight CFD capability has been implemented into the finite-volume solver US3D under the Entry Systems Modeling project. Several simulations of ballistic range experiments have been performed in order to validate the simulation software and methodology. Extension of the software to flight scale trajectories with varying freestream conditions has been carried out. Results show promising ability to predict vehicle behavior as compared to flight. Finally, a multi-body free-flight capability has been developed to generalize the single-body free-flight solver to study multiple bodies in proximal flight.

Brock, Joseph M.↗
NASA NTRS2019

Accuracy, Scalability, and 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 including parallelization, line generation for general unstructured grids, improved discretization scheme, and optimized iterative solver. The paper reports on the new developments to the flow solver and assesses the accuracy, scalability, and efficiency. The USM3D assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Two benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere-cylinder configuration and a transonic flow around the ONERA M6 wing are considered.

Pandya, Mohagna J.↗
NASA NTRSDate not supplied

Verification of Anisotropic Mesh Adaptation for Complex Aerospace Applications

The stabilized finite element solver, FUN3D-SFE, along with the grid mechanics package refine are verified for aerospace applications of laminar and turbulent flow simulations. The current verification exercise represents an extension of previous research using FUN3D-SFEwith adjoint-based mesh adaptation to generate highly anisotropic adapted meshes for inviscid problems. Adaptations are performed using a solution-based approach that controls the Lpnorm of Mach number interpolation error and an adjoint-based approach that controls the error in some output functional. Adaptive results are shown for laminar subsonic flow over a delta wing, laminar subsonic flow over ONERA M6 wing, inviscid supersonic flow over a sonic boom test case, and a turbulent flow over a high lift configuration (JAXA Standard Model). Mesh convergence results are also compared with results from FUN3D-FV (Finite Volume) whenever available. For all test cases considered, FUN3D-SFE gives significantly better accurate results than FUN3D-FV on coarse meshes.

Aravind Balan↗
NASA NTRSDate not supplied

Wall-modeled LES of the NASA Juncture Flow Experiment

We perform wall-modeled large eddy simulation (WMLES) of the NASA Juncture Flow experiment (Kegerise & Neuhart, NASA/TM–2019–20286) using an equilibrium wall model and unstructured finite volume solver, to assess its predictive ability for complex geometries. The flow condition simulated corresponds to 5 degrees angle of attack with a freestream Mach number of 0.189, and a Reynolds number based on the mean crank chord of 2.4 million. All simulations neglect the effects of wind tunnel walls, sting and mast present in the experiment. To assess the performance of WMLES on simpler flows with very coarse grids, we simulate flow in a turbulent channel at a friction Reynolds number, Reτ≈2000, and flow through a square duct at anReτ≈1000. The duct flow results using the equilibrium wall model indicate that the skin friction is inaccurate near the corner for grid topologies and resolutions typically used in WMLES, having implications for more complicated flows such as the juncture flow. For the juncture flow simulation, two different approaches are investigated. The first approach consists of a truncated-domain simulation wherein the inflow plane of the computational domain is placed at≈0.2chord length of the wing, and we prescribe the mean flow from a separate Reynolds-averaged Navier-Stokes (RANS) solution along with synthetic turbulence to initiate realistic unsteadiness in the domain. The second approach involves simulating the entire geometry with trips to trigger transition to turbulence. The truncated- and full-domain simulations contain about 62 and 90 million cells, respectively; with 8-10 points per boundary layer thickness in the wall-parallel directions, 14-20 points in the wall-normal direction and a near-wall viscous spacing (∆n+1) of≈100. Preliminary results are encouraging overall in terms of the prediction of wall pressure, wall skin friction, velocity and stresses; but indicate that further work is required to improve the predictions in the separation bubble and wing-fuselage corner regions. The grid resolution used in this study is still fairly coarse, and the results should be interpreted as work in progress.

Prahladh S Iyer↗
NASA NTRSDate not supplied

Investigation of Improved Wall Heat Flux Behavior Observed When Computing Hypersonic Flow Using a Face-averaged Nodal-gradient Approach On Tetrahedral Grids

We investigate the origin of a significant reduction in computed numerical “noise” in the spatial variation of computed wall heat flux observed when using a node-centered weighted least squares(WLSQ) gradient approach when solving the Reynolds-averaged Navier-Stokes (RANS) equations for hypersonic boundary layer flows on highly stretched tetrahedral grids. This investigation is conducted using a 2nd-order, cell-centered, finite-volume, discretization which we employ to contrast results obtained using two WLSQ gradient approaches: a conventional 3-D, cell-centered, node-neighbor, WLSQ gradient (NN-CCG) approach and a 3-D, face-averaged, node-centered gradient (F-ANG)approach. The F-ANG approach, is shown to produce significantly less “noise” in the spatial variation of wall heat flux, than the NN-CCG approach and we demonstrate the cause of this “noise” via a series of numerical experiments where we compute hypersonic laminar and turbulent boundary layer flows.We further investigate the level of grid resolution required for a tetrahedral grid to approach the fidelity of a grid resolved hexahedral grid solution obtained using an equivalent structured grid method.

Jeffery A. White↗
NASA NTRSDate not supplied

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar↗
NASA NTRSDate not supplied

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar↗
NASA NTRSDate not supplied

Validation and Sensitivity Analyses of Arc-Jet Performance of Woven Thermal Protection Entry Systems

Woven thermal protection materials for the Adaptable, Deployable, Entry and Placement Technology (ADEPT) entry system are simulated using a dual layer formalism in a finite volume implicit framework. The model developed and demonstrated here includes multi-dimensional and time-varying thermal response of the material to arc-jet test conditions. This work further demonstrates the sensitivities of the macro-scale thermal response to the through-thickness and in-plane conductivities of the carbon fabric, and to other experimental measurements such as cloth emissivity and thickness. The modeling can be extended to the ADEPT system that would encompass a single-piece woven heat shield configuration

Pratibha Raghunandan↗
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↗