Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation codes”

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 NTRS2015

Los Alamos RAGE Simulations of the HAIV Mission Concept

The mitigation of potentially hazardous objects (PHOs) can be accomplished by a variety of methods including kinetic impactors, gravity tractors and several nuclear explosion options. Depending on the available lead time prior to Earth impact, non- nuclear options can be very effective at altering a PHOs orbit. However if the warning time is short nuclear options are generally deemed most effective at mitigating the hazard. The NIAC mission concept for a nuclear mission has been presented at several meetings, including the last PDC (2013).We use the adaptive mesh hydrocode RAGE to perform detailed simulations of this Hypervelocity Asteroid Intercept Vehicle (HAIV) mission concept. We use the RAGE code to simulate the crater formation by the kinetic impactor as well as the explosion and energy coupling from the follower nuclear explosive device (NED) timed to detonate below the original surface to enhance the energy coupling. The RAGE code has been well validated for a wide variety of applications. A parametric study will be shown of the energy and momentum transfer to the target 100 m diameter object: 1) the HAIV mission as planned; 2) a surface explosion and 3) a subsurface (contained) explosion; both 2) and 3) use the same source energy as 1).Preliminary RAGE simulations show that the kinetic impactor will carve out a surface crater on the object and the subsequent NED explosion at the bottom of the crater transfers energy and momentum to the target effectively moving it off its Earth crossing orbit. Figure 1 shows the initial (simplified) RAGE 2D setup geometry for this study. Figure 2 shows the crater created by the kinetic impactor and Figure 3 shows the time sequence of the energy transfer to the target by the NED.

Asteroid Mitigation↗
DOE OSTI2024

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗
DOE OSTI2025

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Maximum-likelihood soft-decision decoding of block codes using the A* algorithm

The A* algorithm finds the path in a finite depth binary tree that optimizes a function. Here, it is applied to maximum-likelihood soft-decision decoding of block codes where the function optimized over the codewords is the likelihood function of the received sequence given each codeword. The algorithm considers codewords one bit at a time, making use of the most reliable received symbols first and pursuing only the partially expanded codewords that might be maximally likely. A version of the A* algorithm for maximum-likelihood decoding of block codes has been implemented for block codes up to 64 bits in length. The efficiency of this algorithm makes simulations of codes up to length 64 feasible. This article details the implementation currently in use, compares the decoding complexity with that of exhaustive search and Viterbi decoding algorithms, and presents performance curves obtained with this implementation of the A* algorithm for several codes.

Ekroot, L.↗
NASA NTRS2009

Desktop Application Program to Simulate Cargo-Air-Drop Tests

The DSS Application is a computer program comprising a Windows version of the UNIX-based Decelerator System Simulation (DSS) coupled with an Excel front end. The DSS is an executable code that simulates the dynamics of airdropped cargo from first motion in an aircraft through landing. The bare DSS is difficult to use; the front end makes it easy to use. All inputs to the DSS, control of execution of the DSS, and postprocessing and plotting of outputs are handled in the front end. The front end is graphics-intensive. The Excel software provides the graphical elements without need for additional programming. Categories of input parameters are divided into separate tabbed windows. Pop-up comments describe each parameter. An error-checking software component evaluates combinations of parameters and alerts the user if an error results. Case files can be created from inputs, making it possible to build cases from previous ones. Simulation output is plotted in 16 charts displayed on a separate worksheet, enabling plotting of multiple DSS cases with flight-test data. Variables assigned to each plot can be changed. Selected input parameters can be edited from the plot sheet for quick sensitivity studies.

Cuthbert, Peter↗
NASA NTRS1980

JPL analytical performance prediction capability

The development of a computer code to simulate and predict the performance of large deployable antenna structures is discussed. The computer code is expected to account for all phases of hardware use, i.e., ground handling, boost, deployment, on orbit station keeping, and transfer from LEO to GEO, if applicable. In addition to accounting for static and dynamic loading of the basic structure, the variation in the precision of the reflector surface and the alignment of the feed support structure as a function orbit are scrutinized in order to accommodate electromagnetic analysis. Specific tasks directed towards the development of the computer code are reported and a number of thermal and structural models are discussed.

Chen, R.↗
NASA NTRS1985

Analytical investigation of the dynamics of tethered constellations in Earth orbit, phase 2

This Quarterly Report deals with the deployment maneuver of a single-axis, vertical constellation with three masses. A new, easy to handle, computer code that simulates the two-dimensional dynamics of the constellation has been implemented. This computer code is used for designing control laws for the deployment maneuver that minimizes the acceleration level of the low-g platform during the maneuver.

Lorenzini, E.↗
NASA NTRS1991

Discrete shape sensitivity equations for aerodynamic problems

Intensive research and development in Computational Fluid Dynamics (CFD) has recently produced many powerful CFD codes to simulate complex aerodynamic phenomena. However, in order to enhance the usefulness of these CFD codes for design practice, development of design sensitivity equations compatible to these codes becomes very important. This paper represents a part of such an effort to develop a sensitivity analysis methodology that enables the sensitivity equations to be implemented into existing CFD codes with minimal code modification. The methodology is based upon a preelimination procedure which accounts for consistently linearized boundary conditions. Formulations of both the direct differentiation and the adjoint variable methods will be presented in the paper.

