Search NASA⌕ Search

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.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Adaptive time scaling for accelerating explicit finite element analysis

A method for accelerating an explicit finite element analysis (FEA) simulation of a modeled system or process includes performing an initial iteration of the FEA simulation according to a baseline time interval via an FEA computing network, and calculating a criteria ratio of a predetermined set of scaling criteria for the modeled system or process. The method includes determining a time-scaling factor using the criteria ratio via the FEA computing network as a function of the criteria ratio, and then applying the time-scaling factor to the baseline time interval to generate a scaled time interval. The scaled time interval accelerates simulation time of the FEA simulation. The method includes performing a subsequent iteration of the explicit FEA simulation at the scaled time interval using the FEA computing network. The process continues for subsequent iterations, with the time-scaling factor adapting with each iteration.

Chen, Jian↗
DOE OSTI2017

The use of BMAD in simulating transverse and longitudinal dynamics in RHIC

In accelerator physics, models of a given machine are used to predict the behaviors of the beam, magnets, and radiofrequency cavities. The use of the computational model has become wide spread to ease the development period of the accelerator lattice. There are various programs that are used to create lattices and run simulations of both transverse and longitudinal beam dynamics. The programs include Methodical Accelerator Design(MAD) MAD8, MADX, Zgoubi, Polymorphic Tracking Code (PTC), and many others. In this discussion the BMAD (Baby Methodical Accelerator Design) is presented as an additional tool in creating and simulating accelerator lattices for the study of beam dynamics in the Relativistic Heavy Ion Collider (RHIC).

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

On the utility of GPU accelerated high-order methods for unsteady flow simulations: A comparison with industry-standard tools

First- and second-order accurate numerical methods, implemented for CPUs, underpin the majority of industrial CFD solvers. Whilst this technology has proven very successful at solving steady-state problems via a Reynolds Averaged Navier–Stokes approach, its utility for undertaking scale-resolving simulations of unsteady flows is less clear. High-order methods for unstructured grids and GPU accelerators have been proposed as an enabling technology for unsteady scale-resolving simulations of flow over complex geometries. In this study we systematically compare accuracy and cost of the high-order Flux Reconstruction solver PyFR running on GPUs and the industry-standard solver STAR-CCM+ running on CPUs when applied to a range of unsteady flow problems. Specifically, we perform comparisons of accuracy and cost for isentropic vortex advection (EV), decay of the Taylor–Green vortex (TGV), turbulent flow over a circular cylinder, and turbulent flow over an SD7003 aerofoil. We consider two configurations of STAR-CCM+: a second-order configuration, and a third-order configuration, where the latter was recommended by CD-adapco for more effective computation of unsteady flow problems. Results from both PyFR and STAR-CCM+ demonstrate that third-order schemes can be more accurate than second-order schemes for a given cost e.g. going from second- to third-order, the PyFR simulations of the EV and TGV achieve 75× and 3× error reduction respectively for the same or reduced cost, and STAR-CCM+ simulations of the cylinder recovered wake statistics significantly more accurately for only twice the cost. Moreover, advancing to higher-order schemes on GPUs with PyFR was found to offer even further accuracy vs. cost benefits relative to industry-standard tools.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Study on efficiency of time computation in x-ray imaging simulation base on Monte Carlo algorithm using graphics processing unit

Monte Carlo (MC) is one of the powerful techniques for simulation in x-ray imaging. MC method can simulate the radiation transport within matter with high accuracy and provides a natural way to simulate radiation transport in complex systems. One of the codes based on MC algorithm that are widely used for radiographic images simulation is MC-GPU, a codes developed by Andrea Basal. This study was aimed to investigate the time computation of x-ray imaging simulation in GPU (Graphics Processing Unit) compared to a standard CPU (Central Processing Unit). Furthermore, the effect of physical parameters to the quality of radiographic images and the comparison of image quality resulted from simulation in the GPU and CPU are evaluated in this paper. The simulations were run in CPU which was simulated in serial condition, and in two GPU with 384 cores and 2304 cores. In simulation using GPU, each cores calculates one photon, so, a large number of photon were calculated simultaneously. Results show that the time simulations on GPU were significantly accelerated compared to CPU. The simulations on the 2304 core of GPU were performed about 64 -114 times faster than on CPU, while the simulation on the 384 core of GPU were performed about 20 – 31 times faster than in a single core of CPU. Another result shows that optimum quality of images from the simulation was gained at the history start from 10{sup 8} and the energy from 60 Kev to 90 Kev. Analyzed by statistical approach, the quality of GPU and CPU images are relatively the same.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Fast and self-learning indoor airflow simulation based on in situ adaptive tabulation

