Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid code”

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,279 records · Page 71

An improved supersonic, three-dimensional, external, inviscid flow field code

A numerical procedure was developed to compute the inviscid super/hypersonic flow fields about complex vehicle geometries accurately and efficiently. A second-order accurate finite difference scheme is used to integrate the three-dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine-Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects for equilibrium air are included using curve fits of Mollier charts. This report deals only with modifications to these procedures in four specific areas: inlet mass ingestion, subsonic axial Mach number, improved conformal mappings, and vehicles flying at yaw. In each area both the modifications to the computational procedures and computer code are discussed.

Marconi, F.↗

Laser propagation in a subcritical foam: Subgrid model

Here, we present a subgrid model for laser propagation in a subcritical foam. Our model describes the expansion of laser-irradiated foam elements that are below the resolution of the simulation grid and predicts the plasma conditions that result from burning down the foam. Our model can be included as a module within a larger multiphysics code, and we have implemented it within the code pF3D, which is used for simulating a laser-plasma interaction. The model predicts a reduced propagation velocity for a laser through a subcritical foam compared to simulating that foam as a homogeneous gas. This is attributed to the laser energy that goes into burning down the foam microstructure. We compare our model against experimental data by simulating a 2 mg/cc SiO 2 foam shot performed at the Janus laser facility at the Lawrence Livermore National Laboratory. pF3D simulations with the foam model predict hot ion temperatures. This leads to a reduction in the level of stimulated Brillouin scattering (SBS), bringing the simulated level of SBS into agreement with the data. Intensity fluctuations at the foam front due to laser speckles and refraction result in ion temperature fluctuations when the foam burns down. These drive long-lived electron density fluctuations on scales that are large compared to the pore size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AnisWave2D

AnisWave2D is a 2D finite-difference code for simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.

Toomey, Aoife↗

Software Control Program For Transportable Microgrid State-of-charge Balancing And Frequency Stability Controls

A deterministic state-of-charge (SOC) balancing approach software control code is introduced as an integral secondary management to primary control layer of an islanded small microgrid or nanogrid system made up of multiple grid-forming inverter/battery/solar combination systems, where each set of batteries with each inverter are on independent DC buses (i.e. non-paralleled on the DC sides). A DERMS-level control approach, algorithm and automation controller program was developed to improve coordination and enable microgrid asset compliance and SOC balancing, enabling provision of a system-level power stability support architecture, load support, and asset scalability. The architecture is configured to treat each unit or micro/nano-grid as a node in a microgrid network, allowing for autonomous DERMS control regarding load and SOC balancing and power stability. As the network grows with the addition of units, greater coordination efforts may be required. The ideal small network microgrid ranges from 2-10 inverter/battery units before additional control parameters must be considered in the existing architecture. The control approach focuses on a deterministic state-of-charge analysis as the primary level control process followed by a secondary control loop using a forced frequency-watt droop strategy to conform off-the-shelf components into behaving under a leader-follower configuration. Adopting this control scheme has been shown to allow for a balanced, unit-coordinated microgrid network, enabling stable power flow. The deterministic state-of-charge approach is introduced as an integral primary control layer of an islanded small network microgrid. A standard strategy for SOC balancing is implementing a battery management system (BMS) to control SOC on the DC side. An alternative approach is to determine how to coordinate sending and receiving power on the AC side with multiple units. The latter approach assesses all the integrated units in the microgrid network. Once the individual units are identified, further system data is required to calculate each unit's total kWh, provided information about its capability to supply or consume kWh and availability. The secondary control layer in the multi-layered small network microgrid methodology uses the primary layer’s decision to initiate frequency setpoint changes, initializing the SOC balancing. The secondary control layer considers numerous system-dependent variables to enable a charging and discharging profile based on adjustable frequency setpoints. The combined architecture will result in stable, coordinated power flow enhancing an AC microgrid's functionalities.

