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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Excitation of long electrostatic solitary waves in ion beam neutralization process

Unusually long electrostatic solitary waves (ESWs) are discovered in a particle-in-cell simulation study of the process of ion beam neutralization by electron emission from a filament. These ESWs are long because the density perturbation responses to the potential wells created by the ESWs are very small. Notably, the density perturbation is small because the trapped (positive) and untrapped (negative) electron density perturbations nearly compensate each other because of a non-Maxwellian electron velocity distribution in the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Sensitivity study with a $\textit{D}$ and $\textit{B}$ mesons modular simulation code of heavy flavor $R_{AA}$ and azimuthal anisotropies based on beam energy, initial conditions, hadronization, and suppression mechanisms

Heavy flavor probes provide important information about the in-medium properties of the quark-gluon plasma produced in heavy-ion collisions. In this work, we investigate the effects of 2D + 1 event-by-event fluctuating hydrodynamic backgrounds on the nuclear suppression factor and momentum anisotropies of heavy flavor mesons and nonphotonic electrons. Using the state-of-the-art $\textit{D}$ and $\textit{B}$ mesons modular simulation code (called “DAB-MOD”), we perform a systematic comparison of different transport equations in the same background, including a few energy-loss models—with and without energy-loss fluctuations—and a relativistic Langevin model with different drag parametrizations. We present the resulting $\textit{D}$ and $\textit{B}$ mesons $R_{A A}, v_2, v_3,$ and $v_4$ as well as multiparticle cumulants, in AuAu collisions at $\sqrt{s_{NN}}$ = 200 GeV and PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV and $\sqrt{s_{NN}}$ = 5.02 TeV , and compare them to the available experimental data. The $v_2${4}/$v_2${2} ratio, which is known to be a powerful probe of the initial conditions and flow fluctuations in the soft sector, is also studied in the context of heavy flavor. We also investigate the correlations between the transverse anisotropies of heavy mesons and all charged particles to better understand how heavy quarks couple to the hydrodynamically expanding quark-gluon plasma. Finally, we study the influence that different initial conditions and the implementation of heavy-light quark coalescence has on our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Mutual optical intensity propagation through non-ideal two-dimensional mirrors

The mutual optical intensity (MOI) model is a partially coherent radiation propagation tool that can sequentially simulate beamline optics and provide beam intensity, local degree of coherence and phase distribution at any location along a beamline. This paper extends the MOI model to non-ideal twodimensional (2D) optical systems, such as ellipsoidal and toroidal mirrors with 2D figure errors. Simulation results show that one can tune the trade-off between calculation efficiency and accuracy by varying the number of wavefront elements. The focal spot size of an ellipsoidal mirror calculated with 100 $\times$ 100 elements gives less than 0.4% deviation from that with 250 $\times$ 250 elements, and the computation speed is nearly two orders of magnitude faster. Effects of figure errors on 2D focusing are also demonstrated for a non-ideal ellipsoidal mirror and by comparing the toroidal and ellipsoidal mirrors. Finally, the MOI model is benchmarked against the multi-electron Synchrotron Radiation Workshop (SRW) code showing the model’s high accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Self-Consistent Relativistic Electron Scattering using the Sherlock Scattering Model for X-ray Diagnostics

We present on a new, self-consistent, arbitrary-temperature Romberg integration scheme for modeling electron scattering in materials in a LANL Lagrangian Shock Hydro (LSH) code. Electron beam-target interactions are fundamental to a wide range of scientific and technological applications. When high-energy electron beams hit their target, they may scatter, deposit energy, or ionize the source. These processes govern the behavior and outcomes in nanotechnology manufacturing, electron microscopy, and modern X-ray diagnostics. Simulating these interactions is essential for interpreting experimental results, predicting material responses, and designing efficient tools and experiments. At Los Alamos, this is done using a LSH code, which is a multi-dimension, multi-material, massively parallel, multi-physics code used to simulate applications from asteroid impacts to electron beam interactions. By effectively and efficiently modeling the way that electrons scatter from the beam we can bolster these simulations and more accurately predict experimental outcomes. The model currently implemented in the LSH of interest is based on work by Papp and does not self-consistently preserve momentum in the slightly relativistic regime; here we adopt a model proposed by Braams and Karney and implement a Romberg integration scheme to compute the diffusion tensor. In this paper we will provide background on the Braams-Karney diffusion tensor as well as the Romberg integration scheme we employed to numerically solve for it. We will show that our integration scheme is accurate in solving for the set of scalar potentials used to re-express the diffusion tensor in differential form, and in solving for the diffusion coefficients in the larger LSH code. By using this diffusion tensor rather than the existing Papp one, and numerically integrating it with a Romberg method, we produce much more accurate, self-consistent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Incident beamline design for a modern cold triple axis spectrometer at the High Flux Isotope Reactor

