A self-consistent computer code for particle simulations
A self-consistent computer code for particle simulations which could treat both Landau and cyclotron interactions is discussed.
Source record↗
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.
Minimum Hamiltonian Ascent Trajectory Evaluation (MASTRE) program (update to automatic flight trajectory design, performance prediction, and vehicle sizing for support of Shuttle and Shuttle derived vehicles) engineering manual
The Minimum Hamiltonian Ascent Trajectory Evaluation (MASTRE) program and its predecessors, the ROBOT and the RAGMOP programs, have had a long history of supporting MSFC in the simulation of space boosters for the purpose of performance evaluation. The ROBOT program was used in the simulation of the Saturn 1B and Saturn 5 vehicles in the 1960's and provided the first utilization of the minimum Hamiltonian (or min-H) methodology and the steepest ascent technique to solve the optimum trajectory problem. The advent of the Space Shuttle in the 1970's and its complex airplane design required a redesign of the trajectory simulation code since aerodynamic flight and controllability were required for proper simulation. The RAGMOP program was the first attempt to incorporate the complex equations of the Space Shuttle into an optimization tool by using an optimization method based on steepest ascent techniques (but without the min-H methodology). Development of the complex partial derivatives associated with the Space Shuttle configuration and using techniques from the RAGMOP program, the ROBOT program was redesigned to incorporate these additional complexities. This redesign created the MASTRE program, which was referred to as the Minimum Hamiltonian Ascent Shuttle TRajectory Evaluation program at that time. Unique to this program were first-stage (or booster) nonlinear aerodynamics, upper-stage linear aerodynamics, engine control via moment balance, liquid and solid thrust forces, variable liquid throttling to maintain constant acceleration limits, and a total upgrade of the equations used in the forward and backward integration segments of the program. This modification of the MASTRE code has been used to simulate the new space vehicles associated with the National Launch Systems (NLS). Although not as complicated as the Space Shuttle, the simulation and analysis of the NLS vehicles required additional modifications to the MASTRE program in the areas of providing additional flexibility in the use of the program, allowing additional optimization options, and providing special options for the NLS configuration.
Lyons, J. T.↗
Monte Carlo analysis of dissociation and recombination behind strong shock waves in nitrogen
Computations are presented for the relaxation zone behind strong, 1D shock waves in nitrogen. The analysis is performed with the direct simulation Monte Carlo method (DSMC). The DSMC code is vectorized for efficient use on a supercomputer. The code simulates translational, rotational and vibrational energy exchange and dissociative and recombinative chemical reactions. A model is proposed for the treatment of three body-recombination collisions in the DSMC technique which usually simulates binary collision events. The model improves previous models because it can be employed with a large range of chemical-rate data, does not introduce into the flow field troublesome pairs of atoms which may recombine upon further collision (pseudoparticles) and is compatible with the vectorized code. The computational results are compared with existing experimental data. It is shown that the derivation of chemical-rate coefficients must account for the degree of vibrational nonequilibrium in the flow. A nonequilibrium-chemistry model is employed together with equilibrium-rate data to compute the flow in several different nitrogen shock waves.
Boyd, I. D.↗
Simulating Magnetic Reconnection in ‘Two Ribbon’ Type Solar Flares
Magnetic reconnection is an astrophysical process where neighboring magnetic field lines, facing anti-parallel, are reconfigured. This reconfiguration results in built up magnetic energy being explosively released as it is being converted to plasma kinetic and thermal energies. There are various kinds of simulations used to simulation reconnection; our work begins with Athena++, a magnetohydrodynamic (MHD) simulation code typically used for astrophysical problems, and a reconnection specific code file. The resistive MHD equations are solved with Riemann solvers. There was an initial test run without modifying the code to understand the dynamics of the simulation. We expand on the original reconnection problem file by implementing a radiative cooling term specific to the corona. The radiative cooling is theorized to have an effect on solar coronal plasma and magnetic reconnection dynamics. The cooling term will be tested with various parameters and compared to the case without cooling to study these dynamics. The condensation found in only the with cooling case emphasizes the importance of implementing this feature and will be later tested with a thermal conduction term. We want to determine the parameter regime where non-equilibrium cooling will be important for the reconnection dynamics.
79 ASTRONOMY AND ASTROPHYSICS↗
Boundary conditions which lead to excitation of instabilities in plasma simulations
Two examples of two-dimensional electrostatic particle-code simulations are shown in which one exhibits characteristics of a stable plasma while the other exhibits unstable, long wavelength plasma oscillations. The only difference between the two simulations is a change in the boundary condition on the electrostatic potential. An energy theorem is derived which shows that the rate of change of field and particle energy within a closed volume is related to a surface integral involving the electrostatic potential and the normal component of the electric current. An analytic theory is developed for a one-dimensional plasma to show how boundary effects can excite spurious plasma instabilities. The theory is tested with a series of one-dimensional plasma simulations. Finally, practical considerations on means of avoiding the non-physical instabilities in simulation plasmas are given.
Swift, D. W.↗
Early time interaction of lithium ions with the solar wind in the AMPTE mission
The early time interaction of an artifically injected lithium cloud with the solar wind is simulated with a one-dimensional hybrid code. Simulation results indicate that the lithium cloud presents an obstacle to the solar wind flow, forming a shock-like interaction region. Several notable features are found: (1) The magnetic field is enhanced up to a factor of about 6, followed by a magnetic cavity downstream. (2) Solar wind ions are slowed down inside the lithium cloud, with substantial upstream reflection. (3) Most of the lithium ions gradually pick up the velocity of the solar wind and move downstream. (4) Intense and short-wavelength electric fields exist ahead of the interaction region. (5) Strong electron heating occurs within the lithium cloud. (6) The convection electric field in the solar wind is modulated in the interaction region. The simulation results are in remarkable agreement with in situ spacecraft measurements made during lithium releases in the solar wind by the AMPTE (Active Magnetospheric Particle Tracer Explorers) Program.
Lui, A. T. Y.↗
Lip leakage flow simulation for the Gravity Probe B gas spinup
The lip leakage flow for the Gravity Probe B (GP-B) gas spinup system is investigated using a particle simulation code on the Connection Machine (PSiCM). Particle simulation is employed because the flow conditions are in the transition regime between continuum and free molecule where particle methods are of greatest use. The dominant flow is Couette in nature and the simulation is first validated through comparison to theoretical results for Couette flow in the transition regime. A GP-B type geometry is then simulated and results are presented for two conditions, those corresponding to near the inlet and near the outlet of the spinup channel. Comparison to experiment is not made because experimental data is not yet available.
Dagum, Leonardo↗
Simulations of ICRF Heating for SPARC during First Campaign and Primary Reference-Like Discharge using the Stix Code
High magnetic field tokamaks, like SPARC, rely on ion cyclotron radio frequency heating (ICRF) to reach fusion relevant temperatures. The SPARC tokamak will have 14 ICRF antennas in 7 toroidal locations delivering > 20 MW of power to the plasma. New capabilities with the full wave cold plasma solver, Stix, now allow for resolving the wave-particle resonances using lower order thermal corrections to capture core absorption of Landau damping and ion resonances in devices like SPARC. Favorable comparisons to the TORIC codes give confidence in the single pass absorption of this model to accurately capture the strength of edge interactions of the RF. Using this new dielectric formulation in the Stix code, simulations of the 2D poloidal cross section of SPARC are completed for the first campaign and primary reference-like discharges (PRD-like). A scan of the minority ion concentrations of helium-3 is performed and shows the expected behavior that as the helium-3 decreases the amount of single pass absorption also decreases which is seen in both scenarios. Additionally, both scenarios show only slight differences in single-pass absorption for the range of 3% to 5% helium-3 allowing for more flexibility in experiments. This study also highlights the differences between the first campaign and PRD-like with the first campaign discharges showing much more multi-pass absorption and an effect of confining the wave to a smaller portion of the cross-section due to the fast wave cut-off. This latter result suggests that far-field sheath rectification at the high-field side would be minimal for the first campaign scenario.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Concerning the Perturbations Generated by Flux Transfer Events
The interaction of the solar wind with the Earth's magnetosphere is often highly unsteady. Bursts of magnetic reconnection at multiple locations on the dayside equatorial magnetopause generate flux transfer events or FTEs: twisted ropes of interconnected magnetosheath and magnetospheric magnetic field lines. Once formed, the events move antisunward, displacing and perturbing the ambient media. This talk explores the perturbations predicted by both global numerical simulations and analytical models. Results from global hybrid code simulation confirm the predictions of analytical models indicating that the events generate standing forward slow mode waves as their speeds relative to the magnetosheath flow approach the Alfven velocity. Geometric considerations lead to the conclusion that events generated by component reconnection on the dayside magnetopause move poleward and exhibit strong signatures during intervals of southward IMF orientation, but move towards the flanks and exhibit weak signatures during intervals of northward IMF orientation. Changing event orientations and magnetosheath/magnetospheric magnetic field orientation can enhance the amplitudes of events reaching the flanks. Although the orientations of events on the flanks inferred from multispacecraft timing techniques are consistent with the predictions of the component reconnection model, occurrence patterns versus latitude and IMF orientation require an explanation in terms of both the component and antiparallel reconnection models.
Sibeck, David G.↗
Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields
An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.
Fahrenthold, Eric P.↗
Documentation of the detailed radiation property data for the radiation-ablation code RASLE
This report is a documentation of the necessary radiation property input data for the radiating shock layer simulation code RASLE. The tabulated data are required to simulate systems which are composed of oxygen, nitrogen, carbon, hydrogen, and silicon. These data are needed to compute the flowfield effects of many practical ablative, hypersonic vehicle heat shield materials. A brief outline description is provided for the RASLE code. A more detailed discussion is provided for the RASLE code non-grey gas spectral radiation model. This model is related to the required radiation property data which are tabulated at the end of the report. Other correlations needed for the RASLE simulations are not discussed, since these are automatically included in the program and no input data are required.
Henline, William D.↗
Lidar performance analysis
Section 1 details the theory used to build the lidar model, provides results of using the model to evaluate AEOLUS design instrument designs, and provides snapshots of the visual appearance of the coded model. Appendix A contains a Fortran program to calculate various forms of the refractive index structure function. This program was used to determine the refractive index structure function used in the main lidar simulation code. Appendix B contains a memo on the optimization of the lidar telescope geometry for a line-scan geometry. Appendix C contains the code for the main lidar simulation and brief instruction on running the code. Appendix D contains a Fortran code to calculate the maximum permissible exposure for the eye from the ANSI Z136.1-1992 eye safety standards. Appendix E contains a paper on the eye safety analysis of a space-based coherent lidar presented at the 7th Coherent Laser Radar Applications and Technology Conference, Paris, France, 19-23 July 1993.
Spiers, Gary D.↗
Light Water Reactor Sustainability Program: Use of Time Distributions to Predict Operator Procedure Performance in Dynamic Human Reliability Analysis
The Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) framework affords software capable of conducting human reliability analysis (HRA) using a dynamic approach built around operating procedures (OPs) from nuclear power plants (NPPs). Previous HUNTER reports document the development of this software tool, the coupling of HUNTER to the simulator code, the collection of operator performance data by using simulators to calibrate HUNTER models, and linking HUNTER to probabilistic risk assessment (PRA) software. The present report largely addresses two topics. The first is a new function in HUNTER called the HUNTER Procedure Performance Predictor (P3). HUNTER P3 uses HUNTER’s built in Monte Carlo tools featuring human performance variability to identify potential error traps in procedures. The second topic is time distribution analysis to generate time inputs for dynamic HRA. The current analysis was performed to investigate time distributions for task primitives, which are the minimum task unit of analysis used in dynamic HRA modeling. Using the time distribution data, the elapsed time for human actions in an extended loss of AC power (ELAP) scenario is then investigated. Time data and prediction are essential for modeling procedure performance.
99 GENERAL AND MISCELLANEOUS↗
Accelerating lattice gauge theory studies with Agentic AI
Lattice gauge theory research, with its computationally intensive simulations and complex multi‑stage workflows, is well positioned to benefit from agentic AI systems. We demonstrate how such tools can support key components of lattice gauge theory research, including novel simulation code development using standard LQCD frameworks, HPC job orchestration, simulation data analysis, and expert‑guided tuning of algorithmic parameters such as Hasenbusch mass preconditioning and multigrid solvers. Our results show that agentic AI can reduce manual effort, improve productivity, and accelerate the research cycle while maintaining essential human oversight.
Ayyar, Venkitesh [Fermilab]↗
NASA Glenn Research Center UEET (Ultra-Efficient Engine Technology) Program: Agenda and Abstracts
Topics discussed include: UEET Overview; Technology Benefits; Emissions Overview; P&W Low Emissions Combustor Development; GE Low Emissions Combustor Development; Rolls-Royce Low Emissions Combustor Development; Honeywell Low Emissions Combustor Development; NASA Multipoint LDI Development; Stanford Activities In Concepts for Advanced Gas Turbine Combustors; Large Eddy Simulation (LES) of Gas Turbine Combustion; NASA National Combustion Code Simulations; Materials Overview; Thermal Barrier Coatings for Airfoil Applications; Disk Alloy Development; Turbine Blade Alloy; Ceramic Matrix Composite (CMC) Materials Development; Ceramic Matrix Composite (CMC) Materials Characterization; Environmental Barrier Coatings (EBC) for Ceramic Matrix Composite (CMC) Materials; Ceramic Matrix Composite Vane Rig Testing and Design; Ultra-High Temperature Ceramic (UHTC) Development; Lightweight Structures; NPARC Alliance; Technology Transfer and Commercialization; and Turbomachinery Overview; etc.
Manthey, Lri↗
An Integrated Approach to Modeling Solar Electric Propulsion Vehicles During Long Duration, Near-Earth Orbit Transfers
Recent NASA interest in utilizing solar electronic propulsion (SEP) technology to transfer payloads, e.g. from low-Earth orbit (LEO) to higher energy geostationary-Earth orbit (GEO) or to Earth escape, has necessitated the development of high fidelity SEP vehicle models and simulations. These models and simulations need to be capable of capturing vehicle dynamics and sub-system interactions experienced during the transfer trajectories which are typically accomplished with continuous-burn (potentially interrupted by solar eclipse), long duration "spiral out" maneuvers taking several months or more to complete. This paper presents details of an integrated simulation approach achieved by combining a high fidelity vehicle simulation code with a detailed solar array model. The combined simulation tool gives researchers the functionality to study the integrated effects of various vehicle sub-systems (e.g. vehicle guidance, navigation and control (GN&C), electric propulsion system (EP)) with time varying power production. Results from a simulation model of a vehicle with a 50 kW class SEP system using the integrated tool are presented and compared to the results from another simulation model employing a 50 kW end-of-life (EOL) fixed power level assumption. These models simulate a vehicle under three degree of freedom dynamics (i.e. translational dynamics only) and include the effects of a targeting guidance algorithm (providing a "near optimal" transfer) during a LEO to near Earth escape (C (sub 3) = −2.0 km (sup 2) / sec (sup −2) spiral trajectory. The presented results include the impact of the fully integrated, time-varying solar array model (e.g. cumulative array degradation from traversing the Van Allen belts, impact of solar eclipses on the vehicle and the related temperature responses in the solar arrays due to operating in the Earth's thermal environment, high fidelity array power module, etc.); these are used to assess the impact on vehicle performance (i.e. propellant consumption) and transit times.
vehicle↗
Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on MIMD
A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed.
Massively Parallel Supercomputers↗
Ray-Tracing Simulation and In-Orbit Performance of the ASTRO-H Hard X-Ray Telescope (HXT)
A ray-trace simulation code for the Hard X-ray Telescope (HXT) on board the Hitomi (ASTRO-H) satellite is beingdeveloped. The half power diameter and effective area simulated based on the code are consistent with groundmeasurements within 10%. The HXT observed the pulsar wind nebula G21.5-0.9 for 105 ksec. We confirmed that theencircled energy function and the half power diameter obtained from the data are consistent with the ground measurements.
Matsumoto, H.↗