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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1962

Monte Carlo Calculations of Neutron Number Spectra and Buildup Factors in Infinite Conical Configurations

A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.

CONE↗
NASA NTRS1990

Some aspects of robotics calibration, design and control

The main objective is to introduce techniques in the areas of testing and calibration, design, and control of robotic systems. A statistical technique is described that analyzes a robot's performance and provides quantitative three-dimensional evaluation of its repeatability, accuracy, and linearity. Based on this analysis, a corrective action should be taken to compensate for any existing errors and enhance the robot's overall accuracy and performance. A comparison between robotics simulation software packages that were commercially available (SILMA, IGRIP) and that of Kennedy Space Center (ROBSIM) is also included. These computer codes simulate the kinematics and dynamics patterns of various robot arm geometries to help the design engineer in sizing and building the robot manipulator and control system. A brief discussion on an adaptive control algorithm is provided.

Tawfik, Hazem↗
NASA NTRS1998

LOOPREF: A Fluid Code for the Simulation of Coronal Loops

This report documents the code LOOPREF. LOOPREF is a semi-one dimensional finite element code that is especially well suited to simulate coronal-loop phenomena. It has a full implementation of adaptive mesh refinement (AMR), which is crucial for this type of simulation. The AMR routines are an improved version of AMR1D. LOOPREF's versatility makes is suitable to simulate a wide variety of problems. In addition to efficiently providing very high resolution in rapidly changing regions of the domain, it is equipped to treat loops of variable cross section, any non-linear form of heat conduction, shocks, gravitational effects, and radiative loss.

deFainchtein, Rosalinda↗
NASA NTRS2005

Advanced Combustion Modeling for Complex Turbulent Flows

The next generation of aircraft engines will need to pass stricter efficiency and emission tests. NASA's Ultra-Efficient Engine Technology (UEET) program has set an ambitious goal of 70% reduction of NO(x) emissions and a 15% increase in fuel efficiency of aircraft engines. We will demonstrate the state-of-the-art combustion tools developed a t Stanford's Center for Turbulence Research (CTR) as part of this program. In the last decade, CTR has spear-headed a multi-physics-based combustion modeling program. Key technologies have been transferred to the aerospace industry and are currently being used for engine simulations. In this demo, we will showcase the next-generation combustion modeling tools that integrate a very high level of detailed physics into advanced flow simulation codes. Combustor flows involve multi-phase physics with liquid fuel jet breakup, evaporation, and eventual combustion. Individual components of the simulation are verified against complex test cases and show excellent agreement with experimental data.

Ham, Frank Stanford↗
NASA NTRS2020

InSight Entry, Descent and Landing Pre-Flight Performance Predictions

On November 26, 2018, the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) lander successfully touched down on the surface of Mars. Over its seven-plus year development, NASA Langley Research Center’s (LaRC) Program to Optimize Simulated Trajectories II (POST2) was used to assess the mission’s Entry, Descent and Landing (EDL) vehicle system performance against related requirements across the full range of possible environmental and spacecraft conditions. Much of the simulation code was derived from the Phoenix mission, for which this vehicle is very similar. The InSight six degree-of-freedom simulation included models for Mars atmosphere, gravity and digital elevation maps of the landing location. Additionally, vehicle specific aerodynamic, parachute, engine, navigation sensor, flight software and landing radar models were also included. A set of dispersions for each model, as well as for additional simulation input parameters, were also included in order to provide a statistical, Monte Carlo prediction of the EDL system performance. An overview of the pre-flight performance assessments completed, including the various simulation campaigns used, will be provided. Ultimately, this work was critical in the assessment of readiness for InSight launch. A brief description of the use of this simulation in support of flight operations is also discussed.

Robert W Maddock↗
NASA NTRS1977

NASCAP user's manual

The NASCAP (NASA Charging Analyzer Program) code simulates the charging process for a complex object in either tenuous plasma or ground test environment. Detailed specifications needed to run the code are presented. The object definition section, OBJDEF, allows the test object to be easily defined in the cubic mesh. The test object is composed of conducting sections which may be wholly or partially covered with thin dielectric coatings. The potential section, POTENT, obtains the electrostatic potential in the space surrounding the object. It uses the conjugate gradient method to solve the finite element formulation of Poisson's equation. The CHARGE section of NASCAP treats charge redistribution among the surface cells of the object as well as charging through radiation bombardment. NASCAP has facilities for extensive graphical output, including several types of object display plots, potential contour plots, space charge density contour plots, current density plots, and particle trajectory plots.

Mandell, M. J.↗
NASA NTRSDate not supplied

Computational Icing Analysis on NASA’s SIDRM Geometry to Investigate Collection Efficiency