Fast simulation for stratified indoor airflow distributions is desired for various applications, such as design of advanced indoor environments, emergency management, and coupled annual energy simulation for buildings with stratified air distributions. Reduced order models trained by pre-computed computational fluid dynamics results are fast, but their prediction may be inaccurate when applied for conditions outside the training domain. To overcome this limitation, we propose in this paper a fast and self-learning model based on an in situ adaptive tabulation (ISAT) algorithm, which is trained by a fast fluid dynamics (FFD) model as an example. The idea is that the ISAT will retrieve the solutions from an existing data set if the estimated prediction error is within a pre-defined tolerance. Otherwise, the ISAT will execute the FFD simulation, which is accelerated by running in parallel on a graphics processing unit, for a full-scale simulation. This paper systematically investigates the feasibility of the ISAT for indoor airflow simulations by presenting the ISAT-FFD implementation alongside results related to its overall performance. Using a stratified indoor airflow as an example, we evaluated how the training time of ISAT was impacted by four factors (training methods, error tolerances, number of inputs, and number of outputs). Then we demonstrated that a trained ISAT model can predict the key information for inputs both inside and outside the training domain. The ISAT was able to answer query points both inside and close to training domain using retrieve actions within a time less than 0.001 s for each query. In conclusion, we provided suggestions for using the ISAT for building applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Three-Dimensional Particle-in-Cell Analysis of Mechanical Design Tolerances of the Scorpius Injector With Warp

Here, this article examines the mechanical design tolerances for the Scorpius injector using a 3-D particle-in-cell (PIC) code. In order to meet its objective as a high-fidelity, multipulse flash X-ray source, the accelerator has requirements for the electron beam in order to minimize the spot size of the beam focused onto a foil target to create high fluence X-rays at energies greater than 10 MeV. The requirements for alignment of the cathode, anode beam pipe, focusing solenoids, and transport solenoids specified using Trak and AMBER are independently verified using a 3-D electrostatic/magnetostatic model of the Scorpius injector, including transverse magnetic field errors that arise from mechanical alignment errors. It is predicted that the injector will meet the performance requirements for radiographic operation and has a design margin to correct for potentially unforeseen operating conditions.

42 ENGINEERING↗
DOE OSTI2014

High-gradient two-beam electron accelerator

The main goal for this project was to design, build, and evaluate a detuned-cavity, collinear, two-beam accelerator structure. Testing was to be at the Yale University Beam Physics Laboratory, under terms of a sub-grant from Omega-P to Yale. Facilities available at Yale for this project include a 6-MeV S-band RF gun and associated beam line for forming and transporting a ~1 A drive beam , a 300 kV beam source for use as a test beam, and a full panoply of laboratory infrastructure and test equipment. During the first year of this project, availability and functionality of the 6-MeV drive beam and 300 kV test beam were confirmed, and the beam line was restored to a layout to be used with the two-beam accelerator project. Major efforts during the first year were also focused on computational design and simulation of the accelerator structure itself, on beam dynamics, and on beam transport. Effort during the second year was focussed on building and preparing to test the structure, including extensive cold testing. Detailed results from work under this project have been published in twelve archival journal articles, listed in Section IV of the technical report.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

A Massively Parallel Solver for the Mechanical Harmonic Analysis of Accelerator Cavities

ACE3P is a 3D massively parallel simulation suite that developed at SLAC National Accelerator Laboratory that can perform coupled electromagnetic, thermal and mechanical study. Effectively utilizing supercomputer resources, ACE3P has become a key simulation tool for particle accelerator R and D. A new frequency domain solver to perform mechanical harmonic response analysis of accelerator components is developed within the existing parallel framework. This solver is designed to determine the frequency response of the mechanical system to external harmonic excitations for time-efficient accurate analysis of the large-scale problems. Coupled with the ACE3P electromagnetic modules, this capability complements a set of multi-physics tools for a comprehensive study of microphonics in superconducting accelerating cavities in order to understand the RF response and feedback requirements for the operational reliability of a particle accelerator. (auth)

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

MULTI-SHELL MAGNETIC TWISTERS AS A NEW MECHANISM FOR CORONAL HEATING AND SOLAR WIND ACCELERATION

We perform numerical simulations of impulsively generated Alfvén waves in an isolated photospheric flux tube and explore the propagation of these waves along such magnetic structure that extends from the photosphere, where these waves are triggered, to the solar corona, and we analyze resulting magnetic shells. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical model and specifying the curved magnetic field lines that constitute the magnetic flux tube that is rooted in the solar photosphere. The evolution of the solar atmosphere is described by 3D, ideal MHD equations that are numerically solved by the FLASH code. Our numerical simulations reveal, based on the physical properties of the multi-shell magnetic twisters and the amount of energy and momentum associated with them, that these multi-shell magnetic twisters may be responsible for the observed heating of the lower solar corona and for the formation of solar wind. Moreover, it is likely that the existence of these twisters can be verified by high-resolution observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Harnessing the power of gradient-based simulations for multi-objective optimization in particle accelerators

