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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

On the possibility of footprint compression with one lens in nonlinear accelerator lattice

Electromagnetic interaction of colliding beams along with other nonlinear fields often limits the beams' lifetimes and luminosities. Nonlinearities result in the spread of betatron frequencies (footprint) and, thus, may enhance dynamic diffusion of particles due to high order resonances. One of the possible ways to eliminate nonlinearities and overcome the corresponding difficulties is compensation of nonlinear forces, but, in practice, it is hardly possible to obtain exact linearity of the system. The compensation with a single nonlinear lens cannot cope with distributed nonlinearities, nonlinearities due to parasitic crossings, etc. Here, we present a method to compute parameters of nonlinear element (lens) that eliminates both the footprint and resonance strength without achieving full compensation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗
NASA NTRS2001

Accelerator Test of an Imaging Calorimeter

The Imaging Calorimeter for ACCESS (ICA) utilizes a thin sampling calorimeter concept for direct measurements of high-energy cosmic rays. The ICA design uses arrays of small scintillating fibers to measure the energy and trajectory of the produced cascades. A test instrument has been developed to study the performance of this concept at accelerator energies and for comparison with simulations. Two test exposures have been completed using a CERN test beam. Some results from the accelerator tests are presented.

Christl, Mark J.↗
DOE OSTI2012

Multiphysics Applications of ACE3P

The TEM3P module of ACE3P, a parallel finite-element electromagnetic code suite from SLAC, focuses on the multiphysics simulation capabilities, including thermal and mechanical analysis for accelerator applications. In this pa- per, thermal analysis of coupler feedthroughs to supercon- ducting rf (SRF) cavities will be presented. For the realistic simulation, internal boundary condition is implemented to capture RF heating effects on the surface shared by a di- electric and a conductor. The multiphysics simulation with TEM3P matched the measurement within 0.4%.

43 PARTICLE ACCELERATORS↗
NASA NTRS1977

Study of changes in properties of solar sail materials from radiation exposure

Techniques for monitoring changes in preparation of solar sail materials resulting from space radiation simulation, stressing (e.g., thermal, mechanical) and exposure to terrestrial environments are developed. The properties of interest are: metallic coating deterioration, polymeric film deterioration, interfacial debonding and possible metallic coating diffusion into the polymeric film. Four accelerated tests were devised to simulate the possible degradation processes mentioned above. These four tests are: a thermal shock test to simulate the wide variation of temperature expected in space (260 C to -100 C), a cyclic temperature test to stimulate the 6 minute temperature cycle anticipated in space, a mechanical vibration test to simulate mechanical bonding, folding and handling, and a humidity test to simulate terrestrial environment effects. The techniques for monitoring property changes are: visual and microscopic examination, ellipsometry, surface potential difference (SPD), photoelectron emission (PEE), and water contact angles.

Smith, T.↗
DOE OSTI2018

Pronghorn: A Porous Media Thermal-Hydraulics Core Simulator and its Validation with the SANA Experiments

Pebble bed High Temperature Reactors (HTRs) are characterized by many advantageous design features, such as excellent passive heat removal in accidents, large margins to fuel failure, and online refueling potential. However, a significant challenge in the core modeling of pebble bed reactors is the complex fuel-coolant structure. This paper presents a new porous media simulation code, Pronghorn, that aims to alleviate modeling challenges for pebble bed reactors by providing a fast-running, medium- fidelity core simulator. Pronghorn is intended to accelerate the design and analysis cycle for pebble bed and prismatic HTRs by permitting fast scoping studies and providing boundary conditions for systems-level analysis. Pronghorn is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) using modern software practices and a thorough testing framework. This paper describes the physical models used in Pronghorn and demonstrates Pronghorn’s capability for modeling gas-cooled pebble bed HTRs by presenting simulation results obtained for the German SANA pebble bed decay heat experiments. Within the limitations of the porous media approximation and existing available closure relationships, Pronghorn predicts the SANA experimental pebble temperatures well, expanding the code’s validation base. A brief code-to-code comparison shows a similar level of accuracy with other porous media simulation tools; this enhances Pronghorn’s advantages of an arbitrary equation of state, unstructured mesh capabilities, compressible flow models, high coupling potential to MOOSE fuels performance and systems-level thermal-hydraulics codes, and modern software design.

