Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Beam Simulations”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Cathode to Target Simulations for Scorpius: I. Simulation Codes and Models

The new Scorpius linear induction electron accelerator is under development for multi-pulse flash radiography of large, explosively-driven hydrodynamic experiments. Beam physics from the cathode to the target was examined with computer simulations, including envelope, ray-trace, and particle-in cell (PIC) codes. Beam instabilities investigated included beam breakup (BBU), image displacement, diocotron, parametric envelope, ion hose, and the resistive wall instability. Beam corkscrew motion and emittance growth from beam mismatch were also studied. The results of these simulations is documented in a series of reports. In this report the computer codes and physical models used for these simulations are described. The conclusion of this study is that Scorpius will produce and accelerate a beam with radiographic quality equivalent to the present accelerators at Los Alamos National Laboratory if the same engineering standards and construction details are upheld.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Reconstruction of beam parameters and betatron radiation spectra measured with a Compton spectrometer

The photon flux resulting from high-energy electron beam interactions with high-field systems, such as those found in the upcoming FACET-II experiments at the SLAC National Accelerator Laboratory, yields deep insight into the electron beam’s underlying dynamics during the interaction. However, extracting this information is an intricate process. To demonstrate how to approach this challenge using modern methods, this paper utilizes simulated data that models plasma wakefield acceleration-derived betatron radiation in experiments to determine reliable methods of reconstructing key beam and beam-plasma interaction properties. For betatron radiation measurements, translating the observed 200⁢ keV to 30⁢ MeV photon double-differential energy-angle spectra obtained from an advanced Compton spectrometer requires testing multiple methods to optimize the pipeline from its response to incident electron beam information. The paper compares maximum likelihood estimation and machine learning to refine the translation of photon spectra into precise electron beam metrics, such as spot size, energy, and emittance, enhancing the understanding of beam behavior within these dense, high-field environments. We also introduce machine learning and the expected maximization algorithm to reconstruct the primary photon spectrum, employing a multilayer neural network for regression analysis of the energy and angle spectra. With appropriate modifications, the advanced methods reproduce relevant incident beam parameters with high accuracy, even for beam sizes in the <10 μ⁢m range. This capacity is critical to understanding intense beam propagation and its optimization in plasma.

Beam code development & simulation techniques↗
DOE OSTI2026

Experimental Studies of Nonlinear Integrable Optics with Elliptic Potentials

The stable transport of charged particle beams is the core challenge in designing and operating accelerators. The current paradigm for transverse focusing is based on quadrupole and dipole elements, described as a linear integrable Hamiltonian by Courant and Snyder. The operational limits of high intensity accelerators are often determined by collective instabilities within the beam. Nonlinear elements may be added to suppress certain forms of these collective effects, but restrict the range of stable trajectories. These shortcomings motivate extending to a nonlinear integrable system which can offer the benefits of suppressing collective instabilities without limiting the stable trajectories. Danilov and Nagaitsev proposed a novel nonlinear integrable system with elliptical potentials, which could be implemented as magnetic elements in an accelerator. Such a system has been implemented for practical verification in the integrable optics test accelerator (IOTA), a small storage ring constructed for beam dynamics studies. Electron beam studies in IOTA treat the low-emittance beam as a macroparticle for detailed probing of the expected single particle dynamics. This dissertation details new measurements of the nonlinear integrable system relevant for practical implementation. Turn-by-turn measurements of the kicked beam responses are used for phase space reconstruction and analysis of the predicted nonlinear dynamics. The stability and aperture for various nonlinear configurations are measured with beam losses. Amplitude dependent detuning, a core figure of merit for suppressing instabilities, is measured and compared with high fidelity particle tracking simulations. The synchrotron radiation images of circulating beam allow direct measurements of the topology and lifetime in configurations where nonlinear focusing terms dominate.

