Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerics”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

The Importance of the Cathode Plume and Its Interactions with the Ion Beam in Numerical Simulations of Hall Thrusters

Hall2De is a first-principles, 2-D axisymmetric code that solves the equations of motion for ions, electrons, and neutrals on a magnetic-field-aligned grid. The computational domain downstream of the acceleration channel exit plane is large enough to include self-consistently the cathode boundary. In this paper, we present results from numerical simulations of the H6 laboratory thruster with an internally mounted cathode, with the aim of highlighting the importance of properly accounting for the interactions between the ion beam and cathode plume. The anomalous transport of electrons across magnetic field lines in Hall2De is modelled using an anomalous collision frequency, νanom, yielding νanom approximately equal to omega ce (i.e., the electron cyclotron frequency) in the plume. We first show that restricting the anomalous collision frequency to only regions where the current density of ions is large does not alter the plasma discharge in the Hall thruster as long as the interaction between the ion beam and the cathode plume is captured properly in the computational domain. This implies that the boundary conditions must be placed sufficiently far as to not interfere with the electron transport in this region. These simulation results suggest that electron transport across magnetic field lines occurs largely inside the beam and may be driven by the interactions between beam ions and electrons. A second finding that puts in relevance the importance of including the cathode plume in numerical simulations is on the significance of accounting for the ion acoustic turbulence (IAT), now known to occur in the vicinity of the cathode exit. We have included in the Hall2De simulations a model of the IAT-driven anomalous collision frequency based on Sagdeev's model for saturation of the ion-acoustic instability. This implementation has allowed us to achieve excellent agreement with experimental measurements in the near plume obtained during the operation of the H6 thruster at nominal conditions (300V, 20A) and chamber background pressure of approximately 1.5 x 10(exp -5) Torr. In addition, the numerical results obtained with the latter approach exhibit less sensitivity to background pressure than previous attempts at explaining the features of the plasma properties in the near plume.

Lopez Ortega, Alejandro↗
NASA NTRS2017

Modeling of Electromagnetic Scattering by Discrete and Discretely Heterogeneous Random Media by Using Numerically Exact Solutions of the Maxwell Equations

In this paper, we discuss some aspects of numerical modeling of electromagnetic scattering by discrete random medium by using numerically exact solutions of the macroscopic Maxwell equations. Typical examples of such media are clouds of interstellar dust, clouds of interplanetary dust in the Solar system, dusty atmospheres of comets, particulate planetary rings, clouds in planetary atmospheres, aerosol particles with numerous inclusions and so on. Our study is based on the results of extensive computations of different characteristics of electromagnetic scattering obtained by using the superposition T-matrix method which represents a direct computer solver of the macroscopic Maxwell equations for an arbitrary multisphere configuration. As a result, in particular, we clarify the range of applicability of the low-density theories of radiative transfer and coherent backscattering as well as of widely used effective-medium approximations.

Dlugach, Janna M.↗
NASA NTRS2017

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗
NASA NTRS2018

Numerical Simulation of a Simplified High-Lift CRM Configuration Embedded with Fluidic Actuators

Numerical simulations have been performed for a simplified high-lift configuration that is representative of a modern transport airplane. This configuration includes a leading-edge slat, fuselage, wing, nacelle-pylon and a simple hinged flap. The suction surface of the flap is embedded with multiple rows of fluidic actuators to reduce the extent of reversed flow regions and improve the aerodynamic performance of the configuration with flap in a deployed state. In the current paper, a Lattice Boltzmann Method based high-fidelity computational fluid dynamics (CFD) code, known as PowerFLOW® is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW® code that has been validated for high speed flows is used for the present simulations to accurately represent the transonic flow regimes that are encountered in the flow field generated by the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient active flow control (AFC) systems and actuator arrangement for lift augmentation are emerging based on the parametric studies conducted here prior to wind tunnel tests. These numerical solutions will be compared with experimental data, once such data becomes available.

Vatsa, Veer N.↗
NASA NTRS2015

Numerical Simulations of Shock/Plume Interaction Using Structured Overset Grids

