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,225 records · Page 68

Electron and Proton Acceleration in Trans-relativistic Magnetic Reconnection: Dependence on Plasma Beta and Magnetization

Nonthermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection—where the magnetization σ, defined as the ratio of magnetic energy density to enthalpy density, is ~1—is frequently encountered in such flows. By means of a large suite of two-dimensional particle-in-cell simulations, we investigate electron and proton acceleration in the trans-relativistic regime. Additionally, we focus on the dependence of the electron energy spectrum on σ and the proton β (the ratio of proton thermal pressure to magnetic pressure). We find that the electron spectrum in the reconnection region is nonthermal and can be modeled as a power law. At low β, the slope, p, is independent of β and hardens with increasing σ as $p\simeq 1.8+0.7/\sqrt{\sigma }$. Electrons are primarily accelerated by the nonideal electric field at X-points, either in the initial current layer or in current sheets generated between merging magnetic islands. At higher values of β, the electron power law steepens, and the electron spectrum eventually approaches a Maxwellian distribution for all values of σ. At values of β near β max ≈ 1/4σ, when both electrons and protons are relativistically hot prior to reconnection, the spectra of both species display an additional component at high energies, containing a few percent of particles. Finally, these particles are accelerated via a Fermi-like process by bouncing between the reconnection outflow and a stationary magnetic island

43 PARTICLE ACCELERATORS↗

Parallel Acceleration of Suprathermal Electrons Caused by Whistler‐Mode Hiss Waves

Abstract Suprathermal electrons (~0.1–10 keV) in the inner magnetosphere are usually observed in a 90°‐peaked pitch angle distribution, formed due to the conservation of the first and second adiabatic invariants as they are transported from the plasma sheet. We report a peculiar field‐aligned suprathermal electron (FASE) distribution measured by Van Allen Probes, where parallel fluxes are 1 order of magnitude higher than perpendicular fluxes. Those FASEs are found to be closely correlated with large‐amplitude hiss waves and are observed around the Landau resonant energies. We demonstrate, using quasilinear diffusion simulations, that hiss waves can rapidly accelerate suprathermal electrons through Landau resonance and create the observed FASE population. The proposed mechanism potentially has broad implications for suprathermal electron dynamics as well as whistler mode waves in the Earth's magnetosphere and has been demonstrated in the Jovian magnetosphere.

Li, Jinxing↗

Magnetic field amplification by collisionless shocks in partially ionized plasmas

In this paper, we study the shock structures of collisionless shocks in partially ionized plasmas by means of two-dimensional hybrid simulations, where the shock is a perpendicular shock with shock velocity v{sub sh}≈40 v{sub A}≈1333 km s{sup −1} and the upstream ionization fraction is 0.5. We find that large density fluctuations and large magnetic field fluctuations are generated both in the upstream and downstream regions. In addition, we find that the velocity distribution of downstream hydrogen atoms has three components. The observed shock structures suggest that diffusive shock acceleration can operate at perpendicular shocks propagating into partially ionized plasmas in real three-dimensional systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Development and Validation of a Simulation Testbed for the Intelligent Building Agents Laboratory (IBAL) Using TRNSYS

This paper documents the development and validation of a dynamic primary cooling and thermal storage system simulation testbed. The system simulation testbed, sIBAL, is based on the Intelligent Building Agents Laboratory (IBAL) at the National Institute of Standards and Technology (NIST), which is a research infrastructure and testbed for the development, evaluation, and demonstration of intelligent control algorithms. The sIBAL testbed developed in this project will serve as a virtual twin of the real facility for future control algorithm development. The details of the methodologies used to develop and validate the simulation testbed, which replicates the dynamic behaviors of the primary cooling and ice storage system in the IBAL facility, are presented. The simulation testbed was developed in TRNSYS using built-in component models and MATLAB functions to replicate the two water-cooled chillers, a thermal storage tank, pumps, valves and other components for four different operation modes. Experiments on IBAL components were designed and executed to generate experimental data for model development and verification of the simulation platform. The validation of the simulation results was carried out in two phases: 1) independent component simulations for the chillers and thermal storage tank, and 2) a combined testbed simulation of the entire hydronic system. Comparison of simulation results to the experimental data obtained from the IBAL facility showed errors within 1 °C for the temperatures outputs of both the chiller and the thermal storage model. The error is within an acceptable range for further intelligent control algorithms development. The findings from the study are summarized and presented along with areas where additional research is needed. In addition, data filtering procedures and model refinement measures utilized to improve the accuracy and accelerate the computation time of the simulation are presented.

