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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

UCLA Final Technical Report for the "Community Petascale Project for Accelerator Science and Simulation”.

The UCLA Plasma Simulation Group is a major partner of the “Community Petascale Project for Accelerator Science and Simulation”. This is the final technical report. We include an overall summary, a list of publications, progress for the most recent year, and individual progress reports for each year. We have made tremendous progress during the three years. SciDAC funds have contributed to the development of a large number of skeleton codes that illustrate how to write PIC codes with a hierarchy of parallelism. These codes cover 2D and 3D as well as electrostatic solvers (which are used in beam dynamics codes and quasi-static codes) and electromagnetic solvers (which are used in plasma based accelerator codes). We also used these ideas to develop a GPU enabled version of OSIRIS. SciDAC funds were also contributed to the development of strategies to eliminate the Numerical Cerenkov Instability (NCI) which is an issue when carrying laser wakefield accelerator (LWFA) simulations in a boosted frame and when quantifying the emittance and energy spread of self-injected electron beams. This work included the development of a new code called UPIC-EMMA which is an FFT based electromagnetic PIC code and to new hybrid algorithms in OSIRIS. A new hybrid (PIC in r-z and gridless in φ) algorithm was implemented into OSIRIS. In this algorithm the fields and current are expanded into azimuthal harmonics and the complex amplitude for each harmonic is calculated separately. The contributions from each harmonic are summed and then used to push the particles. This algorithm permits modeling plasma based acceleration with some 3D effects but with the computational load of an 2D r-z PIC code. We developed a rigorously charge conserving current deposit for this algorithm. Very recently, we made progress in combining the speed up from the quasi-3D algorithm with that from the Lorentz boosted frame. SciDAC funds also contributed to the improvement and speed up of the quasi-static PIC code QuickPIC. We have also used our suite of PIC codes to make scientific discovery. Highlights include supporting FACET experiments which achieved the milestones of showing high beam loading and energy transfer efficiency from a drive electron beam to a witness electron beam and the discovery of a self-loading regime a for high gradient acceleration of a positron beam. Both of these experimental milestones were published in Nature together with supporting QuickPIC simulation results. Simulation results from QuickPIC were used on the cover of Nature in one case. We are also making progress on using highly resolved QuickPIC simulations to show that ion motion may not lead to catastrophic emittance growth for tightly focused electron bunches loaded into nonlinear wakefields. This could mean that fully self-consistent beam loading scenarios are possible. This work remains in progress. OSIRIS simulations were used to discover how 200 MeV electron rings are formed in LWFA experiments, on how to generate electrons that have a series of bunches on nanometer scale, and how to transport electron beams from (into) plasma sections into (from) conventional beam optic sections.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

SU-E-T-36: A GPU-Accelerated Monte-Carlo Dose Calculation Platform and Its Application Toward Validating a ViewRay Beam Model

Purpose: To build a fast, accurate and easily-deployable research platform for Monte-Carlo dose calculations. We port the dose calculation engine PENELOPE to C++, and accelerate calculations using GPU acceleration. Simulations of a Co-60 beam model provided by ViewRay demonstrate the capabilities of the platform. Methods: We built software that incorporates a beam model interface, CT-phantom model, GPU-accelerated PENELOPE engine, and GUI front-end. We rewrote the PENELOPE kernel in C++ (from Fortran) and accelerated the code on a GPU. We seamlessly integrated a Co-60 beam model (obtained from ViewRay) into our platform. Simulations of various field sizes and SSDs using a homogeneous water phantom generated PDDs, dose profiles, and output factors that were compared to experiment data. Results: With GPU acceleration using a dated graphics card (Nvidia Tesla C2050), a highly accurate simulation – including 100*100*100 grid, 3×3×3 mm3 voxels, <1% uncertainty, and 4.2×4.2 cm2 field size – runs 24 times faster (20 minutes versus 8 hours) than when parallelizing on 8 threads across a new CPU (Intel i7-4770). Simulated PDDs, profiles and output ratios for the commercial system agree well with experiment data measured using radiographic film or ionization chamber. Based on our analysis, this beam model is precise enough for general applications. Conclusions: Using a beam model for a Co-60 system provided by ViewRay, we evaluate a dose calculation platform that we developed. Comparison to measurements demonstrates the promise of our software for use as a research platform for dose calculations, with applications including quality assurance and treatment plan verification.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations

Small-angle X-ray scattering (SAXS) has become an increasingly popular technique for characterizing the solution ensemble of flexible biomolecules. However, data resulting from SAXS is typically low-dimensional and is therefore difficult to interpret without additional structural knowledge. In theory, molecular dynamics (MD) trajectories can provide this information, but conventional simulations rarely sample the complete ensemble. Here, we demonstrate that accelerated MD simulations can be used to produce higher quality models in shorter time scales than standard simulations, and we present an iterative Bayesian Monte Carlo method that is able to identify multistate ensembles without overfitting. This methodology is applied to several ubiquitin trimers to demonstrate the effect of linkage type on the solution states of the signaling protein. Here, we observe that the linkage site directly affects the solution flexibility of the trimer and theorize that this difference in plasticity contributes to their disparate roles in vivo.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Measurements of Decay and Snapback in Nb₃Sn Accelerator Magnets at Fermilab

In recent years, Fermilab has been executing an inten-sive R&D program on Nb₃Sn accelerator magnets. This program has included dipole and quadrupole models and demonstrators for various programs and projects, including the HL-LHC accelerator upgrade project. A systematic study of the field decay and snapback during the injection portion of a simulated accelerator cycle was executed at the Fermilab Magnet Test Facility. This paper summarizes the recent measurements of the MQXFS1 short quadrupole model and discusses the results of some previously measured Nb₃Sn magnets at CERN

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Accelerating gravitational microlensing simulations using the Xeon Phi coprocessor

Recently Graphics Processing Units (GPUs) have been used to speed up very CPU-intensive gravitational microlensing simulations. In this work, we use the Xeon Phi coprocessor to accelerate such simulations and compare its performance on a microlensing code with that of NVIDIA's GPUs. For the selected set of parameters evaluated in our experiment, we find that the speedup by Intel's Knights Corner coprocessor is comparable to that by NVIDIA's Fermi family of GPUs with compute capability 2.0, but less significant than GPUs with higher compute capabilities such as the Kepler. However, the very recently released second generation Xeon Phi, Knights Landing, is about 5.8 times faster than the Knights Corner, and about 2.9 times faster than the Kepler GPU used in our simulations. Here, we conclude that the Xeon Phi is a very promising alternative to GPUs for modern high performance microlensing simulations. LESS

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Performance Understanding and Analysis for Exascale Data Management Workflows

The goal of the Exascale Data Management Workflows project is to create tools to support techniques for in-situ data analysis and data management. Online data management co-running with simulations accelerates the scientific processes being carried out, provides rapid and timely scientific insights, and can help avoid unnecessary and scientifically invalid simulation computations, particularly when combined with or able to utilize data from experimental instrument and/or from past simulation runs. The University of Oregon team is contributing to the performance instrumentation and analysis aspect of this project. Our focus has been on efficient performance monitoring of data management at scale. The general approach of the MONA project is depicted in Figure 1 below. The figure shows performance monitoring applied to an I/O workflow associated with a scientific simulation, with online measurements captured from this workflow informing methods for workflow reconfiguration and dynamic adjustment. The goal is to maintain suitable levels of Quality of Service for workflow execution, by understanding the under lying causes of workflow performance and behavior. One outcome is workflow performance models able to characterize realistic workflows. Another outcome is performance ’mini-apps’ implementing such workflow behavior. Out of scope for this project are advanced methods for online workflow control.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

PowerCell aboard the EuCROPIS satellite - Results from the first synthetic biology experiment in space

