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,531 records · Page 85

Dual Energy Storage Ring Project

The overall goals of the project were to evaluate the possibility of an electron cooler based on the novel two-energy ring idea, by designing a suitable ring, and by simulating actual ring performance characteristics with the accelerator physics and design program ELEGANT. Ring parameters were first chosen to be consistent with Jefferson Lab’s EIC design. About mid-way through the project, parameters were adjusted to be consistent with a cooler design for the new EIC project at Brookhaven National Lab.

43 PARTICLE ACCELERATORS↗

Project DarkStar: Vision for LLNL in 2030

DarkStar was a Strategic Initiative (FY2021-FY2024) to investigate applications of Artificial Intelligence (AI) and Machine Learning (ML) to scientific problems of complex hydrodynamics, shockwave physics and energetic materials. The research focused on physics and engineering design as a process that can be tremendously accelerated through merging AI with advanced physics simulation on exascale-class platforms, and to experimentally validate this revolutionary new approach through dynamic materials campaigns. A central thread of scientific inquiry was in the application of AI to enable human understanding of how to control hydrodynamic instability (which has impacts to areas such as inertial confinement fusion) via engineering features and time-dependent sources. Motivated by an unfinished line of research started by Dr. Johnny von Neumann, AI-enabled simulation approaches were developed that allowed DarkStar researchers to uncover several ground-breaking discoveries regarding hydrodynamic instability, including how to completely suppress Richtmyer-Meshkov instability (RMI). These S&T discoveries, along with other advances, have shown the way for an entirely new approach to time-dependent problems known as inverse design – the idea that complex systems can be developed directly from a final state that is to be achieved and resolve the initial design via satisfying several constraints simultaneously via AI/ML. Through experimental campaigns conducted across a wide range of facilities in the NNSA complex (the High Explosive Application Facility at LLNL, the Dynamic Compression Sector/Advanced Photon Source at Argonne National Lab, and Special Technologies Laboratory at MSTS) the radical new AI/ML approach to engineering complex material dynamics was verified, establishing a new field of study within the realm of shock physics. As advanced manufacturing capabilities continue to develop, the great importance of inverse design as a means to apply that technology effectively for NNSA missions will feature prominently over this decade. DarkStar has positioned NNSA as a world-leader in this newly emerging cross-disciplinary area of AI methods for advanced physics simulation and pioneered multiple novel approaches that have enabled the broader scientific community. By allowing us to see past the horizon, to 2030 and beyond, DarkStar has illuminated the vast potential of AI/ML to impact a wide range of new national security missions and, consequently, multiple areas of further research have already emerged across the NNSA and DOD complex.

42 ENGINEERING↗

Aerodynamic Characterization of 3D Scanned Wind Turbine Blades Using Experimental and Computational Methods

This study presents an aerodynamic characterization of 3D scanned wind turbine blades using both experimental and computational methods. The research was conducted by Gulf Wind Technology and Sandia National Laboratories. The primary objective was to investigate the aerodynamic impacts of leading-edge manufacturing defects on wind turbine blades. The study utilized the Stratasys NEO 800 3D Printer for high-precision manufacturing and the GWT Accelerator Wind Tunnel for experimental testing. Computational simulations were performed using COMSOL Multiphysics to model the wind tunnel and analyze flow characteristics and OpenFOAM to study the aerodynamic impacts of leading-edge defects. OpenFAST was used to estimate how these defects can lead to revenue losses for wind farm operators as high as 6%. The results demonstrated significant aerodynamic performance variations due to defects, with detailed analysis provided through wind tunnel and CFD data. The findings contribute to the understanding of defect impacts on wind turbine blade performance and offer insights for future design improvements.

17 WIND ENERGY↗

Short-Pulse Photoneutron Production on Beryllium Using the Mercury Pulsed Power X-ray Source

An x-ray beam can produce neutrons using beryllium as a photoneutron converter because of the relatively low 9 Be(ɤ,n) reaction energy threshold. Previous photoneutron production simulations have focused on using a linear accelerator and heavy water to generate photoneutrons for potential medical therapies. This experimental work focused on creating a short, intense neutron pulse using the Mercury pulsed power machine at the U.S. Naval Research Laboratory, which fired a Bremsstrahlung photon beam with an energy endpoint of 4.8 MeV. Several independent detector systems were included in the experiment to measure both the characteristics of the x-ray pulse and resulting photoneutrons. We discovered the neutron pulse follows the time history of the x-ray pulse, enabling neutron pulses as short as 30 ns with less than 10% shot-to-shot variation. The approach offers a high-flux of neutrons with a higher end-point energy than using a heavy water converter. Previous experiments at Mercury have generated fission neutrons using a higher x-ray endpoint and 238 U(n,f). Results indicate that with an optimized diode and beryllium converter target design, neutron production could be further increased.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Neutron structure with (un)polarized deuterons and forward spectator tagging at the Electron-Ion Collider

