Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accelerator Simulation”

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,297 records · Page 72

On the high-energy gamma-ray signature of cosmic-ray sources

Monte Carlo simulations of the gamma-ray emission from hypothetical cosmic-ray sources are performed. Sources which might correspond to acceleration by supernova shocks in 'average' interstellar conditions and deep within giant molecular clouds are considered. The consequences of dropping the common assumption that the cosmic-ray spectrum at the sources is the same as that observed at earth are examined. Spectral effects which can be related to the depth of the material shroud and the population of accelerated particles are explored using these simulations and are described. The results are compared with the COS B catalog of gamma-ray sources, and the implications for the underlying particle populations and source mechanisms are discussed.

Ormes, J. F.↗

The mechanisms of electron heating and acceleration during magnetic reconnection

The heating of electrons in collisionless magnetic reconnection is explored in particle-in-cell simulations with non-zero guide fields so that electrons remain magnetized. In this regime, electric fields parallel to B accelerate particles directly, while those perpendicular to B do so through gradient-B and curvature drifts. The curvature drift drives parallel heating through Fermi reflection, while the gradient B drift changes the perpendicular energy through betatron acceleration. We present simulations in which we evaluate each of these mechanisms in space and time in order to quantify their role in electron heating. For a case with a small guide field (20% of the magnitude of the reconnecting component), the curvature drift is the dominant source of electron heating. However, for a larger guide field (equal to the magnitude of the reconnecting component) electron acceleration by the curvature drift is comparable to that of the parallel electric field. In both cases, the heating by the gradient B drift is negligible in magnitude. It produces net cooling because the conservation of the magnetic moment and the drop of B during reconnection produce a decrease in the perpendicular electron energy. Heating by the curvature drift dominates in the outflow exhausts where bent field lines expand to relax their tension and is therefore distributed over a large area. In contrast, the parallel electric field is localized near X-lines. This suggests that acceleration by parallel electric fields may play a smaller role in large systems where the X-line occupies a vanishing fraction of the system. The curvature drift and the parallel electric field dominate the dynamics and drive parallel heating. A consequence is that the electron energy spectrum becomes extremely anisotropic at late time, which has important implications for quantifying the limits of electron acceleration due to synchrotron emission. An upper limit on electron energy gain that is substantially higher than earlier estimates is obtained by balancing reconnection drive with radiative loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The International Linear Collider: A Global Project

The International Linear Collider (ILC) is now under consideration as the next global project in particle physics. In this report, we review of all aspects of the ILC program: the physics motivation, the accelerator design, the run plan, the proposed detectors, the experimental measurements on the Higgs boson, the top quark, the couplings of the W and Z bosons, and searches for new particles. We review the important role that polarized beams play in the ILC program. The first stage of the ILC is planned to be a Higgs factory at 250 GeV in the centre of mass. Energy upgrades can naturally be implemented based on the concept of a linear collider. We discuss in detail the ILC program of Higgs boson measurements and the expected precision in the determination of Higgs couplings. We compare the ILC capabilities to those of the HL-LHC and to those of other proposed e + e - Higgs factories. We emphasize throughout that the readiness of the accelerator and the estimates of ILC performance are based on detailed simulations backed by extensive RandD and, for the accelerator technology, operational experience.

43 PARTICLE ACCELERATORS↗

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗

An Updated Simulation of the Booster Neutrino Beam

For accelerator neutrino experiments, an accurate prediction of the incoming neutrino flux is crucial for reducing uncertainties for all physics measurements. In this exciting period for the Short-Baseline Neutrino program at Fermilab, with far detector (ICARUS) already operating and the near detector (SBND) nearing operation, an updated flux model for the Booster Neutrino Beam (BNB) is presented. The BNB delivers 8 GeV protons to a beryllium target, subsequently producing neutrinos from the decay of the secondary beam of hadrons. A updated Monte-Carlo simulation of beam production in GEANT4 has been created, allowing predictions to be made for detectors with different baselines, offsets and sizes. This new simulation contains new features, such as a full neutrino ancestry to handle hadron production systematics with more precision, with a view to storing all resulting particles - including neutral mesons - from the proton-Beryllium scatter to allow the study of exotic BSM scenarios. Results are presented, with comparisons to the original flux simulated for the MiniBooNE experiment.

Paton, Josephine↗

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗

EXPERIMENTAL STUDIES OF NONLINEAR INTEGRABLE OPTICS

