SEARCH · Search NASA

Results for “Accelerator Simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2017

WARP

WARP, which can stand for ``Weaving All the Random Particles,'' is a three-dimensional (3D) continuous energy Monte Carlo neutron transport code developed at UC Berkeley to efficiently execute on NVIDIA graphics processing unit (GPU) platforms. WARP accelerates Monte Carlo simulations while preserving the benefits of using the Monte Carlo method, namely, that very few physical and geometrical simplifications are applied. WARP is able to calculate multiplication factors, neutron flux distributions (in both space and energy), and fission source distributions for time-independent neutron transport problems. It can run in both criticality or fixed source modes, but fixed source mode is currently not robust, optimized, or maintained in the newest version. WARP can transport neutrons in unrestricted arrangements of parallelepipeds, hexagonal prisms, cylinders, and spheres. The goal of developing WARP is to investigate algorithms that can grow into a full-featured, continuous energy, Monte Carlo neutron transport code that is accelerated by running on GPUs. The crux of the effort is to make Monte Carlo calculations faster while producing accurate results. Modern supercomputers are commonly being built with GPU coprocessor cards in their nodes to increase their computational efficiency and performance. GPUs execute efficiently on data-parallel problems, but most CPU codes, including those for Monte Carlo neutral particle transport, are predominantly task-parallel. WARP uses a data-parallel neutron transport algorithm to take advantage of the computing power GPUs offer.

