Search NASA⌕ Search

SEARCH · Search NASA

Results for “beam 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 307 records · Page 17

WE-EF-303-06: Feasibility of PET Image-Based On-Line Proton Beam-Range Verification with Simulated Uniform Phantom and Human Brain Studies

Purpose: To study the feasibility of clinical on-line proton beam range verification with PET imaging Methods: We simulated a 179.2-MeV proton beam with 5-mm diameter irradiating a PMMA phantom of human brain size, which was then imaged by a brain PET with 300*300*100-mm{sup 3} FOV and different system sensitivities and spatial resolutions. We calculated the mean and standard deviation of positron activity range (AR) from reconstructed PET images, with respect to different data acquisition times (from 5 sec to 300 sec with 5-sec step). We also developed a technique, “Smoothed Maximum Value (SMV)”, to improve AR measurement under a given dose. Furthermore, we simulated a human brain irradiated by a 110-MeV proton beam of 50-mm diameter with 0.3-Gy dose at Bragg peak and imaged by the above PET system with 40% system sensitivity at the center of FOV and 1.7-mm spatial resolution. Results: MC Simulations on the PMMA phantom showed that, regardless of PET system sensitivities and spatial resolutions, the accuracy and precision of AR were proportional to the reciprocal of the square root of image count if image smoothing was not applied. With image smoothing or SMV method, the accuracy and precision could be substantially improved. For a cylindrical PMMA phantom (200 mm diameter and 290 mm long), the accuracy and precision of AR measurement could reach 1.0 and 1.7 mm, with 100-sec data acquired by the brain PET. The study with a human brain showed it was feasible to achieve sub-millimeter accuracy and precision of AR measurement with acquisition time within 60 sec. Conclusion: This study established the relationship between count statistics and the accuracy and precision of activity-range verification. It showed the feasibility of clinical on-line BR verification with high-performance PET systems and improved AR measurement techniques. Cancer Prevention and Research Institute of Texas grant RP120326, NIH grant R21CA187717, The Cancer Center Support (Core) Grant CA016672 to MD Anderson Cancer Center.

60 APPLIED LIFE SCIENCES↗

DSMC simulations of polarized atomic beam sources including magnetic fields

In recent decades a lot of work has been done to understand and optimize the output of polarized atomic beam sources. However, there are many effects which prevent a complete analytic description of the system. Now, a DSMC (Direct Simulation Monte Carlo ) simulation based on OpenFOAM 1.7.1 using Birds algorithm is set up. So far, the simulated particles have been given spin and a generic utility to include arbitrary magnetic fields has been created. The equation of motion in every timestep is solved by a fourth order Runge-Kutta scheme. In this contribution, we describe promising first results of the simulations. Additionally a tool to measure the collision age has been created. Further effects to be included are recombination on the walls and spin exchange collisions. After that an algorithm will be included to optimize the output of the atomic beam source.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Simulation of supersonic molecular beam injection fueling into H-mode plasmas on EAST using BOUT++

Transport dynamic 2D simulations for supersonic molecular beam injection (SMBI) fueling into H-mode deuterium plasmas on the Experimental Advanced Superconducting Tokamak (EAST) is first simulated using a seven-field two-fluid model in the BOUT++ framework. The SMB is assumed to be injected into plasma from the midplane at the low field side with a fixed width and constant molecular flux. The different densities and injection velocities of SMBI are investigated within the upper single-null geometry of EAST. The simulations indicate that the SMBI has a self-shielding effect on molecules' inward transport into the plasma, and the deposition of SMBI leads to a large increase in plasma density and decrease in plasma temperature. There is a velocity threshold for SMB penetrating the pedestal and depositing at the top of density pedestal. However, the deposition point would be back toward the plasma boundary after SMB arriving at the deepest penetration position. Comparing the different molecular injection velocities and densities, the outcomes demonstrate that the depth of deposition is closely related to the injection velocity rather than the molecular density. The simulated results show good agreement with the EAST experiments by comparing the electron density profiles obtained by simulation and experiment, respectively. These results will be helpful for guiding future experiments, as well as the design of the SMBI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Geometry Effects on Multipole Components and Beam Optics in High-Velocity Multi-Spoke Cavities