State-of-the-art accelerators at energy and intensity frontiers require increasingly bright and powerful particle beams. In conventional linear lattices, intense beams suffer from collective instabilities, resulting in beam losses and maximum beam intensity limits. This thesis presents experimental studies of a novel lattice design concept, the nonlinear integrable optics (NIO), aimed at enhancing beam stability limits with little to no beam dynamics degradation. Single-particle beam dynamics measurements of two NIO devices, the quasi-integrable octupole system and the fully integrable Danilov-Nagaitsev system, were carried out at the purpose-built Integrable Optics Test Accelerator (IOTA) at Fermilab. Their simulation, hardware design, and commissioning process are presented. Extensive model and analysis algorithm development and benchmarking is described. Electron beam data from two scientific runs is analyzed, yielding frequency and phase space dynamics consisten t with m odels. These results demonstrate viability and advantages of the NIO design, providing the groundwork for proton studies in the strong space-charge regime and future integrable accelerators.

Kuklev, Nikita↗

Precipitation of fast ion beams from the plasma sheet boundary layer

This paper presents a model of precipitated fluxes from the PSBL and CPS. Simulation results and data from Aureol-3 spacecraft indicate the presence of velocity dispersed precipitated ion structures (VDIS) at the poleward edge of the auroral oval. These structures are associated with fast ion beams in the PSBL region of the earth's magnetotail, confirming previous experimental results. The simulations also reveal possible substructuring of the VDIS. The bulk of the PSBL population which is not precipitated is very effectively thermalized and quasi-isotropized after multiple interactions with the magnetotail current layer. After each reflection cycle some part of the distribution is precipitated and forms multiple 'echoes' of VDIS. The CPS distributions occurring as a result of scattering, convection, multiple reflections and Fermi acceleration appear isotropic in the simulation model. This paper portrays the important role of the VDIS auroral region medium for complicated and energetically significant processes occurring in different regions of the distant magnetotail.

Ashour-Abdalla, M.↗

Dusty-Plasma Particle Accelerator

A dusty-plasma apparatus is being investigated as means of accelerating nanometer- and micrometer-sized particles. Applications for the dusty-plasma particle accelerators fall into two classes: Simulation of a variety of rapidly moving dust particles and micrometeoroids in outer-space environments that include micrometeoroid streams, comet tails, planetary rings, and nebulae and Deposition or implantation of nanoparticles on substrates for diverse industrial purposes that could include hardening, increasing thermal insulation, altering optical properties, and/or increasing permittivities of substrate materials. Relative to prior apparatuses used for similar applications, dusty-plasma particle accelerators offer such potential advantages as smaller size, lower cost, less complexity, and increased particle flux densities. A dusty-plasma particle accelerator exploits the fact that an isolated particle immersed in plasma acquires a net electric charge that depends on the relative mobilities of electrons and ions. Typically, a particle that is immersed in a low-temperature, partially ionized gas, wherein the average kinetic energy of electrons exceeds that of ions, causes the particle to become negatively charged. The particle can then be accelerated by applying an appropriate electric field. A dusty-plasma particle accelerator (see figure) includes a plasma source such as a radio-frequency induction discharge apparatus containing (1) a shallow cup with a biasable electrode to hold the particles to be accelerated and (2) a holder for the substrate on which the particles are to impinge. Depending on the specific design, a pair of electrostatic-acceleration grids between the substrate and discharge plasma can be used to both collimate and further accelerate particles exiting the particle holder. Once exposed to the discharge plasma, the particles in the cup quickly acquire a negative charge. Application of a negative voltage pulse to the biasable electrode results in the initiation of a low-current, high-voltage cathode spot. Plasma pressure associated with the cathode spot as well as the large voltage drop at the cathode spot accelerates the charged particles toward the substrate. The ultimate kinetic energy attained by particles exiting the particle holder depends in part on the magnitude of the cathode spot sheath potential difference, which is proportional to the magnitude of the voltage pulse, and the on the electric charge on the dust. The magnitude of the voltage pulse can be controlled directly, whereas the particle s electric charge can be controlled indirectly by controlling the operating parameters of the plasma apparatus.

Foster, John E.↗

Ion probe beam experiments and kinetic modeling in a dense plasma focus Z-pinch