An Electron-Ion Collider (EIC) would enable novel measurements of neutron structure through deep-inelastic electron-deuteron scattering with coincidence tagging of the forward-moving spectator proton. The free neutron structure functions can be obtained by extrapolating the measured recoil momentum distributions to the on-shell point. Such measurements provide essential information for the flavor separation of the nucleon parton densities, the nucleon spin decomposition, and precision studies of QCD evolution in the flavor-singlet and nonsinglet sectors. The EMC effect in light nuclei can be elucidated by measuring the recoil momentum dependence of the nuclear modification. A Jefferson Lab LDRD project has supported a study of neutron structure measurements with spectator tagging at an EIC, with Monte Carlo simulation of the physical processes and the accelerator and detector requirements. We present results of a model-independent extraction of the free neutron structure function $F_2^n$, as well as the spin structure function $g_1^n$, through on-shell extrapolation in the spectator momentum.

Park, Kijun↗

STABILITY OF COSMIC-RAY MODIFIED SHOCKS: TWO-FLUID APPROACH

The stability of cosmic-ray modified shocks (CRMSs) is studied by means of numerical simulations. Owing to the nonlinear feedback of cosmic-ray (CR) acceleration, a downstream state of the modified shock can no longer be uniquely determined for given upstream parameters. It is known that up to three distinct solutions exist, which are characterized by CR production efficiency as the 'efficient', 'intermediate', and 'inefficient' branches. The stability of these solutions is investigated by performing direct time-dependent simulations of a two-fluid model. It is found that both the efficient and inefficient branches are stable even against a large-amplitude perturbation, while the intermediate one is always unstable and evolves into the inefficient state as a result of nonlinear time development. This bistable feature is robust in a wide range of parameters and does not depend on the injection model. Fully nonlinear time evolution of a hydrodynamic shock with injection results in the least efficient state in terms of CR production, consistent with the bistable feature. This suggests that the CR production efficiency in supernova remnant shocks may be lower than previously discussed in the framework of the nonlinear shock acceleration theory considering the efficient solution of the CRMS.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the beam halo in negative ion sources by using 3D3V PIC code

The physical mechanism of the formation of the negative ion beam halo and the heat loads of the multi-stage acceleration grids are investigated with the 3D PIC (particle in cell) simulation. The following physical mechanism of the beam halo formation is verified: The beam core and the halo consist of the negative ions extracted from the center and the periphery of the meniscus, respectively. This difference of negative ion extraction location results in a geometrical aberration. Furthermore, it is shown that the heat loads on the first acceleration grid and the second acceleration grid are quantitatively improved compared with those for the 2D PIC simulation result.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Optimized matching conditions for self-guided laser wakefield accelerators

We revisit the matching conditions for self-guided laser pulse propagation in plasma and refine their formulation to maximize the energy of electrons produced via laser wakefield acceleration. Bayesian optimization, combined with particle-in-cell simulations carried out in a quasi-three-dimensional geometry and a Lorentz-boosted frame, is employed. The optimization identifies the maximum electron energy that a self-guided laser wakefield accelerator, driven by a laser of a given energy, can produce, together with the corresponding acceleration distance. Our results further demonstrate that electrons with energies close to the maximum value can be obtained across a relatively wide range of input parameters and without the need for their precise tuning. This provides substantial flexibility for experimental implementation and significantly relaxes the operational constraints associated with self-guided laser wakefield accelerators.

43 PARTICLE ACCELERATORS↗

A NUMERICAL MODEL OF PARSEC-SCALE SSC MORPHOLOGIES AND THEIR RADIO EMISSION

In current models for jets of active galactic nuclei and their emission a shortcoming in the description and understanding of the connection between the largest and smallest scales exists. In this work we present a spatially resolved synchrotron self-Compton model extended to parsec scales, which opens the possibility of probing the connections between the radio and high-energy properties. We simulate an environment that leads to Fermi-I acceleration of leptonic particles and includes the full time dependence of this process. Omitting the restriction of a finite downstream region, we find that the spectral energy distribution produced by the accelerated particles strongly depends on their radial confinement behind the shock. The requirement, for both the restriction of high-energy emission to a small region around the shock and the production of a flat radio spectrum, is an initial linear increase of the radius immediately behind the shock, which then slows down with increasing distance from the shock. A good representation of the data for the blazar Mrk 501 is achieved by a parameterized log function. The prediction for the shape of the radio blob is given by the flux distribution with respect to shock distance.

79 ASTRONOMY AND ASTROPHYSICS↗

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Collisionless Weibel shocks and electron acceleration in gamma-ray bursts