42 ENGINEERING↗
DOE OSTI2015

GEANT4 simulations for beam emittance in a linear collider based on plasma wakefield acceleration

Alternative acceleration technologies are currently under development for cost-effective, robust, compact, and efficient solutions. One such technology is plasma wakefield acceleration, driven by either a charged particle or laser beam. However, the potential issues must be studied in detail. In this paper, the emittance evolution of a witness beam through elastic scattering from gaseous media and under transverse focusing wakefields is studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Commissioning of a Varian Clinac iX 6 MV photon beam using Monte Carlo simulation

Monte Carlo modelling of a linear accelerator is the first and most important step in Monte Carlo dose calculations in radiotherapy. Monte Carlo is considered today to be the most accurate and detailed calculation method in different fields of medical physics. In this research, we developed a photon beam model for Varian Clinac iX 6 MV equipped with MilleniumMLC120 for dose calculation purposes using BEAMnrc/DOSXYZnrc Monte Carlo system based on the underlying EGSnrc particle transport code. Monte Carlo simulation for this commissioning head LINAC divided in two stages are design head Linac model using BEAMnrc, characterize this model using BEAMDP and analyze the difference between simulation and measurement data using DOSXYZnrc. In the first step, to reduce simulation time, a virtual treatment head LINAC was built in two parts (patient-dependent component and patient-independent component). The incident electron energy varied 6.1 MeV, 6.2 MeV and 6.3 MeV, 6.4 MeV, and 6.6 MeV and the FWHM (full width at half maximum) of source is 1 mm. Phase-space file from the virtual model characterized using BEAMDP. The results of MC calculations using DOSXYZnrc in water phantom are percent depth doses (PDDs) and beam profiles at depths 10 cm were compared with measurements. This process has been completed if the dose difference of measured and calculated relative depth-dose data along the central-axis and dose profile at depths 10 cm is ≤ 5%. The effect of beam width on percentage depth doses and beam profiles was studied. Results of the virtual model were in close agreement with measurements in incident energy electron 6.4 MeV. Our results showed that photon beam width could be tuned using large field beam profile at the depth of maximum dose. The Monte Carlo model developed in this study accurately represents the Varian Clinac iX with millennium MLC 120 leaf and can be used for reliable patient dose calculations. In this commissioning process, the good criteria of dose difference in PDD and dose profiles were achieve using incident electron energy 6.4 MeV.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2019

Numerical investigations of background pressure effects and channel erosion in the SPT-140 Hall thruster for the Psyche mission

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

Lenguito, Giovanni↗
DOE OSTI2024

Emulator of PR‐DNS: Accelerating Dynamical Fields With Neural Operators in Particle‐Resolved Direct Numerical Simulation

Abstract Particle‐resolved direct numerical simulations (PR‐DNS) play an increasing role in investigating aerosol‐cloud‐turbulence interactions at the most fundamental level of processes. However, the high computational cost associated with high resolution simulations poses considerable challenges for large domain or long duration simulation using PR‐DNS. To address these issues, here we present an emulator of the complex physics‐based PR‐DNS developed by use of the data‐driven Fourier Neural Operator (FNO) method. The effectiveness of the method is showcased by presenting turbulence and temperature fields in a two‐dimensional space. The results demonstrate high accuracy at various resolutions and the emulator is two orders of magnitude cheaper in terms of computational demand compared to the physics‐based PR‐DNS model. Furthermore, the FNO emulator exhibits strong generalization capabilities for different initial conditions and ultra‐high‐resolution without the need to retrain models. These findings highlight the potential of the FNO method as a promising tool to simulate complex fluid dynamics problems with high accuracy, computational efficiency, and generalization capabilities, enhancing our understanding of the aerosol‐cloud‐precipitation system.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2017

Irradiation Accelerated Corrosion of Stainless Steel and Ferritic-Martensitic Steel in Simulated Primary Water