Synthetic biology is the design and construction of new biological parts and systems not found in nature, a technology that we believe will be a game-changing technology for space exploration. The PowerCell concept is an ecology based on photosynthetic microbes taking advantage of in situ materials and energy to generate, on-demand, useful products (food, fuel, cloths, drugs, etc.) that satisfy specific needs of long-term human presence away from Earth. This first PowerCell experiment is testing this concept as well as several key mechanisms of engineering biology in a 48-well plate microfluidics system. NASA's PowerCell experiment was loaded in the fluidics system at NASA Ames in May 2016, and subsequently integrated on the DLR's compact satellite Eu:CROPIS (Euglena Combined Regenerative Organic food Production In Space), as a secondary payload. Launched into Earth’s orbit from the Vandenberg Air Force Base, California, USA, in December 3, 2019, the Eu:CROPIS platform has provided an artificial gravity through rotational acceleration, simulating gravity ranges compatible to those found on Mars and on the Moon, as well as microgravity in the absence of acceleration. Each set of biological experiments performed at Earth's orbit was compared to ground controls. Over 400 days' worth of mission data showed that all hardware and software worked as planned, and microbial growth was observed in all four gravity regimes, including growth on cyanobacterial cell extract in three out of four gravity regimes. The genetic transformation experiment was successful in microgravity and provided valuable long-term reagent stability data, including antibiotic efficacy. The demonstrated ability of hardy bacterial spores to recover from a long-term stasis period will be fundamental for deep human space missions. Lessons learned from hardware design to experimental operations are leading to improved designs and capabilities for future missions.

Ivan Glaucio Paulino Lima↗
DOE OSTI2016

Stochastic averaging and sensitivity analysis for two scale reaction networks

In the presence of multiscale dynamics in a reaction network, direct simulation methods become inefficient as they can only advance the system on the smallest scale. This work presents stochastic averaging techniques to accelerate computations for obtaining estimates of expected values and sensitivities with respect to the steady state distribution. A two-time-scale formulation is used to establish bounds on the bias induced by the averaging method. Further, this formulation provides a framework to create an accelerated “averaged” version of most single-scale sensitivity estimation methods. In particular, we propose the use of a centered ergodic likelihood ratio method for steady state estimation and show how one can adapt it to accelerated simulations of multiscale systems. Lastly, we develop an adaptive “batch-means” stopping rule for determining when to terminate the micro-equilibration process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Coupling reduced‐order modeling and coarse‐grained CFD‐DEM to accelerate coal gasifier simulation and optimization

Abstract High fidelity three‐dimensional (3‐D) numerical simulations of commercial‐scale coal gasifiers are very time consuming and expensive. This article proposes a reduced‐order modeling approach that uses quasi one‐dimensional (1‐D) CFD–DEM simulation results to serve as an accurate initial condition for the 3‐D simulation, which significantly shortens the physical time from several hours to minutes to achieve steady states. The 3‐D simulation employs a previously validated coarse‐grained CFD‐DEM approach, which further accelerates the simulation. Both 1‐D and 3‐D simulations adopt a particle shrinkage scheme to correctly account for the particle volume change due to chemical reaction and thus provides an accurate bed height. Comparison with available experimental measurements of the bed temperature and synthesis gas composition is used to validate the simulation. The final syngas composition and flowrates are strongly affected by the fuel or gasification agent rate and process operating conditions as expected. This approach is then used to optimize two important operating conditions (air and steam flows) to achieve specific gasifier performance goals for the production of Fischer–Tropsch synthesis gas and production of fuel gas for gas turbines.

Yu, Jia↗
DOE OSTI2020

Subsurface Nitrogen Dissociation Kinetics in Lithium Metal from Metadynamics

The dissociation of molecular nitrogen in lithium is of interest for several promising technologies, such as the catalytic synthesis of ammonia in ambient or mild conditions. In this work we simulate nitrogen dissociation in the lithium BCC (110) surface at ambient and elevated temperatures using density functional theory (DFT) metadynamics simulations. The rate constants at temperatures of 300, 400, and 500 K are calculated by statistical analysis of the reaction time distributions from the accelerated simulations. This approach finds and estimates rate constants for transition pathways out of the initial state; the required input is the stable initial state and a reasonable choice of collective variable. A single collective variable is used in this case: the N–N distance. The results are robust to changes in metadynamics parameters, and the reaction time distributions follow the expected exponential distribution. We show that the metadynamics-derived rate constants are in agreement with results from the conventional harmonic approximation approach using a climbing image nudged elastic band (NEB) transition state search. The reaction barriers from metadynamics and the NEB/harmonic approximation agree to within 0.02–0.04 eV at all temperatures studied. This paper demonstrates that the harmonic approximation provides an accurate description of the rate constants for nitrogen dissociation in lithium metal, even at temperatures near or above the melting point of lithium, lending credence to previous and future theoretical studies using this approximation. Moreover, this work demonstrates a step toward the automated exploration and discovery of reaction mechanisms and associated rate constants for elementary surface-catalyzed reactions using DFT-based metadynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Derivative-free optimization of a rapid-cycling synchrotron

