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,369 records · Page 76

3D simulations of inertial confinement fusion implosions part 2: systematic flow anomalies and impact of low modes on performances in OMEGA experiments

We present the first 3D radiation-hydrodynamic simulations of directly driven inertial confinement fusion implosions with an inline package for polarized crossed-beam energy transfer, which were used to assess the impact of the current distributed polarization rotators (DPRs) on OMEGA as well as other known sources of asymmetry. Applied to OMEGA implosions, the simulations predict bang times with no need for ad hoc multipliers, as well as yields—if you separately account for the impacts of imprint and fuel age. The magnitude of the flow is well reproduced when the low mode sources are large, whereas the modeling of the stalk is thought to be required to match the flow magnitude in the remaining cases. For the cases explored in more detail, polarized cross-beam energy transfer (CBET)—the only known systematic drive asymmetry, brought the results closest to the measured flow vectors. The remaining discrepancies are shown to be stemming from the limited knowledge of the laser pointing modes. For typical current levels of beam mispointing, power imbalance, target offset, and asymmetry caused by polarized CBET, low modes degrade the yield by more than 40%. The current strategy of attempting to compensate the mode-1 asymmetry with a preimposed target offset recovers only about one-third of the losses caused by the low modes due to the dynamic nature of the multiple asymmetries and the presence of low modes other than l = 1. Furthermore, addressing the root causes of the drive asymmetries is apt to be more beneficial. To that end, one possible solution to the specific issue of polarized CBET (10 µm DPRs) is shown to work well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculating the space-charge-limited current density for nonplanar geometries by simulating the charge-free electric field

Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.

