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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experiment Simulation Configurations Approximating DUNE TDR

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long-baseline neutrino oscillation experiment consisting of a high-power, broadband neutrino beam, a highly capable near detector located on site at Fermilab, in Batavia, Illinois, and a massive liquid argon time projection chamber (LArTPC) far detector located at the 4850L of Sanford Underground Research Facility in Lead, South Dakota. The long-baseline physics sensitivity calculations presented in the DUNE Physics TDR, and in a related physics paper, rely upon simulation of the neutrino beam line, simulation of neutrino interactions in the near and far detectors, fully automated event reconstruction and neutrino classification, and detailed implementation of systematic uncertainties. The purpose of this posting is to provide a simplified summary of the simulations that went into this analysis to the community, in order to facilitate phenomenological studies of long-baseline oscillation at DUNE. Simulated neutrino flux files and a GLoBES configuration describing the far detector reconstruction and selection performance are included as ancillary files to this posting. A simple analysis using these configurations in GLoBES produces sensitivity that is similar, but not identical, to the official DUNE sensitivity. DUNE welcomes those interested in performing phenomenological work as members of the collaboration, but also recognizes the benefit of making these configurations readily available to the wider community.

Abi, B↗
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗
DOE OSTI2023

Principal Component Analysis of azimuthal flow in intermediate-energy heavy-ion reactions

Principal Component Analysis (PCA) via Singular Value Decomposition (SVD) of large datasets is an adaptive exploratory method to uncover natural patterns underlying the data. Several recent applications of the PCA-SVD to event-by-event single-particle azimuthal angle distribution matrices in ultra-relativistic heavy-ion collisions at RHIC-LHC energies indicate that the sine and cosine functions chosen a priori in the traditional Fourier analysis are naturally the most optimal basis for azimuthal flow studies according to the data itself. We perform PCA-SVD analyses of mid-central Au+Au collisions at $E$ $beam$ / $A$ =1.23 GeV simulated using an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model to address the following two questions: (1) if the principal components of the covariance matrix of nucleon azimuthal angle distributions in heavy-ion reactions around 1 GeV/nucleon are naturally sine and/or cosine functions and (2) what if any advantages the PCA-SVD may have over the traditional flow analysis using the Fourier expansion for studying the EOS of dense nuclear matter. In conclusion, we find that (1) in none of our analyses the principal components come out naturally as sine and/or cosine functions, (2) while both the eigenvectors and eigenvalues of the covariance matrix are appreciably EOS dependent, the PCA-SVD has no apparent advantage over the traditional Fourier analysis for studying the EOS of dense nuclear matter using the azimuthal collective flow in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗
DOE OSTI2026

Detector alignment for X-ray crystallography using Millepede-II

I describe a method for accurately refining the geometrical parameters of segmented X-ray area detectors on the basis of serial crystallography data, using 'Millepede' – an algorithm created for a very similar problem in high-energy physics. The Millepede method for serial crystallography builds on the approach of Brewster et al. [Acta Cryst. (2018), D74, 877–894], in which the detector parameters are refined simultaneously with the parameters for each individual crystal. This accounts for the mutual dependency between the parameters and thereby avoids the bias and slow convergence problems that have afflicted older approaches in which the deviations between observed and calculated Bragg peak positions were taken directly as the updates for the detector panel positions. The Millepede method uses the special structure of the least-squares normal equations to reduce them to a much smaller form that can be solved very quickly, even compared with the sparse matrix methods used previously. This makes it practical to refine the detector geometry frequently and thereby maintain accurate calibration without specialized alignment campaigns. Tilts of detector panels out of the plane can be reliably refined, as can the overall distance of the detector in the beam direction. With a simulated test case, the new method produced panel shifts within 7% of the correct values with only one iteration, and produced almost exactly correct shifts after a second iteration. A simulated out-of-plane panel rotation was correctly determined to within 0.001°. Applied to experimental data from an X-ray free-electron laser, the method increased the indexable fraction of frames from 30% to 91% in a single iteration, and to 96% after two further iterations. Computing the geometry updates on the basis of 2060 crystals took only 0.819 s on desktop computing hardware, including the time taken to read the required data from disk. The scaling was found to be very close to linear for up to 100 980 sets of crystal parameters, which took only 78.2 s to process under the same conditions. The method has been applied as part of a real-time feedback system at a synchrotron radiation beamline, in which an out-of-plane detector tilt of 0.04° was detected and corrected. Possible further applications are also described here.

Millepede-II↗
DOE OSTI2026

Consist v0.1.0

A Python library for provenance tracking, intelligent caching, and data virtualization in scientific simulation workflows. It automatically records code, configuration, and input data to skip redundant computations and enables querying results across many runs without manual bookkeeping. Designed to support multi-model simulation workflows like the BEAM CORE toolset at LBL, but designed to be extensible to a wide range of research workflows. Combines lineage tracking features as provided by OpenLineage with deterministic hashing like SnakeMake, and adds powerful analysis tools on model outputs.

