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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Non-thermal electron acceleration in low Mach number collisionless shocks. II. Firehose-mediated Fermi acceleration and its dependence on pre-shock conditions

Electron acceleration to non-thermal energies is known to occur in low Mach number (M{sub s} ≲ 5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M{sub s} = 3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (≳ 20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Flux Creep in a Bi-2212 Rutherford Cable for Particle Accelerator Applications

Bi-2212 superconducting cables are being considered for use in the high field magnets needed for the next generation of particle accelerators. Magnetization in these cables and the decay of that magnetization lead to field error and field-error drift, respectively, which need to be compensated. To study this, a segment of the winding pack was extracted from a racetrack coil made from Bi-2212 Rutherford cable. Using a Hall probe measurement technique, we measured the response of the cable's magnetization and its magnetization decay to changes in the applied magnetic field. The effect of adjustments to the cycling of the magnetic field was studied, intended to simulate the preinjection cycles of an accelerator magnet. Three M vs. H loops were constructed by sweeping the magnetic field applied to the sample from 0 to 2.5 T, then to a preinjection field 'x' (where x = 0, 0.25, 0.75 T), and finally up to 1 T. The applied field was then held at 1 T for 1500 s, and the magnetization decay was measured. Finally, the decay was found to vary from 8% to 14% after 1500 s, depending on the preinjection field cycle.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

THE MOST INTENSIVE GAMMA-RAY FLARE OF QUASAR 3C 279 WITH THE SECOND-ORDER FERMI ACCELERATION

The very short and bright flare of 3C 279 detected with Fermi-Large Area Telescope in 2013 December is tested by a model with stochastic electron acceleration by turbulences. Our time-dependent simulation shows that the very hard spectrum and asymmetric light curve are successfully reproduced by changing only the magnetic field from the value in the steady period. The maximum energy of electrons drastically grows with the decrease of the magnetic field, which yields a hard photon spectrum as observed. Rapid cooling due to the inverse-Compton scattering with the external photons reproduces the decaying feature of the light curve.The inferred energy density of the magnetic field is much less than the electron and photon energy densities. The low magnetic field and short variability timescale are unfavorable for the jet acceleration model from the gradual Poynting flux dissipation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2016

Numerical investigation of initial condition effects on Rayleigh-Taylor instability with acceleration reversals

The influence of initial conditions on miscible incompressible baroclinically driven Rayleigh-Taylor instability undergoing nonuniform acceleration is explored computationally using an implicit large eddy simulation (ILES) technique. We consider here the particular case of evolution during multiple reversals of acceleration direction, where the flow is alternately statically stable or unstable. In the unstable phase, the flow is driven by the baroclinic release of potential energy, whereas in the stable phase, work is done against the density stratification with the energy exchange taking place by wavelike mechanisms. These dynamics are fundamentally different; here, we track the evolution of volume-averaged turbulent statistics that are most sensitive to changes in the distribution of spectral power and bandwidth of the initial conditions as the flow alternates between dynamical regimes due to acceleration reversal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Celeritas: Accelerating Geant4 with GPUs

Celeritas [1] is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications (specifically, High Lumi- nosity Large Hadron Collider (HL-LHC) simulation) on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a GPU-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.3, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling integration of Celeritas into high energy physics (HEP) experimental frameworks such as CMSSW. On the Summit supercomputer, Celeritas performs EM physics between 6 and 32 faster using the machine’s Nvidia GPUs compared to using only CPUs. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.8–2.3× on GPU and 1.2× on CPU.

Johnson, Seth R.↗
DOE OSTI2023

Polarized electron injector for positron production at CEBAF

As a part of the effort to expand the capabilities of CE-BAF 12 GeV (Continuous Electron Beam Accelerator Facility) at Jefferson Lab, the addition of a polarized positron source is considered. This capability would provide acceleration of high duty-cycle polarized posi-trons, with spin >60% polarization, through the same main CEBAF accelerator machine with appropriate mag-net field reversals and linac phasing to the four CEBAF experimental halls. To produce this positron beam, a high average current (3-10 mA) highly polarized electron beam with energy of 100 ? 150 MeV is required at the positron source target. The focus of this paper is the de-sign of that polarized electron beam injector. We will describe the production and delivery of a >3 mA highly polarized electron beam. We will discuss different aspects of the design, the photocathode gun, beam dynamics simulation results, spin manipulation, bunching and accelerating process and final electron beam parameters.

Kazimi, R.↗
DOE OSTI2014

Comparative study of beam losses and heat loads reduction methods in MITICA beam source

In negative ion electrostatic accelerators a considerable fraction of extracted ions is lost by collision processes causing efficiency loss and heat deposition over the components. Stripping is proportional to the local density of gas, which is steadily injected in the plasma source; its pumping from the extraction and acceleration stages is a key functionality for the prototype of the ITER Neutral Beam Injector, and it can be simulated with the 3D code AVOCADO. Different geometric solutions were tested aiming at the reduction of the gas density. The parameter space considered is limited by constraints given by optics, aiming, voltage holding, beam uniformity, and mechanical feasibility. The guidelines of the optimization process are presented together with the proposed solutions and the results of numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transformational Regional-Scale Earthquake Simulations with the DOE EarthQuake SIMulation Exascale Framework

Earthquakes present worldwide risk to economic and human safety. The 2023 earthquakes in Turkiye provided a reminder of the potential for catastrophic consequences with 50,700 deaths and 15.7 million people affected. The ability to predict ground motions and infrastructure damage for earthquakes continues to be a challenging problem for scientists and engineers. Until now, estimates of ground motions have been performed empirically by looking at sparse data from past earthquakes. This approach can provide statistical information on intensity amplitudes but cannot inform site-specific ground motions essential to developing the most effective resilience. Interest has grown in large-scale computational models to simulate earthquakes at regional scale. The U.S. Department of Energy EarthQuake SIMulation (EQSIM) framework was developed for regional-scale earthquake simulations at unprecedented fidelity, taking advantage of emerging GPU-accelerated systems. This article describes the EQSIM workflow and demonstrates regional-scale simulations with the new computational capability available to scientists in their quest to mitigate future disasters.

58 GEOSCIENCES↗
DOE OSTI2019

Determining the Dominant Acceleration Mechanism during Relativistic Magnetic Reconnection in Large-scale Systems

While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, in this work we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, as it opens up the possibility of modeling the energetic spectra without resolving microscopic diffusion regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Direct laser acceleration of electrons in the plasma bubble by tightly focused laser pulses

We present an analytical theory that reveals the importance of the longitudinal laser electric field in the course of the resonant acceleration of relativistic electrons by a tightly confined laser beam. It is shown that this laser field component always counteracts the transverse one and effectively decreases the final energy gain of electrons via the direct laser acceleration (DLA) mechanism. This effect is demonstrated by carrying out particle-in-cell simulations of the DLA of the electrons injected into the accelerating phase of the plasma wake. It is shown that the electron energy gain from the wakefield is substantially compensated by the quasiresonant energy loss to the longitudinal laser field component. The analytically obtained scalings and estimates are in good agreement with the results of the numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Numerical simulations used for a validity check on the laser induced photo-detachment diagnostic method in electronegative plasmas

Laser photo-detachment is used as a method to measure or determine the negative ion density and temperature in electronegative plasmas. In essence, the method consists of producing an electropositive channel (negative ion free region) via pulsed laser photo-detachment within an electronegative plasma bulk. Electrostatic probes placed in this channel measure the change in the electron density. A second pulse might be used to track the negative ion recovery. From this, the negative ion density and temperature can be determined. We study the formation and relaxation of the electropositive channel via a two-dimensional Particle-In-Cell/Mote Carlo collision model. The simulation is mainly carried out in a Hydrogen plasma with an electronegativity of α = 1, with a parametric study for α up to 20. The temporal and spatial evolution of the plasma potential and the electron densities shows the formation of a double layer (DL) confining the photo-detached electrons within the electropositive channel. This DL evolves into two fronts that move in the opposite directions inside and outside of the laser spot region. As a consequence, within the laser spot region, the background and photo-detached electron energy distribution function relaxes/thermalizes via collisionless effects such as Fermi acceleration and Landau damping. Moreover, the simulations show that collisional effects and the DL electric field strength might play a non-negligible role in the negative ion recovery within the laser spot region, leading to a two-temperature negative ion distribution. The latter result might have important effects in the determination of the negative ion density and temperature from laser photo detachment diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Analysis of rapid increase in the plasma density during the ramp-up phase in a radio frequency negative ion source by large-scale particle simulation

Numerical simulations become useful for the developing RF-ICP (Radio Frequency Inductively Coupled Plasma) negative ion sources. We are developing and parallelizing a two-dimensional three velocity electromagnetic Particle-In-Cell code. The result shows rapid increase in the electron density during the density ramp-up phase. A radial electric field due to the space charge is produced with increase in the electron density and the electron transport in the radial direction is suppressed. As a result, electrons stay for a long period in the region where the inductive electric field is strong, and this leads efficient electron acceleration and a rapid increasing of the electron density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Force field preconditioned ab initio structure relaxation method

We present a general method to accelerate ab initio atomic structural relaxations. In this method, the conjugate gradient ab initio relaxation is preconditioned by a force field relaxation. This force field is constructed on-the-flight with the information of the ab initio forces and the current atomic configuration. At each ab initio relaxation step, the so-constructed force field is used to prerelax the system, with the relaxation direction as the preconditioned direction for the ab initio conjugated gradient method. The force field model and its parameters are rather general making this method applicable for a wide range of systems commonly used in material simulations. We report more than 80 different systems have been tested representing different cases in material simulations. Across the board, we find accelerations mostly over a factor of 2 and for some large systems with a factor of 6 ~ 9 . There is no case where the relaxation becomes worse. The code and lib for this method are provided, which can be used as a plugin in a standard ab initio atomic relaxation code.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Beam dynamics study of a 30 MeV electron linear accelerator to drive a neutron source

An experimental neutron facility based on 32 MeV/18.47 kW electron linac has been studied by means of PARMELA simulation code. Beam dynamics study for a traveling wave constant gradient electron accelerator is carried out to reach the preferential operation parameters (E = 30 MeV, P = 18 kW, dE/E < 12.47% for 99% particles). The whole linac comprises mainly E-gun, pre-buncher, buncher, and 2 accelerating columns. A disk-loaded, on-axis-coupled, 2π/3-mode type accelerating rf cavity is considered for this linac. After numerous optimizations of linac parameters, 32 MeV beam energy is obtained at the end of the linac. As high electron energy is required to produce acceptable neutron flux. The final neutron flux is estimated to be 5 × 10{sup 11} n/cm{sup 2}/s/mA. Future development will be the real design of a 30 MeV electron linac based on S band traveling wave.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

9-D polarized proton transport in the MEIC figure 8 collider ring - first steps

Spin tracking studies in the MEIC figure-8 collider ion ring are presented, based on a very preliminary design of the lattice. They provide numerical illustrations of some of the aspects of the figure-8 concept, including spin-rotator based spin control, and lay out the path towards a complete spin tracking simulation of a figure-8 ring.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

AB INITIO PULSAR MAGNETOSPHERE: THE ROLE OF GENERAL RELATIVITY

It has recently been demonstrated that self-consistent particle-in-cell simulations of low-obliquity pulsar magnetospheres in flat spacetime show weak particle acceleration and no pair production near the poles. We investigate the validity of this conclusion in a more realistic spacetime geometry via general-relativistic particle-in-cell simulations of the aligned pulsar magnetosphere with pair formation. We find that the addition of the frame-dragging effect makes the local current density along the magnetic field larger than the Goldreich–Julian value, which leads to unscreened parallel electric fields and the ignition of a pair cascade. When pair production is active, we observe field oscillations in the open field bundle, which could be related to pulsar radio emission. We conclude that general-relativistic effects are essential for the existence of the pulsar mechanism in low-obliquity rotators.

79 ASTRONOMY AND ASTROPHYSICS↗