The Z-pinch phase of a dense plasma focus (DPF) emits multiple-MeV ions in a ∼cm length. The mechanisms through which these physically simple devices generate such high energy beams in a relatively short distance are not fully understood. We are exploring the origins of these large gradients using measurements of an ion probe beam injected into a DPF during the pinch phase and the first kinetic simulations of a DPF Z-pinch. To probe the accelerating fields in our table top experiment, we inject a 4 MeV deuteron beam along the z-axis and then sample the beam energy distribution after it passes through the pinch region. Using this technique, we have directly measured for the first time the acceleration of an injected ion beam. Our particle-in-cell simulations have been benchmarked on both a kJ-scale DPF and a MJ-scale DPF. They have reproduced experimentally measured neutron yields as well as ion beams and EM oscillations which fluid simulations do not exhibit. Direct comparisons between the experiment and simulations enhance our understanding of these plasmas and provide predictive design capability for accelerator and neutron source applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hubble Space Telescope Pointing Performance Due to Micro-Dynamic Disturbances from the NICMOS Cryogenic Cooler

The Vibration Admittance Test (VET) was performed to measure the emitted disturbances of the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) Cryogenic Cooler (NCC) in preparation for NCC installation onboard the Hubble Space Telescope (HST) during Servicing Mission 3B (SM3B). Details of the VET ground-test are described, including facility characteristics, sensor complement and configuration, NCC suspension, and background noise measurements. Kinematic equations used to compute NCC mass center displacements and accelerations from raw measurements are presented, and dynamic equations of motion for the NCC VET system are developed and verified using modal test data. A MIMO linear frequency-domain analysis method is used to compute NCC-induced loads and HST boresight jitter from VET measurements. These results are verified by a nonlinear time-domain analysis approach using a high-fidelity structural dynamics and pointing control simulation for HST. NCC emitted acceleration levels not exceeding 35 micro-g rms were measured in the VET and analysis methods herein predict 3.1 milli-areseconds rms jitter for HST on-orbit. Because the NCC is predicted to become the predominant disturbance source for HST, VET results indicate that HST will continue to meet the 7 milli-arcsecond pointing stability mission requirement in the post-SM3B era.

Clapp, Brian R.↗

RADIATION DOSE SIMULATIONS FOR JEFFERSON LAB’S PERMANENT MAGNET RESILIENCY LDRD STUDY

In late 2023, Thomas Jefferson National Accelerator Facil- ity (Jefferson Lab) funded a Laboratory Directed Research and Development (LDRD) grant dedicated to investigating the impact of radiation on permanent magnet materials. This research initiative is specifically geared towards assessing materials slated for use in the CEBAF energy upgrade. The experimental approach involves strategically placing perma- nent magnet samples throughout the accelerator, exposing them to varying radiation doses. The simulation code BD- SIM is used to first validate the data and then to simulate the effects on future higher energy passes to study the degra- dation effects on the permanent magnets. In this paper we present the progress of that work.

Deitrick, K.↗

Accelerating transients with NekRS: GPU overlapping domain implementation and multi-rate timestepping

The simulation of nuclear transients using Computational Fluid Dynamics (CFD) presents significant computational challenges due to the inherent complexity and the wide separation in temporal scales between various flow physical phenomena. These disparities lead to high computational costs, often making the simulation of transients impractical without advanced techniques. Consequently, multiple research initiatives are being pursued by the NEAMS thermal-hydraulic area, some driven by academic institutions and some by national laboratories. Overall, they are exploring novel methods to make transient simulations more feasible and efficient. This report delves into recent advancements within the CFD code NekRS, specifically those achieved in Fiscal Year 2024 under the CONNECT effort, aimed at improving the performance and feasibility of transient simulations. The first major advancement involves the porting of NekRS to Aurora, one of the Department of Energy’s (DOE) most powerful supercomputers. Additionally, the report discusses the implementation of an overlapping domain capability within NekRS. This novel GPU-accelerated capability allows different spatial regions of the domain to be solved independently, enhancing the code’s efficiency, particularly when running large-scale simulations in complex domains. The scalability of this approach is demonstrated, highlighting its potential to transform how transients are approached in CFD simulations. Lastly, the report focuses on how this overlapping domain capability specifically accelerates transient simulations through multi-rate timestepping. By decoupling different regions and facilitating faster computations, this method offers a promising pathway to making nuclear transient simulations more computationally feasible, addressing one of the critical bottlenecks in the field. Together, these advancements represent a significant leap forward in transient simulation technology, bringing closer the possibility of handling highly complex nuclear scenarios with greater efficiency and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fermilab Booster Transition Crossing Simulations and Beam Studies

The Fermilab Booster accelerates beam from 400 MeV to 8 GeV at 15 Hz. In the PIP (Proton Improvement Plan) era, it is required that Booster deliver 4.2 x $10^{12}$ protons per pulse to extraction. One of the obstacles for providing quality beam to the users is the longitudinal quadrupole oscillation that the beam suffers from right after transition. Although this oscillation is well taken care of with quadrupole dampers, it is important to understand the source of these oscillations in light of the PIP II requirements that require 6.5 x $10^{12}$ protons per pulse at extraction. This paper explores the results from machine studies, computer simulations and solutions to prevent the quadrupole oscillations after transition.