cooling systems↗

HPSim, Version 1.0

HPSim is a GPU-accelerated online multi-particle beam dynamics simulation tool for ion linacs. It was originally developed for use on the Los Alamos 800-MeV proton linac. It is a “z-code” that contains typical linac beam transport elements. The linac RF-gap transformation utilizes transit-time-factors to calculate the beam acceleration therein. The space-charge effects are computed using the 2D SCHEFF (Space CHarge EFFect) algorithm, which calculates the radial and longitudinal space charge forces for cylindrically symmetric beam distributions. Other space- charge routines to be incorporated include the 3D PICNIC and a 3D Poisson solver. HPSim can simulate beam dynamics in drift tube linacs (DTLs) and coupled cavity linacs (CCLs). Elliptical superconducting cavity (SC) structures will also be incorporated into the code. The computational core of the code is written in C++ and accelerated using the NVIDIA CUDA technology. Users access the core code, which is wrapped in Python/C APIs, via Pythons scripts that enable ease-of-use and automation of the simulations. The overall linac description including the EPICS PV machine control parameters is kept in an SQLite database that also contains calibration and conversion factors required to transform the machine set points into model values used in the simulation.

Pang, Xiaoying↗

HPSim, Version 1.0

HPSim is a GPU-accelerated online multi-particle beam dynamics simulation tool for ion linacs. It was originally developed for use on the Los Alamos 800-MeV proton linac. It is a “z-code” that contains typical linac beam transport elements. The linac RF-gap transformation utilizes transit-time-factors to calculate the beam acceleration therein. The space-charge effects are computed using the 2D SCHEFF (Space CHarge EFFect) algorithm, which calculates the radial and longitudinal space charge forces for cylindrically symmetric beam distributions. Other space- charge routines to be incorporated include the 3D PICNIC and a 3D Poisson solver. HPSim can simulate beam dynamics in drift tube linacs (DTLs) and coupled cavity linacs (CCLs). Elliptical superconducting cavity (SC) structures will also be incorporated into the code. The computational core of the code is written in C++ and accelerated using the NVIDIA CUDA technology. Users access the core code, which is wrapped in Python/C APIs, via Pythons scripts that enable ease-of-use and automation of the simulations. The overall linac description including the EPICS PV machine control parameters is kept in an SQLite database that also contains calibration and conversion factors required to transform the machine set points into model values used in the simulation.

Pang, Xiaoying↗

A Simulation of the Fermilab Main Injector Cavities

The Fermilab Main Injector accelerating cavities have sparking issues when they are run at voltages higher than those required by the PIP-II project. This is a problem Fermilab is working on as planning begins for the next upgrade to the accelerator complex. One of the methods being used to address the issue is the development of a CST Microwave Studio simulation to accurately model the PIP-II dual power amplifier cavities and identify which part(s) of the cavity is causing sparking to develop. The model will also be used to determine if changes to the cavity geometry may allow the cavity to be used at higher voltages before sparking occurs

