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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Galactic Cosmic Ray Simulator at the NASA Space Radiation Laboratory

The external Galactic Cosmic Ray (GCR) spectrum is significantly modified when it passes through spacecraft shielding and astronauts. One approach for simulating the GCR space radiation environment is to attempt to reproduce the unmodified, external GCR spectrum at a ground based accelerator. A possibly better approach would use the modified, shielded tissue spectrum, to select accelerator beams impinging on biological targets. NASA plans for implementation of a GCR simulator at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory will be discussed.

Norbury, John W.↗
DOE OSTI2022

Improved quantitation of SIMS depth profile measurements of niobium via sample holder design improvements and characterization of grain orientation effects

The importance of SIMS analyses for “N-doped” impurity alloyed niobium and other surface-alloyed materials continues to increase. A major hurdle is the uncertainty of instrument calibration due to changes in sample height either from sample surface topography or from the sample holder itself. The CAMECA sample holder design allows for many types of samples to be analyzed. However, a drawback is that the holder faceplate can bend, contributing an uncertainty in the relative sensitivity factor (RSF) used to quantify the SIMS results. Here, we describe an improved sample holder having a reinforced faceplate, which prevents deflection and reduces uncertainty. Simulations show that the new design significantly reduces deflection from 10 μm to 5 nm. Sample measurements show a reduction of RSF uncertainty from this source from 4.1% to 0.95%. Grain orientation has long been suspected to affect RSF measurement as well. A bicrystal implant standard, consisting of randomly oriented and [001] grains, was successively rotated 15° between analyses. It was observed that 20% of the analyses performed on the randomly oriented grain exhibited anomalously high RSF values as well as slow sputter rates. These features were associated with the changing grain normal orientation with respect to the primary Cs + beam. The grain orientation associated with the rise in RSF was simulated and determined to be the [101] crystallographic plane, thus indicating that ion channeling was responsible for the significantly increased RSF. Focused ion beam analysis confirmed slower sputter rates for the cardinal crystallographic orientations, indicating that ion channeling occurred for each.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Simulating Magnetized Electron Cooling for EIC with JSPEC

We present a possible electron cooling configuration for the proposed Electron Ion Collider (EIC) facility, developed using a Nelder-Mead Simplex optimization procedure built into JSPEC, an electron cooling code developed at Jefferson Lab. We show the time evolution of the emittance of the ion beam in the presence of this cooler evaluated assuming the ion distribution remains Gaussian. We also show that bi-gaussian distributions emerge in simulations of ion macro-particles. We show how intra-beam scattering can be treated with a core-tail model in simulations of ion macro-particles. The Sirepo/JSPEC* and Sirepo/Jupyter** apps will be presented, with instructions enabling the community to reproduce our simulations.

Coleman, S. J.↗
NASA NTRS2008

An Approach to Assess Delamination Propagation Simulation Capabilities in Commercial Finite Element Codes

An approach for assessing the delamination propagation simulation capabilities in commercial finite element codes is presented and demonstrated. For this investigation, the Double Cantilever Beam (DCB) specimen and the Single Leg Bending (SLB) specimen were chosen for full three-dimensional finite element simulations. First, benchmark results were created for both specimens. Second, starting from an initially straight front, the delamination was allowed to propagate. The load-displacement relationship and the total strain energy obtained from the propagation analysis results and the benchmark results were compared and good agreements could be achieved by selecting the appropriate input parameters. Selecting the appropriate input parameters, however, was not straightforward and often required an iterative procedure. Qualitatively, the delamination front computed for the DCB specimen did not take the shape of a curved front as expected. However, the analysis of the SLB specimen yielded a curved front as was expected from the distribution of the energy release rate and the failure index across the width of the specimen. Overall, the results are encouraging but further assessment on a structural level is required.

Krueger, Ronald↗
DOE OSTI2021

Numerical simulation of NSLS-II fast orbit feedback system

We discuss recent work on the construction of a numerical simulation of the fast orbit feedback system at NSLS-II. The simulation operates in the time domain and also includes the spatial domain, i.e. all beam position monitors (BPMs) and correctors. It can accept inputs of the real beam orbit and the measured orbit response matrix. It can also add errors to every stage of the calculation. We present methods to verify the simulation results by comparing the simulation results with experimental data collected at NSLS-II. The results are in very good agreement. The effects of errors from BPMs and correctors, as well as the amplitude of excitation, on the feedback performance, are also explored. This simulation can be used to predict the behavior and performance of the system for a future upgrade.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Improved Treatment of Host-galaxy Correlations in Cosmological Analyses with Type Ia Supernovae

