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,009 records · Page 56

Radiation from laser-microplasma-waveguide interactions in the ultra-intense regime

When a high-contrast ultra-relativistic (>10{sup 20} W/cm{sup 2}) laser beam enters a micro-sized plasma waveguide, the pulse energy is coupled into waveguide modes, which significantly modifies the interaction between the electrons and electromagnetic wave. Electrons pulled out from the walls of the waveguide form a dense helical bunch inside the channel and are efficiently accelerated by the transverse magnetic modes to hundreds of MeV. The asymmetry in the transverse electric and magnetic fields drives strong oscillations, which lead to the emission of bright, well-collimated, hard X-rays. In this paper, we present our study on the underlying physics in the aforementioned process using 3D particle-in-cell simulations. The mechanism of electron acceleration and the dependence of radiation properties on different laser plasma parameters are addressed. An analytic model and basic scalings for X-ray emission are also presented by considering the lowest optical modes in the waveguide, which is adequate to describe the basic phenomenon. In addition, the effects of high-order modes as well as laser polarization are also qualitatively discussed. The considered X-ray source has promising features, potentially making it a competitive candidate for a future tabletop synchrotron source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling laser-wakefield accelerators using the time-averaged ponderomotive approximation in a Lorentz boosted frame

Future, high-fidelity simulations of multi-GeV-class laser Wakefield accelerators (LWFAs) will need to model the propagation of high-intensity laser drivers over meter-scale plasmas with high spatial and temporal resolutions, thus requiring high amounts of computational resources. Various techniques have been devised over the years to reduce the computational cost of such simulations, including the time-averaged ponderomotive approximation, and the use of the Lorentz boosted frame technique. In this paper we discuss the combination of these two computational techniques, highlighting the resulting significant reduction in the computational cost of LWFA simulations and the limitations of this approach. The combination of the two techniques can potentially become essential for the modeling of a multi-TeV, LWFA-based collider.

Laser Wakefield Acceleration↗

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗

Advances in the simulation of space-charge effects

The nonlinear space-charge effects are an important topic in high intensity accelerator beam dynamics and have been extensively studied using macroparticle tracking simulations. Here, we report on recent advances in the simulation of space-charge effects using a symplectic multiparticle space-charge model. The transverse space-charge limit was explored for a periodic focusing and defocusing lattice. The artificial numerical emittance growth in the macroparticle space-charge simulation was analyzed using a one-dimensional model.

47 OTHER INSTRUMENTATION↗

Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm

Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.

43 PARTICLE ACCELERATORS↗

Pandemic drugs at pandemic speed: infrastructure for accelerating COVID-19 drug discovery with hybrid machine learning- and physics-based simulations on high-performance computers

The race to meet the challenges of the global pandemic has served as a reminder that the existing drug discovery process is expensive, inefficient and slow. There is a major bottleneck screening the vast number of potential small molecules to shortlist lead compounds for antiviral drug development. New opportunities to accelerate drug discovery lie at the interface between machine learning methods, in this case, developed for linear accelerators, and physics-based methods. The two in silico methods, each have their own advantages and limitations which, interestingly, complement each other. Here, we present an innovative infrastructural development that combines both approaches to accelerate drug discovery. The scale of the potential resulting workflow is such that it is dependent on supercomputing to achieve extremely high throughput. We have demonstrated the viability of this workflow for the study of inhibitors for four COVID-19 target proteins and our ability to perform the required large-scale calculations to identify lead antiviral compounds through repurposing on a variety of supercomputers.

97 MATHEMATICS AND COMPUTING↗

Supervised paralled-in-time algorithm for long-time Lagrangian simulations of stochastic dynamics: Application to hydrodynamics

Lagrangian particle methods based on detailed atomic and molecular and nanoscale systems. However, the maximum time step is limited by the smallest oscillation period of the fastest atomic motion, rendering long-time simulations very expensive. To resolve this bottleneck, we propose a supervised parallel-in-time algorithm for stochastic dynamics to accelerate long-time Lagrangian particle simulations. Our method is inspired by bottom-up coarse-graining projections that yield mean-field hydrodynamic behavior in the continuum limit.

Multiscale modeling, Lagrangian method, Parallel-i↗

Understanding thermal and mechanical effects on lithium plating in lithium-ion batteries

Lithium plating is one of the main concerns for cell durability, as cyclable lithium and cell capacity are decreased during the process. Thermal and mechanical loads can influence cell behavior, including lithium plating driving forces, as suggested by the Arrhenius equation and the Bruggeman relation, respectively. Hence, a 1D and a 3D model are employed to study the thermal and mechanical effects on lithium plating in a lithium-ion battery with a high areal capacity (4 mAh/cm 2 ). The 1D model assumes uniform mechanical and temperature profile, while the 3D model is capable of capturing the nonuniformity of stress and temperature in the cell. Further, the effects of mechanical deformation and different temperatures on the process of lithium plating are investigated in this study. With the help of the 1D charging simulations, it is observed that lower temperatures increase the risk of lithium plating and the presence of mechanical deformations can accelerate its initiation. Furthermore, 3D charging simulations indicated that folds and boundaries of jellyrolls are more susceptible to lithium plating due to the higher stress concentration and lower temperature, respectively.

25 ENERGY STORAGE↗

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of a Green Parabolic Trough Collector

This report describes in detail a research project that assesses the potential of a parabolic trough solar collector with a lumber construction to meet ambitious performance, cost, and durability goals. Instead of exploring lumber as a general question for solar collectors, the project focuses on the design optimization, laboratory testing, outdoor testing, and economic modeling of a particular scale of parabolic trough, namely, the 6m aperture and 12m long module that has become the standard in the utility-scale power plant market. As will be appreciated in the report, it is also a scale that matches well the capabilities of lumber as a structural material and the commonly available cross sections. The design and its analysis extend far beyond traditional considerations for wooden structures in the built environment, and includes a computationally efficient coupled structural-optical simulation and the empirical characterization of joint stiffness at small deformations. Laboratory testing is used to characterize cyclic joint behavior, accelerated weathering, and to validate the simulation used as the basis for building a full-scale prototype. The construction and outdoor testing of this prototype is described. A total installed cost model of the collector is presented and the impact of the cost is seen in the context of a full plant techno-economic model. The results of this model in combination with the experience with the first prototype suggest that with improvements to the design and assembly process, this collector type has the potential to lower the levelized cost of heat for industrial solar steam plants used for process heating or cogeneration.

14 SOLAR ENERGY↗

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↗

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↗

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.↗