Austenitic stainless steel, and ferritic-martensitic steel were irradiated with a proton beam while simultaneously exposed to simulated PWR primary water. Samples were exposed for 4-72 hr in 320°C water with 3 wppm hydrogen while simultaneously irradiated at surface dose rates from 400-4000 kGy/s (4x10-7 to 7x10-6 dpa/s). After exposure, samples were characterized with TEM, EDS, and Raman spectroscopy to determine the effects of radiation on the growth and stability of the oxide films. Samples which were irradiated during exposure had thinner inner oxides compared to unirradiated samples. Irradiated samples also had hematite on their surfaces, while unirradiated samples did not, indicating an increase in corrosion potential due to irradiation. Sample regions which were not irradiated but were exposed to the flow of irradiated water were similar to directly irradiated regions, suggesting that water radiolysis is the mechanism primarily responsible for the observed effects.

Stainless Steel, ferritic-martensitic steel, radia↗
DOE OSTI2022

Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration

The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the δƒ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K with the δƒ method, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes have been achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Parameter Sweep and Optimization of Loosely Coupled Simulations Using the DAKOTA Toolkit

The increasing availability of large scale computing capabilities has accelerated the development of high-fidelity coupled simulations. Such simulations typically involve the integration of models that implement various aspects of the complex phenomena under investigation. Coupled simulations are playing an integral role in fields such as climate modeling, earth systems modeling, rocket simulations, computational chemistry, fusion research, and many other computational fields. Model coupling provides scientists with systematic ways to virtually explore the physical, mathematical, and computational aspects of the problem. Such exploration is rarely done using a single execution of a simulation, but rather by aggregating the results from many simulation runs that, together, serve to bring to light novel knowledge about the system under investigation. Furthermore, it is often the case (particularly in engineering disciplines) that the study of the underlying system takes the form of an optimization regime, where the control parameter space is explored to optimize an objective functions that captures system realizability, cost, performance, or a combination thereof. Novel and flexible frameworks that facilitate the integration of the disparate models into a holistic simulation are used to perform this research, while making efficient use of the available computational resources. In this paper, we describe the integration of the DAKOTA optimization and parameter sweep toolkit with the Integrated Plasma Simulator (IPS), a component-based framework for loosely coupled simulations. The integration allows DAKOTA to exploit the internal task and resource management of the IPS to dynamically instantiate simulation instances within a single IPS instance, allowing for greater control over the trade-off between efficiency of resource utilization and time to completion. We present a case study showing the use of the combined DAKOTA-IPS system to aid in the design of a lithium ion battery (LIB) cell, by studying a coupled system involving the electrochemistry and ion transport at the lower length scales and thermal energy transport at the device scales. The DAKOTA-IPS system provides a flexible tool for use in optimization and parameter sweep studies involving loosely coupled simulations that is suitable for use in situations where changes to the constituent components in the coupled simulation are impractical due to intellectual property or code heritage issues.

IPS CAEBAT Coupled Battery Simulations↗
DOE OSTI2025

Application of a temporal multiscale method for efficient simulation of degradation in PEM Water Electrolysis under dynamic operating conditions

Hydrogen is emerging as a vital energy carrier, driven by the need to reduce carbon emissions. Proton Electrolyte Membrane Water Electrolysis (PEMWE) enables hydrogen production under fluctuating renewable power conditions but requires improved understanding and stability of the anode catalyst layer under dynamic operating conditions, especially with low noble metal loadings. Long-term degradation experiments are both time-consuming and costly; therefore, a systematic, model-aided approach is essential. In the present work, a temporal multiscale method is applied to reduce the computational effort of simulating long-term degradation processes in PEMWE, with an exemplary focus on catalyst dissolution. A mechanistic model incorporating the oxygen evolution reaction, catalyst dissolution, and hydrogen permeation from the cathode to the anode was hypothesized and implemented. In this way, the local periodicity of transport and reaction processes in dynamic PEMWE operation, which influence the gradual degradation of the catalyst layer, is captured. The temporal multiscale method significantly reduces the computational effort of simulation, decreasing processing time from hours to mere minutes. This efficiency gain is attributed to the limited evolution of Slow-Scale variables during each period of time P of the Fast-Scale variables. Consequently, simulation is required only until local periodicity is achieved within each Slow-Scale time step. Hence, the fully resolved dynamic problem is decoupled into these two scales, employing a heterogeneous multiscale technique. The developed approach effectively accelerates parameter estimation and predictive simulations, supporting systematic modeling of PEMWE degradation under dynamic conditions.