Velocity-of-light, multi-spoke cavities are being proposed to accelerate electrons in a compact light-source. There are strict requirements on the beam quality which require that the linac have only small non-uniformities in the accelerating field. Beam dynamics simulations have uncovered varying levels of focusing and defocusing in the proposed cavities, which is dependent on the geometry of the spoke in the vicinity of the beam path. Here we present results for the influence different spoke geometries have on the multipole components of the accelerating field and how these components, in turn, impact the simulated beam properties.

43 PARTICLE ACCELERATORS↗

Energy Production Demonstrator and Material Testing Station optimization for Megawatt proton beams

A simulation study of the Energy Production Demonstrator (EPD) concept – a solid heavy metal target irradiated by GeV-range intense proton beams and producing more energy than consuming – is carried out. Neutron production, fission, energy deposition, energy gain, testing volume and helium production are simulated using the MARS15 code for tungsten, thorium, and natural uranium targets in the proton energy range 0.5–120 GeV. This study shows that the proton energy range of 2–4 GeV is optimal for both a nat U EPD and the tungsten-based testing station for proton accelerator facilities. Simulation-based conservative estimates not including breeding and fission of plutonium suggest the proton beam current sufficient to produce 1 GW of thermal output power with the nat U EPD while supplying a relatively small fraction of that power to operate the accelerator. The thermal analysis has been performed and has shown that the EPD with a parallel proton beam feeding the core has a potential problem due to a possible core meltdown. In conclusion, a scheme has been proposed with a beam steering on the outer surface of the target; its thermal analysis indicates that slicing the target in shorter parts and maintaining a temperature of 25 °C on their surfaces would to avoid the target meltdown.

43 PARTICLE ACCELERATORS↗

Exact wave solver for nonparaxial laser beam propagation

Simulations of inertial confinement fusion (ICF) experiments require high-fidelity models for laser beam propagation in a nonuniform plasma with varying index of refraction. We describe a new numerical wave solver that is applicable to centimeter-scale length plasmas encountered in indirect drive ICF applications. The one-way Helmholtz equation (OHE) generalizes the time-harmonic paraxial wave equation to large angles. Here, we present a methodology to numerically evaluate the exact solution to the OHE. This solution is computed by analytically advancing eigenfunctions of the one-way Helmholtz operator along a propagation direction and is applicable to any given index of a refraction profile. We compare our exact method with a commonly used approximate split-step technique for solving the OHE. As a test problem, we consider nonparaxial propagation of Gaussian and speckled beams in a plasma density channel with internal reflection. We find that the split-step approach incurs significant errors compared to the exact solution computed using the novel algorithm.

Belyaev, Mikhail A. (ORCID:0000000224908887)↗

Simulations of pulsed electron beam injection during active experiments

Two-dimensional electromagnetic particle simulations are used to investigate the characteristics of the spacecraft charging, particle acceleration, and wave emissions during beam injection, turnoff and subsequent pulsing of the beam. It is shown that, during beam injection, the beam current is neutralized by a spatially separate return current region extending several tens of meters from the beam region, with the currents being closed across the field lines by the perpendicular acceleration of ambient plasma ions into the beam region. After beam turnoff, this current system reverses after a time lag of about an ion plasma period. The current reversal is accompanied by prolonged electron collection by the spacecraft near the beam region, preferential ion collection by sections of the spacecraft magnetically connected to the initial return current regions, and the creation of hot plasma extending well into the return current regions. Because of the time lag, the currents induced in the plasma during periodic beam injection can be strongly modified from the imposed beam current.

Winglee, R. M.↗

Theory and simulations of cross-beam energy transfer between speckled laser beams