Wieland, John [Michigan U.] (ORCID:000000028971852↗
DOE OSTI2023

Preliminary Design of a LAMP DTL (Update)

Present report is a continuation of the previously published report LA-UR-23-26600 from 19-Jun-2023 and incorporates the mentioned report as its part. All the additional information is contained in the part III of this report, the DTL Transverse Focusing Scheme. A preliminary Drift Tube Linac (DTL) layout was designed to create an algorithm for developing the conceptual design of the proposed LANSCE Accelerator Modernization Project (LAMP) final section: the proton linear accelerator from 3 MeV to 100 MeV. Preceding reports describe the proposed layouts of the LEBT, the RFQ, and the MEBT subsections of the linac. Initial estimates of the needed RF power for the DTL are also presented in the report. Present report describes the initial layout of the DTL and the first longitudinal beam dynamics (BD) results of the simulations in simplified models. The part III of this report includes the transverse focusing scheme and transverse BD in details.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Molten pool dynamics and humping suppression in high-speed laser welding via tailored beam configurations

High-speed laser welding is essential for increasing the production rate of fuel cell fabrication. However, when the welding speed exceeds a critical limit, humping occurs and reduces the weld quality. In this study, two tailored beam configurations, including an adjustable ring mode and a dual-beam configuration, were employed to suppress humping. Computational fluid dynamics simulations were performed to elucidate the underlying suppression mechanisms. Here, the results show that, in the adjustable ring mode, humping mitigation arises from a reduced backward cross-sectional melt flow rate and a more stable molten pool. In the dual-beam configuration, humping suppression is attributed to the deceleration of melt flow, the conduction-mode behavior of the trailing beam, and the widening of the molten pool induced by the trailing laser. Furthermore, because the dual-beam configuration directly modifies the trailing molten pool dynamics, it achieves more effective humping suppression, extending the welding speed limit to 1.50 m/s, compared with 1.00 m/s for the adjustable ring mode.

08 HYDROGEN↗
DOE OSTI2023

Energy Characterization of the UMD Electron Accelerator

Multiple different groups at Los Alamos National Laboratory have a cited interest in understanding the effects of ionizing radiation in materials. Of especial interest to the ISR division is understanding the effects in materials that are used in the space environment, where radiation is omnipresent and impossible to fully shield from. Development is under way to simulate the effects of radiation damage in detectors and electronics, primarily with particle transport codes such as Geant4 and RAMPART, a code that provides a user friendly interface to Geant4. However, these simulations need benchmarking to experimental tests, some of which have been performed at a variety of small electron accelerator facilities. One of these facilities is the Linear Electron Accelerator at the University of Maryland Radiation Facilities. This facility can provide electron beams on the order of 100 mA delivered in tight macropulse bunches, which allows for experiments where the behavior of electronics is monitored as a function of pulse number. While the beam rates per macropulse of this accelerator is monitored during an experiment via a current pick-off connected to a large graphite beamstop, the energy distribution of this electron beam is not measured during the experiment and has not been well characterized. The energy distribution of the beam is pertinent for these types of experiments, since dose is directly related to energy deposited in a medium. Further, these experiments are performed in open air and sample placement in the room can dictate what dose rate they receive per macropulse. We have created simulated dosemaps for the UMD vault room, but these simulations assume a given energy of the beam, both for calculating the dose and to simulate how much the beam diverges in the air. There have been two attempts now to characterize this accelerator beam energy with a plate type spectrometer, including the most recent experiment earlier this year to quantify the beam energy distribution. These results will also be compared to a dosimetry experiment performed by collaborators at the UMD facility, the details of which are described in Reference. Understanding the electron beam energy is important for future experiments, because the energy directly impacts the dose and radiation damage given to a sample and is important for properly simulating the electron beam when comparing simulation to experiment.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Impact of beam–beam effects on absolute luminosity calibrations at the CERN Large Hadron Collider

At the Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (vdM) method are affected by the mutual electromagnetic interaction of the two beams. The colliding bunches experience relative orbit shifts, as well as optical distortions akin to the dynamic-β effect, that both depend on the transverse beam separation and must therefore be corrected for when deriving the absolute luminosity scale. In the vdM regime, the beam–beam parameter is small enough that the orbit shift can be calculated analytically. The dynamic-β corrections to the luminometer calibrations, however, had until the end of Run 2 been estimated in the linear approximation only. In this report, the influence of beam–beam effects on the vdM-based luminosity scale is quantified, together with the associated systematic uncertainties, by means of simulations that fully take into account the non-linearity of the beam–beam force, as well as the resulting non-Gaussian distortions of the transverse beam distributions. Two independent multiparticle simulations, one limited to the weak-strong approximation and one that models strong-strong effects in a self-consistent manner, are found in excellent agreement; both predict a percent-level shift of the absolute pp-luminosity values with respect to those assumed until recently in the physics publications of the LHC experiments. These results also provide guidance regarding further studies aimed at reducing the beam–beam-related systematic uncertainty on beam–beam corrections to absolute luminosity calibrations by the van der Meer method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multiphysics simulation of recent experiments on alkali‐silica reaction expansion in reinforced concrete members

Alkali‐silica reaction (ASR) is an important degradation process that causes volumetric expansion and damage in concrete, and is affected significantly by the local temperature, moisture and stress conditions that often vary across the regions of a structure. Numerical simulation is essential to predict the progression and effects of ASR on the performance of structures. Because of the interactions between thermal and moisture transport and mechanical deformation, it is important for numerical models to represent all these physical phenomena and the coupling between them. Simulations of ASR in reinforced concrete (RC) structures are further complicated by the need to capture interactions between concrete and embedded reinforcing bars. Here, this paper describes the implementation of a scalable, coupled‐physics ASR model for simulating RC structures and assesses the ability of that model to predict ASR‐induced expansion in recent laboratory tests on RC block and beam specimens. These laboratory tests and the simulation approach were selected because of their applicability to RC structural‐scale simulations. This validation study helps builds confidence the ability of this approach to model ASR expansion in large, complex RC structures, which is a current high‐priority need.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Drift-cyclotron loss-cone instability in 3-D simulations of a sloshing-ion simple mirror

The kinetic stability of collisionless, sloshing beam-ion (45° pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over 1 – 10 μs, drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool ∼1keV ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Virtual Tokamak for Test and Development of Plasma Control Applied to NSTX-U

Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Vertical Septum Magnet Design for the APS Upgrade

The vertical injection scheme proposed for the APS Upgrade (APS-U) Project requires a challenging septum magnet that must meet stringent beam physics, magnetic field leakage, and vacuum requirements. The current iteration of this magnet design includes an enlarged stored-beam chamber aperture of 9 mm x 12 mm and a reduction of the septum thickness to 1.5 mm. The enlarged aperture accommodates a non-evaporable getter (NEG)-coated stored beam chamber to better achieve the required vacuum. A prototype septum magnet has been built and measurements confirm the cancellation of a peak leakage field even though the value is six times larger than the design. The leakage field measured at the upstream (US) end cancels the downstream (DS) end as was expected by design. The measured and simulated leakage field and the stored beam trajectories are reported.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Generation of High-Charge Magnetized Electron Beams Consistent With JLEIC Electron Cooling Requirements

The proposed Jefferson Lab Electron-Ion Collider (JLEIC), currently under design, relies on electron cooling in order to achieve the desired luminosity. This includes an electron beam with >55 Mev, 3.2 nC bunches that cools hadron beams with energies up to 100 GeV. To enhance the cooling, the electron beam must be magnetized with a specific eigen-emittance partition. This paper explores the use of the Fermilab Accelerator Science and Technology (FAST) facility to demonstrate the generation of an electron beam with parameters consistent with those required in the JLEIC high-energy cooler. We demonstrate via simulations the generation of the required electron-beam parameters and perform a preliminary experiment to validate FAST capabilities to produce such beams.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

New Simulations for Ion-Production and Back-Bombardment in GaAs Photo-guns

GaAs-based DC high voltage photo-guns used at accelerators with extensive user programs must exhibit long photocathode operating lifetime. Achieving this goal represents a significant challenge for proposed high average current facilities that must operate at tens of milliamperes or more. Specifically, the operating lifetime is dominated by ion back-bombardment of the photocathode from ionized residual gas. While numerous experiments have been performed to characterize the operating lifetime under various conditions, detailed simulations of the ion back-bombardment mechanism that explains these experiments are lacking. Recently, a new user routine was implemented using the code General Particle Tracer (GPT) to simulate electron impact ionization of residual beam line gas and simultaneously track the incident electron, the ejected electron, and the newly formed ion. This new routine was benchmarked against analytical calculations and then applied to experiments performed at the CEBAF injector at the Thomas Jefferson National Accelerator Facility to study the effectiveness of limiting ions from entering the cathode-anode gap using a positively biased anode. These simulations were performed using detailed 3D field maps produced with CST Microwave Studio describing the photo-gun electrostatics. Discussion of the experiment and the application of this new GPT routine to model the experiments will be presented at the workshop.

Yoskowitz, Joshua↗
DOE OSTI2021

Lattice Design of the Beam Transfer Line (BTL) from PIP-II LINAC to the Booster at Fermilab

PIP-II beam transfer line (BTL) to transport the beam from PIP-II Linac to the Booster ring at Fermilab. The latest design eliminates rolling the dipoles in the beam line to cross over the Tevatron tunnel. Also re-designed is the lattice in the region of the Booster Injection to meet the request of the civil construction needs and accommodate the constrains of the Booster injection request. A beam line to the beam absorber (BAL) is designed based on the request from the results of Mars simulations and ANASYS calculation of the absorber. Simulations with dipole and quadrupole field errors for the Beam Transport Line (BTL) to the Booster, which provides the specifications for all the magnets and Power supplies, will be presented too.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Single Shot Characterization of High Transformer Ratio Wakefields in Nonlinear Plasma Acceleration

We report that plasma wakefields can enable very high accelerating gradients for frontier high energy particle accelerators, in excess of 10 GeV/m. To overcome limits on single stage acceleration, specially shaped drive beams can be used in both linear and nonlinear plasma wakefield accelerators (PWFA), to increase the transformer ratio, implying that the drive beam deceleration is minimized relative to acceleration obtained in the wake. In this Letter, we report the results of a nonlinear PWFA, high transformer ratio experiment using high-charge, longitudinally asymmetric drive beams in a plasma cell. An emittance exchange process is used to generate variable drive current profiles, in conjunction with a long (multiple plasma wavelength) witness beam. The witness beam is energy modulated by the wakefield, yielding a response that contains detailed spectral information in a single-shot measurement. Using these methods, we generate a variety of beam profiles and characterize the wakefields, directly observing transformer ratios up to R=7.8. Lastly, a spectrally based reconstruction technique, validated by 3D particle-in-cell simulations, is introduced to obtain the drive beam current profile from the decelerating wake data.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Absorption of charged particles in perfectly matched layers by optimal damping of the deposited current

Perfectly matched layers (PMLs) are widely used in particle-in-cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. Here, in order to mitigate these artifacts, we introduce a PML algorithm whereby the current deposited by the macroparticles in the PML is damped by an analytically derived optimal coefficient. The benefits of this algorithm are illustrated in practical simulations. In particular, it is shown that this algorithm is well suited for particles exiting the box in near-normal incidence, in the sense that the fields behave as if the exiting particle is propagating in an infinite vacuum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