08 HYDROGEN↗
DOE OSTI2026

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗
DOE OSTI2017

Preferential Heating and Acceleration of Heavy Ions in Impulsive Solar Flares

We simulate decaying turbulence in a homogeneous pair plasma using a three-dimensional electromagnetic particle-in-cell method. A uniform background magnetic field permeates the plasma such that the magnetic pressure is three times larger than the thermal pressure and the turbulence is generated by counter-propagating shear Alfvén waves. The energy predominately cascades transverse to the background magnetic field, rendering the turbulence anisotropic at smaller scales. We simultaneously move several ion species of varying charge to mass ratios in our simulation and show that the particles of smaller charge to mass ratios are heated and accelerated to non-thermal energies at a faster rate. This is in accordance with the enhancement of heavy ions and a non-thermal tail in their energy spectrum observed in the impulsive solar flares. We further show that the heavy ions are energized mostly in the direction perpendicular to the background magnetic field, with a rate consistent with our analytical estimate of the rate of heating due to cyclotron resonance with the Alfvén waves, of which a large fraction is due to obliquely propagating waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Advances in Multi-Physics Modeling for Parallel Finite Element Code Suite ACE3P

SLAC developed ACE3P code suite is a comprehensive set of parallel finite-element codes for multi-physics modeling of particle accelerators. Running on massively parallel computer platforms for high fidelity and high accuracy simulation, ACE3P enables rapid virtual prototyping of accelerator and RF component design, optimization and analysis. Furthermore, recent advances of ACE3P have been focused on the implementation of advanced numerical algorithms, enhancement of multi-physics capabilities, and improvement of code performance on state-of-the-art high performance computing (HPC) platforms for large-scale RF modeling in accelerator applications and beyond. In this paper, we will present a nonlinear electromagnetic (EM) eigensolver for damping calculation of resonant modes for cavities coupled with waveguides, the thermal solver for accurate evaluation of heat loads in superconducting RF (SRF) cavity power coupler, the elastic solver for investigation of SRF cavity cryomodule (CM) bowing due to thermal effect, the frequency domain mechanical solver for studying SRF cavity mechanical oscillation modes, the particle tracking code for dark current simulation in SRF CM using a hybrid MPI+OpenMP parallel programing model, the time domain solver for high fidelity modeling of EM wave propagation for medical device in human body.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Three-dimensional simulations of the non-thermal broadband emission from young supernova remnants including efficient particle acceleration

Supernova remnants are believed to be major contributors to Galactic cosmic rays. In this paper, we explore how the non-thermal emission from young remnants can be used to probe the production of energetic particles at the shock (both protons and electrons). Our model couples hydrodynamic simulations of a supernova remnant with a kinetic treatment of particle acceleration. We include two important back-reaction loops upstream of the shock: energetic particles can (1) modify the flow structure and (2) amplify the magnetic field. As the latter process is not fully understood, we use different limit cases that encompass a wide range of possibilities. We follow the history of the shock dynamics and of the particle transport downstream of the shock, which allows us to compute the non-thermal emission from the remnant at any given age. We do this in three dimensions, in order to generate projected maps that can be compared with observations. We observe that completely different recipes for the magnetic field can lead to similar modifications of the shock structure, although to very different configurations of the field and particles. We show how this affects the emission patterns in different energy bands, from radio to X-rays and γ-rays. High magnetic fields (>100 μG) directly impact the synchrotron emission from electrons, by restricting their emission to thin rims, and indirectly impact the inverse Compton emission from electrons and also the pion decay emission from protons, mostly by shifting their cut-off energies to respectively lower and higher energies.

79 ASTRONOMY AND ASTROPHYSICS↗