Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Corkscrew Motion of an Electron Beam due to Coherent Variations in Accelerating Potentials

Corkscrew motion results from the interaction of fluctuations of beam electron energy with accidental magnetic dipoles caused by misalignment of the beam transport solenoids. Corkscrew is a serious concern for high-current linear induction accelerators (LIA). A simple scaling law for corkscrew amplitude derived from a theory based on a constant-energy beam coasting through a uniform magnetic field has often been used to assess LIA vulnerability to this effect. We use a beam dynamics code to verify that this scaling also holds for an accelerated beam in a non-uniform magnetic field, as in a real accelerator. Results of simulations with this code are strikingly similar to measurements on one of the LIAs at Los Alamos National Laboratory.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Particle Acceleration during Magnetic Reconnection in a Low-beta Plasma

Magnetic reconnection is a primary mechanism for particle energization in space and astrophysical plasmas. By carrying out two-dimensional (2D) fully kinetic simulations, we study particle acceleration during magnetic reconnection in plasmas with different plasma β (the ratio between the thermal pressure and the magnetic pressure). For the high-β cases, we do not observe significant particle acceleration. In the low-β regime (β<0.1), we find that reconnection is efficient at energizing both electrons and ions. While the distribution of accelerated particles integrated over the whole simulation box appears highly non-thermal, it is actually the superposition of a series of distributions in different sectors of a 2D magnetic island. Each of those distributions has only a small non-thermal component compared with its thermal core. By tracking a large number of particles, we show that particles get energized in X-line regions, contracting magnetic islands, and magnetic island coalescence regions. We obtain the particle energization rate j⋅E by averaging over particle drift motions and find that it agrees well with the particle kinetic energy change. We quantify the contribution of curvature drift, gradient drift, polarization drift, magnetization, non-gyrotropic effect, and parallel electric field in different acceleration regions. We find that the major energization is due to particle curvature drift along the motional electric field. The other particle motions contribute less but may become important in different acceleration regions. The highly efficient particle energization in low-β plasmas may help us understand the strong particle energization in solar flares and accretion disk coronae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

PONDEROMOTIVE ACCELERATION IN CORONAL LOOPS

Ponderomotive acceleration has been asserted to be a cause of the first ionization potential (FIP) effect, the well-known enhancement in abundance by a factor of 3–4 over photospheric values of elements in the solar corona with FIP less than about 10 eV. It is shown here by means of numerical simulations that ponderomotive acceleration occurs in solar coronal loops, with the appropriate magnitude and direction, as a “by-product” of coronal heating. The numerical simulations are performed with the HYPERION code, which solves the fully compressible three-dimensional magnetohydrodynamic equations including nonlinear thermal conduction and optically thin radiation. Numerical simulations of coronal loops with an axial magnetic field from 0.005 to 0.02 T and lengths from 25,000 to 75,000 km are presented. In the simulations the footpoints of the axial loop magnetic field are convected by random, large-scale motions. There is a continuous formation and dissipation of field-aligned current sheets, which act to heat the loop. As a consequence of coronal magnetic reconnection, small-scale, high-speed jets form. The familiar vortex quadrupoles form at reconnection sites. Between the magnetic footpoints and the corona the reconnection flow merges with the boundary flow. It is in this region that the ponderomotive acceleration occurs. Mirroring the character of the coronal reconnection, the ponderomotive acceleration is also found to be intermittent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Assessment of the Pan-Arctic Accelerated Rate of Sea Ice Decline in CMIP6 Historical Simulations

