Search NASA⌕ Search

SEARCH · Search NASA

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,657 records · Page 92

Efficient Parallelization of Irregular Applications on GPU Architectures

With the enlarging computation capacity of general Graphics Processing Units (GPUs), leveraging GPUs to accelerate parallel applications has become a critical topic in academia and industry. However, a wide range of irregular applications with the computation-/memory-intensive nature cannot easily achieve high GPU utilization. The challenges mainly involve the following aspects: first, data dependence leads to coarse-grained kernel and inefficient parallelism; second, heavy GPU memory usage may cause frequent memory evictions and extra overhead of I/O; third, specific computation patterns produce memory redundancies; last, workload balance and data reusability conjunctly benefit the overall performance, but there may exist a dynamic trade-off between them. Targeting these challenges, this dissertation proposes multiple optimizations to accelerate two real-world applications: many-body correlation functions to simulate nuclear physics in a large-scale scientific system; the other is the eALS-based matrix factorization recommendation system. To accelerate the calculations of many-body correlation functions, this dissertation presents three frameworks in GPU memory management and multi-GPU scheduling. Firstly, an optimized systematic GPU memory management framework, MemHC, utilizes a series of new memory reduction designs in GPU memory allocation, CPU/GPU communications, and GPU memory oversubscription. Secondly, an enhanced multi-GPU scheduling framework, MICCO, particularly by taking both data dimension (e.g., data reuse and data eviction) and computation dimension into account. MICCO designs a heuristic scheduling algorithm and a machine learning-based regression model to generate the optimal settings of a proposed new concept to manage the trade-off. Thirdly, a locality-aware multi-GPU scheduling framework. This scheduler leverages pipeline batch generation with a looking-ahead strategy by building local dependency graphs for memory transfer reduction and better data reuse, achieving up to 79.92% memory cost reduction and 1.67x speedup. To parallelize the eALS-based recommendation system, this dissertation proposes an efficient CPU/GPU heterogeneous recommendation system, HEALS. HEALS employs newly designed architecture-adaptive data formats to achieve load balance and good data locality on CPU and GPU. To mitigate the data dependence, HEALS presents a CPU/GPU collaboration model for both task parallelism and data parallelism with multiple kernel computation optimizations. In summary, this dissertation efficiently accelerates two typical irregular applications on GPUs by building four frameworks, including CPU/GPU collaboration, GPU memory management, and multi-GPU scheduling.

Wang, Qihan↗

Note: The design of thin gap chamber simulation signal source based on field programmable gate array

The Thin Gap Chamber (TGC) is an important part of ATLAS detector and LHC accelerator. Targeting the feature of the output signal of TGC detector, we have designed a simulation signal source. The core of the design is based on field programmable gate array, randomly outputting 256-channel simulation signals. The signal is generated by true random number generator. The source of randomness originates from the timing jitter in ring oscillators. The experimental results show that the random number is uniform in histogram, and the whole system has high reliability.

97 MATHEMATICS AND COMPUTING↗

eSHIVA in IOTA: Experiment on Stochasticity in High Intensity Variable Accelerators