Abstract Particle accelerator operation requires simultaneous optimization of multiple objectives. Multi-objective optimization (MOO) is particularly challenging due to trade-offs between the objectives. Evolutionary algorithms, such as genetic algorithms (GAs), have been leveraged for many optimization problems, however, they do not apply to complex control problems by design. This paper demonstrates the power of differentiability for solving MOO problems in particle accelerators using a deep differentiable reinforcement learning (DDRL) algorithm. We compare the DDRL algorithm with model-free reinforcement learning (MFRL), GA, and Bayesian optimization (BO) for simultaneous optimization of heat load and trip rates in the continuous electron beam accelerator facility. The underlying problem enforces strict constraints on both individual states and actions as well as cumulative (global) constraints on energy requirements of the beam. Using historical accelerator data, we develop a physics-based surrogate model which is differentiable and allows for back-propagation of gradients. The results are evaluated in the form of a Pareto-front with two objectives. We show that the DDRL outperforms MFRL, BO, and GA on high dimensional problems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Kinetic Simulations of the Kruskal–Schwarzschild Instability in Accelerating Striped Outflows: Dynamics and Energy Dissipation

Astrophysical relativistic outflows are launched as Poynting-flux dominated, yet the mechanism governing efficient magnetic dissipation, which powers the observed emission, is still poorly understood. We study magnetic energy dissipation in relativistic "striped" jets, which host current sheets separating magnetically dominated regions with opposite field polarity. The effective gravity force g in the rest frame of accelerating jets drives the Kruskal–Schwarzschild instability (KSI), a magnetic analog of the Rayleigh–Taylor instability. By means of 2D and 3D particle-in-cell simulations, we study the linear and nonlinear evolution of the KSI. The linear stage is well described by linear stability analysis. The nonlinear stages of the KSI generate thin (skin-depth-thick) current layers, with length comparable to the dominant KSI wavelength. There, the relativistic drift-kink mode and the tearing mode drive efficient magnetic dissipation. The dissipation rate can be cast as an increase in the effective width Δ eff of the dissipative region, which follows dΔ eff /dt ≃ 0.05 $\sqrt{Δ_{eff} g}$. Our results have important implications for the location of the dissipation region in gamma-ray burst and active galactic nuclei jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Ionization injection of highly-charged copper ions for laser driven acceleration from ultra-thin foils

Laser-driven ion acceleration is often analyzed assuming that ionization reaches a steady state early in the interaction of the laser pulse with the target. This assumption breaks down for materials of high atomic number for which the ionization occurs concurrently with the acceleration process. Using particle-in-cell simulations, we have examined acceleration and simultaneous field ionization of copper ions in ultra-thin targets (20-150 nm thick) irradiated by a laser pulse with intensity 1 × 10 21 W/cm 2 . At this intensity, the laser pulse drives strong electric fields at the rear side of the target that can ionize Cu to charge states with valence L-shell or full K-shell. The highly-charged ions are produced only in a very localized region due to a significant gap between the M- and L-shells' ionization potentials and can be accelerated by strong, forward-directed sections of the field. Such an "ionization injection" leads to well-pronounced bunches of energetic, highly-charged ions. We also find that for the thinnest target (20 nm) a push by the laser further increases the ion energy gain. Thus, the field ionization, concurrent with the acceleration, offers a promising mechanism for the production of energetic, high-charge ion bunches.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

MOOSE: A Modular Platform for Fission and Fusion Multiphysics

The Multiphysics Object-Oriented Simulation Environment (MOOSE) Framework, as well as MOOSE-based simulation tools, have accelerated the development of fission energy and advanced reactor technologies through the United States Department of Energy, Office of Nuclear Science, Nuclear Energy Advanced Modeling & Simulation (NEAMS) Program. MOOSE contains a complete platform of multiphysics simulation capabilities, capable of running on massively parallel systems, and is developed in an open-source manner with great attention paid to high-quality software quality assurance practices. This overall approach could greatly benefit the fusion energy community, which requires rapid design iteration and improvement in order to facilitate the successful development of fusion as an alternative energy source to fossil fuels. In the first half of this talk, applications of MOOSE and MOOSE-based tools for advanced reactor designs will be showcased, as well as MOOSE ecosystem infrastructure (such as the NEAMS Virtual Test Bed) that enables and accelerates fission reactor design. In the second half, a discussion of how the MOOSE approach to modeling and simulation is currently being applied internationally in fusion energy research and development at the United Kingdom Atomic Energy Authority will be discussed, and ongoing/future domestic research efforts will be highlighted.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Radiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration

Here, we use the public code EBHLIGHT to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

High-level support for simulations in astro- and elementary particle physics

To fully exploit the potential of modern supereomputing systems, when performing particle and astrophysics simulations, the special user support unite called simulation laboratory “Elementary- and Astro-Particle Physics” (Sim-Lab E&A Particle) is established at Steinbuch Centre for Computing of the Karlsruhe Institute for Technology. SimLab is providing advanced support to scientific groups in developing the simulation software and porting them into up-to-date super-computing infrastructures. Some details of how the SimLab is governed and cooperation with developers and users of codes are presented in this work in progress paper including examples about code THiSMPI for simulation of particle accelerations in supernovae shock fronts and the KB3D code for solving quantum kinematic equations that could be used in simulations of heavy ion collisions or laser excited semiconductor systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

SU-E-T-244: Designing Low-Z Targets To Enhance Surface Dose: A Monte Carlo Simulation

Purpose: Recent developments in The Varian Truebeam linac platform allows for the introduction of low-Z targets into the beam line for the imaging purposes. We have proposed using a low-Z target for radiation therapy purposes to enhance the surface dose during radiation treatment. The target arm of the Varian Truebeam accelerator consists of multiple targets with are linearly translated into the beam line. We have designed two Low-Z targets made of carbon: 1) a step target consisting of three steps of 15%, 30% and 60% CSDA range for 2.5 MeV electrons Figure 1a; 2) and a ramp target, an incline plane 2cm long with thicknesses ranging from 0% to 60% CSDA range, Figure 1b. The purpose of this work will determine the spectral characteristics of these target designs and determine if they have practical clinical applications for enhancing surface dose. Methods: To calculate the spectral characteristics of these targets, a standard Monte Carlo model of a Varian Clinac accelerator was used. Simulations were performed with a carbon step target, and a carbon ramp target, located at the same position as the electron foil in the rotating carousel. Simulations were carried out using a 2.5 MeV electron beam. Results: The step target design produced spectral characteristics which were similar to spectral model using a single disk target of the same thickness. The ramp target provides a means to have positional variation of the spectral components of the beam, however, the electron component as 60% CSDA us much broader than the step target. Conclusion: The carbon step-target provides a spectral distribution which is similar to a carbon disk of comparable thickness. The spectral distribution from the ramp-target can be modified as a function of position to provide a wide range of low energy electrons for surface dose enhancement.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Advanced Modeling of Conventional Particle Accelerators

SciDAC-5 goals: Deliver particle accelerator and beam simulations tools that go beyond the current state of the art, up to the realization of virtual twins of particle accelerators, enabling design and modeling of particle accelerators at unprecedented speed, levels of accuracy, and realism; and apply these tools to key accelerator facilities relevant to DOE HEP (such as PIP-II/DUNE, FACET-II).

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

LOFAR, VLA, and Chandra observations of the Toothbrush Galaxy Cluster

We present deep LOFAR observations between 120 {181 MHz of the `Toothbrush' (RX J0603.3+4214), a cluster that contains one of the brightest radio relic sources known. Our LOFAR observations exploit a new and novel calibration scheme to probe 10 times deeper than any previous study in this relatively unexplored part of the spectrum. The LOFAR observations, when combined with VLA, GMRT, and Chandra X-ray data, provide new information about the nature of cluster merger shocks and their role in re-accelerating relativistic particles. We derive a spectral index of α = -0:8±0:1 at the northern edge of the main radio relic, steepening towards the south to α ≈ -2. The spectral index of the radio halo is remarkably uniform (α = -1:16, with an intrinsic scatter of ≤ 0:04). The observed radio relic spectral index gives a Mach number of M = 2:8 +0:5 -0:3 , assuming diffusive shock acceleration (DSA). However, the gas density jump at the northern edge of the large radio relic implies a much weaker shock (M≈1:2, with an upper limit ofM≈1:5). The discrepancy between the Mach numbers calculated from the radio and X-rays can be explained if either (i) the relic traces a complex shock surface along the line of sight, or (ii) if the radio relic emission is produced by a re-accelerated population of fossil particles from a radio galaxy. Our results highlight the need for additional theoretical work and numerical simulations of particle acceleration and re-acceleration at cluster merger shocks.

79 ASTRONOMY AND ASTROPHYSICS↗