A study of collisionless external shocks in gamma-ray bursts is presented. The shock structure, electromagnetic field, and process of electron acceleration are assessed by performing a self-consistent 3D particle-in-cell simulation. In accordance with hydrodynamic shock systems, the shock consists of a reverse shock (RS) and forward shock separated by a contact discontinuity. The development and structure are controlled by the ion Weibel instability. The ion filaments are sources of strong transverse electromagnetic fields at both sides of the double shock structure over a length of 30–100 ion skin depths. Electrons are heated up to a maximum energy ϵ{sub ele}≈√(ϵ{sub b}), where ϵ is the energy normalized to the total incoming energy. Jet electrons are trapped in the RS transition region due to the presence of an ambipolar electric field and reflection by the strong transverse magnetic fields in the shocked region. In a process similar to shock surfing acceleration for ions, electrons experience drift motion and acceleration by ion filament transverse electric fields in the plane perpendicular to the shock propagation direction. Ultimately, accelerated jet electrons are convected back into the upstream.

79 ASTRONOMY AND ASTROPHYSICS↗

High contrast ion acceleration at intensities exceeding 10{sup 21} Wcm{sup −2}

Ion acceleration from short pulse laser interactions at intensities of 2×10{sup 21}Wcm{sup −2} was studied experimentally under a wide variety of parameters, including laser contrast, incidence angle, and target thickness. Trends in maximum proton energy were observed, as well as evidence of improvement in the acceleration gradients by using dual plasma mirrors over traditional pulse cleaning techniques. Extremely high efficiency acceleration gradients were produced, accelerating both the contaminant layer and high charge state ions from the bulk of the target. Two dimensional particle-in-cell simulations enabled the study of the influence of scale length on submicron targets, where hydrodynamic expansion affects the rear surface as well as the front. Experimental evidence of larger electric fields for sharp density plasmas is observed in simulation results as well for such targets, where target ions are accelerated without the need for contaminant removal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Piecewise empirical model (PEM) of resistive memory for pulsed analog and neuromorphic applications

With the end of Dennard scaling and the ever-increasing need for more efficient, faster computation, resistive switching devices (ReRAM), often referred to as memristors, are a promising candidate for next generation computer hardware. These devices show particular promise for use in an analog neuromorphic computing accelerator as they can be tuned to multiple states and be updated like the weights in neuromorphic algorithms. Modeling a ReRAM-based neuromorphic computing accelerator requires a compact model capable of correctly simulating the small weight update behavior associated with neuromorphic training. These small updates have a nonlinear dependence on the initial state, which has a significant impact on neural network training. Consequently, we propose an empirically derived general purpose model that can accurately capture the nonlinearity of an arbitrary two-terminal device to match pulse measurements important for neuromorphic computing applications. By defining the state of the device to be proportional to its current, the model parameters can be extracted from a series of voltages pulses that mimic the behavior of a device in an analog neuromorphic computing accelerator. This allows for a general, accurate, and intuitive compact circuit model that is applicable to different resistance switching device technologies. In this work we explain the details of the model, implement the model in the circuit simulator Xyce, and give an example of its usage to model a specific Ta/TaO x device.

97 MATHEMATICS AND COMPUTING↗

Numerical and Analytical Analysis of Thin Laser-Driven Flyer Plates

Computational and analytical investigations have been performed for laser-driven flyer plate experiments and compared with data. These have resulted in a more fundamental understanding of the physical processes in the laser absorption, and acceleration of the flyers by the expansion of laser produced plasma. The analytical work describes the dependences of final foil velocity, while the computational simulations permit a detailed look at energy coupling and foil acceleration issues.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Design of a high-power, high-gain, 2nd harmonic, 22.848 GHz gyroklystron

In this paper we consider the design of a four-cavity, high-gain K-band gyroklystron experiment for high gradient structure testing. The frequency doubling gyroklystron utilizes a beam voltage of 500 kV and a beam current of 200 A from a magnetron injection gun (MIG) originally designed for a lower-frequency device. The microwave circuit features input and gain cavities in the circular TE{sub 011} mode and penultimate and output cavities that operate at the second harmonic in the TE{sub 021} mode. We investigate the MIG performance and study the behavior of the circuit for different values of perpendicular to parallel velocity ratio (α= V{sub ⊥}/ V{sub z}). This microwave tube is expected to be able to produce at least 20 MW of power in 1μs pulses at a repetition rate of at least 120 Hz. A maximum efficiency of 26% and a large signal gain of 58 dB under zero-drive stable conditions were simulated for a velocity ratio equal to 1.35.

43 PARTICLE ACCELERATORS↗

Giga-electronvolt electrons due to a transition from laser wakefield acceleration to plasma wakefield acceleration

We show through experiments that a transition from laser wakefield acceleration (LWFA) regime to a plasma wakefield acceleration (PWFA) regime can drive electrons up to energies close to the GeV level. Initially, the acceleration mechanism is dominated by the bubble created by the laser in the nonlinear regime of LWFA, leading to an injection of a large number of electrons. After propagation beyond the depletion length, leading to a depletion of the laser pulse, whose transverse ponderomotive force is not able to sustain the bubble anymore, the high energy dense bunch of electrons propagating inside bubble will drive its own wakefield by a PWFA regime. This wakefield will be able to trap and accelerate a population of electrons up to the GeV level during this second stage. Three dimensional particle-in-cell simulations support this analysis and confirm the scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Recent developments in the PySCF program package

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. Here, we explain the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