Abstract Arctic sea ice loss in response to a warming climate is assessed in 42 models participating in phase 6 of the Coupled Model Intercomparison Project (CMIP6). Sea ice observations show a significant acceleration in the rate of decline commencing near the turn of the twenty-first century. It is our assertion that state-of-the-art climate models should qualitatively reflect this accelerated trend within the limitations of internal variability and observational uncertainty. Our analysis shows that individual CMIP6 simulations of sea ice depict a wide range of model spread on biases and anomaly trends both across models and among their ensemble members. While the CMIP6 multimodel mean captures the observed sea ice area (SIA) decline relatively well, an individual model’s ability to represent the acceleration in sea ice decline remains a challenge. Seventeen (40%) out of 42 CMIP6 models and 37 (13%) out of the total 286 ensemble members reasonably capture the observed trends and acceleration in SIA decline. In addition, a larger ensemble size appears to increase the odds for a model to include at least one ensemble member skillfully representing the accelerated SIA trends. Simulations of sea ice volume (SIV) show much larger spread and uncertainty than SIA; however, due to limited availability of sea ice thickness data, these are not as well constrained by observations. Finally, we find that models with more ocean heat transport simulate larger sea ice declines, which suggests an emergent constraint in CMIP6 ensembles. This relationship points to the need for better understanding and modeling of ice–ocean interactions, especially with respect to frazil ice growth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

AlphaBuilding ResCommunity: A multi-agent virtual testbed for community-level load coordination

Training and validating algorithms in a simulation testbed can accelerate research and applications of optimal control of residential loads to improve energy flexibility and grid resilience. We developed an open-source simulation environment, AlphaBuilding ResCommunity, that can be used to train and validate algorithms to control a single thermostatically controlled load (TCL) or coordinate a group of TCLs. We used reduced-order models to simulate the thermodynamics of TCLs, and the parameter values were determined from the connected smart thermostat data of real households. The environment was built upon the standardized OpenAI Gym interface. Ancillary functions, such as retrieving the parameters and weather forecasts, are provided to facilitate control strategies that require predictive information. Compared with existing efforts, AlphaBuilding ResCommunity has three advantages: (1) more realistic model settings because the parameter values are identified from actual household operating data, and modelling and measurement uncertainty are considered; (2) passive thermal storage control; and (3) ease of use due to a simple software dependency and standardized interface. We demonstrated the applications of the environment by implementing a Kalman Filter and Model Predictive Control on a single TCL and a Priority-Stack-Based Control and Alternating Direction Method of Multipliers to coordinate multiple TCLs for load tracking.

Wang, Z↗
DOE OSTI2020

Confining Liquids inside Carbon Nanotubes: Accelerated Molecular Dynamics with Spliced, Soft-Core Potentials and Simulated Annealing

Understanding emergent phenomena of fluids under physical confinement requires the development of advanced tools for rapid and accurate simulation of their physiochemical properties. Simulating liquid molecules commensurate in size with the nanoscale enclosures that confine them is a key challenge. In this work, we demonstrate an accelerated molecular dynamics simulation technique that combines soft-core potentials (SCP) and simulated annealing (SA) to analyze confined liquids. This integrated SCP/SA method relies on a new spliced soft-core potential (SSCP), which enables tunable accuracy with respect to the target hard-core potential (HCP). SCP/SA enables the packing of enclosures with bulk material in a controlled, thermodynamically consistent manner. The enhanced SSCP accuracy is a critical feature of SCP/SA, enabling a smooth transition between the SCP and the HCP at a desired SCP hardness. We applied SCP/SA to the problem of filling a carbon nanotube (CNT) in periodic boundary conditions with a popular ionic liquid (IL), 1- butyl-3-methylimidazolium hexafluorophosphate $[\text{BMIM}^+][\text{PF}^–_6]$. We performed a series of triplicate simulations on systems with varying CNT diameter and charge to demonstrate SCP/SA’s versatility. Beyond this IL/CNT system, the SCP/SA simulation framework has a broad range of potential applications, not limited to nanoscale enclosures and interfaces, including both solid-state and biological systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deep quantum circuit simulations of low-energy nuclear states

Numerical simulation is an important method for verifying the quantum circuits used to simulate low-energy nuclear states. However, real-world applications of quantum computing for nuclear theory often generate deep quantum circuits that place demanding memory and processing requirements on conventional simulation methods. Here, we present advances in high-performance numerical simulations of deep quantum circuits to efficiently verify the accuracy of low-energy nuclear physics applications. Our approach employs novel methods for accelerating the numerical simulation including management of simulated mid-circuit measurements to verify projection based state preparation circuits. In this study, we test these methods across a variety of high-performance computing systems and our results show that circuits up to 21 qubits and more than 115,000,000 gates can be efficiently simulated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Simulation of sugar kelp ( Saccharina latissima ) breeding guided by practices to accelerate genetic gains