Stevenson, Susanna M. [Fermilab] (ORCID:0009000891↗

Computationally Guided Design of Multiple Impurities Tolerant Cathode for Solid Oxide Fuel Cell Applications

Solid oxide fuel cells (SOFCs) have been recognized as promising candidates in the demand of reducing the dependence on fossil fuels, which is one of the solutions to the threat of global warming. It was reported that the direct chemical-to-electrical energy conversion in SOFCs comes with negligible emission of greenhouse gases and up to 60% working efficiency. SOFC has unique advantages over traditional power generation technologies and other types of fuel cells due to the use of relatively inexpensive materials, low sensitivity to poisoning impurities in fuels, relatively high conversion efficiency, and fuel flexibility (like carbon-based fuels). However, during the SOFC operation, the cathode is exposed to intrinsic and extrinsic airborne contaminants at high temperatures (600–900°C), leading to their reaction with electrode materials, irreversible chemical and structural changes, and electrochemical performance degradation. of which CO2, SO2, Cr6+ and H2O are known to be concerns. Thus, developing a comprehensive understanding of the multiple impurities poisoning of SOFC cathodes in the presence of single or multiple impurities is necessary to address the long-term degradation of the SOFC cathode systems at operating conditions. An integrated computational and experimental approach has been utilized in the current work aiming at the comprehensive investigation of the cathode materials under single or multiple impurities poisoning effects, understanding the multiple impurities poisoning mechanism(s), and improving the long-term durability of the SOFC cathodes to develop better performing impurity-tolerant cathode materials. We first examined the agreement between the simulation approach and the experimental observations at various atmospheric and temperature conditions to confirm the reliability of CALPHAD for poisoning simulations in these cathodes. Then, using this approach, we investigated the accelerated testing protocol, which has become the standard method for conducting impurity poisoning experiments involving these cathode materials. We found our simulation approach is able to predict in which systems, i.e., cathode material and treatment environment, the accelerated testing protocols reflect actual operating conditions. Finally, in-depth simulations were done for these three candidate cathodes under different experimental conditions and suggested directions toward alternative cathode materials that have superior impurity tolerance.

36 MATERIALS SCIENCE↗

GPU acceleration and performance of the particle-beam-dynamics code Elegant

Elegant is an accelerator physics and particle-beam dynamics code widely used for modeling and design of a variety of high-energy particle accelerators and accelerator-based systems. We discuss a recently developed version of the code that can take advantage of CUDA-enabled graphics processing units (GPUs) to achieve significantly improved performance for a large class of simulations that are important in practice. The GPU version is largely defined by a framework that simplifies implementations of the fundamental kernel types that are used by Elegant: particle operations, reductions, particle loss, histograms, array convolutions and random number generation. Accelerated performance on the Titan Cray XK-7 supercomputer is approximately 6–10 times better with the GPU than all the CPU cores associated with the same node count. In addition to performance, the maintainability of the GPU-accelerated version of the code was considered a key design objective. Accuracy with respect to the CPU implementation is also a core consideration. Finally, four different methods are used to ensure that the accelerated code faithfully reproduces the CPU results.

43 PARTICLE ACCELERATORS↗

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗

Generation of gamma radiation by a subterawatt ultrashort laser pulse: optimisation of preplasma and pulse duration

We report an experimental and numerical study of the acceleration of electrons in a plasma interacting with a subterawatt laser pulse (intensity of ∼3 × 10{sup 18} W cm{sup −2} at a pulse duration of 50 fs). A preplasma layer on the surface of a molybdenum target is formed by an additional laser pulse with a duration of 8 ns and an intensity of ∼2 × 10{sup 12} W cm{sup −2}. It is shown that an increase in the laser pulse duration to 1700 fs at a constant energy (and a proportional decrease in intensity) leads to an increase in the yield of bremsstrahlung γ-radiation by more than an order of magnitude when the nanosecond pulse is ahead of the femtosecond one by 15 − 25 ns. Interferometry data and results of diagnostics of optical and γ-radiation of a plasma demonstrate that the collisional ionisation of atoms by electrons oscillating in the field of such a laser pulse plays an essential role in the formation of electron density profile. The sensitivity of the described effect to the level of amplified spontaneous emission is determined, despite the nanosecond pulse impact. Numerical simulations show that at a large pulse duration, the acceleration of electrons is stipulated by the beaking of plasma waves excited during stimulated Raman scattering of laser radiation. (paper)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

End-to-end simulation of bunch merging for a muon collider

Muon accelerator beams are commonly produced indirectly through pion decay by interaction of a charged particle beam with a target. Efficient muon capture requires the muons to be first phase-rotated by rf cavities into a train of 21 bunches with much reduced energy spread. Since luminosity is proportional to the square of the number of muons per bunch, it is crucial for a Muon Collider to use relatively few bunches with many muons per bunch. In this paper we will describe a bunch merging scheme that should achieve this goal. We present for the first time a complete end-to-end simulation of a 6D bunch merger for a Muon Collider. The 21 bunches arising from the phase-rotator, after some initial cooling, are merged in longitudinal phase space into seven bunches, which then go through seven paths with different lengths and reach the final collecting "funnel" at the same time. The final single bunch has a transverse and a longitudinal emittance that matches well with the subsequent 6D rectilinear cooling scheme.

43 PARTICLE ACCELERATORS↗

Enhanced proton acceleration in an applied longitudinal magnetic field

Using two-dimensional particle-in-cell simulations, we examine how an externally applied strong magnetic field impacts proton acceleration in laser-irradiated solid-density targets. We find that a kT-level external magnetic field can sufficiently inhibit transverse transport of hot electrons in a flat laser-irradiated target. While the electron heating by the laser remains mostly unaffected, the reduced electron transport during proton acceleration leads to an enhancement of maximum proton energies and the overall number of energetic protons. The resulting proton beam is much better collimated compared to a beam generated without applying a kT-level magnetic field. A factor of three enhancement of the laser energy conversion efficiency into multi-MeV protons is another effect of the magnetic field. The required kT-level magnetic fields are becoming feasible due to a significant progress that has been made in generating magnetic fields with laser-driven coils using ns-long laser pulses. The possibility of improving characteristics of laser-driven proton beams using such fields is a strong motivation for further development of laser-driven magnetic field capabilities.

43 PARTICLE ACCELERATORS↗

Measurement and Analysis of Beam Phase-Space Distributions for the Fermilab Muon Campus Accelerator Complex

The Muon g-2 experiment at Fermilab is tasked with measuring the muon’s anomalous magnetic moment with high precision. Since the experiment requires large amounts of muons, it is imperative to systematically study the behavior of the beam along the transport line. Unfortunately, the available diagnostics only provide beam information in X-Y space. For a complete evaluation, information of the phase-space is required. This paper demonstrates a technique to measure the beam phase-direction distribution by using a set of beam profiles. First, we establish the theoretical framework that describes the principle of the technique. Next, we apply the technique at four different locations along the accelerator delivery line. Finally, we compare our findings to predictions from tracking simulations.

43 PARTICLE ACCELERATORS↗

Plasma acceleration using a radio frequency self-bias effect

In this work plasma acceleration using a RF self-bias effect is experimentally studied. The experiments are conducted using a novel plasma accelerator system, called Neptune, consisting of an inductively coupled plasma source and a RF-biased set of grids. The plasma accelerator can operate in a steady state mode, producing a plasma flow with separately controlled plasma flux and velocity without any magnetic configuration. The operating pressure at the source output is as low as 0.2 mTorr and can further be decreased. The ion and electron flows are investigated by measuring the ion and electron energy distribution functions both space resolved and with different orientations with respect to the flow direction. It is found that the flow of electrons from the source is highly anisotropic and directed along the ion flow and this global flow of accelerated plasma is well localized in the plasma transport chamber. The maximum flux is about 7.5·10{sup 15} ions s{sup −1} m{sup −2} (at standard conditions) on the axis and decreasing to almost zero at a radial distances of more than 15 cm from the flow axis. Varying the RF acceleration voltage in the range 20–350 V, the plasma flow velocity can be changed between 10 and 35 km/s. The system is prospective for different technology such as space propulsion and surface modification and also interesting for fundamental studies for space-related plasma simulations and investigation of the dynamo effect using accelerated rotating plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