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

Simplified Predictive Models for CO 2 Sequestration Performance Assessment: Research Topical Report on Task #4 - Reduced-Order Method (ROM) Based Models

Reduced-order models provide a means for greatly accelerating the detailed simulations that will be required to manage CO 2 storage operations. In this work, we investigate the use of one such method, POD-TPWL, which has previously been shown to be effective in oil reservoir simulation problems. This method combines trajectory piecewise linearization (TPWL), in which the solution to a new (test) problem is represented through a linearization around the solution to a previously-simulated (training) problem, with proper orthogonal decomposition (POD), which enables solution states to be expressed in terms of a relatively small number of parameters. We describe the application of POD-TPWL for CO 2 -water systems simulated using a compositional procedure. Stanford’s Automatic Differentiation-based General Purpose Research Simulator (AD-GPRS) performs the full-order training simulations and provides the output (derivative matrices and system states) required by the POD-TPWL method. A new POD-TPWL capability introduced in this work is the use of horizontal injection wells that operate under rate (rather than bottom-hole pressure) control. Simulation results are presented for CO 2 injection into a synthetic aquifer and into a simplified model of the Mount Simon formation. Test cases involve the use of time-varying well controls that differ from those used in training runs. Results of reasonable accuracy are consistently achieved for relevant well quantities. Runtime speedups of around a factor of 370 relative to full- order AD-GPRS simulations are achieved, though the preprocessing needed for POD-TPWL model construction corresponds to the computational requirements for about 2.3 full-order simulation runs. A preliminary treatment for POD-TPWL modeling in which test cases differ from training runs in terms of geological parameters (rather than well controls) is also presented. Results in this case involve only small differences between training and test runs, though they do demonstrate that the approach is able to capture basic solution trends. The impact of some of the detailed numerical treatments within the POD-TPWL formulation is considered in an Appendix.

Mishra, Srikanta↗

A detailed examination of laser-ion acceleration mechanisms in the relativistic transparency regime using tracers

Here we present a particle-in-cell study of linearly polarized laser-ion acceleration systems, in which we use both two-dimensional (2D) and three-dimensional (3D) simulations to characterize the ion acceleration mechanisms in targets which become transparent to the laser pulse during irradiation. First, we perform a target length scan to optimize the peak ion energies in both 2D and 3D, and the predictive capabilities of 2D simulations are discussed. Tracer analysis allows us to isolate the acceleration into stages of target normal sheath acceleration (TNSA), hole boring (HB), and break-out afterburner (BOA) acceleration, which vary in effectiveness based on the simulation parameters. The thinnest targets reveal that enhanced TNSA is responsible for accelerating the most energetic ions, whereas the thickest targets have ions undergoing successive phases of HB and TNSA (in 2D) or BOA and TNSA (in 3D); HB is not observed to be a dominant acceleration mechanism in the 3D simulations. It is in the intermediate optimal regime, both when the laser breaks through the target with appreciable amplitude and when there is enough plasma to form a sustained high density flow, that BOA is most effective and is responsible for the most energetic ions. Eliminating the transverse laser spot size effects by performing a plane wave simulation, we can isolate with greater confidence the underlying physics behind the ion dynamics we observe. Specifically, supplemented by wavelet and FFT analyses, we match the post-transparency BOA acceleration with a wave-particle resonance with a high-amplitude low-frequency electrostatic wave of increasing phase velocity, consistent with that predicted by the Buneman instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Symplectic machine learning model for fast simulation of space-charge effects

Symplectic simulation of space-charge effects is crucial for the design and operation of high-intensity particle accelerators. Traditional methods for simulating these effects are often computationally expensive, resulting in significant overhead. In this work, we introduce a generative model based on a U-Net architecture within a generative adversarial network framework to efficiently simulate space-charge effects. The model is trained to predict the transverse multiparticle space-charge Hamiltonian, which can be physically computed using a gridless spectral method. The one-step symplectic transverse transfer map for the particles is then obtained by differentiating the predicted Hamiltonian. Benchmarking results demonstrate that this generative model achieves an order of magnitude higher computational efficiency compared to the spectral method, providing a highly efficient alternative for simulating space-charge effects with a large number of particles. By maintaining symplecticity, the model effectively preserves the phase-space structure and mitigates nonphysical errors in long-term simulations. This model has been integrated into jutrack, a novel autodifferentiable accelerator modeling code developed in the julia programming language.