Myers, KurtS [Idaho National Laboratory (INL), Ida↗

Dayflow-PR: High-Resolution Streamflow Reanalysis for Puerto Rico, Version 1.0

This dataset presents a high-resolution historical streamflow reanalysis for NHDPlusV2 stream reaches across Puerto Rico (PR) spanning 1950 - 2019. The reanalysis is generated using the calibrated VIC-RAPID hydrologic modeling framework at the Hydrologic Unit Code Sub-basin (HUC08) scale, forced with sub-daily and daily meteorological forcings from Daymet. Runoff is simulated on 1- and 6-km grids, and the resulting total runoff is routed through the NHDPlusV2 river network using the RAPID routing model to produce Naturalized Streamflow Reanalysis. Where complete observational records are available over 1980 - 2019, streamflows are assimilated (substituted) and subsequently routed downstream through the river network to produce Assimilated Streamflow Reanalysis. The dataset includes streamflow outputs from eight distinct hydrologic modeling configurations along with key performanc evaluation metrics at daily and monthly scales, supporting a wide range of water resource applications. This dataset is derived to support the Non-Powered Dam Assessment, as well as 9505 Secure Water Assessment projects for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, refer to Ghimire et al. (2023), Kao et al. (2024), and Ghimire et al. (2025).

13 HYDRO ENERGY↗

NRAP-Open-IAM: Generic Aquifer Component Development and Testing

The Generic Aquifer Model calculates the concentrations of dissolved salt and dissolved CO 2 surrounding a leaking legacy well. The Generic Aquifer model can also estimate the size of an “impact plume” where concentration changes exceed user-specified thresholds. The model is a component of NRAP-Open-IAM, an open-source Integrated Assessment Model (IAM) developed by the National Risk Assessment Partnership (NRAP) to perform risk assessment for geologic CO 2 storage. The input parameters were selected to cover a wide range of groundwater aquifers and leakage rates. The generic aquifer model was developed using a generative adversarial deep learning network, trained using a large synthetic dataset of STOMP multiphase flow simulations. The deep learning model predictions of dissolved salt and dissolved CO 2 in the aquifer compare well to the original STOMP simulation results. The extent of aquifer impacted by leaking CO 2 or brine is calculated using a user-defined mass fraction threshold. The aquifer impact volumes calculated based on STOMP simulation results compare well to those calculated based on the deep learning model. In a provided python script, gridded observation results from the generic aquifer component of NRAP-Open-IAM are converted to HDF5 format files for monitoring design with the DREAM code.

54 ENVIRONMENTAL SCIENCES↗

Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. II. Test Particles in Electromagnetic Fields and GRMHD

Observations of compact objects, in the form of radiation spectra, gravitational waves by LIGO/Virgo, and direct imaging with the Event Horizon Telescope, are currently the main information sources on plasma physics in extreme gravity. Modeling such physical phenomena Requires numerical methods that allow for the simulation of microscopic plasma dynamics in the presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018), we presented a detailed study of numerical techniques for the integration of free geodesic motion. In this paper, we extend the study by introducing electromagnetic forces in the simulation of charged particles in curved spacetimes. We extend the Hamiltonian energy-conserving method presented in Bacchini et al. to include the Lorentz force, and we test its performance compared to that of standard explicit Runge–Kutta and implicit midpoint rule schemes against analytic solutions. Then, we show the application of the numerical schemes to the integration of test particle trajectories in general relativistic magnetohydrodynamic (GRMHD) simulations by modifying the algorithms to handle grid-based electromagnetic fields. We test this approach by simulating ensembles of charged particles in a static GRMHD configuration obtained with the black hole accretion code (BHAC).

79 ASTRONOMY AND ASTROPHYSICS↗

Crossed-grid charge locator

Circuit locates center of cloud of charge on wire grid to within 6.5 micrometers. Wires in vicinity of charge cloud develop voltages that are processed by priority encoders to develop coarse and fine position codes. Device is used with microchannel plate amplifier in X-ray photon detectors, electron microscopes, and closed-circuit television.

Harrison, D. C.↗

An enhanced version of an implicit code for the Euler equations

A two-dimensional implicit finite-difference code is applied to the inviscid Euler equations to compute transonic flow past airfoils in order to provide well-documented standard test cases for the general user community. The code is an improved version of Steger's 1976 implicit code. Enhancements include the use of up-wind differencing in supersonic regions before shocks and a variable time step to accelerate convergence. An airfoil grid generation routine based on algebraic techniques is employed. The grids are clustered near shocks to improve resolution. Computed results are compared with other numerical results from the literature.