Abstract Though Saccharina japonica cultivation has been established for many decades in East Asian countries, the domestication process of sugar kelp (Saccharina latissima) in the Northeast United States is still at its infancy. In this study, by using data from our breeding experience, we will demonstrate how obstacles for accelerated genetic gain can be assessed using simulation approaches that inform resource allocation decisions. Thus far, we have used 140 wild sporophytes that were sampled in 2018 from the northern Gulf of Maine to southern New England. From these sporophytes, we sampled gametophytes and made and evaluated over 600 progeny sporophytes from crosses among the gametophytes in 2019 and 2020. The biphasic life cycle of kelp gives a great advantage in selective breeding as we can potentially select both on the sporophytes and gametophytes. However, several obstacles exist, such as the amount of time it takes to complete a breeding cycle, the number of gametophytes that can be maintained in the laboratory, and whether positive selection can be conducted on farm-tested sporophytes. Using the Gulf of Maine population characteristics for heritability and effective population size, we simulated a founder population of 1,000 individuals and evaluated the impact of overcoming these obstacles on rate of genetic gain. Our results showed that key factors to improve current genetic gain rely mainly on our ability to induce reproduction of the best farm-tested sporophytes, and to accelerate the clonal vegetative growth of released gametophytes so that enough gametophyte biomass is ready for making crosses by the next growing season. Overcoming these challenges could improve rates of genetic gain more than 2-fold. Future research should focus on conditions favorable for inducing spring reproduction, and on increasing the amount of gametophyte tissue available in time to make fall crosses in the same year.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Accelerating Convergence of Free Energy Computations with Hamiltonian Simulated Annealing of Solvent (HSAS)

Coupling between binding of a ligand to a receptor and the displacement of a number of bound water molecules is a common event in molecular recognition processes. When the binding site is deeply buried and the exchange of water molecules with the bulk region is difficult to sample, the convergence and accuracy in free energy calculations can be severely compromised. Traditionally, Grand Canonical Monte Carlo (GCMC) based methods have been used to accelerate equilibration of water-at the expense, however, of lengthy trials before a molecular dynamics (MD) simulation. In this paper, a user-friendly and cost-efficient method, Hamiltonian simulated annealing of solvent in combination with lambda-exchange of free energy perturbation (FEP) is proposed to accelerate the sampling of water molecules in free energy calculations. As an illustrative example with reliable data from previous GCMC simulations, absolute binding affinity of camphor to cytochrome P450 was calculated. The simulated hydration state change in the buried binding pocket quantitatively agrees with GCMC simulations. Finally, it is shown that the new protocol significantly accelerates sampling of water in a buried binding pocket and the convergence of free energy, with negligible setup and computing costs compared to GCMC methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

The acceleration of electrons at collisionless shocks moving through a turbulent magnetic field

We perform a numerical-simulation study of the acceleration of electrons at shocks that propagate through a prespecified, kinematically defined turbulent magnetic field. The turbulence consists of broadband magnetic fluctuations that are embedded in the plasma and cover a range of wavelengths, the smallest of which is larger than the gyroradii of electrons that are initially injected into the system. We find that when the variance of the turbulent component of the upstream magnetic field is sufficiently large—σ{sup 2} ∼ 10 B{sub 0}{sup 2}, where B{sub 0} is the strength of the background magnetic field—electrons can be efficiently accelerated at a collisionless shock regardless of the orientation of the mean upstream magnetic field relative to the shock-normal direction. Since the local angle between the incident magnetic-field vector and the shock-normal vector can be quite large, electrons can be accelerated through shock-drift acceleration at the shock front. In the upstream region, electrons are mirrored back to the shock front leading to multiple shock encounters. Eventually the accelerated electrons are energetic enough that their gyroradii are of the same order as the wavelength of waves that are included in our description of the turbulent magnetic field. Our results are consistent with recent in situ observations at Saturn’s bow shock. This study may help us to understand the acceleration of electrons at shocks in space and astrophysical systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Study of magnetic reconnection at low- β using laser-powered capacitor coils