In this paper, we introduce a realistic modeling of laser beams smoothed by random phase plates to investigate cross-beam energy transfer (CBET). We compute the total field of the smoothed laser beams, which allows us to perform accurate simulations and to build a model accounting for the realistic structure of the speckles. This lets us confirm the results already discussed in a previous work given by Oudin et al. [Phys. Plasmas 29, 112112 (2022)]. When CBET is induced by a plasma flow, plane wave models can predict the exchange between smoothed laser beams. By contrast, when the exchange is caused by a wavelength shift between the laser fields of each beam, plane wave models overestimate the energy transfer, even for strongly Landau-damped ion acoustic waves. We also show that the laser angular spread, which results from its focusing, affects the resonance between the ponderomotive force induced by the laser fields and the ion acoustic wave. It increases the resonance width, which weakly depends on the Landau damping rate when there is a wavelength shift and no flow. In addition, the power transfer obtained with or without a phase plate is essentially the same in the perturbative regime and for small crossing angles.

Laser plasma interactions↗

Breakdown simulations in a focused microwave beam within the simplified model

The simplified model is proposed to simulate numerically air breakdown in a focused microwave beam. The model is 1D from the mathematical point of view, but it takes into account the spatial non-uniformity of microwave field amplitude along the beam axis. The simulations are completed for different frequencies and different focal lengths of microwave beams. The results demonstrate complicated regimes of the breakdown evolution which represents a series of repeated ionization waves. These waves start at the focal point and propagate towards incident microwave radiation. The ionization wave parameters vary during propagation. At relatively low frequencies, the propagation regime of subsequent waves can also change qualitatively. Each next ionization wave is less pronounced than the previous one, and the breakdown evolution approaches the steady state with relatively small plasma density. The ionization wave parameters are sensitive to the weak source of external ionization, but the steady state is independent on such a source. As the beam focal length decreases, the stationary plasma density increases and the onset of the steady state occurs faster.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations relevant to the beam instability in the foreshock

The results presently obtained from two-dimensional simulations of the reactive instability for Maxwellian beams and cutoff distributions are noted to be consistent with recent suggestions that electrons backstreaming into earth's foreshock have steep-sided cutoff distributions, which are initially unstable to the reactive instability, and that the back-reaction to the wave growth causes the instability to pass into its kinetic phase. It is demonstrated that the reactive instability is a bunching instability, and that the reactive instability saturates and passes over into the kinetic phase by particle trapping.

Cairns, I. H.↗

Space Charge Simulations of High Intensity Proton Beams in the AGS Booster