Beam code development & simulation techniques↗

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control of target-normal-sheath-accelerated protons from a guiding cone

It is demonstrated through particle-in-cell simulations that target-normal-sheath-accelerated protons can be well controlled by using a guiding cone. Compared to a conventional planar target, both the collimation and number density of proton beams are substantially improved, giving a high-quality proton beam which maintained for a longer distance without degradation. The effect is attributed to the radial electric field resulting from the charge due to the hot target electrons propagating along the cone surface. This electric field can effectively suppress the spatial spread of the protons after the expansion of the hot electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle acceleration in laser-driven magnetic reconnection

Particle acceleration induced by magnetic reconnection is thought to be a promising candidate for producing the nonthermal emissions associated with explosive phenomena such as solar flares, pulsar wind nebulae, and jets from active galactic nuclei. Laboratory experiments can play an important role in the study of the detailed microphysics of magnetic reconnection and the dominant particle acceleration mechanisms. We have used two- and three-dimensional particle-in-cell simulations to study particle acceleration in high Lundquist number reconnection regimes associated with laser-driven plasma experiments. For current experimental conditions, we show that nonthermal electrons can be accelerated to energies more than an order of magnitude larger than the initial thermal energy. The nonthermal electrons gain their energy mainly from the reconnection electric field near the X points, and particle injection into the reconnection layer and escape from the finite system establish a distribution of energies that resembles a power-law spectrum. Energetic electrons can also become trapped inside the plasmoids that form in the current layer and gain additional energy from the electric field arising from the motion of the plasmoid. We compare simulations for finite and infinite periodic systems to demonstrate the importance of particle escape on the shape of the spectrum. Based on our findings, we provide an analytical estimate of the maximum electron energy and threshold condition for observing suprathermal electron acceleration in terms of experimentally tunable parameters. We also discuss experimental signatures, including the angular distribution of the accelerated particles, and construct synthetic detector spectra. Finally, these results open the way for novel experimental studies of particle acceleration induced by reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

INJECTION TO RAPID DIFFUSIVE SHOCK ACCELERATION AT PERPENDICULAR SHOCKS IN PARTIALLY IONIZED PLASMAS

We present a three-dimensional hybrid simulation of a collisionless perpendicular shock in a partially ionized plasma for the first time. In this simulation, the shock velocity and upstream ionization fraction are v {sub sh} ≈ 1333 km s{sup −1} and f {sub i} ∼ 0.5, which are typical values for isolated young supernova remnants (SNRs) in the interstellar medium. We confirm previous two-dimensional simulation results showing that downstream hydrogen atoms leak into the upstream region and are accelerated by the pickup process in the upstream region, and large magnetic field fluctuations are generated both in the upstream and downstream regions. In addition, we find that the magnetic field fluctuations have three-dimensional structures and the leaking hydrogen atoms are injected into the diffusive shock acceleration (DSA) at the perpendicular shock after the pickup process. The observed DSA can be interpreted as shock drift acceleration with scattering. In this simulation, particles are accelerated to v ∼ 100 v {sub sh} ∼ 0.3 c within ∼100 gyroperiods. The acceleration timescale is faster than that of DSA in parallel shocks. Our simulation results suggest that SNRs can accelerate cosmic rays to 10{sup 15.5} eV (the knee) during the Sedov phase.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling of ion acceleration through drift and diffusion at interplanetary shocks

A test particle simulation designed to model ion acceleration through drift and diffusion at interplanetary shocks is described. The technique consists of integrating along exact particle orbits in a system where the angle between the shock normal and mean upstream magnetic field, the level of magnetic fluctuations, and the energy of injected particles can assume a range of values. The technique makes it possible to study time-dependent shock acceleration under conditions not amenable to analytical techniques. To illustrate the capability of the numerical model, proton acceleration was considered under conditions appropriate for interplanetary shocks at 1 AU, including large-amplitude transverse magnetic fluctuations derived from power spectra of both ambient and shock-associated MHD waves.