Needell, Zachary [Lawrence Berkeley National Labor↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Direct Acceleration of an Electron Beam with a Radially Polarized Long-Wave Infrared Laser

Direct laser acceleration with radially polarized lasers is an intriguing variant of laser-based particle acceleration that has the potential of offering GeV/cm-level energy while avoiding the instabilities and complex beam dynamics associated with plasma wakefield accelerators. A major limiting factor is the difficulty of generating high-power radially polarized beams. In this paper, we propose the use of CO2-based long-wave infrared (LWIR) lasers as a driver for direct laser acceleration, as the polarization insensitivity of the gain medium allows a radially polarized beam to be amplified. Additionally, the larger waist sizes, Rayleigh lengths, and pulse lengths associated with the long wavelength could improve the injection efficiency of the electron beam. By comparing acceleration simulations using a near-infrared laser and an LWIR laser, we show that the injection efficiency is indeed improved by up to an order of magnitude with the longer wavelength. Furthermore, we show that even sub-TW peak powers with an LWIR laser can provide MeV-level energy gains. Thus, radially polarized LWIR lasers show significant promise as a driver of a direct laser-driven demonstration accelerator.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Neutron Calibration Studies for the SuperCDMS SNOLAB Experiment

My time at Fermilab has been dedicated to advancing the field of dark matter detection through my contributions to the SuperCDMS collaboration. The first project involved the development and testing of a delivery system for a Cf-252 neutron source to calibrate the SuperCDMS silicon and germanium detectors. Subsequently, I focused on parametrizing the detectors' instrumental response by optimizing a plastic scintillator backing array to measure the scattering angle of a collimated neutron beam. The Monte Carlo simulation I designed provided insights into the detection efficiency of the array, and the best-fit design I arrived at for the scattering counters will be implemented into the final array. Additionally, I assembled a muon detector with legacy hardware to enhance the ability of the calibration experiment to reduce dark current and false signals. As my work nears completion, I am confident that these contributions will significantly enhance the precision and reliability of the SuperCDMS experiment, pushing the boundaries of our understanding of dark matter and contributing to the broader scientific community's quest to unravel the mysteries of the universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

MAIA: A new detector concept for a 10 TeV muon collider

Muon colliders offer a compelling opportunity to explore the TeV scale and conduct precision tests of the Standard Model, all within a relatively compact geographical footprint. This paper introduces a new detector concept, MAIA (Muon Accelerator Instrumented Apparatus), optimized for $\sqrt{s}=10$ TeV $\mu\mu$ collisions. The detector features an all-silicon tracker immersed in a 5T solenoid field. High-granularity silicon-tungsten and iron-scintillator calorimeters surrounding the solenoid capture high-energy electronic and hadronic showers, respectively, and support particle-flow reconstruction. The outermost subsystem comprises an air-gap muon spectrometer, which enables standalone track reconstruction for high-momentum muons. The performance of the MAIA detector is evaluated in terms of differential particle reconstruction efficiencies and resolutions. Beam-induced background (BIB) simulations generated in FLUKA are overlaid with single particle gun samples to assess detector reconstruction capabilities under realistic experimental conditions. Even with BIB, reconstruction efficiencies exceed 95% for energetic tracks, photons, and neutrons in the central region of the detector. This paper outlines promising avenues of future work, including forward region optimization and opportunities for enhanced flavor/boosted object tagging, and addresses the technological assumptions needed to achieve the desired detector performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Increased spatial coherence length from an asymmetric crystal reflection at grazing exit

Coherent X-ray imaging is an active field at synchrotron sources. The images rely on the available coherent flux over a limited field of view. At many synchrotron beamlines a double-crystal monochromator (DCM) is employed in a standard nondispersive arrangement. For coherent diffraction imaging it is advantageous to increase the available field of view by increasing the spatial coherence length (SCL) of a beam exiting such a DCM. Here, Talbot interferometry data together with ray-tracing simulations for a (+ − − +) four-reflection experimental arrangement are presented, wherein the first two reflections are in the DCM and the final fourth reflection is asymmetric at grazing exit. Analyses of the interferometry data combined with the simulations show that compared with the beam exiting the DCM a gain of 76% in the SCL was achieved, albeit with a factor of 20 reduction in flux density, which may not be a severe penalty at a synchrotron beamline. Previous efforts reported in the literature to increase the SCL that employed asymmetric crystal diffraction at grazing incidence are also discussed. A much reduced SCL is found presently in simulations wherein the same asymmetric crystal is set for grazing incidence instead of grazing exit. In addition, the present study is compared and contrasted with two other means of increasing the SCL. These are (i) focusing the beam onto an aperture to act as a secondary source, and (ii) allowing the beam to propagate in vacuum an additional distance along the beamline.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Enhanced Beam Diagnostics with Existing BPPMs via GPU-powered Multi-Particle Simulation

This research aims to utilize the multi-particle code, High-Performance Simulator (HPSim), to realistically model the Side-Coupled-Cavity Linac (CCL) lattice of the LANSCE accelerator. This new model would allow us to predict the beam’s bunch length (the longitudinal spread), which is unavailable for individual accelerating modules or only accessible at the end of the linac. However, a correct bunch length is critical for the high-energy beam transport after the CCL. Its impact would be most significant in the Proton Storage Ring (PSR), where we should be able to reduce losses for the circulating beam. The PSR is scheduled to have a 25% current increase for the neutron spallation target upgrade at the Lujan Center. A highly bunched beam would be necessary to reduce the particle losses and lower the radiation levels produced from the ring. A realistic HPSim model with >1M macro-particles can help tackle the beam losses at the sub-percent level. This new work would also create a realistic surrogate model for future machine learning projects.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Plasma expansion and relativistic filamentation in intense laser-irradiated cone targets

We report that compound parabolic concentrator (CPC) cone targets have been shown to produce increased MeV photons on the NIF-ARC by 10× over flat targets. Multiple x-ray frames can potentially be generated by firing the NIF-ARC's beamlets into distinct cone targets at few nanosecond relative delays. This requires that the cone targets with delayed beams are not degraded by their proximity to previous targets. One concern is that the spatial wings of a beam fired into one target can fall on neighboring targets, producing a preformed plasma that may interfere with laser light reaching the tip of the cone. In this work, 3D hydra simulations of realistic targets and beam parameters show that hundreds of micrometer scale length preplasmas are produced in cones within 1 mm of the laser spot. 2D particle-in-cell simulations of the intense main pulse in this preplasma indicate a density threshold for the onset of relativistic filamentation in our conditions. Applying our modeling approach to a NIF-ARC shot with an intentional 15 J prepulse yields good agreement with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Lost Muons Analysis

The Muon g-2 experiment at FermiLab represents a significant scientific endeavor aimed at measuring the anomalous magnetic moment of the muon, $a_µ$, with a precision of 140 parts per billion (ppb). During my internship at FermiLab, I focused on the simulation and analysis of Lost Muons, a critical aspect to ensure the accuracy of the experimental results. I conducted a detailed study of the beam dynamics using the existing simulation, I generated a new simulation sample enrich of events with lost muons, I analyzed it and compared the results I obtained with those from the Run-2 dataset acquired by the experiment. In this report I will briefly describe my work, discuss the results and suggest future improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Latest Updates on JSPEC - An IBS and Electron Cooling Simulation Program

JSPEC (JLab Simulation Package on Electron Cooling) is an open-source C++ program developed at Jefferson Lab, which simulates the evolution of the ion beam under the influence of both IBS and electron cooling effects. In this paper, we will report the latest updates to JSPEC. Firstly, we have added theoretical and numerical models that simulate the effect of the electron beam dispersion on non-magnetized cooling. Secondly, the cooler can now be treated as an element with length, rather than a thin lens. This change will impact the modelling of the ions and the electrons in the cooling rate calculation for both magnetized and non-magnetized cooling. Numeri-cal results will be provided to demonstrate the perfor-mance of the new models.

Zhang, H.↗
DOE OSTI2023

Dose-Weighted Radiographic Source-Spot Size

The spatial resolution of a flash radiograph is determined by the time-integrated size of the radiation source. Since the radiation pulse includes the beam head and tail, which may be substantially larger than the beam body, the radiation source is larger than the focused flattop of the electron beam pulse. Moreover, the radiation dose rate is weighted by both beam current and electron energy, and it is the radiation dose that forms the radiograph image. Therefore, in order to estimate radiographic resolution from computer simulations of a focused electron beam, one should use a dose weighted moment analysis of the radiation source distribution. This approach is used to show that, in the absence of beam effects other than space-charge and emittance, radiographic resolution for Scorpius exceeds requirements, regardless of reasonable timing options for current and gap-voltage pulses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Inversion of dynamical Bragg intensities to complex structure factors by iterated projections. For Ultramic. 2020. ("Pico" Festschrift, May 2021)

We discuss a method for recovering complex structure factors from many simultaneously excited Bragg beam intensities is described. The method is applied to simulated transmission electron diffraction data over a wide range of crystal thickness and beam energies. The method is based on iterated projections between structure and scattering matrices, which are related by a matrix unitary transformation, exponential, which we invert. The algorithm removes multiple-scattering perturbations from diffraction data and might be extended to other fields, including X-ray and neutron diffraction and cryo-electron microscopy. Because coherent multiple scattering involves interference between Bragg beams, the method also solves the phase problem. Unlike dynamical inversion from electron microscope images or ptychography data, the method, which starts with Bragg beam intensities, provides complex structure factors unaffected by focusing errors or resolution limitations imposed by lenses. We provide inversions from simulated data with 441 simultaneously excited Bragg beams over a range of thickness and beam energy. We discuss the retrieval of chirality information from enantiomorphs, the efficient incorporation of symmetry information using the irreducible representation of the group of structure matrices, and the effect of HOLZ lines to provide three-dimensional information.

97 MATHEMATICS AND COMPUTING↗