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 1,117 records · Page 62

Ion injection simulations of quasi-parallel shock re-formation

One-dimensional hybrid simulations are used to investigate the process of quasi-parallel shock reformation and to examine the coupling of a beam of ions reflected at the shock to the incoming solar wind. A simple simulation configuration is constructed that makes it possible to control the properties of the background plasma and of the reflected ions. The length and time scales for the coupling of the reflected ions to the background plasma are investigated as functions of the upstream magnetic field direction, beam density, and beam temperature. The coupling length and time scales are found to vary systematically with the upstream magnetic field direction. The coupling occurs at roughly the time and location where the injected ions become deflected transverse to the shock normal direction.

Onsager, T. G.↗

Measurements of particle emission from discharge sites in Teflon irradiated by high energy electron beams

Anomalous behavior of synchronous orbit satellites manifested by overall degradation of system performance and reduced operating life is associated with electrical discharges resulting from differential charging of the spacecraft surface by fluxes of high energy electrons. During a laboratory simulation silver-backed Teflon samples have been irradiated by electron beams having energies in the range 16-26 keV. Charged particles emitted from the resultant electrical discharges have been measured with a biased Faraday cup and retarding potential analyser. Measurements indicate the presence of two distinct fluxes of particles, the first being an early pulse (0-600ns) of high energy (about 7keV) electrons, while the second is a late pulse (1-5 microseconds) of low energy electrons (less than 1eV) and ions (70eV) leaving the discharge site as a quasi plasma. Calculations indicate an electrostatic field as the dominant accelerating mechanism for charged particles.

Hazelton, R. C.↗

Simulation Validation for the DUNE Near Detector

DUNE will study neutrino oscillations as they are beamed from Fermilab in Batavia, Illinois, 1300 km away to Sanford Underground Research Facility in South Dakota. However, when preparing for an experiment of this scale, it is crucial that we run simulations to validate the geometry and overall design of the project. This poster is comprised of my work during a 2026 SULI internship at Fermilab, involving reading the simulated data and assessing its validity.