Here, we develop and solve a constrained optimization model to identify an integrable optics rapid-cycling synchrotron lattice design that performs well in several capacities. Our model encodes the design criteria into 78 linear and nonlinear constraints, as well as a single nonsmooth objective, where the objective and some constraints are defined from the output of Synergia, an accelerator simulator. We detail the difficulties of the 23-dimensional simulation-constrained decision space and establish that the space is nonempty. We use a derivative-free manifold sampling algorithm to account for structured nondifferentiability in the objective function. Our numerical results quantify the dependence of solutions on constraint parameters and the effect of the form of objective function.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Technical Report for Bayesian Optimization and Reinforcement Learning for Beam Polarization Increase in the BNL Hadron Injectors

This project developed and evaluated physics-informed Bayesian learning and machine learning (ML)-based optimization methods for improving beam polarization preservation in the BNL hadron injector chain. The work focused on uncertainty-aware digital twin modeling, Bayesian calibration of accelerator simulations using beam measurements, and data-efficient optimization strategies including Bayesian optimization and reinforcement learning. These methods were applied to injector tuning and RF control problems in realistic accelerator settings to support improved operational robustness and readiness for RHIC operations and future Electron–Ion Collider facilities. No subject inventions were disclosed under this award.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

High Performance Computing Modeling Advances Accelerator Science for High-Energy Physics

The development and optimization of particle accelerators are essential for advancing our understanding of the properties of matter, energy, space, and time. Particle accelerators are complex devices whose behavior involves many physical effects on multiple scales. Therefore, advanced computational tools utilizing high-performance computing are essential for accurately modeling them. In the past decade, the US Department of Energy's SciDAC program has produced accelerator-modeling tools that have been employed to tackle some of the most difficult accelerator science problems. The authors discuss the Synergia framework and its applications to high-intensity particle accelerator physics. Synergia is an accelerator simulation package capable of handling the entire spectrum of beam dynamics simulations. Our authors present Synergia's design principles and its performance on HPC platforms.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Curvature of the spectral energy distributions of blazars

In this paper, spectral energy distributions (SED) of both synchrotron and inverse Compton (IC) components of a sample of Fermi bright blazars are fitted by a log-parabolic law. The second-degree term of the log parabola measures the curvature of an SED. We find a statistically significant correlation between the synchrotron peak frequency and its curvature. This result is in agreement with the theoretical prediction and confirms previous studies that dealt with a single source with observations at various epochs or a small sample. If a broken power law is employed to fit the SED, the difference between the two spectral indices (i.e., |α{sub 2} – α{sub 1}|) can be considered a 'surrogate' of the SED curvature. We collect data from the literature and find a correlation between the synchrotron peak frequency and the spectral difference. We do not find a significant correlation between the IC peak frequency and its curvature, which may be caused by a complicated seed photon field. It is also found that the synchrotron curvatures are on average larger than those of IC curvatures, and there is no correlation between these two parameters. As suggested by previous works, both the log-parabolic law of the SED and the above correlation can be explained by statistical and/or stochastic particle accelerations. Based on a comparison of the slops of the correlation, our result seems to favor stochastic acceleration mechanisms and emission processes. Additional evidence, including SED modeling, particle acceleration simulation, and comparisons between some predictions and empirical relations/correlations, also seems to support the idea that the electron energy distribution (and/or synchrotron SED) may be log-parabolic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Environmental exposure to TiO 2 nanomaterials incorporated in building material