Computational simulations using structured overset grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting oblique shock/plume interaction effects to near-field sonic boom signatures. Standard second-order accurate as well as higher-resolution numerical discretizations are utilized and compared in the study. The numerical approach is compared with supersonic wind-tunnel data for three cases. The cases include an empty wind-tunnel at the operating conditions, an isolated shockgenerating diamond wedge within the tunnel, and a nozzle with diamond wedge configuration at five different nozzle pressure ratios. Solution sensitivity to numerical discretization is analyzed. Favorable comparisons between the computational results and experimental data of near-field pressure signatures are obtained. A simple prediction method for plume induced shock deflection is developed and results are compared with the CFD data.

Housman, Jeffrey A.↗
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 NTRS2017

A new concept of stability in orbit propagation, useful for quantifying numerical errors

We present the concept of topological stability in the numerical propagation of orbits, and show how it results in a useful new method for measuring the global numerical error of an orbit propagation. The concept applies to any problem in orbital dynamics. Moreover, it can be extended to any three-dimensional system of di erential equations of second order. In order to assess the topological stability of a given integration a special metric is introduced, which can be used to estimate the numerical errors robustly. The method is particularly well suited for dealing with strongly perturbed and chaotic systems. The construction is based on the constraint imposed by the Hopf map that supports the Kustaanheimo-Stiefel transformation. Generic concepts of stability are translated to KS space.

Pelaez, Jesus↗
NASA NTRSDate not supplied

Numerical Simulation of the High-Speed Leg of the National Transonic Facility

Numerical simulations for the flow inside the high-speed leg of the National Transonic Facility was conducted. The NASA Tetrahedral Unstructured Software System (TetrUSS) with its USM3D_ME solver was used to perform the numerical simulations. USM3D_ME is developed and maintained by NASA Langley Research Center. Simulations were conducted for three configurations: empty tunnel, body of revolution installed, and NASA Common Research Model installed in the test section. Simulations were performed for a test section Mach number of 0.7 and 0.85 and a corresponding Reynolds number of 8 million per ft. A controller was developed that automated dynamic outflow boundary and streamlined the process of running multiple simulations. The use of a dynamic outflow boundary was the key parameter to drive simulation to the desired tunnel conditions. The numerical simulations captured expected flow features in NTF test section and in the tunnel plenum. The simulations depicted that separation of the flow inside the diffuser appears to be asymmetric and more extensive toward the top and bottom walls. Analysis of the flow field were conducted, and the computed drag coefficient compared to wind tunnel data.

NTF↗
NASA NTRSDate not supplied

Assessment of Numerical and Modeling Errors of RANS based Transition Models for Low-Reynolds Numbers 2-D Flows

In this paper we report the outcome of selected workshops organized as part of the NATO Applied Vehicle Technology (AVT)-313 activity Incompressible Laminar-to-Turbulent Flow Transition Study that focused on assessing the numerical and modeling accuracy of the γ−Reθ and γ transition models coupled to the k−ω Shear-Stress Transport (SST) two-equation eddy-viscosity model. Three different test cases involving nominally 2D flow configurations were selected: flow over a flat plate with two different levels of turbulence intensity at the inlet; flow around the Eppler 387 foil at a Reynolds number of 3×10^5 and angles of attack of 1 deg. and 7 deg. flow around the NACA 0015 foil at a Reynolds number of 1.8×10^5 and angles of attack of 5 deg. and10 deg. The flat plate flow conditions correspond to natural and by-pass transition, whereas the other two test cases include laminar separation bubbles that lead to separation-induced transition. For each test case, the selected quantities of interest include both integral and local flow quantities. Geometrically similar grids with a wide range of grid refinement ratios were generated for each of the test cases to allow the estimation of numerical uncertainties for all quantities of interest selected for this study. Several RANS flow solvers were used, employing common grids with the same boundary conditions and mathematical models. Therefore, it is possible to analyze the consistency of the results, i.e., to check if the intervals defined by the different numerical solutions with their respective uncertainties overlap with each other. Modeling errors can also be addressed for the selected flow quantities that have experimental data available. However, the experimental information available in these cases is not sufficient to guarantee that experiments and simulations are performed with the same settings. Nonetheless, the available experimental data is sufficient to guarantee that modeling errors are significantly reduced with the use of the transition models when compared to simulations performed using only the k−ω SST model.

CFD Modeling↗
NASA NTRSDate not supplied

Graphical and Numerical Description of Strain-Gage Balance Interactions