Improving the use of Type Ia supernovae (SNe Ia) as standard candles requires a better approach to incorporate the relationship between SNe Ia and the properties of their host galaxies. Using a spectroscopically confirmed sample of ~1600 SNe Ia, we develop the first empirical model of underlying populations for SNe Ia light-curve properties that includes their dependence on host-galaxy stellar mass; we find a significant correlation between stretch population and stellar mass (99.9% confidence) and a weaker correlation between color and stellar mass (90% confidence). These populations are important inputs to simulations that are used to model selection effects and correct distance biases within the BEAMS with Bias Correction (BBC) framework. Here we improve BBC to also account for SNe Ia-host correlations, and we validate this technique on simulated data samples. Here, we recover the input relationship between SNe Ia luminosity and host-galaxy stellar mass (the mass step, γ) with a bias of 0.004 ±0.001 mag, which is a factor of 5 improvement over previous methods that have a γ bias of ~0.02 ± 0.001 mag. We adapt BBC for a novel dust-based model of intrinsic brightness variations, which results in a greatly reduced mass step for data (γ = 0.017 ± 0.008) and for simulations (γ = 0.006 ± 0.007). Analyzing simulated SNe Ia, the biases on the dark energy equation of state, w, vary from Δw = 0.006(5) to 0.010(5) with our new BBC method; these biases are significantly smaller than the 0.02(5) w bias using previous BBC methods that ignore SNe Ia-host correlations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Focussing Protons from a Kilojoule Laser for Intense Beam Heating using Proximal Target Structures

Proton beams driven by chirped pulse amplified lasers have multi-picosecond duration and can isochorically and volumetrically heat material samples, potentially providing an approach for creating samples of warm dense matter with conditions not present on Earth. Envisioned on a larger scale, they could heat fusion fuel to achieve ignition. We have shown in an experiment that a kilojoule-class, multi-picosecond short pulse laser is particularly effective for heating materials. The proton beam can be focussed via target design to achieve exceptionally high flux, important for the applications mentioned. The laser irradiated spherically curved diamond-like-carbon targets with intensity 4 × 10 18 W/cm 2 , producing proton beams with 3 MeV slope temperature. A Cu witness foil was positioned behind the curved target, and the gap between was either empty or spanned with a structure. With a structured target, the total emission of Cu Kα fluorescence was increased 18 fold and the emission profile was consistent with a tightly focussed beam. Transverse proton radiography probed the target with ps order temporal and 10 μm spatial resolution, revealing the fast-acting focussing electric field. Complementary particle-in-cell simulations show how the structures funnel protons to the tight focus. Here, the beam of protons and neutralizing electrons induce the bright Kα emission observed and heat the Cu to 100 eV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Real-Time Inference For MI/RR Deblending

The Fermilab Main Injector (MI) and Recycler Ring (RR) share a common beam loss monitor (BLM) system, making loss events difficult to attribute to their source machine when beam is present in both simultaneously. The Real-time Edge AI for Distributed Systems (READS) project addresses this by deblending BLM readings in real time using machine learning (ML). The current FPGA based implementation meets the sub-3 ms latency requirement but carries a resource intensive hls4ml development cycle, motivating exploration of GPU based deployment. This paper characterizes inference latency on an NVIDIA Jetson Orin Nano and introduces a packet organization scheme for assembling synchronized event frames from seven distributed BLM DAQ streams. Using a Python based DAQ simulation with injected timing jitter in place of unavailable live beam data, the pipeline achieved an average end to end latency of 0.456 ms (σ = 0.122 ms) across 167,000 test frames, comfortably meeting the timing constraint. Early outliers were attributed to TensorRT warm-up rather than steady state limitations, suggesting GPU based inference is a viable alternative to the existing FPGA implementation.

Yu, Kellen [Cornell U.]↗
DOE OSTI2026

Exploration of Real Time Inference for MI-RR Deblending on GPU/TPU Systems

The Fermilab Main Injector (MI) and Recycler Ring (RR) share a common beam loss monitor (BLM) system, making loss events difficult to attribute to their source machine when beam is present in both simultaneously. The Real-time Edge AI for Distributed Systems (READS) project addresses this by deblending BLM readings in real time using machine learning (ML). The current FPGA based implementation meets the sub-3 ms latency requirement but carries a resource intensive hls4ml development cycle, motivating exploration of GPU based deployment. This paper characterizes inference latency on an NVIDIA Jetson Orin Nano and introduces a packet organization scheme for assembling synchronized event frames from seven distributed BLM DAQ streams. Using a Python based DAQ simulation with injected timing jitter in place of unavailable live beam data, the pipeline achieved an average end to end latency of 0.456 ms (σ = 0.122 ms) across 167,000 test frames, comfortably meeting the timing constraint. Early outliers were attributed to TensorRT warm-up rather than steady state limitations, suggesting GPU based inference is a viable alternative to the existing FPGA implementation.

Yu, Kellen [Fermilab; Cornell U.]↗
DOE OSTI2023