Pulliam, T. H.↗

Adaptive finite difference methods in fluid dynamics

An adaptive method to solve partial differential equations in fluid mechanics is presented. The approach requires internal boundary conditions that must be conservative, data structures for keeping track of several layers of fine grid patches, error estimation, and heuristics for automatic grid generation. In practical calculations gains in computer efficiency up to 10 over nonadaptive methods are observed. The whole procedure takes 3000 lines of FORTRAN code.

Berger, Marsha J.↗

BRYNTRN: A baryon transport computer code, computation procedures and data base

The development is described of an interaction data base and a numerical solution to the transport of baryons through the arbitrary shield material based on a straight ahead approximation of the Boltzmann equation. The code is most accurate for continuous energy boundary values but gives reasonable results for discrete spectra at the boundary with even a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H2O).

Wilson, John W.↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the 'multiple grids' or 'zonal grids' approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, M.↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-state Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, Monchai↗

Thermal History and Mantle Dynamics of Venus

One objective of this research proposal is to develop a 3-D thermal history model for Venus. The basis of our study is a finite-element computer model to simulate thermal convection of fluids with highly temperature- and pressure-dependent viscosities in a three-dimensional spherical shell. A three-dimensional model for thermal history studies is necessary for the following reasons. To study planetary thermal evolution, one needs to consider global heat budgets of a planet throughout its evolution history. Hence, three-dimensional models are necessary. This is in contrasts to studies of some local phenomena or local structures where models of lower dimensions may be sufficient. There are different approaches to treat three-dimensional thermal convection problems. Each approach has its own advantages and disadvantages. Therefore, the choice of the various approaches is subjective and dependent on the problem addressed. In our case, we are interested in the effects of viscosities that are highly temperature dependent and that their magnitudes within the computing domain can vary over many orders of magnitude. In order to resolve the rapid change of viscosities, small grid spacings are often necessary. To optimize the amount of computing, variable grids become desirable. Thus, the finite-element numerical approach is chosen for its ability to place grid elements of different sizes over the complete computational domain. For this research proposal, we did not start from scratch and develop the finite element codes from the beginning. Instead, we adopted a finite-element model developed by Baumgardner, a collaborator of this research proposal, for three-dimensional thermal convection with constant viscosity. Over the duration supported by this research proposal, a significant amount of advancements have been accomplished.

Hsui, Albert T.↗

Efficiency and Accuracy of Time-Accurate Turbulent Navier-Stokes Computations

The accuracy and efficiency of two types of subiterations in both explicit and implicit Navier-Stokes codes are explored for unsteady laminar circular-cylinder flow and unsteady turbulent flow over an 18-percent-thick circular-arc (biconvex) airfoil. Grid and time-step studies are used to assess the numerical accuracy of the methods. Nonsubiterative time-stepping schemes and schemes with physical time subiterations are subject to time-step limitations in practice that are removed by pseudo time sub-iterations. Computations for the circular-arc airfoil indicate that a one-equation turbulence model predicts the unsteady separated flow better than an algebraic turbulence model; also, the hysteresis with Mach number of the self-excited unsteadiness due to shock and boundary-layer separation is well predicted.

Rumsey, Christopher L.↗

NASA Tech Briefs, May 2007