Bergmann, RyanM. [University of California - Berke↗

DOE OSTI2017

WARP

Bergmann, Ryan M.↗

DOE OSTI2014

Ion pre-acceleration in fully self-consistent particle-in-cell simulations of supercritical perpendicular reforming shocks in multiple ion species plasmas

Supernova remnant and heliopause termination shock plasmas may contain significant populations of minority heavy ions, with relative number densities n{sub α}/n{sub i} up to 50%. Preliminary kinetic simulations of collisionless shocks in these environments showed that the reformation cycle and acceleration mechanisms at quasi-perpendicular shocks can depend on the value of n{sub α}/n{sub i} . Shock reformation unfolds on ion spatio-temporal scales, requiring fully kinetic simulations of particle dynamics, together with the self-consistent electric and magnetic fields. This paper presents the first set of particle-in-cell simulations for two ion species, protons (n{sub p} ) and α-particles (n{sub α}), with differing mass and charge-to-mass ratios, that spans the entire range of n{sub α}/n{sub i} from 0% to 100%. The interplay between the differing gyro length scales and timescales of the ion species is crucial to the time-evolving phenomenology of the shocks, the downstream turbulence, and the particle acceleration at different n{sub α}/n{sub i} . We show how the overall energization changes with n{sub α}/n{sub i} , and relate this to the processes individual ions undergo in the shock region and in the downstream turbulence, and to the power spectra of magnetic field fluctuations. The crossover between shocks dominated by the respective ion species happens when n{sub α}/n{sub i} = 25%, and minority ion energization is strongest in this regime. Energization of the majority ion species scales with injection energy. The power spectrum of the downstream turbulence includes peaks at sequential ion cyclotron harmonics, suggestive of ion ring-beam collective instability.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Small-angle X-ray scattering characteristics of mouse brain: Planar imaging measurements and tomographic imaging simulations

Small-angle x-ray scattering (SAXS) imaging can differentiate tissue types based on their nanoscale molecular structure. However, characterization of the coherent scattering cross-section profile of relevant tissues is needed to optimally design SAXS imaging techniques for a variety of biomedical applications. Reported measured nervous tissue x-ray scattering cross sections under a synchrotron source have had limited agreement. We report a set of x-ray cross-section measurements obtained from planar SAXS imaging of 1 mm thick mouse brain (APP/PS1 wild-type) coronal slices using an 8 keV laboratory x-ray source. Two characteristic peaks were found at 0.96 and 1.60 nm -1 attributed to myelin. The peak intensities varied by location in the slice. We found that regions of gray matter, white matter, and corpus callosum could be segmented by their increasing intensities of myelin peaks respectively. Measured small-angle x-ray scattering cross sections were then used to define brain tissue scattering properties in a GPU-accelerated Monte Carlo simulation of SAXS computed tomography (CT) using a higher monochromatic x-ray energy (20 keV) to study design trade-offs for noninvasive in vivo SAXS imaging on a small-animal head including radiation dose, signal-to-noise ratio (SNR), and the effect of skull presence on the previous two metrics. Simulation results show the estimated total dose to the mouse head for a single SAXS-CT slice was 149.4 mGy. The pixel SNR was approximately 30.8 for white matter material whether or not a skull was present. In this early-stage proof-of-principle work, we have demonstrated our brain cross-section data and simulation tools can be used to assess optimal instrument parameters for dedicated small-animal SAXS-CT prototypes.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2000

Simulated Space Environmental Testing on Thin Films

An exploratory program has been conducted, to irradiate some mature commercial and some experimental polymer films with radiation simulating certain Earth orbits, and to obtain data about the response of each test film's reflective and tensile properties. Protocols to conduct optimized tests were considered and developed to a "prototype" level during this program. Fifteen polymer film specimens were arranged on a specially designed test fixture. The fixture featured controlled exposure areas, and protected the ends of the samples for later gripping in tensile tests. The fixture featured controlled exposure areas, and protected the ends of the samples for later gripping in tensile tests. The fixture containing the films was installed in a clean vacuum chamber where protons, electrons and solar ultraviolet (UV) radiation could simultaneously irradiate the specimens. Near realtime UV rates were used, whereas proton and electron rates were accelerated appreciably to simulate 5 years in orbit during a two month test. Periodically, the spectral reflectance of each film was measured in situ. After the end of the irradiation, final reflectance measurements were made in situ, and solar absorptance values were derived for each specimen. These samples were then measured in air for thermal emittance and for tensile strength. Most specimens withstood the irradiation intact, but with reduced reflectance (increased solar absorptance). Thermal emittance changed slightly in several materials, as did their tensile strength and elongation at break. Conclusions are drawn about the performance of the films. Simulated testing to an expected 5 year dose of electrons and protons consistent with those expected at L2 and 0.98 AU orbits and 100 equivalent solar hours exposure.

Russell, Dennis A.↗

DOE OSTI2021

Reduced bandwidth Compton photons from a laser-plasma accelerator using tailored plasma channels

It has been demonstrated experimentally that laser plasma accelerators can produce multi-100 MeV electron bunches with a few percent energy spread, and from these electrons, multi-MeV quasi-monoenergetic photons have been demonstrated based on Compton up-scattering from a counter-propagating laser. This offers the potential of a high-quality, narrow-bandwidth, compact, photon source with broad application. The bandwidth of the resulting photons depends directly on the distribution of the electron bunch and is limited, in particular, by the bunch divergence (i.e., the spread in transverse velocity angle). At the same time, the ability to decelerate electrons after scattering is important to source deployment. We describe a series of plasma structures that expand and then collimate the electron bunch, reducing its divergence and thus reducing the bandwidth of the scattered photons while enabling both high performance scattering and deceleration. Furthermore, these plasma structures are demonstrated in simulations of the accelerator system, showing the potential to reach few-percent photon spread which is important for applications using nuclear resonance fluorescence.

42 ENGINEERING↗

DOE OSTI2024

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗

DOE OSTI2023

Quantum computing for fusion energy science applications

This is a review of recent research exploring and extending present-day quantum computing capabilities for fusion energy science applications. We begin with a brief tutorial on both ideal and open quantum dynamics, universal quantum computation, and quantum algorithms. Then, we explore the topic of using quantum computers to simulate both linear and nonlinear dynamics in greater detail. Because quantum computers can only efficiently perform linear operations on the quantum state, it is challenging to perform nonlinear operations that are generically required to describe the nonlinear differential equations of interest. In this work, we extend previous results on embedding nonlinear systems within linear systems by explicitly deriving the connection between the Koopman evolution operator, the Perron–Frobenius evolution operator, and the Koopman–von Neumann evolution (KvN) operator. We also explicitly derive the connection between the Koopman and Carleman approaches to embedding. Extension of the KvN framework to the complex-analytic setting relevant to Carleman embedding, and the proof that different choices of complex analytic reproducing kernel Hilbert spaces depend on the choice of Hilbert space metric are covered in the appendixes. Finally, we conclude with a review of recent quantum hardware implementations of algorithms on present-day quantum hardware platforms that may one day be accelerated through Hamiltonian simulation. We discuss the simulation of toy models of wave–particle interactions through the simulation of quantum maps and of wave–wave interactions important in nonlinear plasma dynamics.

Joseph, I. (ORCID:0000000255400840)↗

NASA NTRS2025

Rheology of Lunar Regolith Simulant Under Varying Gravitational Conditions

Understanding structure-property-process relationship aka rheology of regolith in varying gravity conditions is critical for exploration and future missions. In this work, a framework for studying rheology of lunar regolith simulant in varying gravity conditions (Terrestrial, Lunar and Martian) is described theoretically/numerically, where the driving force of flow of simulants was gravitational acceleration. In the analysis particle-particle interaction forces were included. It was found that the dynamic behavior of granular material is extremely sensitive to external conditions by virtue of the myriad of forces present between particle grains. The results obtained were validated against results obtained in earth and Lunar gravity conditions. The theoretical/numerical and experimental results showed that the complex interaction of these forces can drastically change the dynamics of the material, which is not captured by standard design rules, generally used in industry, for variable gravity applications.

Simulation↗

DOE OSTI2018

Enhanced ion acceleration in transition from opaque to transparent plasmas

Using particle-in-cell simulations, we investigate ion acceleration in the interaction of high intensity lasers with plasmas which transition from opaque to transparent during the interaction process. We show that the highest ion energies are achieved when the laser traverses the target around the peak intensity and re-heats the electron population responsible for the plasma expansion, enhancing the corresponding sheath electric field. This process can lead to an increase of up to 2x in ion energy when compared with the standard Target Normal Sheath Acceleration in opaque targets under the same laser conditions. A theoretical model is developed to predict the optimal target areal density as a function of laser intensity and pulse duration. A systematic parametric scan for a wide range of target densities and thicknesses is performed in 1D, 2D and 3D and shown consistent with the theory and with recent experimental results. Thus, these results open the way for a better optimization of the ion energy in future laser–solid experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

MuSIKAL: Multiphysics Simulations and Knowledge Discovery through AI/ML Technologies

Under the MuSiKAL project, we developed a framework for a coastal digital twin (DT) platform capable of integrating diverse data resources, configuring multiscale model simulations, performing SciML‐accelerated predictions, with applications primarily driven by storm surge and heavily rainfall events impacting the Gulf Coast of the U.S.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Pulsed-focusing recirculating linacs for muon acceleration

Since the muon has a short lifetime, fast acceleration is essential for high-energy applications such as muon colliders, Higgs factories, or neutrino factories. The best one can do is to make a linear accelerator with the highest possible accelerating gradient to make the accelerating time as short as possible. However, the cost of such a single linear accelerator is prohibitively large due to expensive power sources, cavities, tunnels, and related infrastructure. As was demonstrated in the Thomas Jefferson Accelerator Facility (Jefferson Lab) Continuous Electron Beam Accelerator Facility (CEBAF), an elegant solution to reduce cost is to use magnetic return arcs to recirculate the beam through the accelerating RF cavities many times, where they gain energy on each pass. In such a Recirculating Linear Accelerator (RLA), the magnetic focusing strength diminishes as the beam energy increases in a conventional linac that has constant strength quadrupoles. After some number of passes the focusing strength is insufficient to keep the beam from going unstable and being lost. In this project, the use of fast pulsed quadrupoles in the linac sections was considered for stronger focusing as a function of time to allow more successive passes of a muon beam in a recirculating linear accelerator. In one simulation, it was shown that the number of passes could be increased from 8 to 12 using pulsed magnet designs that have been developed and tested. This could reduce the cost of linac sections of a muon RLA by 8/12, where more improvement is still possible. The expense of a greater number of passes and corresponding number of return arcs was also addressed in this project by exploring the use of ramped or FFAG-style magnets in the return arcs. A better solution, invented in this project, is to use combined-function dipole-quadrupole magnets to simultaneously transport two beams of different energies through one magnet string to reduce costs of return arcs by almost a factor of two. A patent application was filed for this invention and a detailed report published in Physical Review Special Topics. A scaled model using an electron beam was developed and proposed to test the concept of a dog bone RLA with combined-function return arcs. The efforts supported by this grant were reported in a series of contributions to particle accelerator conferences that are reproduced in the appendices and summarized in the body of this report.

43 PARTICLE ACCELERATORS↗

DOE OSTI2019

Influence of Longer Range Transfer of Vapor Interface Modified Caging Constraints on the Spatially Heterogeneous Dynamics of Glass-Forming Liquids

Based on the elastically collective nonlinear Langevin equation (ECNLE) theory of bulk relaxation in glass-forming liquids, and our recent ideas of how interface-nucleated modification of caging constraints are spatially transferred into the interior of a thick film, we present a force-based theory for dynamical gradients in thick films with one vapor interface that includes collective elasticity effects. Quantitative applications to the foundational hard sphere fluid and polymer melts of diverse fragilities are presented. In this work, we predict a roughly exponential spatial variation of the total activation barrier which is nonperturbatively modified to similar to 10 particle diameters into the film. This leads, to leading order, to the prediction of a reduced alpha relaxation time gradient of a double-exponential form characterized by a nearly constant spatial decay length, in qualitative accord with simulations. The relaxation acceleration at the surface grows exponentially with increasing packing fraction or decreasing temperature. The rate of increase with cooling of the alpha time strongly weakens upon approaching the vapor interface, with the top two layers remaining liquid-like down to similar to ~80-85% of the bulk glass transition temperature. These spatially heterogeneous changes of the temperature variation of the alpha time result in a large and long-range gradient of the local glass transition temperature. An average interfacial layer thickness relevant to ensemble-averaged dielectric and other experiments is also computed. It is predicted to be rather large at the bulk glass transition temperature, decreasing roughly linearly with heating. Remarkably, to leading order the ratio of this layer-averaged interfacial relaxation time to its bulk analogue is invariant to chemistry, volume fraction, and temperature and of modest absolute magnitude. Spatially inhomogeneous power law decoupling of the alpha relaxation time from its bulk value is predicted, with an effective exponent that decays to zero with distance from the free surface in a nearly exponential manner, trends which are in qualitative accord with recent simulations. Finally, this behavior and the double-exponential alpha time gradient are related and can be viewed as consequences of an effective quasi-universal factorization of the total barrier in films into the product of its bulk-temperature-dependent value times a function solely of location in the film.

36 MATERIALS SCIENCE↗

NASA NTRS1998

[Activities of Psychology Dept., California Univ.]

We have completed two studies during the grant period, with manuscripts published or ready for submission for publication: (1) Dual adaptation and adaptive generalization in the human vestibuloocular reflex and (2) Frequency vs. acceleration specificity in human VOR adaptation. In the 1st study two studies examined the possibility that rotational VOR plasticity is subject to dual adaptation and adaptive generalization. Subjects in the experimental condition were exposed to an altered visual-vestibular environment for about four minutes every day for five consecutive days. The working hours between these testing sessions constituted re-exposure to the normal visual environment. Thus, subjects were repeatedly adapting and re-adapting to both environments which is a condition designed to produce dual adaptation. In each training session a measure of baseline VOR gain was obtained (in the dark). A small laser spot (the only visual stimulus) was systematically moved in the same direction as the subject's head, but by half the angle of rotation (target/head gain = 0.5). This resulted in adaptation values relativized to the non-adapted gain of each subject. These values were then analyzed using an analysis of variance with day and session (within a day) as factors. In the 2nd study human VOR adaption has been assumed to be frequency specific, despite the fact that the semicircular canals are simulated by rotational acceleration and not frequency per se.

Bridgeman, Bruce↗

DOE OSTI2016

Multigroup Radiation-Hydrodynamics with a High-Order, Low-Order Method

Recent efforts at Los Alamos National Laboratory to develop a moment-based, scale-bridging [or high-order (HO)–low-order (LO)] algorithm for solving large varieties of the transport (kinetic) systems have shown promising results. A part of our ongoing effort is incorporating this methodology into the framework of the Eulerian Applications Project to achieve algorithmic acceleration of radiationhydrodynamics simulations in production software. By starting from the thermal radiative transfer equations with a simple material-motion correction, we derive a discretely consistent energy balance equation (LO equation). We demonstrate that the corresponding LO system for the Monte Carlo HO solver is closely related to the original LO system without material-motion corrections. We test the implementation on a radiative shock problem and show consistency between the energy densities and temperatures in the HO and LO solutions as well as agreement with the semianalytic solution. We also test the approach on a more challenging two-dimensional problem and demonstrate accuracy enhancements and algorithmic speedups. This paper extends a recent conference paper by including multigroup effects.

42 ENGINEERING↗

NASA NTRS1963

Effects of gust-induced and maneuvering acceleration stress on pilot-vehicle performance

Piloted motion simulator study of gust induced and maneuvering acceleration stress on pilot performance in low level mission

COCKPIT SIMULATOR↗

DOE OSTI2018

The Roles of Fluid Compression and Shear in Electron Energization during Magnetic Reconnection

Particle acceleration in space and astrophysical reconnection sites is an important unsolved problem in studies of magnetic reconnection. Earlier kinetic simulations have identified several acceleration mechanisms that are associated with particle drift motions. Here, we show that, for sufficiently large systems, the energization processes due to particle drift motions can be described as fluid compression and shear, and that the shear energization is proportional to the pressure anisotropy of energetic particles. By analyzing results from fully kinetic simulations, we show that the compression energization dominates the acceleration of high-energy particles in reconnection with a weak guide field, and the compression and shear effects are comparable when the guide field is 50% of the reconnecting component. Spatial distributions of those energization effects reveal that reconnection exhausts, contracting islands, and island-merging regions are the three most important regions for compression and shear acceleration. This study connects particle energization by particle guiding-center drift motions with that due to background fluid motions, as in the energetic particle transport theory. It provides foundations for building particle transport models for large-scale reconnection acceleration such as those in solar flares.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Start-to-End Simulation of the Drive-Beam Longitudinal Dynamics for Beam-Driven Wakefield Acceleration

Collinear beam-driven wakefield acceleration (WFA) relies on shaped driver beam to provide higher accelerating gradient at a smaller cost and physical footprint. This acceleration scheme is currently envisioned to accelerate electron beams capable of driving free-electron laser *. Start-to-end simulation of drive-bunch beam dynamics is crucial for the evaluation of the design of accelerators built upon WFA. We report the start-to-end longitudinal beam dynamics simulations of an accelerator beamline capable of producing high charge drive beam. The generated wakefield when it passes through a corrugated waveguide results in a transformer ratio of 5. This paper especially discusses the challenges and criteria associated with the generation of temporally-shaped driver beam, including the beam formation in the photoinjector, and the influence of energy chirp control on beam transport stability.

43 PARTICLE ACCELERATORS↗