Hou, Gene J.-W.↗
NASA NTRS2007

Simulating Space Capsule Water Landing with Explicit Finite Element Method

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Wang, John T.↗
DOE OSTI2025

Simulation of an inductively coupled plasma with a two-dimensional Darwin particle-in-cell code

A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

An RL10A-3-3A rocket engine model using the rocket engine transient simulator (ROCETS) software

Steady-state and transient computer models of the RL10A-3-3A rocket engine have been created using the Rocket Engine Transient Simulation (ROCETS) code. These models were created for several purposes. The RL10 engine is a critical component of past, present, and future space missions; the model will give NASA an in-house capability to simulate the performance of the engine under various operating conditions and mission profiles. The RL10 simulation activity is also an opportunity to further validate the ROCETS program. The ROCETS code is an important tool for modeling rocket engine systems at NASA Lewis. ROCETS provides a modular and general framework for simulating the steady-state and transient behavior of any desired propulsion system. Although the ROCETS code is being used in a number of different analysis and design projects within NASA, it has not been extensively validated for any system using actual test data. The RL10A-3-3A has a ten year history of test and flight applications; it should provide sufficient data to validate the ROCETS program capability. The ROCETS models of the RL10 system were created using design information provided by Pratt & Whitney, the engine manufacturer. These models are in the process of being validated using test-stand and flight data. This paper includes a brief description of the models and comparison of preliminary simulation output against flight and test-stand data.

Binder, Michael↗
NASA NTRS1994

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗
NASA NTRS2003

Further Development and Implementation of Implicit Time Marching in the CAA Code

The fellowship research project continued last-year work on implementing implicit time marching concepts in the Broadband Aeroacoustic System Simulator (BASS) code. This code is being developed at NASA Glenn for analysis of unsteady flow and sources of noise in propulsion systems, including jet noise and fan noise.

Golubev, Vladimir V.↗
NASA NTRS2005

Simulation of gaseous diffusion in partially saturated porous media under variable gravity with lattice Boltzmann methods

Liquid distributions in unsaturated porous media under different gravitational accelerations and corresponding macroscopic gaseous diffusion coefficients were investigated to enhance understanding of plant growth conditions in microgravity. We used a single-component, multiphase lattice Boltzmann code to simulate liquid configurations in two-dimensional porous media at varying water contents for different gravity conditions and measured gas diffusion through the media using a multicomponent lattice Boltzmann code. The relative diffusion coefficients (D rel) for simulations with and without gravity as functions of air-filled porosity were in good agreement with measured data and established models. We found significant differences in liquid configuration in porous media, leading to reductions in D rel of up to 25% under zero gravity. The study highlights potential applications of the lattice Boltzmann method for rapid and cost-effective evaluation of alternative plant growth media designs under variable gravity.

NASA Center JSC↗
NASA NTRS1984

Simulation of transonic separated airfoil flow by finite difference viscous-inviscid interaction

A finite difference viscous inviscid interaction program has been developed for simulating the separated transonic flow about lifting airfoils, including the wake. In contrast to most interaction programs, this code combines a finite difference boundary layer algorithm with the inviscid program. The recently developed finite difference boundary layer code efficiently simulates attached and reversed compressible boundary layer and wake flows. New viscous inviscid interaction algorithms were also developed to couple the boundary layer code with the inviscid transonic full potential program. Transonic cases with shock induced and trailing edge separation are computed and compared with experimental and Navier-Stokes results.

Vandalsem, W. R.↗
DOE OSTI2026

LANL T1 t25_mlgpu

The CPU allocation was primarily used for nucleosyn thesis calculations using LANL’s PRISM code, radiative transfer simulations using LANL’s SuperNu code, and parameter inference of the radiative transfer simulations. Below, I document all relevant publications where Inti tutional Computing (IC) resources were used in relation to this allocation and their scientific impacts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

Monte Carlo Simulation of a Knudsen Effusion Mass Spectrometer Sampling System

Knudsen flow is easily simulated with a Monte Carlo method. In this study we develop a Visual Basic for Excel (VBA) code to simulate the molecular beam from a vaporizing solid. The system at NASA Glenn uses the restricted collimation method of Chatillon and colleagues, which consists of two apertures between the effusion cell and the ionizer. The diameter of the first aperture is smaller than the diameter of the effusion cell orifice, so the ionizer effectively only sees inside the effusion cell. The code is able to calculate the transmission coefficient through the cell orifice, through the cell orifice and the first aperture, and through the cell orifice and first and second aperatures. Calculated transmission coefficients through the cell orifice are compared to tabulated values to validate the code. Then transmission coefficients are calculated through the cell orifice and both apertures to the ionizer. This allows the geometry (aperture spacing and diameters) of the sampling system the sampling system to be optimized. Calculated transmission factors are also compared to literature values calculated via an analytic method.

Radke, Michael J.↗