Computer simulation studies have been performed to understand the beam behavior and to explore intensity limitations of proton beams in the AGS Booster at higher beam intensities. During the 100 GeV polarized proton operations of RHIC Run 2024, sPHENIX operated in modes with a crossing angle at collisions in order to mitigate beam-beam effects. Three different running modes were employed: (a) sPHENIX operated with a negative (-2 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. (b) sPHENIX was brought into collisions with 0 mrad first. Then STAR was brought into collisions after the beam-beam parameter from sPHENIX reduced to below $10 \times 10^{-3}$. (c) sPHENIX operated with a positive (+1.5 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. The collisions with a crossing angle of up to $\pm 2$ mrad, as in running modes (a) and (c), lead to large Piwinski angle in the new sPHENIX detector, which reduces luminosity if other parameters are unchanged. There are two ways to compensate the reduction in luminosity: squeeze $\beta^{*}$ if there is sufficient dynamic aperture, or increase the injected beam intensity. The first part of polarized proton operation during RHIC Run 2024 was dedicated to increasing the intensity. Different configurations were tested with crossing angle and lattice adjustments on RHIC. At the same time, new injector configurations were developed and tested in an effort to push for both higher intensity and better quality of the beam injected into RHIC. When the beam intensity is increased, space charge is a concern particularly in the lower energy stages of acceleration, such as during the injection and the early part of the Booster cycle, which could become a dominant effect in limiting the intensity of the beam that can be delivered to RHIC.

43 PARTICLE ACCELERATORS↗

Start-to-End Simulation of the Drive-Beam Longitudinal Dynamics for Beam-Driven Wakefield Acceleration

Collinear beam-driven wakefield acceleration (WFA) relies on shaped driver beam to provide higher accelerating gradient at a smaller cost and physical footprint. This acceleration scheme is currently envisioned to accelerate electron beams capable of driving free-electron laser *. Start-to-end simulation of drive-bunch beam dynamics is crucial for the evaluation of the design of accelerators built upon WFA. We report the start-to-end longitudinal beam dynamics simulations of an accelerator beamline capable of producing high charge drive beam. The generated wakefield when it passes through a corrugated waveguide results in a transformer ratio of 5. This paper especially discusses the challenges and criteria associated with the generation of temporally-shaped driver beam, including the beam formation in the photoinjector, and the influence of energy chirp control on beam transport stability.

43 PARTICLE ACCELERATORS↗

Aligning the Thomson scattering and charge exchange recombination diagnostics using neutral beam emission at DIII-D

This work addresses discrepancies in the alignment of the H-mode pedestal profiles of the electron and ion properties in the DIII-D tokamak as measured by Thomson Scattering (TS) and Charge Exchange Recombination Spectroscopy (CER) diagnostics. While the alignment of these profiles is key for accurate studies of tokamak physics and plasma confinement, misalignments can occur due to inaccuracies, such as in magnetic equilibrium reconstructions required to map measurements in different poloidal and toroidal locations. Both FIDASIM, an established simulation package, and a simplified collisional radiative model are used to simulate neutral beam state densities and neutral beam emission. Simulated neutral beam emissions are calculated based on shifted TS profiles and compared to beam emission measurements from the Main Ion CER system to determine the best shift for aligning TS with CER. This analysis is performed on various DIII-D discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron-beam-ion-source (EBIS) modeling progress at FAR-TECH, Inc

FAR-TECH, Inc. has been developing a numerical modeling tool for Electron-Beam-Ion-Sources (EBISs). The tool consists of two codes. One is the Particle-Beam-Gun-Simulation (PBGUNS) code to simulate a steady state electron beam and the other is the EBIS-Particle-In-Cell (EBIS-PIC) code to simulate ion charge breeding with the electron beam. PBGUNS, a 2D (r,z) electron gun and ion source simulation code, has been extended for efficient modeling of EBISs and the work was presented previously. EBIS-PIC is a space charge self-consistent PIC code and is written to simulate charge breeding in an axisymmetric 2D (r,z) device allowing for full three-dimensional ion dynamics. This 2D code has been successfully benchmarked with Test-EBIS measurements at Brookhaven National Laboratory. For long timescale (< tens of ms) ion charge breeding, the 2D EBIS-PIC simulations take a long computational time making the simulation less practical. Most of the EBIS charge breeding, however, may be modeled in 1D (r) as the axial dependence of the ion dynamics may be ignored in the trap. Where 1D approximations are valid, simulations of charge breeding in an EBIS over long time scales become possible, using EBIS-PIC together with PBGUNS. Initial 1D results are presented. The significance of the magnetic field to ion dynamics, ion cooling effects due to collisions with neutral gas, and the role of Coulomb collisions are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using kernel-based statistical distance to study the dynamics of charged particle beams in particle-based simulation codes

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. Here, the focus is on kernel-based methods such as maximum mean discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

47 OTHER INSTRUMENTATION↗

An Updated Simulation of the Booster Neutrino Beam

For accelerator neutrino experiments, an accurate prediction of the incoming neutrino flux is crucial for reducing uncertainties for all physics measurements. In this exciting period for the Short-Baseline Neutrino program at Fermilab, with far detector (ICARUS) already operating and the near detector (SBND) nearing operation, an updated flux model for the Booster Neutrino Beam (BNB) is presented. The BNB delivers 8 GeV protons to a beryllium target, subsequently producing neutrinos from the decay of the secondary beam of hadrons. A updated Monte-Carlo simulation of beam production in GEANT4 has been created, allowing predictions to be made for detectors with different baselines, offsets and sizes. This new simulation contains new features, such as a full neutrino ancestry to handle hadron production systematics with more precision, with a view to storing all resulting particles - including neutral mesons - from the proton-Beryllium scatter to allow the study of exotic BSM scenarios. Results are presented, with comparisons to the original flux simulated for the MiniBooNE experiment.

Paton, Josephine↗