Bhat, C. M.↗

Validation of Finite Element Crash Test Dummy Models for Predicting Orion Crew Member Injuries During a Simulated Vehicle Landing

A series of crash tests were conducted with dummies during simulated Orion crew module landings at the Wright-Patterson Air Force Base. These tests consisted of several crew configurations with and without astronaut suits. Some test results were collected and are presented. In addition, finite element models of the tests were developed and are presented. The finite element models were validated using the experimental data, and the test responses were compared with the computed results. Occupant crash data, such as forces, moments, and accelerations, were collected from the simulations and compared with injury criteria to assess occupant survivability and injury. Some of the injury criteria published in the literature is summarized for completeness. These criteria were used to determine potential injury during crew impact events.

Tabiei, Al↗

Technical Note: Motion-perturbation method applied to dosimetry of dynamic MLC target tracking—A proof-of-concept

Purpose: Previous studies show that dose to a moving target can be estimated using 4D measurement-guided dose reconstruction based on a process called virtual motion simulation, or VMS. A potential extension of VMS is to estimate dose during dynamic multileaf collimator (MLC)-tracking treatments. The authors introduce a modified VMS method and quantify its performance as proof-of-concept for tracking applications. Methods: Direct measurements with a moving biplanar diode array were used to verify accuracy of the VMS dose estimates. A tracking environment for variably sized circular MLC apertures was simulated by sending preprogrammed control points to the MLC while simultaneously moving the accelerator treatment table. Sensitivity of the method to simulated tracking latency (0–700 ms) was also studied. Potential applicability of VMS to fast changing beam apertures was evaluated by modeling, based on the demonstrated dependence of the cumulative dose on the temporal dose gradient. Results: When physical and virtual latencies were matched, the agreement rates (2% global/2 mm gamma) between the VMS and the biplanar dosimeter were above 96%. When compared to their own reference dose (0 induced latency), the agreement rates for VMS and biplanar array track closely up to 200 ms of induced latency with 10% low-dose cutoff threshold and 300 ms with 50% cutoff. Time-resolved measurements suggest that even in the modulated beams, the error in the cumulative dose introduced by the 200 ms VMS time resolution is not likely to exceed 0.5%. Conclusions: Based on current results and prior benchmarks of VMS accuracy, the authors postulate that this approach should be applicable to any MLC-tracking treatments where leaf speeds do not exceed those of the current Varian accelerators.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Beam-Dynamics Analysis of Long-Range Wakefield Effects on the SCRF Cavities at the Fast Facility

Long-range wakefields in superconducting RF (SCRF) cavities create complicated effects on beam dynamics in SCRF-based FEL beamlines. The driving bunch excites effectively an infinite number of structure modes (including HOMs) which oscillate within the SCRF cavity. Couplers with loads are used to damp the HOMs. However, these HOMs can persist for long periods of time in superconducting structures, which leads to long-range wakefields. Clear understanding of the long-range wakefield effects is a critical element for risk mitigation of future SCRF accelerators such as XFEL at DESY, LCLS-II XFEL, and MaRIE XFEL. We are currently developing numerical tools for simulating long-range wakefields in SCRF accelerators and plan to experimentally verify the tools by measuring these wakefields at the Fermilab Accelerator Science and Technology (FAST) facility. This paper previews the experimental conditions at the FAST 50 MeV beamline based on the simulation results.

43 PARTICLE ACCELERATORS↗

Modeling laser-driven ion acceleration with deep learning

Developments in machine learning promise to ameliorate some of the challenges of modeling complex physical systems through neural-network-based surrogate models. High-intensity, short-pulse lasers can be used to accelerate ions to mega-electronvolt energies, but to model such interactions requires computationally expensive techniques such as particle-in-cell simulations. Multilayer neural networks allow one to take a relatively sparse ensemble of simulations and generate a surrogate model that can be used to rapidly search the parameter space of interest. In this study, we created an ensemble of over 1,000 simulations modeling laser-driven ion acceleration and developed a surrogate to study the resulting parameter space. A neural-network-based approach allows for rapid feature discovery not possible for traditional parameter scans given the computational cost. A notable observation made during this study was the dependence of ion energy on the pre-plasma gradient length scale. While this methodology harbors great promise for ion acceleration, it has ready application to all topics in which large-scale parameter scans are restricted by significant computational cost or relatively large, but sparse, domains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