Computational icing analysis results were compared to experimental icing tunnel data including aerothermal (e.g. dry air) and supercooled water droplet rime-ice conditions from tests conducted in early 2022 at the NASA Icing Research Tunnel (IRT). The Simulated Intercompressor Duct Research Model (SIDRM) test article was used in this study, and its geometry represents the inter-compressor duct region of a turbofan engine. The test article’s purpose is to study the physics of supercooled water icing and ice crystal icing. This study compared three different icing codes: FENSAP-ICE (Eulerian approach), LEWICE3D (Lagrangian approach), and GlennICE (Lagrangian approach). All three icing codes were conducted on SIDRM’s complex body flow-field and compared to different experimental supercooled water rime runs. The test article instrumentation (pressure taps, thermocouples, etc.) and 3D laser scans of final ice shapes were used to compare against the different icing code simulations. The overall objectives are to understand how the icing codes handle capturing collection efficiency on the complex test article’s unheated surfaces. In the aerothermal cases, pressure tap readings matched the CFD results, but dry air CFD underpredicted thermocouple readings. Collection efficiency results from all three icing codes matched well together on the main body leading edge, main body slope, and the strut leading edges. All three icing codes underpredicted the final ice shape using a single-shot constant ice density approach, with more difficulty coming from the strut leading edge ice shape due to the swept wing like flow field. NASA’s overall goal for this effort is to develop computational icing tools to assist in the design and certification of engines for flight in icing conditions.

Supercooled icing, SIDRM, aerodynamics, computatio↗
NASA NTRS1992

Numerical simulations of the electrodynamic interactions between the Tethered-Satellite-System and space plasma

The first Tethered-Satellite-System (TSS-1), scheduled for a flight in late 1992, is expected to provide relevant information related to the concept of generating an emf in a 20-km-long (or longer) conducting wire. This paper presents numerical simulations of the electrodynamic interactions between the TSS system and space plasma, using a 2D and 3D models of the system. The 2D case code simulates the motion of a long cylinder past a plasma, which is composed of electrons and H(+) ions. The system is solved by allowing the plasma to flow past the cylinder with an imposed magnetic field. The more complex 3D case is considered to study the dynamics in great detail. Results of 2D simulation show that the interaction of a satellite with plasma flowing perpendicularly to the magnetic field results in an enhancement in the current collection.

Vashi, Bharat I.↗
NASA NTRS2017

Neutron Angular Scatter Effects in 3DHZETRN: Quasi-Elastic

The current 3DHZETRN code has a detailed three dimensional (3D) treatment of neutron transport based on a forward/isotropic assumption and has been compared to Monte Carlo (MC) simulation codes in various geometries. In most cases, it has been found that 3DHZETRN agrees with the MC codes to the extent they agree with each other. However, a recent study of neutron leakage from finite geometries revealed that further improvements to the 3DHZETRN formalism are needed. In the present report, angular scattering corrections to the neutron fluence are provided in an attempt to improve fluence estimates from a uniform sphere. It is found that further developments in the nuclear production models are required to fully evaluate the impact of transport model updates. A model for the quasi-elastic neutron production spectra is therefore developed and implemented into 3DHZETRN.

Wilson, John W.↗
NASA NTRS2012

Calculation of the Frequency Distribution of the Energy Deposition in DNA Volumes by Heavy Ions

Radiation quality effects are largely determined by energy deposition in small volumes of characteristic sizes less than 10 nm representative of short-segments of DNA, the DNA nucleosome, or molecules initiating oxidative stress in the nucleus, mitochondria, or extra-cellular matrix. On this scale, qualitatively distinct types of molecular damage are possible for high linear energy transfer (LET) radiation such as heavy ions compared to low LET radiation. Unique types of DNA lesions or oxidative damages are the likely outcome of the energy deposition. The frequency distribution for energy imparted to 1-20 nm targets per unit dose or particle fluence is a useful descriptor and can be evaluated as a function of impact parameter from an ions track. In this work, the simulation of 1-Gy irradiation of a cubic volume of 5 micron by: 1) 450 (1)H(+) ions, 300 MeV; 2) 10 (12)C(6+) ions, 290 MeV/amu and 3) (56)Fe(26+) ions, 1000 MeV/amu was done with the Monte-Carlo simulation code RITRACKS. Cylindrical targets are generated in the irradiated volume, with random orientation. The frequency distribution curves of the energy deposited in the targets is obtained. For small targets (i.e. <25 nm size), the probability of an ion to hit a target is very small; therefore a large number of tracks and targets as well as a large number of histories are necessary to obtain statistically significant results. This simulation is very time-consuming and is difficult to perform by using the original version of RITRACKS. Consequently, the code RITRACKS was adapted to use multiple CPU on a workstation or on a computer cluster. To validate the simulation results, similar calculations were performed using targets with fixed position and orientation, for which experimental data are available [5]. Since the probability of single- and double-strand breaks in DNA as function of energy deposited is well know, the results that were obtained can be used to estimate the yield of DSB, and can be extended to include other targeted or non-target effects.