Nonlinear integrable optics, the basis of IOTA, is a very special form of Hamiltonian system, called completely integrable. The consequences are many, most importantly the absence of fixed points that are in one-to-one correspondence with (overlapping) resonances. This, in turn, does not limit the dynamic aperture and allows large amplitude dependent tune shifts without adverse effects. In fact, it helps significantly in alleviating issues with collective instabilities. We developed and tested different numerical methods to give a more complete picture of the dynamics of integrable and near-integrable nonlinear Hamiltonian systems consisting of motion of charged particles in electromagnetic fields. Comparison of the results allowed us to discern real phenomena from numerical artifacts and come up with conclusions that will facilitate pushing the intensity frontier in beam dynamics even further for the benefit of the accelerator community as a whole and society as the consequence of the science and technological advancements it will make possible. We concluded that, while novel structure preserving integrators would be ideal, symplectic methods under mild restrictions offer adequate invariant preservation. We especially highlighted the very good performance of the Stormer-Verlet method at low accuracy and a modified Lobatto method at high accuracy for the Integrable Optics Test Accelerator (IOTA). However, we showed that symplectic integrators are limited by numerically induced resonances, which make interpretation of some simulation results challenging. These limitations and workarounds have been addressed.The nonlinear integrable optics test accelerator IOTA at Fermilab -- a very special experimental accelerator designed to investigate nonlinear optics under severe space-charge conditions – has recently completed its construction and initial commissioning of the ring has begun. We report first, under 2 (a) below: the progress on the nonlinear tune-shift measurements for various values of the nonlinear strength potential and their comparison with theoretical expectations, thus benchmarking the nonlinear characteristics of the ring; then, under 2 (b) below: the progress on simulation and design of a gas-jet monitor to measure beam profile and in particular its diffusive halo, expected to develop during stochastic diffusion of large amplitude particles under nonlinear forces of the magnetic lattice and Coulomb space-charge. Since the development of a proton beam source in IOTA is delayed and to be completed only in the future pending financial support, the gas-jet monitor is planned for testing with real proton beams elsewhere at Fermilab (e.g. in the Main Injector or PIP-II source possibly). Once proton beams are available in IOTA, the understanding of the nonlinear properties of IOTA combined with the developed gas-jet monitor, will enable a thorough experimental analysis of nonlinear stochastic diffusion in IOTA nonlinear ring, and boosted by proper theoretical understanding of the same.

43 PARTICLE ACCELERATORS↗

Integrable Optics Design Principles for Beam Halo Suppression in Accelerator Rings at the Intensity Frontier

The project was dedicated to quantifying the possibilities and limitations of the integrable optics, a fundamentally new approach to high-intensity accelerator ring design, using the parallel Synergia framework for simulations. The following topics were studied: a) the effects of space charge forces, field errors, magnet misalignments, and other sources of parametric resonance; b) the effect of controlled nonlinearities, such as sextupoles for chromaticity control; c) the effect of longitudinal dynamics driven by finite bunch length and rf cavities; and d) reexamination of relevant lattices in present accelerators to explore potential benefits of these new design principles.

43 PARTICLE ACCELERATORS↗

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗

Anderson Acceleration for Fixed-Point Iterations

The purpose of this grant was to support research on acceleration methods for fixed-point iterations, with applications to computational frameworks and simulation problems that are of interest to DOE.

43 PARTICLE ACCELERATORS↗

INTERSTELLAR PICKUP ION ACCELERATION IN THE TURBULENT MAGNETIC FIELD AT THE SOLAR WIND TERMINATION SHOCK USING A FOCUSED TRANSPORT APPROACH

We study the physics of locally born interstellar pickup proton acceleration at the nearly perpendicular solar wind termination shock (SWTS) in the presence of a random magnetic field spiral angle using a focused transport model. Guided by Voyager 2 observations, the spiral angle is modeled with a q -Gaussian distribution. The spiral angle fluctuations, which are used to generate the perpendicular diffusion of pickup protons across the SWTS, play a key role in enabling efficient injection and rapid diffusive shock acceleration (DSA) when these particles follow field lines. Our simulations suggest that variation of both the shape ( q -value) and the standard deviation ( σ -value) of the q -Gaussian distribution significantly affect the injection speed, pitch-angle anisotropy, radial distribution, and the efficiency of the DSA of pickup protons at the SWTS. For example, increasing q and especially reducing σ enhances the DSA rate.

79 ASTRONOMY AND ASTROPHYSICS↗

New constraints on warm dark matter from the Lyman- α forest power spectrum

The forest of Lyman-α absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable—the 1D flux power spectrum (FPS)—should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard Λ⁢CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the graphics processing unit-accelerated code cholla to study the evolution of the Lyman-α forest under a wide range of physically motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts 4.0≲z≲5.2 [E. Boera et al., Revealing reionization with the thermal history of the intergalactic medium: New constraints from the Ly⁢α flux power spectrum, Astrophys. J. 872, 101 (2019)] yields a lower-limit m WDM >3.1 keV (95% C.L.) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of m WDM =4.5 keV. In conclusion, we find that the suppression of the FPS from free-streaming saturates at k≳0.1 s km -1 because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of a Laser-Electron-Interaction Region for Compact High Repetition-Rate Mono-Energetic Laser Based Gamma Sources