Vershaw, Andre [Unlisted, US, IL] (ORCID:000900084↗

Particle-in-cell simulation of laser wakefield accelerators with oblique lasers in quasicylindrical geometry

In the studies of optical injections in laser or beam-driven wakefield accelerators, there is a frequent demand for using obliquely propagating assistant lasers in particle-in-cell simulations. In conventional methods, this is only possible in either two- or three-dimensional Cartesian geometries, which have the drawbacks of either lack of fidelity for the actual situation or requiring a huge amount of computational resources. In this work, we develop a new method that uses an expression-defined oblique laser field to simulate such situations in quasicylindrical geometry particle-in-cell simulations, having the advantages of both maintaining good fidelity and saving computational resources. As an example, we use this method in the scissor-cross ionization injection scheme for the optimization of the injected beam quality. This method is widely applicable to particle injections with assistant lasers in wakefield accelerators, as long as the assistant lasers only influence the injected particles during the injection process. Published by the American Physical Society 2025

Ma, Minghao (ORCID:0009000602317775)↗

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mixed Effect Modeling of Dose and Linear Energy Transfer Correlations With Brain Image Changes After Intensity Modulated Proton Therapy for Skull Base Head and Neck Cancer

Intensity modulated proton therapy (IMPT) could yield high linear energy transfer (LET) in critical structures and increased biological effect. For head and neck cancers at the skull base this could potentially result in radiation-associated brain image change (RAIC). The purpose of the current study was to investigate voxel-wise dose and LET correlations with RAIC after IMPT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Spatially dependent modeling and simulation of runaway electron mitigation in DIII-D

New simulations with the Kinetic Orbit Runaway electron (RE) Code (KORC) show RE deconfinement losses to the wall during plasma scrape off are the primary current dissipation mechanism in DIII-D experiments with high-Z impurity injection, and not collisional slowing down. The majority of simulations also exhibit an increase in the RE beam energy due to acceleration by the induced toroidal electric field, even while the RE beam current is decreasing. In this study, KORC integrates RE orbits using the relativistic guiding center equations of motion and incorporates time-sequenced, experimental reconstructions of the magnetic and electric fields and line integrated electron density to construct spatiotemporal models of electron and partially ionized impurity transport in the companion plasma. Comparisons of experimental current evolution and KORC results demonstrate the importance of including Coulomb collisions with partially ionized impurity physics, initial RE energy, pitch angle, and spatial distributions, and spatiotemporal electron and partially ionized impurity transport. The research presented here provides an initial quantification of the efficacy of RE mitigation via injected impurities and identification of the critical role played by loss of confinement due to plasma scrape off on the inner wall as compared to the relatively slow collisional damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of laser ablation of silicon as a function of pulse length at constant fluence via time-resolved x-ray spectroscopy

We investigate the ablation of silicon as a function of laser pulse length at a constant fluence using time-resolved x-ray spectroscopy data obtained from OMEGA EP experiments at the University of Rochester's Laboratory for Laser Energetics. Our targets consisted of three-layer planar structures composed of Si (50 μm), Cu (25 μm), and SiO 2 (500 μm) layers. The Si layer was irradiated by a 351-nm laser with varying pulse widths of 250 ps, 500 ps, 1 ns, and 10 ns while maintaining a constant fluence of ~27.9 kJ/cm 2 . Electron temperatures and densities of the ablated plasma were determined by analyzing the time-resolved x-ray spectroscopy data through a comparison of experimental measurements with synthetic results obtained from Si atomic calculations in a steady state and non-local thermodynamic equilibrium. These calculations were computed using PrismSPECT. Additionally, radiation-hydrodynamics simulations with FLASH are used to generate simulated plasma-density and plasma-temperature profiles, which are then compared with the experimental measurements. Our analyses reveal that increasing the laser pulse length at a constant fluence results in a decrease in electron temperatures and densities. Furthermore, the longer pulses with lower intensities lead to deeper ablation regions before reaching the peak ablation but lower ionization balances in the silicon layer. Furthermore, these findings emphasize the critical role of laser pulse length in plasma ablation and shock generation for laser-impulse studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The FRB-searching Pipeline of the Tianlai Cylinder Pathfinder Array

This paper presents the design, calibration, and survey strategy of the Fast Radio Burst (FRB) digital backend and its real-time data processing pipeline employed in the Tianlai Cylinder Pathfinder Array. The array, consisting of three parallel cylindrical reflectors and equipped with 96 dual-polarization feeds, is a radio interferometer array designed for conducting drift scans of the northern celestial semi-sphere. The FRB digital backend enables the formation of 96 digital beams, effectively covering an area of approximately 40 square degrees with the 3 dB beam. Our pipeline demonstrates the capability to conduct an automatic search of FRBs, detecting at quasi-real-time and classifying FRB candidates automatically. The current FRB searching pipeline has an overall recall rate of 88%. During the commissioning phase, we successfully detected signals emitted by four well-known pulsars: PSR B0329+54, B2021+51, B0823+26, and B2020+28. We report the first discovery of an FRB by our array, designated as FRB 20220414A. We also investigate the optimal arrangement for the digitally formed beams to achieve maximum detection rate by numerical simulation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Startup regime of high-efficiency tapering-enhanced FEL oscillator

In this paper, we present a design of a high-efficiency high-gain free-electron laser oscillator based on the use of a strongly tapered undulator for extracting energy from high-brightness electron beams. We provide an analytical model of the setup followed by numerical simulations for lasing at the wavelength of 13.5 nm. We discuss the optimization of the system in steady state and the conditions necessary for the pass-per-pass buildup of the power from shot noise level. We propose the use of fast phase shifters as a way to accelerate the buildup. Finally, we present time-dependent simulations of the oscillator and discuss the role of spectral filtering. The optimized working point yields a total energy conversion efficiency from the electron beam to output radiation above 1% at the wavelength of 13.5 nm.

Beam dynamics↗

Demonstration of a ring-FEL as an EUV lithography tool

In this work, we present the required structure and function of a ring-FEL as a radiation source for extreme ultraviolet radiation lithography (EUVL). A 100 m-long straight section that conducts an extremely low emittance beam from a fourth-generation storage ring can increase the average power at 13.5 nm wavelength to up to 1 kW without degrading the beam in the rest of the ring. Here, simulation results for a ring-FEL as a EUVL source are described.

07 ISOTOPE AND RADIATION SOURCES↗

Mbtrack2 - Application on EIC 5GE V Electron Ring Reverse Phase Configuration

In this paper, we present a report on simulations conducted using the Mbtrack2 code for the 5GeV EIC Electron Storage Ring (ESR) RF system with a reversed phasing scheme. The Mbtrack2 code [1] has been thoroughly benchmarked, as described in references [2] and [3]. The simulation results presented in this paper include the beam loading effect but do not incorporate a feedback system.

43 PARTICLE ACCELERATORS↗

ACHILLES-GENIE Interface for Neutrino Simulations, ADRIANO2 Tile Prototype for High-Granularity Dual-Readout Calorimetry

ICARUS (Imaging Cosmic And Rare Underground Signals) is a liquid argon time projection chamber (LArTPC) detector that pursues the sterile neutrino, which relies on accurate simulations of neutrino-argon interactions. REDTOP (Rare Eta Decays To Observe new Physics) is a proposed low-energy, high-intensity meson factory designed to explore rare $\eta$/$\eta'$ meson decays and probe physics beyond the Standard Model. As a next-generation experiment, this requires both accurate simulations and innovative detector technologies. This project contributes to both ICARUS, from a simulation perspective, and REDTOP, from both a simulation and detection perspective, through the event generation of lepton-nucleon interactions and the physical enhancement of the calorimeter technology within the REDTOP detector. We developed an interface between ACHILLES (A CHIcago Land Lepton Event Simulator), a theory-driven lepton-level event generator, and GENIE, a robust event generator framework used for neutrino physics. By incorporating the precise theoretical cross-section calculations of ACHILLES into the experimental realism of GENIE, the interface allows for improved accuracy of neutrino-nucleon simulations, which can be adapted for the proton beam specifications of the REDTOP meson factory as well as for the ICARUS experiment. In parallel, we developed an improved prototype for the ADRIANO2 (A Dual Readout Integrally Active Non-segmented Option) dual-readout calorimeter tiles for the REDTOP detector. To improve the efficiency of the lead-glass tiles trapping Cherenkov light for energy reconstruction and particle identification, we optimized the application of a highly reflective coating. Through viscosity and thickness control, masking, and a custom spray technique, we refined the coating process to reduce surface defects and improve light yield. Together, these efforts strengthen the ICARUS neutrino program and REDTOP's capability of detecting rare decay events.

Visser, Erin [Michigan State U.]↗

Multibeam antenna study, Phase 2

An antenna system concept for point to point communication between any two points in the Continental United States was developed using a synchronous orbit satellite. The design and fabrication of a two dimensional artificial dielectric lens contained in a parallel plate waveguide system is presented. The experimental pattern and voltage standing wave ratio results are listed from tests with four different lenses: a plano-convex natural dielectric (Teflon) lens, a plano-convex artificial dielectric lens, a two-point corrected Teflon lens, and a two-point corrected artificial dielectric lens. A computer simulated design of an operational two lens, 24 beam configuration is also presented. The computer design presents overlays of each beam onto the Continental United States from a fixed synchronous orbit satellite with a definition of these areas where beam-to-beam isolation is not sufficient to permit the reuse of frequency.

Byrnes, P. J.↗

Two-dimensiosnal electron beam charging model for polymer films

A two dimensional model was developed to describe the charging of thin polymer films exposed to a uniform mon-energetic electron beam. The study was motivated by observed anomalous behavior of geosynchronous satellites which was attributed to electrical discharges associated with the differential charging of satellite surfaces of magnetospheric electrons. Electric fields both internal and external to the irradiated specimen were calculated at steady state in order to identify regions of high electrical stress. Particular emphasis was placed on evaluating the charging characteristics near the material's edge. The model was used to identify and quantify the effects of some of the experimental parameters notably: beam energy; beam angle of incidence; beam current density; material thickness; and material width. Simulations of the following situations were also conducted: positive or negative precharging over part of the surface; a central gap in the material; and a discontinuity in the material's thickness.

Reeves, R. D.↗

On-line structural parameter identification

Algorithms are presented for on-line parameter identification of structural dynamic systems. As an example, they are used to calculate the parameters of a modal model of a flexible beam. The algorithms are tested using hardware consisting of a 12 ft. beam with four voice coil actuators and nine noncontacting displacement sensors. They are programmed in a CDC Cyber 175 digital computer which provides input command signals for the actuators, reads the sensor data, and processes the algorithm to calculate consistent estimates of the modal parameters of the beam. Experimental results are compared with those of simulation analysis.

Thau, F. E.↗

Scattering and collapse of Langmuir waves driven by a weak electron beam

Wave-wave interactions are examined, particularly for significance in relation to Langmuir turbulence associated with type III solar radio emissions. The Zakharov equations are considered as two coupled nonlinear partial differential equations for the low-frequency density and the Langmuir envelope, and the electron and ion temperatures are regarded as equal. The Zhakharov equations are solved in one dimension, and display no collapse of a driven broadband pump, indicating that scattering instabilities are important in the saturation of a beam-plasma instability. A two-dimensional numerical simulation is developed to interpret the role of parametric instabilities, and it is concluded that the distortion of a wave packet is due to scattering instability, and not a change in the physics of direct collapse due to an artificial numerical damping. Further discussions are presented for comparisons of models of Langmuir turbulence in type III bursts.

Hafizi, B.↗