Plante, Ianik↗
NASA NTRS1988

Overview of numerical codes developed for predicted electrothermal deicing of aircraft blades

An overview of the deicing computer codes that have been developed at the University of Toledo under sponsorship of the NASA-Lewis Research Center is presented. These codes simulate the transient heat conduction and phase change occurring in an electrothermal deicier pad that has an arbitrary accreted ice shape on its surface. The codes are one-dimensional rectangular, two-dimensional rectangular, and two-dimensional with a coordinate transformation to model the true blade geometry. All modifications relating to the thermal physics of the deicing problem that have been incorporated into the codes will be discussed. Recent results of reformulating the codes using different numerical methods to increase program efficiency are described. In particular, this reformulation has enabled a more comprehensive two-dimensional code to run in much less CPU time than the original version. The code predictions are compared with experimental data obtained in the NASA-Lewis Icing Research Tunnel with a UH1H blade fitted with a B. F. Goodrich electrothermal deicer pad. Both continuous and cyclic heater firing cases are considered. The major objective in this comparison is to illustrate which codes give acceptable results in different regions of the airfoil for different heater firing sequences.

Keith, Theo G.↗
DOE OSTI2026

A Comparison of Three Neodymium Atomic Data Sets for Kilonova Modeling

We examine the impact of input neodymium (Nd) atomic data on the light curves and spectra of kilonovae (KNe), probing the sensitivity of kilonova observables to the atomic physics of this important lanthanide element. We use the SuperNu Monte Carlo radiative transfer code, simulating a simple semianalytic 1D kilonova (KN) with a pure Nd atmosphere, fixing the radiative transfer method while using input atomic data generated by three different codes: the LANL suite of atomic physics codes, HULLAC, and Autostructure. We see that the choice of atomic data significantly shapes the resulting light curves and spectra. Peak bolometric luminosities differ by a ratio of nearly 1.5 between HULLAC/Autostructure and LANL data sets. Moreover, we observe significant near- to mid-IR differences in the structure of the spectra. We specifically attribute these differences to the choice of atomic data for neutral Nd I. Many of the results here have been adapted from a presentation at “Radiative Transfer and Atomic Physics of Kilonovae” in Stockholm, 2023. We additionally present a LANL data set with energies calibrated to available values in the NIST Atomic Spectra Database, and demonstrate that this calibration also significantly affects IR spectral structure at late time. The substantial differences in KN observables that arise from tuning the atomic data of just one lanthanide element highlight the special attention that must be paid to atomic physics uncertainties when modeling KNe, from AT2017gfo to beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Electron acceleration at nearly perpendicular collisionless shocks. I - One-dimensional simulations without electron scale fluctuations

Under certain conditions electrons can be reflected and effectively energized at quasi-perpendicular shocks. This process is most prominent close to the point where the upstream magnetic field is tangent to the curved shock. A theoretical explanation of the underlying physical mechanism has been proposed which assumes conservation of magnetic moment and a static, simplified shock profile are performed. Test particle calculations of the electron reflection process in order to examine the results of the theoretical analysis without imposing these restrictive conditions. A one-dimensional hybrid simulation code generates the characteristic field variations across the shock. Special emphasis is placed on the spatial and temporal length scales involved in the mirroring process. The simulation results agree generally well with the predictions from adiabatic theory. The effects of the cross-shock potential and unsteadiness are quantified, and the influence of field fluctuations on the reflection process is discussed.

Krauss-Varban, D.↗
NASA NTRS2001

Formation and Early Evolution of Solar and Extra-Solar Giant Planets

This project investigates the origin of giant planets, both in the Solar System and around other stars. It is assumed that the planets form by the core accretion process: small solid particles in a disk surrounding a young star gradually coagulate into objects of a few kilometers in size, known as planetesimals, which then accumulate into solid protoplanetary cores. Once the cores have become large enough, they are able to attract gas from the surrounding disk to form the deep gaseous envelope of the giant planet. Our code simulates giant planet growth in a spherical approximation, and it has been quite successful in addressing a number of basic planetary properties. Further improvements to the code have been made to achieve a more realistic understanding of planetary formation. The computations of the models were based on an earlier version of our code and were stopped at the onset of runaway gas accretion. Now, improved boundary conditions have been incorporated into the code to allow for hydrodynamic inflow of gas and to handle the late stages of evolution when the planet evolves at constant mass. These changes were made to the version of the code that uses a constant accretion rate and to the version that uses a self-consistent method for calculating both the solid and gas accretion rates. The equation of state has been updated to incorporate the detailed tables of Saumon, Chabrier, and Van Horn. The opacities were updated to include the results of Alexander and Ferguson. The outer boundary conditions were modified. During the accretion phase when the planet's radius is between the accretion radius and the tidal radius, we set the outer boundary at a 'modified' accretion radius, which is the point where thermal energy is enough to bring gas to the edge of the Hill sphere.