Decker, R. B.↗

Accelerating laser ray tracing in high fidelity physics simulations of laser melting using squeeze U-net

Laser melting is a core component of the ongoing industrial revolution, dubbed Industry 4.0, as lasers facilitate fast and precise melting and fusion in advanced manufacturing. There is a strong need to optimize the laser process using simulations. However, this has proven challenging as high fidelity simulations are needed for predictive modeling and this is currently prohibitively expensive even when run on hundreds of processors on high performance computers. The challenge is capturing complex physics of laser material interaction, fluid dynamics, thermal physics and material phase transformations at various length and time scales. To close this technological gap, we modified a squeeze U-net to accelerate the laser ray tracing component of such high fidelity models by ~4x–40x while preserving the core physics principle of conservation of energy with 97% accuracy. This approach enables the accurate modeling of global laser energy absorption as a function of local surface temperatures and complex surface topologies, which govern the reflection directions and energy losses of laser rays upon interacting with the material surface.

Computer science↗

Start-to-end beam dynamics simulation of double triangular current profile generation in Argonne Wakefield Accelerator

Double triangular current profile (DT) gives a high transformer ratio which is the determining factor of the performance of collinear wakefield accelerator. This current profile can be generated using the emittance exchange (EEX) beam line. Argonne Wakefield Accelerator (AWA) facility plans to generate DT using the EEX beam line. We conducted start-to-end simulation for the AWA beam line using PARMELA code. Also, we discuss requirements of beam parameters for the generation of DT.

43 PARTICLE ACCELERATORS↗

Vortex paradigm for shock-accelerated density-stratified interfaces

Numerical simulations are reported that capture the observations of Haas and Sturtevant (1985) when a planar shock accelerates an inclined planar interface between two media. A vortex paradigm is presented for interpreting the evolution of shock-accelerated density-stratified interfaces beyond early times. The simulations illustrate the paradigm and are in good agreement with the results of Haas and Sturtevant. Color images of space-time diagrams of the one-space integrated vorticity function are given to show the primary and secondary features of the emerging vortex structures.

Hawley, John F.↗

Global adjoint tomography: first-generation model

We present the first-generation global tomographic model constructed based on adjoint tomography, an iterative full-waveform inversion technique. Synthetic seismograms were calculated using GPU-accelerated spectral-element simulations of global seismic wave propagation, accommodating effects due to 3-D anelastic crust & mantle structure, topography & bathymetry, the ocean load, ellipticity, rotation, and self-gravitation. Fréchet derivatives were calculated in 3-D anelastic models based on an adjoint-state method. The simulations were performed on the Cray XK7 named ‘Titan’, a computer with 18 688 GPU accelerators housed at Oak Ridge National Laboratory. The transversely isotropic global model is the result of 15 tomographic iterations, which systematically reduced differences between observed and simulated three-component seismograms. Our starting model combined 3-D mantle model S362ANI with 3-D crustal model Crust2.0. We simultaneously inverted for structure in the crust and mantle, thereby eliminating the need for widely used ‘crustal corrections’. We used data from 253 earthquakes in the magnitude range 5.8 ≤ M w ≤ 7.0. We started inversions by combining ~30 s body-wave data with ~60 s surface-wave data. The shortest period of the surface waves was gradually decreased, and in the last three iterations we combined ~17 s body waves with ~45 s surface waves. We started using 180 min long seismograms after the 12th iteration and assimilated minor- and major-arc body and surface waves. The 15th iteration model features enhancements of well-known slabs, an enhanced image of the Samoa/Tahiti plume, as well as various other plumes and hotspots, such as Caroline, Galapagos, Yellowstone and Erebus. Furthermore, we see clear improvements in slab resolution along the Hellenic and Japan Arcs, as well as subduction along the East of Scotia Plate, which does not exist in the starting model. Point-spread function tests demonstrate that we are approaching the resolution of continental-scale studies in some areas, for example, underneath Yellowstone. Here, this is a consequence of our multiscale smoothing strategy in which we define our smoothing operator as a function of the approximate Hessian kernel, thereby smoothing gradients less wherever we have good ray coverage, such as underneath North America.