ITER Toroidal Interferometer and Polarimeter (TIP) beam refraction in 3D density profiles

Calculations of the expected refraction of the five ITER Toroidal Interferometer and Polarimeter (TIP) chords are presented for a range of conditions including high-density axisymmetric, ELMing H-mode cases and shattered pellet disruption mitigation cases. The calculations are carried out with a newly developed ray tracing code capable of following TIP’s 10.59 μm laser beams through arbitrary 3D electron density profiles. Using JOREK simulations of ELM density perturbations in a 15MA ITER baseline plasma, it is shown that refraction from ELMs is expected to be negligible. It is found, however, that TIP interferometers will be able to clearly resolve the line-integrated density perturbation from ELMs and contribute to the ITER measurement requirement “14. H-mode, ELMs and L-H mode transition indicator”. Calculations of the expected refraction in a NIMROD simulated shattered pellet disruption mitigation scenario with peak local electron densities of n e = 3.7 x 10 21 m -3 (max line-averaged densities of n avg = 1.3 x 10 21 m -3 ) also show tolerable refraction and it is likely that most, if not all chords, would avoid signal loss. In conclusion, for both axisymmetric and structured plasmas with line-averaged densities in the mid to upper 10 21 m -3 range, values likely present only during disruption mitigation, refraction becomes significant and could limit the ability of TIP to make reliable density measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Commissioning a time-gated camera for fast neutron beamline spatial-energy characterization at LANSCE-WNR spallation source

An energy-resolved fast neutron beam imaging diagnostic has been successfully commissioned at the Weapons Neutron Research (WNR) spallation source within the Los Alamos Neutron Science Center (LANSCE) facility. This diagnostic replaces the existing analog phosphor image plates, which integrate across all neutron energies, as well as other particles, with a near-real-time energy-sensitive imaging capability. The system uses a fast plastic scintillator coupled with an intensified CCD camera. Specifically, the Teledyne Pi-MAX4 camera is coupled with either a 4 mm thick Eljen (EJ) 204 or 228 plastic scintillator. These scintillators are most sensitive to the fast neutrons (0.8-800 MeV) directly from the spallation source rather than low energy background radiation. Experimentally, these plastic scintillators were shown to have sufficiently fast decay to differentiate the bright gamma flash from the spallation neutrons. The spatial resolution is dominated by neutron beam divergence, with minimal additional contributions from scatter and light divergence. The system successfully resolved changes in neutron beam characteristics caused by intentional proton steering variations. Additionally, simulations of scintillator light yield as a function of thickness conducted using PHITS (with Scinful-QMD package) found that increasing scintillator thickness from 4 mm to 6 or 8 mm could potentially increase brightness ~ 3x. This may be explored if there is a need to reduce image acquisition time from several minutes to under one minute.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Directed Gas-Phase Formation of Azulene (C 10 H 8 ): Unraveling the Bottom-Up Chemistry of Saddle-Shaped Aromatics

The azulene (C 10 H 8 ) molecule, the simplest polycyclic aromatic hydrocarbon (PAH) carrying a fused seven- and five-membered ring, is regarded as a fundamental molecular building block of saddle-shaped carbonaceous nanostructures such as curved nanographenes in the interstellar medium. However, an understanding of the underlying gas-phase formation mechanisms of this nonbenzenoid 10π-Hückel aromatic molecule under low-temperature conditions is in its infancy. Here, by merging crossed molecular beam experiments with electronic structure calculations and molecular dynamics simulations, our investigations unravel an unconventional low-temperature, barrierless route to azulene via the reaction of the simplest organic radical, methylidyne (CH), with indene (C 9 H 8 ) through ring expansion. This reaction might represent the initial step toward to the formation of saddle-shaped PAHs with seven-membered ring moieties in hydrocarbon-rich cold molecular clouds such as the Taurus Molecular Cloud-1 (TMC-1). These findings challenge conventional wisdom that molecular mass growth processes to nonplanar PAHs, especially those containing seven-membered rings, operate only at elevated pressure and high-temperature conditions, thus affording a versatile low-temperature route to contorted aromatics in our galaxy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Ultra-short pulse laser acceleration of protons to 80 MeV from cryogenic hydrogen jets tailored to near-critical density

Laser plasma-based particle accelerators attract great interest in fields where conventional accelerators reach limits based on size, cost or beam parameters. Despite the fact that particle in cell simulations have predicted several advantageous ion acceleration schemes, laser accelerators have not yet reached their full potential in producing simultaneous high-radiation doses at high particle energies. The most stringent limitation is the lack of a suitable high-repetition rate target that also provides a high degree of control of the plasma conditions required to access these advanced regimes. Here, we demonstrate that the interaction of petawatt-class laser pulses with a pre-formed micrometer-sized cryogenic hydrogen jet plasma overcomes these limitations enabling tailored density scans from the solid to the underdense regime. Our proof-of-concept experiment demonstrates that the near-critical plasma density profile produces proton energies of up to 80 MeV. Based on hydrodynamic and three-dimensional particle in cell simulations, transition between different acceleration schemes are shown, suggesting enhanced proton acceleration at the relativistic transparency front for the optimal case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Measuring characteristic differences between high- and low-performing discharges on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF

A dense plasma focus (DPF) is a compact coaxial plasma gun, which completes its discharge as a Z-pinch, producing short (<100 ns) pulses of ions, x rays, and/or neutrons. Lawrence Livermore National Laboratory recently constructed and began operating a new device, the MJOLNIR (MegaJOuLe Neutron Imaging Radiography) DPF, which is designed for single-pulse flash neutron radiography. This device has achieved neutron yields of up to 4.1×10 11 neutrons/pulse at 3.3 MA peak current, and higher-current commissioning is under way. Like most DPFs, MJOLNIR exhibits variable yields in some configurations. Here, we present evidence of the role of parasitic current paths within the gun in stochastically influencing the yield. First through “conditioning shots,” where new hardware has been introduced, we show that increased run-down and run-in speeds correlate with higher yields. These observations are consistent with current being delivered to the electrodes but not to the main plasma sheath, degrading the implosion-driving force. Once nominal conditions are established, we correlate low-performing discharges with smaller current dip and associated voltage spike for a fixed machine configuration. A snow-plow model is able to recreate small-magnitude current dips through the introduction of a parasitic current path, and particle-in-cell simulations establish how parasitic current paths lower the ion beam energy available to produce neutrons. Finally, we observe an increased likelihood of shots with low yield and smaller current dip with increasing fill pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Macroscopic approach to the radar echo scatter from high-energy particle cascades

To probe the cosmic particle flux at the highest energies, large volumes of dense material like ice have to be monitored. This can be achieved by exploiting the radio signal. In this work, we provide a macroscopic model to predict the radar echo signatures found when a radio signal is reflected from a cosmic-ray or neutrino-induced particle cascade propagating in a dense medium like ice. Its macroscopic nature allows for an energy independent run-time, taking less than 10 s for simulating a single scatter event. As a first application, we discuss basic signal properties and simulate the expected signal for the T-576 beam-test experiment at the Stanford Linear Accelerator Center. Finally, we find good signal strength agreement with the only observed radar echo from a high-energy particle cascade to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Precision Measurement of Neutrino Oscillation Parameters with 10 Years of Data from the NOvA Experiment

This Letter reports measurements of muon-neutrino disappearance and electron-neutrino appearance and the corresponding antineutrino processes between the two NOvA detectors in the NuMI neutrino beam. These measurements use a dataset with double the neutrino mode beam exposure that was previously analyzed, along with improved simulation and analysis techniques. A joint fit to these samples in the three-flavor paradigm results in the most precise single-experiment constraint on the atmospheric neutrino mass-splitting, $Δm^2_{32}= 2.431^{+0.036}_{-0.034} (-2.479^{+0.036}_{-0.036}) \times 10^{-3}$~eV$^2$ if the mass ordering is Normal (Inverted). In both orderings, a region close to maximal mixing with $\sin^2θ_{23}=0.55_{+0.06}^{-0.02}$ is preferred. The NOvA data show a mild preference for the Normal mass ordering with a Bayes factor of 2.4 (corresponding to 70% of the posterior probability), indicating that the Normal ordering is 2.4 times more probable than the Inverted ordering. When incorporating a 2D $Δm^2_{32}\textrm{--}\sin^2 2θ_{13}$ constraint based on Daya Bay data, this preference strengthens to a Bayes factor of 6.6 (87\%).

FOS: Physical sciences↗
DOE OSTI2021

ACCELERATOR CHALLENGES OF THE LHeC PROJECT

The LHeC project studies the design of a future electron-proton collider at CERN that will run in parallel to the standard LHC operation. For this purpose, the existing LHC storage ring will be combined with an Energy Recovery Linac (ERL), to accelerate electrons up to kinetic energy of 50 GeV. This concept - also applicable to the FCC-eh collider and studied at the PERLE project as prototype version - allows a peak luminosity of 10³⁴ cm⁻² s⁻¹. A sophisticated design of the RF structures, linacs, arcs, and interaction region is required. The electrons are accelerated and, after the interaction point, their energy is recovered through the same RF structures. While this energy recovery concept is a very promising approach, severe challenges are set by the layout of the interaction region, the beam separation concept and the design of the linac and arc lattice for the highest possible momentum acceptance. Emittance control and beam-beam effect of both, electron and proton beams, have been studied in front-to-end simulations and will be presented. We summarise the design principles of the ERL, the optimization of the arc lattice, and the main parameters of the project.

Holzer, B. J.↗