Topics include: Noise-Canceling Helmet Audio System; Program Analyzes Spacecraft/Ground Radio Links; Two-Way Communication Using RFID Equipment and Techniques; Six-Message Electromechanical Display System; Scanning Terahertz Heterodyne Imaging Systems; Master Clock and Time-Signal-Distribution System; Synchronous Phase-Resolving Flash Range Imaging; Integrated Radial Probe Transition From MMIC to Waveguide; Bar-Code System for a Microbiological Laboratory; MMIC Amplifier Produces Gain of 10 dB at 235 GHz; Mapping Nearby Terrain in 3D by Use of a Grid of Laser Spots; Digital Beam Deflectors Based Partly on Liquid Crystals; Narrow-Band WGM Optical Filters With Tunable FSRs; Better Finite-Element Analysis of Composite Shell Structures; Computing Spacecraft-Pointing Vectors for Limb Tracking; Enhanced Master Controller Unit Tester; Rover Graphical Simulator; Increasing Durability of Flame-Sprayed Strain Gauges; Multifunctional, High-Temperature Nanocomposites; Multilayer Impregnated Fibrous Thermal Insulation Tiles; Radiation-Shielding Polymer/Soil Composites; Film/Adhesive Processing Module for Fiber-Placement Processing of Composites; Fabrication of Submillimeter Axisymmetric Optical Components; Electrochemical Disposal of Hydrazines in Water; Statistical Model of Evaporating Multicomponent Fuel Drops; Resistively Heated SiC Nozzle for Generating Molecular Beams; Compact Packaging of Photonic Millimeter-Wave Receiver; Diffractive Combiner of Single-Mode Pump Laser-Diode Beams; Wide-Band, High-Quantum-Efficiency Photodetector; A Robustly Stabilizing Model Predictive Control Algorithm; Modeling Evaporation of Drops of Different Kerosenes; Development of Vapor-Phase Catalytic Ammonia Removal System; Several Developments in Space Tethers; Design Concept for a Nuclear Reactor-Powered Mars Rover; Formation-Initialization Algorithm for N Spacecraft; and DNSs of Multicomponent Gaseous and Drop-Laden Mixing Layers Achieving Transition to Turbulence.

Source record↗

A Versatile Simulation Framework for Elastodynamic Modeling of Structural Health Monitoring

Structural health monitoring (SHM) has the capacity to reduce failure by detecting damage during service life, by periodic, automated monitoring. Guided Wave (GW) Ultrasound is a common SHM approach for aerospace structures. Modelling the physics of GW SHM systems provides a route for understanding system dependencies, capabilities and limitations as damage evolves during service life. Such a toolset can strengthen the understanding of the connection between GW SHM results and the true material state. The most useful modelling tools are those that provide versatile solutions with respect to the simulated component geometry and computational grid connectivity. This work details a versatile application programming interface (API) for the elastodynamic finite integration technique for modelling GW SHM of metals. The custom code implementation, EFIT-CompCell, allows for the modelling of diverse geometries by automatically balancing the message passing interface parallelization layout. The user provides the basic parameters of the simulation and the software automatically performs an initial balancing based on anticipated computational loads, and establishes the CPU communication patterns for any geometry. This work describes the programming philosophy and code structure used to create EFIT-CompCell and compares its performance and capacity to simulation tools that are more specialized for specific architectures. Results are presented for a simulation of GW SHM of an aluminum fuselage section being tested by the FAA. The simulation consists of 733M voxels which took approximately 70 hours to complete 25000 time steps using 40 Intel Xeon E5-4650v2 Ivy Bridge processor cores.

Gregory, Elizabeth D.↗

CFD Simulations of the Space Launch System Ascent Aerodynamics and Booster Separation

This paper presents details of Computational Fluid Dynamic modeling of the Space Launch System during ascent. The primary focus of the paper is the flow simulation of the vehicle during ascent using the Overflow Navier-Stokes code. Computations of 739 first-stage flight conditions covering a range of Mach numbers, angles of attack, and roll angles were computed. The overset grid system contained 375 million grid points, and over 28 mil- lion CPU hours were used in the simulations. The simulations were run on the Pleiades supercomputer at the NASA Advanced Supercomputer Center at Ames Research Center. The data products from this work include integrated line-loads, surface pressure coefficients, venting pressures, and protuberance air-loads. Detailed comparisons were made of the aerodynamic performance predicted by Overflow and the wind-tunnel derived aero- dynamic database. A small number of the cases were run with two different turbulence models and with two differencing schemes. These results were used to quantify the sensitivity to the choice of the turbulence model and to the differencing scheme. The paper also introduces an effort to use the inviscid, unstructured Cartesian solver Cart3D to compute the aerodynamics during booster separation. Adaptive mesh refinement is being used to enable accurate simulations of sixteen booster-separation-motor plumes. The use of this tool is explored in preparation for building a booster-separation aerodynamic database.

Ascent Aerodynamics↗