Inverse-Compton scattering (ICS) is a well-established technique capable of producing monochromatic radiation over a large range of the electromagnetic spectrum. Most recently, ICS has attracted attention, owing to its compactness and possible coupling with novel acceleration methods to yield compact X- and 𝛄-ray sources. Under this award, we carried out a systematic optimization of the ICS and demonstrate via numerical simulation that combining ICS with a conventional superconducting linac could produce Watt-level average-power 𝛄 rays with high spectral brightness. The study is performed with idealized laser and electron-beam distributions and confirmed with the realistic distributions simulated from start-to-end numerical modeling of the electron accelerator. The studies led to the design and engineering design of an ICS interaction point for the FAST facility along with some beam measurement. Due to the cancellation of the program in Year 3 long with delay in the commissioning of the FAST facility, we were not able to demonstrate the generation of 𝛄 rays, we, however, procured some of the hardware, installed some of the required infrastructures and plan to eventually generate 𝛄 ray pending future interest on the topics and availability of funding.

43 PARTICLE ACCELERATORS↗

Ocean and Sea Ice and their Interactions around Greenland and the West Antarctic Peninsula in Forced Fine-Resolution Global Simulations

Anomalously warm ocean waters have been implicated as the cause of accelerated ablation along the margin of the Greenland ice sheet, as well as the cause of increasing basal melt rate, mass loss, and grounding line retreat of many Antarctic ice shelves. The overarching science objective of this project is to investigate how the interplay of regional processes and decadal changes in local and remote forcing impact the delivery and end member composition of waters over the continental shelves of Greenland and Antarctica. To address this goal, an existing global eddy-active ocean general circulation model (OGCM), coupled to a thermodynamic/dynamic ice model and forced with atmospheric reanalysis fluxes from 1948-2009, was analyzed to understand these processes around the continental margins of Greenland and Antarctica. This simulation did not contain a representation of ice sheet melt or resolve mesoscale eddies over these continental margins, so two new coupled OGCM/sea-ice models were constructed and analyzed to address the roles of small scale features and processes (< 5km) and freshwater flux surpluses.

54 ENVIRONMENTAL SCIENCES↗

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗

A Faster-Than-Real-Time Framework for Reliability-Oriented Simulation of PV Inverters

Physics-of-Failure (PoF) based reliability assessment for photovoltaic (PV) inverters requires long-duration electrical and electrothermal stress histories, yet generating such stress histories with high-fidelity switching models over year long mission profiles is computationally prohibitive. Conventional methods either sacrifice modeling fidelity for speed or require runtimes that are impractical for design iteration and uncertainty studies. To address this bottleneck, this paper presents a High-Performance Computing (HPC) based simulation frame work for faster-than-real-time reliability-oriented simulation. The proposed framework integrates the Average-to-Switching (A2S) method with parallel computing techniques to accelerate switching-level waveform reconstruction. We further introduce optimization strategies, including cluster merging and sensitivity based mission profile screening, to reduce the computational burden. Evaluated using real-world mission profile inputs and a MATLAB/Simulink switching-model reference, the framework reduces the simulation time for a one-year mission from an intractable multi-year duration to approximately 7.3 minutes while maintaining low waveform error. This acceleration provides a practical reliability-oriented simulation engine that can be coupled with component-specific aging models for subsequent PV inverter PoF assessment.

High-performance Computing↗

Polarization simulations in the RHIC run 15 lattice

RHIC polarized proton Run 15 uses a new acceleration ramp optics, compared to RHIC Run 13 and earlier runs, in relation with electron-lens beam-beam compensation developments. The new optics induces different strengths in the depolarizing snake resonance sequence, from injection to top energy. As a consequence, polarization transport along the new ramp has been investigated, based on spin tracking simulations. Sample results are reported and discussed.

43 PARTICLE ACCELERATORS↗

Application of recoil-imaging time projection chambers to directional neutron background measurements in the SuperKEKB accelerator tunnel