A new approach for the graphical and numerical description of balance interactions is presented. The approach uses data from single-component loads as input. This choice has the two advantages. First, the number of applied loads is at the minimum needed for interactions to be observed. In addition, the loads used for the description of interactions can easily be repeated at different sites. Output differences relative to the outputs of the zero load point of a load series are used for the description of interactions. Similarly, load differences relative to the loads of the zero load point of a load series are used for the description of the loads. Interactions are plotted versus the load differences for each load component while omitting the outputs of the primary gage of the chosen load component. The resulting plots have a characteristic star pattern as interactions are zero at zero load. Numerical estimates of the slopes of the interactions can be reverse-engineered from the load prediction equations of a balance if calibration data is examined. The slopes are the off-diagonal coefficients of the inverse of the matrix that has the coefficients of the linear terms of the fitted loads if the Non-Iterative Method is used for the analysis. Similarly, the slopes are the off-diagonal coefficients of the inverse of the matrix that is the non-iterative part of the primary load iteration equation if the Iterative Method is used for the analysis. Data from the calibration of a force balance is processed to illustrate the proposed graphical and numerical description of interactions.

wind tunnel test↗
NASA NTRS2022

A Framework for Deep Learning Emulation of Numerical Models With a Case Study in Satellite Remote Sensing

Numerical models based on physics represent the state of the art in Earth system modeling and comprise our best tools for generating insights and predictions. Despite rapid growth in computational power, the perceived need for higher model resolutions overwhelms the latest generation computers, reducing the ability of modelers to generate simulations for understanding parameter sensitivities and characterizing variability and uncertainty. Thus, surrogate models are often developed to capture the essential attributes of the full-blown numerical models. Recent successes of machine learning methods, especially deep learning (DL), across many disciplines offer the possibility that complex nonlinear connectionist representations may be able to capture the underlying complex structures and nonlinear processes in Earth systems. A difficult test for DL-based emulation, which refers to function approximation of numerical models, is to understand whether they can be comparable to traditional forms of surrogate models in terms of computational efficiency while simultaneously reproducing model results in a credible manner. A DL emulation that passes this test may be expected to perform even better than simple models with respect to capturing complex processes and spatiotemporal dependencies. Here, we examine, with a case study in satellite-based remote sensing, the hypothesis that DL approaches can credibly represent the simulations from a surrogate model with comparable computational efficiency. Our results are encouraging in that the DL emulation reproduces the results with acceptable accuracy and often even faster performance. We discuss the broader implications of our results in light of the pace of improvements in high-performance implementations of DL and the growing desire for higher resolution simulations in the Earth sciences.

Bayesian Deep Learning↗
NASA NTRSDate not supplied

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗
NASA NTRSDate not supplied

Updated Numerical Analysis Benchmarks for Meteoroid Relevant Materials

Proposal for update of numerical analysis benchmark for meteoroid relevant materials. - Two recently performed shots are proposed to be numerical analysis benchmarks for numerical simulations of impacts of high-density meteoroids (Al 2 O 3 surrogate) and low-density meteoroids (Nylon surrogate). - A pair of general Whipple shields have been studied: - Bumper and rear walls are the same material and thickness between shields - Separation is 4.5 cm for Al 2 O 3 and 1.5 cm for Nylon - Information gathered includes high speed (1 MHz) shadowgraphs of debris cloud, bumper hole size and rear wall hole area.

Orbital Debris↗
NASA NTRS2023

Numerical Investigation of C L,max Prediction on the NASA High-Lift Common Research Model

An assessment of a Hybrid RANS/LES (HRLES) approach for C L,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both free air and wind tunnel configurations of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES is shown to address some of the known shortcomings with RANS and prevent inboard and outboard flow separation particularly in the region of C L,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly underperformed RANS. For the in-tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.

TTT↗
NASA NTRSDate not supplied

Numerical Simulation of No-Vent Chill and Fill of a Large Liquid Hydrogen Tank