Magnetic reconnection is a ubiquitous fundamental process in space and astrophysical plasmas that rapidly converts magnetic energy into some combination of flow energy, thermal energy, and non-thermal energetic particles. Over the past decade, a new experimental platform has been developed to study magnetic reconnection using strong coil currents powered by high-power lasers at low plasma beta, typical conditions under which reconnection is energetically important in space and astrophysics. KJ-class lasers were used to drive parallel currents to reconnect MG-level magnetic fields in a quasi-axisymmetric geometry, similar to the magnetic reconnection experiment or MRX, and thus this platform is named micro-MRX. This presentation summarizes two major findings from micro-MRX: direct measurement of accelerated electrons and observation of ion acoustic waves during anti-parallel reconnection. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that direct acceleration by the out-of-plane reconnection electric field is at work. Furthermore, a sudden onset of ion acoustic bursts has been measured by collective Thomson scattering in the exhaust of magnetic reconnection, followed by electron acoustic bursts with electron heating and bulk acceleration. These results demonstrate that the micro-MRX platform offers a novel and unique approach to study magnetic reconnection in the laboratory in addition to the capabilities provided by traditional magnetized plasma experiments such as MRX and the upcoming Facility for Laboratory Reconnection experiments (FLARE). Future prospects to study other particle acceleration mechanisms and ion acoustic waves from magnetic reconnection are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Electromagnetic Transient (EMT) Simulation Algorithm for Evaluation of Photovoltaic (PV) Generation Systems

Use of inverter-based resources facilitating renewable energy resources such as photovoltaic (PV) generation is increasing rapidly with decreasing costs and reduced emissions associated. To accommodate such rapid growth of inverter-based resources like PV systems, electromagnetic transient (EMT) simulation models of both PV systems and grids are required to analyze the interaction of PVs in the grid (like the post-event analysis). In addition, the EMT simulation would help with the planning of future power grid with a large number of PVs as well as other inverter-based distributed generation systems. In this paper, the EMT simulation models of PV systems and grids are developed based on the differential algebraic equations (DAEs) representing their EMT dynamics. Furthermore, advanced simulation algorithms including numerical stiffness-based hybrid discretization, DAE clustering and aggregation, multi-order integration, and matrix splitting approaches are applied to accelerate the EMT simulation. The proposed algorithm was applied to 125 PV inverters within 52-bus medium-voltage (MV) distribution grid.

Choi, Jongchan↗
DOE OSTI2022

Beam Breakup Simulations for the Scorpius Flash-Radiography Accelerator

Beam breakup (BBU) in linear induction accelerators (LIAs) used for flash radiography is problematic because the high-frequency beam motion can blur the source spot, thereby degrading resolution. Amplification of BBU depends directly on details of cell design and is suppressed in operation by the solenoidal magnetic field focusing on the electron beam. Therefore, much effort has gone into design of the cells and magnetic focusing for the new Scorpius radiography LIA. Here, in this article, we use computer simulations to demonstrate that BBU in Scorpius should be no more than in present radiography LIAs at the Los Alamos National Laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A generative artificial intelligence framework for long-time plasma turbulence simulations

Generative deep learning techniques are employed in a novel framework for the construction of surrogate models capturing the spatiotemporal dynamics of 2D plasma turbulence. The proposed Generative Artificial Intelligence Turbulence (GAIT) framework enables the acceleration of turbulence simulations for long-time transport studies. GAIT leverages a convolutional variational auto-encoder and a recurrent neural network to generate new turbulence data from existing simulations, extending the time horizon of transport studies with minimal computational cost. The application of the GAIT framework to plasma turbulence using the Hasegawa–Wakatani (HW) model is presented, evaluating its performance via various analyses. Very good agreement is found between the GAIT and the HW models in the spatiotemporal Fourier and Proper Orthogonal Decomposition spectra, the flow topology characterized by the Okubo–Weiss parameter, and the time autocorrelation function of turbulent fluctuations. Excellent agreement has also been obtained in the probability distribution function of particle displacements and the effective turbulent diffusivity. In-depth analyses of the latent space of turbulent states, choice of hyperparameters and alternative deep learning models for the time prediction are presented. Our results highlight the potential of Artificial Intelligence-based surrogate models to overcome the computational challenges in turbulence simulation, which can be extended to other situations such as geophysical fluid dynamics.