Gaseous time projection chambers (TPCs) with high readout segmentation are capable of reconstructing detailed 3D ionization distributions of nuclear recoils resulting from elastic neutron scattering. Here, using a system of six compact TPCs with pixel ASIC readout, filled with a 70:30 mixture of He:CO 2 gas, we analyze the first directional measurements of beam-induced neutron backgrounds in the tunnel regions surrounding the Belle II detector at the SuperKEKB e + e – collider. With the use of 3D recoil tracking, we show that these TPCs are capable of maintaining nearly 100% nuclear recoil purity to reconstructed ionization energies (E reco ) as low as 5 keV ee . Using a large sample of Monte-Carlo (MC)-simulated 4 He, 12 C, and 16 O recoil tracks, we find consistency between predicted and measured recoil energy spectra in five of the six TPCs, providing useful validation of the neutron production mechanisms modeled in simulation. Restricting this sample to 4 He recoil tracks with E reco > 40 keV ee , we further demonstrate axial angular resolutions within 8° and we introduce a procedure that under suitable conditions, correctly assigns the vector direction to 91% of these simulated 4 He recoils. Applying this procedure to assign vector directions to measured 4 He recoil tracks, we observe consistency between the angular distributions of observed and simulated recoils, providing first experimental evidence of localized neutron “hotspots” in the accelerator tunnel. Observed rates of nuclear recoils in these TPCs suggest that simulation overestimates the neutron flux from these hotspots. Despite this, we estimate these hotspots to produce the majority of neutron backgrounds in the accelerator tunnel at SuperKEKB’s target luminosity of 6.3 x 10 35 cm –2 s –1 , making them important regions to continue to monitor.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-F-T-496: An Investigation of Two Novel Devices for Testing Linac Clearance During CT Simulation

Purpose: This work’s objective is to determine the efficacy of two newly patented devices termed the Mor-O rings, Mark1 and Mark2, developed to predict collisions between the patient and a linac at the time of CT simulation. Methods: Mark1: A ring with an inner radius equal to the distance between the isocenter and the nearest portion of a linac head(diso) was made. This is mounted to a stand that allows vertical repositioning. The ring is placed around the patient on the CT table and aligned with isocenter. The patient is moved through the ring. If the ring touches the patient, a collision is predicted. To test this device, predicted collisions were marked on a phantom. The phantom was then repositioned on the linac table where the collisions were verified. Mark2: An arc with the radius diso was created with a re-locatable half-linac head wire-frame. The Mark2 is positioned in the same way as the Mark1 but can additionally mimic couch and gantry angulations. The Mark2 was tested with a volunteer using multiple couch, gantry and isocenter positions. The volunteer was then repositioned on the linac table to verify the angles of collisions. Results: Mark1: One isocenter out of ten showed negative clearance (0.9mm) on the linac table. All other collisions were predicted with 1 to 36mm of additional clearance. Mark2: All collisions were prevented with an additional 1.3 to 14.8 degrees of clearance. Conclusion: The Mark1 prototype is able to predict all collisions when no couch angulation is used. This device takes less than a minute to setup and is simple to use. The Mark2, when testing beam geometries used for noncoplanar SBRT, was able to prevent all collisions with 1.3 to 14.8 degrees additional clearance. Improvements in construction for both devices could increase accuracy and usability. Andrew Morrow owns Morrow Physics, LLC and Brian Massingill owns Spur Physics, LLC. We are both listed on the patent for the devices investigated in this work (patent 9220922). Our current employer, Baylor Scott and White, is also on the patent.

60 APPLIED LIFE SCIENCES↗

TH-C-18A-11: Investigating the Minimum Scan Parameters Required to Generate Free-Breathing Fast-Helical CT Scans Without Motion-Artifacts