Nanomaterials are increasingly being used to improve the properties and functions of common building materials. A new type of self-cleaning cement incorporating TiO 2 nanomaterials (TiO 2 -NMs) with photocatalytic properties is now marketed. This promising cement might provide air pollution-reducing properties but its environmental impact must be validated. During cement use and aging, an altered surface layer is formed that exhibits increased porosity. The surface layer thickness alteration and porosity increase with the cement degradation rate. The hardened cement paste leaching behavior has been fully documented, but the fate of incorporated TiO 2 -NMs and their state during/after potential release is currently unknown. In this study, photocatalytic cement pastes with increasing initial porosity were leached at a lab-scale to produce a range of degradation rates concerning the altered layer porosity and thickness. No dissolved Ti was released during leaching, only particulate TiO 2 -NM release was detected. The extent of release from this batch test simulating accelerated worst-case scenario was limited and ranged from 18.7 ± 2.1 to 33.5 ± 5.1 mg of Ti/m 2 of cement after 168 h of leaching. TiO 2 -NMs released into neutral aquatic media (simulate pH of surface water) were not associated or coated by cement minerals. The TiO 2 -NM release mechanism is suspected to start from freeing of TiO 2 -NMs in the altered layer pore network due to partial cement paste dissolution followed by diffusion into the bulk pore solution to the surface. The extent of TiO 2 -NM release was not solely related to the cement degradation rate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Supercomputing Pipelines Search for Therapeutics Against COVID-19

The urgent search for drugs to combat SARS-CoV-2 has included the use of supercomputers. The use of general-purpose graphical processing units (GPUs), massive parallelism, and new software for high-performance computing (HPC) has allowed researchers to search the vast chemical space of potential drugs faster than ever before. We developed a new drug discovery pipeline using the Summit supercomputer at Oak Ridge National Laboratory to help pioneer this effort, with new platforms that incorporate GPU-accelerated simulation and allow for the virtual screening of billions of potential drug compounds in days compared to weeks or months for their ability to inhibit SARS-COV-2 proteins. Here, this effort will accelerate the process of developing drugs to combat the current COVID-19 pandemic and other diseases.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1995

Monte Carlo simulations of particle acceleration at oblique shocks: Including cross-field diffusion

The Monte Carlo technique of simulating diffusive particle acceleration at shocks has made spectral predictions that compare extremely well with particle distributions observed at the quasi-parallel region of the earth's bow shock. The current extension of this work to compare simulation predictions with particle spectra at oblique interplanetary shocks has required the inclusion of significant cross-field diffusion (strong scattering) in the simulation technique, since oblique shocks are intrinsically inefficient in the limit of weak scattering. In this paper, we present results from the method we have developed for the inclusion of cross-field diffusion in our simulations, namely model predictions of particle spectra downstream of oblique subluminal shocks. While the high-energy spectral index is independent of the shock obliquity and the strength of the scattering, the latter is observed to profoundly influence the efficiency of injection of cosmic rays into the acceleration process.

Baring, M. G.↗
DOE OSTI2019

DiaMonD: An Integrated Multifaceted Approach to Mathematics at the Interfaces of Data, Models, and Decisions

The major goal of the project was to invent, analyze, implement, and apply new algorithms to speed up the computational simulation of complex systems in settings in which multiple simulations are needed. The main efforts we undertook are as follows: 1. Develop novel reduced order methods to reduce the costs of simulations, especially in settings (e.g., UQ, inference, optimization and control) for which numerous simulation are needed and also for nonlocal problems (e.g., fractional derivative, anomalous diffusion, nonlocal mechanics) for which lack of sparsity results in more expensive full-order simulation than encountered in partial differential equations settings; rigorous analyses as well as computational tests were carried out, 2. Develop multifidelity methods based on the usage of a variety of models having different fidelities and different costs to accelerate simulations again in settings where numerous simulations are needed or are especially expensive; rigorous analyses as well as computational tests were carried out, 3. Develop ensemble methods to reduced the costs of simulations again in settings where numerous simulations are needed or are especially expensive; rigorous analyses as well as computational tests were carried out, and 4. Analyze a model for Richardson pair dispersion in turbulent flows and develop discretization methods for the efficient and accurate approximations of the solution of the model; rigorous analyses as well as computational tests were carried out.

97 MATHEMATICS AND COMPUTING↗