58 GEOSCIENCES↗

Soft Classification of Diffractive Interactions at the LHC

Multivariate machine learning techniques provide an alternative to the rapidity gap method for event-by-event identification and classification of diffraction in hadron-hadron collisions. Traditionally, such methods assign each event exclusively to a single class producing classification errors in overlap regions of data space. As an alternative to this so called hard classification approach, we propose estimating posterior probabilities of each diffractive class and using these estimates to weigh event contributions to physical observables. It is shown with a Monte Carlo study that such a soft classification scheme is able to reproduce observables such as multiplicity distributions and relative event rates with a much higher accuracy than hard classification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental longitudinal emittance manipulation using laser-based photoionization in the Fermilab Linac

A series of simulations and beam studies were conducted at Fermilab’s linear accelerator to evaluate the effectiveness of longitudinal emittance control via laser-induced photoionization. While similar laser techniques have been employed at Fermilab to enhance injection and extraction efficiency into the Booster, the work presented here focuses on extending these methods to bunch-by-bunch manipulation. This approach utilizes fine-scale correction of the H- bunches’ longitudinal spatial distribution. In theory, loosely confined particles in longitudinal phase space contribute to emittance growth during acceleration. By selectively removing these outlying particles through laser scraping (H- + γ → H + e-), this growth can be reduced. This report presents experimental results from both symmetric and asymmetric longitudinal scraping of H⁻ bunches in the Fermilab linac, which were subsequently injected into Booster, and evaluates the broader applicability of this method for future high-intensity accelerator operations.

Landon, Parker [Boston U.; Fermilab]↗

Fatigue-test acceleration with flight-by-flight loading and heating to simulate supersonic-transport operation

Possibilities for reducing fatigue-test time for supersonic-transport materials and structures were studied in tests with simulated flight-by-flight loading. In order to determine whether short-time tests were feasible, the results of accelerated tests (2 sec per flight) were compared with the results of real-time tests (96 min per flight). The effects of design mean stress, the stress range for ground-air-ground cycles, simulated thermal stress, the number of stress cycles in each flight, and salt corrosion were studied. The flight-by-flight stress sequences were applied to notched sheet specimens of Ti-8Al-1Mo-1V and Ti-6Al-4V titanium alloys. A linear cumulative-damage analysis accounted for large changes in stress range of the simulated flights but did not account for the differences between real-time and accelerated tests. The fatigue lives from accelerated tests were generally within a factor of two of the lives from real-time tests; thus, within the scope of the investigation, accelerated testing seems feasible.

Imig, L. A.↗

Design, simulation and conditioning of the fundamental power couplers for BNL SRF gun

The 704 MHz SRF gun for the BNL Energy Recovery Linac (ERL) prototype uses two fundamental power couplers (FPCs) to deliver up to 1 MW of CW RF power to the half-cell cavity. To prepare the couplers for high-power RF service and process multipacting, the FPCs should be conditioned prior to installation into the gun cryomodule. A room-temperature test stand was configured for conditioning FPCs in full reflection regime with varied phase of the reflecting wave. The FPCs have been conditioned up to 250 kW in pulse mode and 125 kW in CW mode. The multipacting simulations were carried out with Track3P code developed at SLAC. The simulations matched the experimental results very well. This paper presents the FPC RF and thermal design, multipacting simulations and conditioning of the BNL gun FPCs.

43 PARTICLE ACCELERATORS↗

Simulation of Electron Beam Dynamics in the 22 MeV Accelerator for a Coherent Electron Cooling Proof of Principle Experiment

Coherent electron cooling (CeC) offers a potential new method of cooling hadron beams in colliders such as the Relativistic Heavy Ion Collider (RHIC) or the future electron ion collider eRHIC. A 22 MeV linear accelerator is currently being built as part of a proof of principle experiment for CeC at Brookhaven National Laboratory (BNL). In this thesis we present a simulation of electron beam dynamics including space charge in the 22 MeV CeC proof of principle experiment using the program ASTRA (A Space charge TRacking Algorithm).

43 PARTICLE ACCELERATORS↗