Wright, Jack K. [Purdue Univ., West Lafayette, IN ↗

Applying Linear Optics from Closed Orbits Modulation for finding beam-based alignment of harmonic sextupoles

A fast and accurate beam-based alignment (BBA) method for harmonic sextupoles has been devel oped at NSLS-II using Linear Optics from Closed Orbit Modulation (LOCOM). The approach excites the beam with simultaneous sine-wave signals at two fast correctors, chosen with an appropriate phase advance to span the full betatron phase space. This strategy suppresses systematic errors from hysteresis, while additional errors from orbit drift and power-supply calibration are minimized by the short measurement time (a few minutes) and reliance solely on beam-based current-to-field conver sion of the sextupoles, with hysteresis explicitly included. Simulations indicate that Linear Optics from Closed Orbits (LOCO) combined with 0.5 mm local orbit bumps can resolve relative sextupole field offsets (∆k₂) with precision better than 10% of k₂, reflecting to beam-based alignment (BBA) accuracy finer than 50 µm. Moreover, employing machine-learning-optimized local orbit bumps en hances sextupole-induced quadrupole signals in a deterministic manner and maintains orbit stability under large sextupole strength variations (±40%). The proposed method is experimentally validated through proof-of-principle measurements at NSLS-II, demonstrating its potential as a fast, precise, and robust tool for harmonic sextupole alignment.

43 PARTICLE ACCELERATORS↗

Injection of electron beams into two laser wakefields and generation of electron rings

Throughout this work, mutual injection of electron beams into two laser plasma wakefields was observed experimentally when driving laser pulses interfered in plasma at a small crossing angle and were slightly relatively delayed, approximately by one pulse duration. Particle-in-cell simulations revealed that the mutual injection was sensitive to the spatial overlap of the laser pulses, which therefore could be used to control the mutual injection. The dual synchronized, femtosecond electron beams are potentially useful for pump-probe experiments in ultrafast science. In addition, out-of-axis ring-shaped electron beams were detected in both experiments and simulations.

47 OTHER INSTRUMENTATION↗

Status of the Measurement of Proton Scattering on Carbon Nuclei in EMPHATIC for Neutrino Flux Uncertainty Reduction

In long-baseline neutrino oscillation experiments, Monte Carlo (MC) simulations based on hadron interactions and decays are used to predict the neutrino flux. The 10%-level systematic uncertainty of the predicted neutrino fluxes from these simulations is dominated by uncertainties in hadron interaction cross sections due to limited hadron scattering data. EMPHATIC aims to reduce the neutrino flux uncertainty by providing additional data. Using a table-top-sized spectrometer located at the Fermilab Test Beam Facility (FTBF), its physics program includes precise measurements of hadron scattering and production cross sections at various beam momenta and target species that are relevant for GeV-scale neutrino production. Using simulation, we have developed a simple single-track reconstruction algorithm that has a momentum resolution of 3-4\%. We will demonstrate the progress in developing one of EMPHATIC’s first track reconstruction algorithms – an important step in making a new single-track forward scattering measurement (p + C $\rightarrow$ p + C at several beam momenta) using Phase 1 data collected between 2022 and 2023.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The scaling of relativistic double-year widths - Poisson-Vlasov solutions and particle-in-cell simulations

The study of relativistic plasma double layers is described through the solution of the one-dimensional, unmagnetized, steady-state Poisson-Vlasov equations and by means of one-dimensional, unmagnetized, particle-in-cell simulations. The thickness vs potential-drop scaling law is extended to relativistic potential drops and relativistic plasma temperatures. The transition in the scaling law for 'strong' double layers suggested by analytical two-beam models by Carlqvist (1982) is confirmed, and causality problems of standard double-layer simulation techniques applied to relativistic plasma systems are discussed.

Sulkanen, Martin E.↗

3D Simulations of Ion Thruster Accelerator Grid Erosion Accounting for Charge Exchange Ion Space Charge

Accelerator (accel) grid sputtering by ions formed through charge-exchange (CEX) reactions between beam ions and residual neutral gas is a critical life-limiting mechanism for gridded ion thrusters. The three-dimensional ion optics code CEX3D is designed to simulate this grid erosion for a single beamlet, with a particular emphasis on non-axisymmetric features such as the "pits and grooves" erosion commonly observed on the accel grid downstream face in two-grid thrusters. The treatment of CEX ions in the code was recently upgraded with a new particle-in-cell (PIC) module to account for the influence of these ions' space charge on the electrostatic potential downstream of the grids. In order to achieve reasonable computation times while resolving the Debye length near the grids and avoiding gross violations of the Courant-Friedrichs-Lewy (CFL) condition, macroparticle velocities in the PIC calculation are limited through a rescaling procedure that preserves ion trajectories and space charge density. The code accounts for beam divergence, finite momentum transfer in CEX collisions, and radial losses of CEX ions from the beam; these effects are important for determining the CEX ion flux to the accel grid because the calculated potential downstream of the grids can become very flat. The upgraded code has been used to simulate operation of NASA's Evolutionary Xenon Thruster (NEXT) during the 51 kHr Long Duration Test - a selection of results is presented and compared with experimental data.

Chaplin, Vernon H.↗

Molecular beam facility for studying mass spectrometer performance

An apparatus which produces neutral gas beams in the velocity range from thermal to 6 km/s makes possible studies simulating the motion of instruments through tenuous atmospheres. Use of the apparatus in studying the response of an open-source mass spectrometer is examined. Experiments performed include responses when the instrument is in a normal mode of operation, when the ion source potentials are adjusted to reject striking gas particles, when the beam is restricted so that portions of the ion source are struck, and when the beam strikes at angles other than normal incidence. Results indicate specular radiation of particles out of the source and stagnation ratios close to unity for normal operation. Retarding potential studies dropped the sensitivity by a factor of about 10, and as angles of attack are varied, the effect was found to depend upon beam velocity and the way the ions are initially accelerated.

Ballenthin, J. O.↗

6D Simulations of PIP-II Booster Injection

The PIP-II superconducting linac will deliver 2 mA average H⁻ beam current at 800 MeV to the existing Booster synchrotron over a period of 0.55 ms (285 turns). As a result, the injected beam power will quadruple to 17 kW. Safe operation at the increased beam power implies careful attention to the origin, magnitude, and distribution of both controlled and uncontrolled losses. Uncontrolled losses are due to neutral ions in excited states stripped in downstream magnets and large angle scattered protons from parasitic foil hits. The relative magnitudes of these loss mechanisms is used to determine the optimal foil thickness. A transverse painting scheme involving closed orbit motion will be used to mitigate space charge effects and minimize parasitic foil hits. Using a detailed full 6D simulation of the injection process, we compute large angle scattering losses and compare results to back of the envelope estimates. We investigate possible impact of space charge on the emittance and beam distribution both during and at the conclusion of the injection period.

43 PARTICLE ACCELERATORS↗

Proton beam dynamics in bare IOTA with intense space-charge

We are commissioning a 2.5-MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MADX simulation codes. Our findings reveal agreement within a factor of 2 in core phase-space density predictions up to 100 synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.

43 PARTICLE ACCELERATORS↗

Laboratory and On-Sky Validation of the Shaped Pupil Coronagraph's Sensitivity to Low-Order Aberrations with Active Wavefront Control

We present early laboratory simulations and extensive on-sky tests validating of the performance of a shaped pupil coronagraph (SPC) behind an extreme-AO (Adaptive Optics) corrected beam of the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system. In tests with the SCExAO internal source/wavefront error simulator, the normalized intensity profile for the SPC degrades more slowly than for the Lyot coronagraph as low-order aberrations reduce the Strehl ratio (S.R.) from extremely high values (S.R. approximately 0.93-0.99) to those characteristic of current ground-based extreme AO systems (S.R. approximately 0.74-0.93) and then slightly lower values down to S.R. approximately 0.57. On-sky SCExAO data taken with the SPC and other coronagraphs for brown dwarf/planet-hosting stars HD 1160 and HR 8799 provide further evidence for the SPC's robustness to low-order aberrations. From H-band Strehl ratios of 80 percent to 70 percent, the Lyot coronagraph's performance versus that of the SPC may degrade even faster on sky than is seen in our internal source simulations. The 5-sigma contrast also degrades faster (by a factor of two) for the Lyot than the SPC. The SPC we use was designed as a technology demonstrator only, with a contrast floor, throughput, and outer working angle poorly matched for SCExAO's current AO performance and poorly tuned for imaging the HR 8799 planets. Nevertheless, we detect HR 8799 cde with SCExAO/CHARIS using the SPC in broadband mode, where the signal-to-noise ratio for planet e is within 30 percent of that obtained using the vortex coronagraph. The shaped-pupil coronagraph is a promising design demonstrated to be robust in the presence of low-order aberrations and may be well-suited for future ground and space-based direct imaging observations, especially those focused on follow-up exoplanet characterization and technology demonstration of deep contrast within well-defined regions of the image plane.

Currie, Thayne↗

Proton Beam Dynamics in Bare IOTA with Intense Space-Charge

We are commissioning a 2.5~MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MAD-X simulation codes. Our findings reveal agreement within 30\% in core phase-space density predictions up to 100~synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.

43 PARTICLE ACCELERATORS↗

Theoretical studies of a molecular beam generator

The following is a proposed baseline model that is being developed for the simulation of hydrodynamic generator, which can be converted at a later date to magnetohydrodynamic MHD thruster by adding the necessary electric and magnetic fields. The following development will include the electric and magnetic terms, however, the initial computer program will not include these terms. The analysis that follows is a one species, single temperature model constructed over a domain that is defined by an enclosed region.

Heinbockel, John H.↗

Analysis and Simulation for a Spotlight-Mode Aircraft SAR in Circular Flight Path

A wide azimuth beam SAR can offer higher resolution or wider azimuth viewing angle; two factors that help better characterize the backscattering property of targets for various science applications. One disadvantage of wide beam SAR is that a much higher pulse repetition frequency (PRF) is usually required since PRF is proportional to the radar beam angle. This problem can be resolved using a spotlight-mode concept: steering a narrow beam SAR to a fixed spot on the ground. The drawback of a spotlight-mode SAR is its limited coverage. A conventional spotlight-mode SAR operates along a straight line path as shown in Figure 1. It can be shown that spotlight-mode SAR that follows a straight line path has difficulty in achieving the ultimate resolution of lambda/4. It also cannot utilize the full 180 degree of azimuth viewing angle that can be attained only when the synthetic aperture length approaches infinity.

Jin, Michael Y.↗

Computational Modeling, Simulation, and Potential Applications of Optical Stochastic Cooling

With the rising demand for intense particles beams, much research is being conducted in the area of particle beam cooling. One of these methods, called Stochastic Cooling (SC) (developed at CERN in the 1970's), delivered a feedback method to improve the quality and lifetime of circulating proton beams by reducing their 6D phase-space and has been widely implemented in a number of hadron machines. However, traditional stochastic cooling schemes are limited by the bandwidth of microwave frequency systems. Optical Stochastic Cooling (OSC) is a promising extension of the stochastic cooling beam cooling technique. OSC instead uses optical wavelengths which allows for improved control and increased cooling but creates its own technical challenges. This lays out work conducted toward the experimental demonstration of OSC at Fermilab's Integrable Optics Test Accelerator (IOTA) storage ring. This includes the design and characterization of parts of the optical delay system, the development and validation of a high-fidelity computational model of the OSC process, and the investigations into possible applications of the OSC mechanism to advanced beam manipulations.

43 PARTICLE ACCELERATORS↗

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗

Impact of flow-induced beam deflection on beam propagation in ignition scale hohlraums

Experiments examining the amount of specular reflection (or “glint”) within hohlraums containing different gas fill densities have recently been performed. Simulations of these experiments are presented that show using a single flux limiter cannot explain the decrease in glinted power with increasing gas fill density. The hypothesis that flow-induced beam deflection alters laser absorption is presented. A model is proposed that can be implemented into a ray tracing description of the laser commonly used in radiation hydrodynamic codes. It is shown that simulations using this model capture the trend with gas fill density improving agreement with measurements. This formulation is then applied to an ensemble of laser-driven inertial confinement fusion experiments performed at the National Ignition Facility. The proposed model shows little impact on the total x-ray drive on the capsule but a large impact on the resulting implosion symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