Artificial intelligence↗
DOE OSTI2015

WE-G-BRD-09: Novel MRI Compatible Electron Accelerator for MRI-Linac Radiotherapy

Purpose: MRI guided radiotherapy is a rapidly growing field; however current linacs are not designed to operate in MRI fringe fields. As such, current MRI- Linac systems require magnetic shielding, impairing MR image quality and system flexibility. Here, we present a bespoke electron accelerator concept with robust operation in in-line magnetic fields. Methods: For in-line MRI-Linac systems, electron gun performance is the major constraint on accelerator performance. To overcome this, we propose placing a cathode directly within the first accelerating cavity. Such a configuration is used extensively in high energy particle physics, but not previously for radiotherapy. Benchmarked computational modelling (CST, Darmstadt, Germany) was employed to design and assess a 5.5 cell side coupled accelerator with a temperature limited thermionic cathode in the first accelerating cell. This simulation was coupled to magnetic fields from a 1T MRI model to assess robustness in magnetic fields for Source to Isocenter Distance between 1 and 2 meters. Performance was compared to a conventional electron gun based system in the same magnetic field. Results: A temperature limited cathode (work function 1.8eV, temperature 1245K, emission constant 60A/K/cm{sup 2}) will emit a mean current density of 24mA/mm{sup 2} (Richardson’s Law). We modeled a circular cathode with radius 2mm and mean current 300mA. Capture efficiency of the device was 43%, resulting in target current of 130 mA. The electron beam had a FWHM of 0.2mm, and mean energy of 5.9MeV (interquartile spread of 0.1MeV). Such an electron beam is suitable for radiotherapy, comparing favourably to conventional systems. This model was robust to operation the MRI fringe field, with a maximum current loss of 6% compared to 85% for the conventional system. Conclusion: The bespoke electron accelerator is robust to operation in in-line magnetic fields. This will enable MRI-Linacs with no accelerator magnetic shielding, and minimise painstaking optimisation of the MRI fringe field. This work was supported by US (NIH) and Australian (NHMRC & Cancer Institute NSW) government research funding. In addition, I would like to thank cancer institute NSW and the Ingham Institute for scholarship support.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A gas-dynamical approach to radiation pressure acceleration

The study of high intensity ion beams driven by high power pulsed lasers is an active field of research. Of particular interest is the radiation pressure acceleration, for which simulations predict narrow band ion energies up to GeV. We derive a laser-piston model by applying techniques for non-relativistic gas-dynamics. The model reveals a laser intensity limit, below which sufficient laser-piston acceleration is impossible. The relation between target thickness and piston velocity as a function of the laser pulse length yields an approximation for the permissible target thickness. We performed one-dimensional Particle-In-Cell simulations to confirm the predictions of the analytical model. These simulations also reveal the importance of electromagnetic energy transport. We find that this energy transport limits the achievable compression and rarefies the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The design of the electron beam dump unit of Turkish Accelerator Center (TAC)

The required simulations of the electron beam interactions for the design of electron beam dump unit for an accelerator which will operate to get two Infra-Red Free Electron Lasers (IR-FEL) covering the range of 3-250 microns is presented in this work. Simulations have been carried out to understand the interactions of a bulk of specially shaped of four different and widely used materials for the dump materials for a 77 pC, 40 MeV, 13 MHz repetition rate e-beam. In the simulation studies dump materials are chosen to absorb the 99% of the beam energy and to restrict the radio-isotope production in the bulk of the dump. A Lead shielding also designed around the dump core to prevent the leakage out of the all the emitted secondary radiations, e.g., neutrons, photons. The necessary dump material requirements, for the overall design considerations and the possible radiation originated effects on the dump unit, are discussed and presented.

07 ISOTOPE AND RADIATION SOURCES↗