A modern cold triple axis spectrometer is being planned for the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. Here, we describe the design of an incident beamline that will put a flux of ~ 10 8 $\frac {n}{cm^{2}s}$ on a sample with an area of 2 cm X 2 cm. It takes current physical constraints at HFIR into account and it can accommodate both single and multiplexed analyzer-detector secondary spectrometers and large superconducting magnets. The proposed incident beamline includes a multi-channel guide with horizontal focusing, a neutron velocity selector, components to facilitate an incident beam polarization option, and a double-focusing pyrolytic graphite monochromator. Here, this work describes the process of optimizing the guide system and monochromator and summarizes the expected performance of the incident beamline for non-polarized operation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Divertor detachment characterization in negative triangularity discharges in DIII-D via 2D edge-plasma transport modeling

Edge fluid modeling of the first divertor-plasma detachment experiments in negative triangularity (NT) discharges in the DIII-D tokamak is presented using the 2D multifluid edge transport code UEDGE, including cross-field particle drifts. Experiments were performed where the lower single-null magnetic equilibrium had a strong NT (δ≈−0.5), that is, where the magnetic X-point is at a larger major radius than the core magnetic axis. Divertor-plasma detachment was induced by increasing the core plasma density in DIII-D via intrinsic gas puffing. Here density scans are performed with UEDGE to reach a detached plasma and to quantitatively recover the experimental rollover of the ion saturation current on the outer divertor target plate. The simulations cover experiments with both signs of the toroidal magnetic field, B T , where the ion magnetic Grad-B drifts are directed into (forward B T ) and out of (reverse B T ) the divertor region. Consistent with experiments with neutral beam power injection, the NT simulations reproduce: 40% higher density is needed to reach detachment onset with forward B T compared with reverse B T , and the absence of deep detachment is found with reverse B T . Similarly, comparison between Ohmic discharges in NT and positive triangularity (PT) shaping confirms that a substantially higher density is needed to achieve detachment in NT than in PT, with NT requiring an line-average density of at least the Greenwald density or higher. Simulation results suggest that higher densities are needed to reach detachment in negative compared to PT because these discharges have a shorter midplane-to-target distance along the total magnetic field B, a shorter outer divertor poloidal leg length (0.06 m vs 0.2 m), and reduced radial transport [near-scrape-off layer (SOL) D ⊥ /χ ⊥ =0.3/0.5 vs D ⊥ /χ ⊥ =1.0/1.0, all in (m 2 /s)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Lateral Beam Shifts and Depolarization Upon Oblique Reflection from Dielectric Mirrors

Abstract Dielectric mirrors comprising thin‐film multilayers are widely used in optical experiments because they can achieve substantially higher reflectance compared to metal mirrors. Here, potential problems are investigated that can arise when dielectric mirrors are used at oblique incidence, in particular for focused beams. It is found that light beams reflected from dielectric mirrors can experience lateral beam shifts, beam‐shape distortion, and depolarization, and these effects have a strong dependence on wavelength, incident angle, and incident polarization. Because vendors of dielectric mirrors typically do not share the particular layer structure of their products, several dielectric‐mirror stacks are designed and simulated, and then the lateral beam shift from two commercial dielectric mirrors and one coated metal mirror is also measured. This paper brings awareness of the tradeoffs between dielectric mirrors and front‐surface metal mirrors in certain optics experiments, and it is suggested that vendors of dielectric mirrors provide information about beam shifts, distortion, and depolarization when their products are used at oblique incidence.

Physics↗
DOE OSTI2025

Development of an FGPA-Based Cavity Simulator for Testing RF Controls

LLRF is used to precisely control the amplitude and phase of the RF field in cavities. Often times, access to test the control algorithms with RF equipment, especially in the presence of beam, is limited or beyond reach. In such cases, testing must be done through computer modeling or simulations. Computer modeling is often too slow and difficult to interface with the LLRF hardware. Analog or digital cavity simulators are preferred as they allow for interaction with the LLRF controls platform in real-time, and compared to their analog counterparts, FPGA-based digital cavity simulators allow for a more adjustable and sophisticated implementation. The newly developed FPGA-based cavity simulator includes the cavity electrical model, the cavity mechanical model including Lorentz Force Detuning and microphonics, an amplifier model which can simulate real amplifier nonlinearities, and a beam model. The simulator has been validated using measurements from BNL’s CeC 704 MHz 5-cell SRF cryomodule.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Understanding Type Ia Supernova Distance Biases by Simulating Spectral Variations

In the next decade, transient searches from the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope will increase the sample of known Type Ia supernovae (SNe Ia) from ~10 3 to 10 5 . With this reduction of statistical uncertainties on cosmological measurements, new methods are needed to reduce systematic uncertainties. Here, characterizing the underlying spectroscopic evolution of SN Ia remains a major systematic uncertainty in current cosmological analyses, motivating a new simulation tool for the next era of SN Ia cosmology: Build Your Own Spectral Energy Distribution (byosed). byosed is used within the SNANA framework to simulate light curves by applying spectral variations to model SEDs, enabling flexible testing of possible systematic shifts in SN Ia distance measurements. We test the framework by comparing a nominal Roman SN Ia survey simulation using a baseline SED model to simulations using SEDs perturbed with byosed, and investigating the impact of ignoring specific SED features in the analysis. These features include semiempirical models of two possible, predicted relationships: between SN ejecta velocity and light-curve observables, and a redshift-dependent relationship between SN Hubble residuals and host-galaxy mass. We analyze each byosed simulation using the SALT2 and BEAMS with Bias Corrections framework, and estimate changes in the measured value of the dark-energy equation-of-state parameter, w. We find a difference of Δw = –0.023 for SN velocity and Δw = 0.021 for redshift-evolving host mass when compared to simulations without these features. By using byosed for SN Ia cosmology simulations, future analyses (e.g., the Rubin and Roman SN Ia samples) will have greater flexibility to constrain or reduce such SN Ia modeling uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Intense gamma-ray source based on focused electron beams from a laser wakefield accelerator

Laser wakefield accelerators generate ultrashort electron bunches with the capability to produce γ-rays. Here, we produce focused laser wakefield acceleration electron beams using three quadrupole magnets. Electron beams are then focused into a 3 mm lead converter to generate intense, focused bremsstrahlung γ beams. Experimental results demonstrate the generation and propagation of focused γ beams to a best focus spot size of 2.3 ± 0.1 × 2.7 ± 0.2 mm 2 using a copper stack calorimeter. Monte Carlo simulations conducted using GEANT4 are in good agreement with experimental results and enable detailed examination of γ-ray generation. Simulations indicate that the focused γ beams contained 2.6 × 10 9 photons in the range of 100 keV to 33 MeV with an average energy of 6.4 MeV. A γ-ray intensity of 7 × 10 10 W/cm 2 was estimated from simulations. In conclusion, the generation of focused bremsstrahlung γ-ray sources can have important applications in medical imaging applications and laboratory astrophysics experiments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

The Simons Observatory: Beam Characterization for the Small Aperture Telescopes

Abstract We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a mapmaker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the mapmaker. We test our reconstruction performance for four different frequency bands against various algorithmic parameters, atmospheric conditions, and input beams. We additionally show the reconstruction quality as a function of the number of available observations and investigate how different calibration strategies affect the beam uncertainty. For all of the cases considered, we find good agreement between the fitted results and the input beam model within an ∼1.5% error for a multipole range ℓ = 30–700 and an ∼0.5% error for a multipole range ℓ = 50–200. We conclude by using a harmonic-domain component separation algorithm to verify that the beam reconstruction errors and biases observed in our analysis do not significantly bias the Simons Observatory r -measurement

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Self-driven ion deflectometry measurements using MeV fusion-driven protons and accelerated deuterons in the deuterated hybrid x-pinch on the MAIZE LTD generator

Abstract We report on the results of point-projection ion deflectometry measurements from a mid-size university z-pinch experiment. A 1 MA 8 kJ LTD generator at the University of Michigan (called MAIZE) drove a hybrid x-pinch (HXP) with a deuterated polyethylene fiber load to produce a point-like source of MeV ions for backlighting. In these experiments, 2.7 MeV protons were generated by DD beam-target fusion reactions. Due to the kinematics of beam-target fusion, the proton energies were down-shifted from the more standard 3.02 MeV proton energy that is released from the center-of-mass rest frame of a DD reaction. In addition to the 2.7 MeV protons, strongly anisotropic beams of 3 MeV accelerated deuterons were detected by ion diagnostics placed at a radial distance of 90 mm from the x-pinch. Numerical reconstruction of experimental data generated by deflected hydrogen ion trajectories evaluated the total current in the vacuum load region. Numerical ion-tracking simulations show that accelerated deuteron beams exited the ion source region at large angles with respect to the pinch current direction.

Physics↗
DOE OSTI2022

Multiphysics Simulation of the Thermal Response of a Nanofibrous Target in a High-Intensity Beam

Nanofibrous structures are of high interest to the fields of engineering and materials science, and investigation of their properties as well as discovery of novel applications for them both constitute lively areas of research. A very promising application of nanofiber mats lies in the field of accelerator technology: beam targets made from nanofiber mats offer a solution to the problem of advancing the "intensity frontier"–-the limit on the beam intensities that can be realized in fixed target experiments and neutrino production facilities. However, testing has shown that the survivability of these nanofiber targets depends strongly on their manufacturing parameters, such as the packing density of fibers. In this work, we will use multiphysics simulations to perform a thermal study on how nanofiber targets react to high intensity beams, so that the dependency of the targets’ lifetime on their construction parameters can be better understood.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Recent Beam Stability Analysis for the EIC

This document summarizes work done in WBS 6.02.02 during FY23. Additional details can be found in the EIC overview paper presented at IPAC’23 [1], and the references therein. Progress has been made in the ESR and HSR lattice design, dynamic aperture optimization, and the RCS lattice design. The impedance and Collective effects are progressing and include the impedance optimization process of the vacuum systems for RCS, HSR, and ESR; collective effects studies; collective effects and beam-beam interaction; coupled bunch instabilities and the crab cavities; low-level RF feedback system design and beam-ion instability. The reversed phasing RF system has been numerically studied for the ESR to mitigate Robinson instability and demonstrate reliable stable beam operation. Various codes, including C++, SPACE [2], ELEGANT [3], and Mbtrack2 [4] have been employed to benchmark the results. The simulation of HSR bunch splitting with beam loading has been performed at 275 GeV energy. The evaluation of Beam Position Monitors aimed to validate their expected performance, with the primary objective being to verify their accuracy. Calculations of Electron Polarization in the ESR, RCS, and HSR are showing good progress. In the RCS, preliminary studies indicate excellent polarization transmission over intrinsic spin resonances, achieving over 90% transmission with improved performance compared to the previous lattice. Preliminary simulations in the ESR indicate encouraging results in minimizing depolarization and improving equilibrium polarization.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer

An imaging neutral particle analyzer (INPA) provides energy and radially resolved measurements of the confined fast-ion population ranging from the high-field side to the edge on the midplane of the DIII-D tokamak. In recent experiments, it was used to diagnose fast-ion flow in the INPA-interrogated phase-space driven by multiple, marginally unstable Alfvén eigenmodes (AEs). Here, the key features of this measured fast-ion flow are: (I) a fast-ion flow from q min and the injection energy (81 keV) towards lower energies and plasma periphery.(II) A flow from the same location towards higher energies and the plasma core, (III) a phase-space 'hole' at the injected energy and plasma core and (IV) a pile-up at the plasma core at lower energies (~60 keV). Ad hoc energetic particle diffusivity modelling of TRANSP significantly deviates from the observation. Comparably, a reduced modelling, i.e. a combination of NOVA-K and ASCOT5 code with the measured mode structure and amplitude, generally reproduce some key features of the observed phase-space flow, but largely failed to interpret fast ion depletion near the plasma axis. At last, self-consistent, first-principle multi-phase hybrid simulations that include realistic neutral beam injection and collisions are able to reproduce most features of the time-resolved phase-space flow. During consecutive hybrid phases, an RSAE consistent with the experiment grows and saturates, redistributing the injected fast ions. The resulting synthetic INPA images are in good agreement with the measurement near the injection energy. The simulations track the fast-ion redistribution within the INPA range, confirming that the measured fast-ion flow follows streamlines defined by the intersection of phase-space surfaces of constant magnetic moment μ and constant E' = nE + ωP φ , where n and ω are the instability toroidal mode number and frequency, and E and P φ the ion energy and toroidal canonical momentum. Nonperturbative effects are required to reproduce the depletion of fast ions near the magnetic axis at the injection energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

3D simulations of inertial confinement fusion implosions part 1: inline modeling of polarized cross beam energy transfer and subsequent drive anomalies on OMEGA and NIF

Inertial confinement fusion experiments are sensitive to cross-beam energy transfer (CBET), a laser-plasma instability that redistributes laser energy in the coronal plasma through self-generated ion acoustic wave (IAW) gratings. The detailed CBET coupling depends on the polarization state of the crossing wavefields. CBET itself can also scramble the beam polarizations by inducing ellipticity through the IAW grating, and rotating the seed polarization toward that of the pump. We develop a ray-based model that describes the polarized CBET coupling and that is compatible with the framework of 3D inline radiative hydrodynamics simulations. The model is implemented in the ASTER/IFRIIT code and verified against an academic test case and an offline polarized CBET post-processor. It is then applied to the detailed configuration of the distributed polarization rotator system on OMEGA, where results highlight how polarized CBET induces significant low modes in the collisional absorption source term. Finally, the modeling is applied to a simple indirect-drive configuration, comparing CBET calculations with 96 unpolarized or polarized beams with 24 unpolarized quads. It is shown that these cases produce similar power amplification per cone of beams grouped with similar polar angles. However, the 96 beam geometry itself is found to reduce azimuthal variations in quad power after the interaction and favors beams with larger polar angles within the cones, an effect that is amplified by the polarized CBET. Application of the model to inline calculations of OMEGA implosions are presented in a companion paper.

3D inline modeling↗
DOE OSTI2021

A Chemical Dynamics Study of the Reaction of the Methylidyne Radical (CH, X 2 Π) with Dimethylacetylene (CH 3 CCCH 3 , X 1 A 1g )

The gas-phase reaction of the methylidyne (CH; X 2 Π) radical with dimethylacetylene (CH 3 CCCH 3 ; X 1 A 1g ) was studied at a collision energy of 20.6 kJ mol -1 under single collision conditions with experimental results merged with ab initio calculations of the potential energy surface (PES) and ab initio molecule dynamics (AIMD) simulations. Here, the crossed molecular beams experiment reveals that the reaction proceeds barrierless via indirect scattering dynamics through long-lived C 5 H 7 reaction intermediate(s) ultimately dissociating to C 5 H 6 isomers along with atomic hydrogen with atomic hydrogen predominantly released from the methyl groups as verified by replacing the methylidyne with the D1-methylidyne reactant. AIMD simulations reveal that the reaction dynamics are statistical leading predominantly to p28 (1-methyl-3-methylene-cyclopropene, 13 %) and p8 (1-penten-3-yne, 81 %) plus atomic hydrogen with a significant amount of available energy being channeled into the internal excitation of the polyatomic reaction products.

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

PARMELA simulation of Line-A beamline

This is a short re-work of the previously accomplished simulations, using PARMELA code package. This exercise was done primarily for increase of fidelity and re-formulating the given task for a different beam dynamics (BD) code. The performed simulations showed very good agreement with the abovementioned simulations and confirm the previously produced results with the BEAMPATH code.

43 PARTICLE ACCELERATORS↗