To enable the design of future in-space cryogenic propellant transfer vehicles, high accuracy models of the cryogenic propellant transfer process are desired. For in-space refueling applications, accurate simulations are desired of the propellant tank chilldown and no-vent fill (NVF) process in order to estimate the amount of propellant required, to determine the maximum sustainable flow rate to fill the tank, and to establish a timeline for the filling process. NVF is a complex process that involves numerous two-phase heat and mass transfer phenomena, such as boiling at the wall and evaporation and condensation at the liquid/vapor interface. In recent years, progress has been made in developing a numerical modeling technique using the Generalized Fluid System Simulation Program (GFSSP) that has been validated against several small-scale liquid nitrogen NVF tests. This paper presents GFSSP numerical model validation against three liquid hydrogen NVF tests on a large-scale propellant tank. Comparing the GFSSP model to the data for pressure, fill level, and wall temperature, the mean absolute percentage error is below 20% across the range of simulated conditions, indicating very good predictive agreement for this large-scale hydrogen NVF test series.

Alok Majumdar↗
NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

Tip clearance effects in centrifugal compressors have been extensively investigated to understand the losses associated with the flow in the impeller tip clearance gap between the rotating blades and the stationary shroud. In Part 3 of this multipart investigation, experimental data and validated numerical simulations from Parts 1 and 2 were used to analyze the effects of the size of the tip clearance gap on the High Efficiency Centrifugal Compressor performance and aerodynamics. Four tip clearance gaps ranging from 0.012-in to 0.030-in (2.0% to 4.9% of the impeller exit blade height) were considered at both design and off-design operating conditions. The total pressure ratio and efficiency of the stage are found to decay linearly with increasing tip gaps. The sensitivity of the impeller performance to the tip gap was found to vary with the rotational speed of the compressor. Spanwise surveys of flow angle, total pressure, and total temperature collected at the impeller exit at design speed are used to validate numerical simulations at each experimental tip gap condition. Numerical simulations show that increased turbulence kinetic energy near the shroud at larger tip gaps leads to spanwise mixing of high entropy fluid near the impeller trailing edge which decreases the useful work input by the impeller. The data presented have been made available to the public in the HECC Data Archive located at https://storage.googleapis.com/hecc-data/NASA-HECC-Data-Archive.zip.

centrifugal compressor↗
NASA NTRSDate not supplied

Numerical Simulation of No-Vent Chill and Fill of a Large Liquid Hydrogen Tank

To enable the design of future in-space cryogenic propellant transfer vehicles, high accuracy models of the cryogenic propellant transfer process are desired. For in-space refueling applications, accurate simulations are desired of the propellant tank chilldown and no-vent fill (NVF) process in order to estimate the amount of propellant required, to determine the maximum sustainable flow rate to fill the tank, and to establish a timeline for the filling process. NVF is a complex process that involves numerous two-phase heat and mass transfer phenomena, such as boiling at the wall and evaporation and condensation at the liquid/vapor interface. In recent years, progress has been made in developing a numerical modeling technique using the Generalized Fluid System Simulation Program (GFSSP) that has been validated against several small-scale liquid nitrogen NVF tests. This paper presents GFSSP numerical model validation against three liquid hydrogen NVF tests on a large-scale propellant tank. Comparing the GFSSP model to the data for pressure, fill level, and wall temperature, the mean absolute percentage error is below 20% across the range of simulated conditions, indicating very good predictive agreement for this large-scale hydrogen NVF test series.

Alok Majumdar↗
NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

Tip clearance effects in centrifugal compressors have been extensively investigated to understand the losses associated with the flow in the impeller tip clearance gap between the rotating blades and the stationary shroud. In Part 3 of this multipart investigation, experimental data and validated numerical simulations from Parts 1 and 2 were used to analyze the effects of the size of the tip clearance gap on the High Efficiency Centrifugal Compressor performance and aerodynamics. Four tip clearance gaps ranging from 0.012-in to 0.030-in (2.0% to 4.9% of the impeller exit blade height) were considered at both design and off-design operating conditions. The total pressure ratio and efficiency of the stage are found to decay linearly with increasing tip gaps. The sensitivity of the impeller performance to the tip gap was found to vary with the rotational speed of the compressor. Spanwise surveys of flow angle, total pressure, and total temperature collected at the impeller exit at design speed are used to validate numerical simulations at each experimental tip gap condition. Numerical simulations show that increased turbulence kinetic energy near the shroud at larger tip gaps leads to spanwise mixing of high entropy fluid near the impeller trailing edge which decreases the useful work input by the impeller. The data presented have been made available to the public in the HECC Data Archive located at https://storage.googleapis.com/hecc-data/NASA-HECC-Data-Archive.zip.

centrifugal compressor↗