Bodenheimer, P. H.↗
NASA NTRS2004

Ascent, Stage Separation and Glideback Performance of a Partially Reusable Small Launch Vehicle

An integrated analysis is presented for ascent, stage separation and glide back performance of a small, partially reusable launch vehicle sized for a payload of about 330 lbs to a 150 nm polar orbit. The altitude margin was used a performance metric for the glideback performance. Aerodynamic databases for each of these three phases of flight were developed using a combination of engineering level code, free stream and proximity wind tunnel test data and Euler CFD results. The ascent and glideback trajectories were generated using POST and the stage separation simulation was done using the in-house software Sep-Sim as a front end to the commercially available multi-body dynamic simulation code ADAMS. The payload to the designated polar orbit was optimized subject to the constraints imposed by stage separation and adequate performance reserve for the glideback booster in addition to the usual ascent trajectory constraints.

Pamadi, Bandu N.↗
NASA NTRS1987

Computational chemistry

With the advent of supercomputers, modern computational chemistry algorithms and codes, a powerful tool was created to help fill NASA's continuing need for information on the properties of matter in hostile or unusual environments. Computational resources provided under the National Aerodynamics Simulator (NAS) program were a cornerstone for recent advancements in this field. Properties of gases, materials, and their interactions can be determined from solutions of the governing equations. In the case of gases, for example, radiative transition probabilites per particle, bond-dissociation energies, and rates of simple chemical reactions can be determined computationally as reliably as from experiment. The data are proving to be quite valuable in providing inputs to real-gas flow simulation codes used to compute aerothermodynamic loads on NASA's aeroassist orbital transfer vehicles and a host of problems related to the National Aerospace Plane Program. Although more approximate, similar solutions can be obtained for ensembles of atoms simulating small particles of materials with and without the presence of gases. Computational chemistry has application in studying catalysis, properties of polymers, all of interest to various NASA missions, including those previously mentioned. In addition to discussing these applications of computational chemistry within NASA, the governing equations and the need for supercomputers for their solution is outlined.

Arnold, J. O.↗
NASA NTRS2018

Optimization of Elastodynamic Finite Integration Technique on Intel Xeon Phi Knights Landing Processors

This work describes the development and optimization of an implementation of an isotropic elastodynamic finite integration technique (EFIT) code for parallelized computation on Intel Knights Landing (KNL) hardware. EFIT is a numerical approach resulting in standard staggered-grid finite difference equations for the elastodynamic equations of motion to simulate bulk waves is solids. The computationally efficient simulation of elastodynamic wave propagation and interactions in aerospace materials is of high-interest in the fields of nondestructive evaluation (NDE) and structural health monitoring (SHM). Ultrasonic inspection uses an ultrasonic signal, generated at the surface of the material/structure via use of a piezoelectric transducer, to propagate sound waves into the material where it interacts with any existing defects, as well as with structural boundaries and any material inhomogeneity. Reflections from defects and boundaries are then measured by a transducer. Realistic ultrasound simulation tools can significantly aid the development and optimization of inspection techniques and can assist in the interpretation of experimental data. The optimization of an elastodynamics simulation code for the KNL Many Integrated Core processor was performed. The optimization focused on data locality and vectorization. Results show that tiling of the data to exploit the cache behavior and allow for significant utilization of the KNL hardware. The MPI implementation allows for a scalable implementation enabling large problems to be simulated. The model results were validated against theoretical dispersion curves to within 2% of the group velocity, and within 0.5% of the phase velocity of the A0 mode. Aggressive use of tiling, threading, and vectorization techniques allowed for dramatically improved time to solution.

William C. Schneck↗
DOE OSTI2026

Kinetic Model of HoxEFU reduction by NADH [SWR-26-087]

This repository is used to release code generated for manuscripts on the Photosynthetic Energy Transduction core program. This code simulates the reduction of HoxEFU by NADH. The electro transfer rate constants for the simulation are specified in the .csv files. The two .csv files correspond tot he two models described in Dawson et al. Cell. Rep. Phys. Sci. 2026. The code utilizes a chemical master equation, a set of differential equations, defining the time evolution of the oxidation and reduction kinetics of NAD+, NADH, a FMN flavin, and a set of iron sulfur clusters. The kinetics of HoxEFU reduction by NADH are evaluated by numerical integration of the chemical master equation using a variable-time-step Runge-Kutta algorithm.

Dahl, Peter [National Laboratory of the Rockies (N↗