Purpose: A recently proposed 4D-CT protocol uses deformable registration of free-breathing fast-helical CT scans to generate a breathing motion model. In order to allow accurate registration, free-breathing images are required to be free of doubling-artifacts, which arise when tissue motion is greater than scan speed. This work identifies the minimum scanner parameters required to successfully generate free-breathing fast-helical scans without doubling-artifacts. Methods: 10 patients were imaged under free breathing conditions 25 times in alternating directions with a 64-slice CT scanner using a low dose fast helical protocol. A high temporal resolution (0.1s) 4D-CT was generated using a patient specific motion model and patient breathing waveforms, and used as the input for a scanner simulation. Forward projections were calculated using helical cone-beam geometry (800 projections per rotation) and a GPU accelerated reconstruction algorithm was implemented. Various CT scanner detector widths and rotation times were simulated, and verified using a motion phantom. Doubling-artifacts were quantified in patient images using structural similarity maps to determine the similarity between axial slices. Results: Increasing amounts of doubling-artifacts were observed with increasing rotation times > 0.2s for 16×1mm slice scan geometry. No significant increase in doubling artifacts was observed for 64×1mm slice scan geometry up to 1.0s rotation time although blurring artifacts were observed >0.6s. Using a 16×1mm slice scan geometry, a rotation time of less than 0.3s (53mm/s scan speed) would be required to produce images of similar quality to a 64×1mm slice scan geometry. Conclusion: The current generation of 16 slice CT scanners, which are present in most Radiation Oncology departments, are not capable of generating free-breathing sorting-artifact-free images in the majority of patients. The next generation of CT scanners should be capable of at least 53mm/s scan speed in order to use a fast-helical 4D-CT protocol to generate a motion-artifact free 4D-CT. NIH R01CA096679.

60 APPLIED LIFE SCIENCES↗

Speculation and replication in temperature accelerated dynamics

Accelerated Molecular Dynamics (AMD) is a class of MD-based algorithms for the long-time scale simulation of atomistic systems that are characterized by rare-event transitions. Temperature-Accelerated Dynamics (TAD), a traditional AMD approach, hastens state-to-state transitions by performing MD at an elevated temperature. Recently, Speculatively-Parallel TAD (SpecTAD) was introduced, allowing the TAD procedure to exploit parallel computing systems by concurrently executing in a dynamically generated list of speculative future states. Although speculation can be very powerful, it is not always the most efficient use of parallel resources. In this paper, we compare the performance of speculative parallelism with a replica-based technique, similar to the Parallel Replica Dynamics method. A hybrid SpecTAD approach is also presented, in which each speculation process is further accelerated by a local set of replicas. Finally and overall, this work motivates the use of hybrid parallelism whenever possible, as some combination of speculation and replication is typically most efficient.

36 MATERIALS SCIENCE↗

Emittance preservation for the electron arm in a single PWFA-LC stage using quasi-adiabatic plasma density ramp matching sections

Plasma-based acceleration (PBA) is being considered for a next generation linear collider (LC). In some PBA-LC designs for the electron arm, the extreme beam parameters are expected to trigger background ion motion within the witness beam, which can lead to longitudinally varying nonlinear focusing forces and result in an unacceptable emittance growth of the beam. To mitigate this, we propose to use quasi-adiabatic plasma density ramps as matching sections at the entrance and exit of each stage. We match the witness electron beam to the low density plasma entrance, where the beam initially has a large matched spot size so the ion motion effects are relatively small. As the beam propagates in the plasma density upramp, it is quasi-adiabatically focused, and its distribution maintains a non-Gaussian equilibrium distribution in each longitudinal slice throughout the process, even when severe ion collapse has occurred. This only causes small amounts of slice emittance growth. The phase mixing between slices with different betatron frequencies leads to additional projected emittance growth within the acceleration stage. A density downramp at the exit of an acceleration section can eliminate much of the slice and projected emittance growth as the beam and ion motion adiabatically defocuses and decreases, respectively. Simulation results from QuickPIC with Azimuthal Decomposition show that within a single acceleration stage with a 25 GeV energy gain, this concept can limit the projected emittance growth to only ∼2% for a 25 GeV, 100 nm emittance witness beam and ∼20% for a 100 GeV, 100 nm normalized emittance witness beam. The trade-off between the adiabaticity of the plasma density ramp and the initial ion motion at the entrance for a given length of the plasma density ramp is also discussed.

Physics↗