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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Parallel sorting algorithm classification: is manual instrumentation necessary?

Understanding parallel algorithms is crucial for accelerating scientific simulations on complex, distributed memory, high-performance computers. Modern algorithm classification approaches learn semantics directly from source code to differentiate between algorithms, however, accessing source code is not always possible. We can learn about parallel algorithms from observing their performance, as programs running the same algorithms and using the same hardware should exhibit similar performance characteristics. We present an approach to learn algorithm classes from parallel performance data directly in order to classify algorithms without access to the source code. We extend previous work to enable classifying parallel sorting algorithms using automatic instrumentation instead of requiring manual region annotations in the source code. In this work, we design and demonstrate a study for classification of parallel sorting algorithms using parallel performance data collected from automatic instrumentation, and evaluate the performance of our new methodology on classification. We leverage Caliper to collect the performance data, Thicket for our exploratory data analysis (EDA), and PyTorch and Scikit-learn to evaluate the effectiveness of random forests, support vector machines (SVMs), decision trees, neural networks, and logistic regressions on parallel performance data. Additionally, we study noise in parallel performance data, whether the removal of noise and pre-processing of the data is necessary to accurately classify parallel sorting algorithms, and determine the effectiveness of features created from performance data. In conclusion, we demonstrate classification accuracy for these five different models of up to 97.7% across four different parallel algorithm classes.

Algorithm Classification↗
DOE OSTI2018

Additive manufacturing of functionally graded transition joints between ferritic and austenitic alloys

Dissimilar metal joints between ferritic and austenitic alloys are susceptible to premature failure due to diffusive carbon loss from the ferritic alloy driven by abrupt changes in carbon chemical potential. Compositional grading of transition joints fabricated using laser-based directed energy deposition additive manufacturing offers a means for limiting carbon diffusion. In this report we fabricate functionally graded transition joints between a ferritic and austenitic alloy, characterize spatial variations of chemical composition, microstructure and hardness, and test their effectiveness to limit carbon loss from the ferritic alloy. Microstructural studies and carbon potential variations in the functionally graded material showed that the length of the joints can be shorter, and there is no benefit to continue compositional grading once the microstructure becomes fully austenitic. Since dissimilar joints have an expected lifetime of several decades, long service times were simulated through accelerated heat treatment experiments at elevated temperatures for both a dissimilar metal weld and a functionally graded transition joint. While the dissimilar weld showed pronounced carbon loss from the ferritic side, there was insignificant change in the carbon concentration profile in the functionally graded joint indicating effectiveness of the graded joints to perform under service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A fast and accurate physics-informed neural network reduced order model with shallow masked autoencoder

Traditional linear subspace reduced order models (LS-ROMs) are able to accelerate physical simulations in which the intrinsic solution space falls into a subspace with a small dimension, i.e., the solution space has a small Kolmogorov n-width. However, for physical phenomena not of this type, e.g., any advection-dominated flow phenomena such as in traffic flow, atmospheric flows, and air flow over vehicles, a low-dimensional linear subspace poorly approximates the solution. To address cases such as these, we have developed a fast and accurate physics-informed neural network ROM, namely nonlinear manifold ROM (NM-ROM), which can better approximate high-fidelity model solutions with a smaller latent space dimension than the LS-ROMs. Our method takes advantage of the existing numerical methods that are used to solve the corresponding full order models. The efficiency is achieved by developing a hyper-reduction technique in the context of the NM-ROM. Numerical results show that neural networks can learn a more efficient latent space representation on advection-dominated data from 1D and 2D Burgers' equations. A speedup of up to 2.6 for 1D Burgers' and a speedup of 11.7 for 2D Burgers' equations are achieved with an appropriate treatment of the nonlinear terms through a hyper-reduction technique. Lastly, a posteriori error bounds for the NM-ROMs are derived that take account of the hyper-reduced operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Irradiation of ultrasonic sensors and adhesive couplants for application in light water reactor primary loop piping and components

The Electric Power Research Institute (EPRI) Nuclear Sector and US Department of Energy Light Water Reactor Sustainability Program are committed to engaging in research and development endeavors to address materials aging issues specific to long term operation of light water power reactors. To this effect, EPRI launched an industry initiative to develop nondestructive evaluation systems for online monitoring of existing cracks in light water reactor primary coolant loop piping and components. One of the goals of this initiative is to develop a sensor system (or systems) that can determine nondestructively if cracks are growing or arrested and, in the case of the former, to characterize their growth rates. A missing component of this initiative is an experimental assessment of how sensors and adhesive couplants will perform in service when exposed to chronic energetic neutron radiation, particularly at the primary coolant loop hot and cold leg dissimilar metal welds, which join the primary loop piping to the reactor pressure vessel and reside in the vicinity of the reactor core. The objective of this experimental study was to determine how ultrasonic transducers and adhesive couplants perform when exposed to irradiation in a test reactor to simulate and accelerate in-service exposure. Further, to achieve this objective, the signal stability of piezoelectric transducers and performance of adhesive couplants as a function of accumulated fast neutron fluence were characterized by collecting ultrasonic data in-situ during irradiation. Of particular interest were the ultrasonic signal quality and time decay of the amplitude of acoustic reflections as a function of fast neutron fluence. The results of the study showed that, of the 8 transducer/substrate sample assemblies tested, only 3 generated usable ultrasonic signals through the conclusion of the irradiation campaign. It was found that high temperature epoxy tends to ultrasonically couple the sensors to the substrates better than three types of refractory ceramic cements studied, as is supported by post irradiation examination. The results obtained through this experimental study will be utilized in the achievement of the overall goal of development of a sensor system to perform online monitoring of primary loop components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

SMASH – A new hadronic transport approach

Microscopic transport approaches are the tool to describe the non-equilibrium evolution in low energy collisions as well as in the late dilute stages of high-energy collisions. Here, a newly developed hadronic transport approach, SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) is introduced. The overall bulk dynamics in low energy heavy ion collisions is shown including the excitation function of elliptic flow employing several equations of state. The implications of this new approach for dilepton production are discussed and preliminary results for afterburner calculations at the highest RHIC energy are presented and compared to previous UrQMD results. A detailed understanding of a hadron gas with vacuum properties is required to establish the baseline for the exploration of the transition to the quark-gluon plasma in heavy ion collisions at high net baryon densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Application of Wavelet-Based Methods for Accelerating Multi-Time-Scale Simulation of Bistable Heterogeneous Catalysis

Here, we report results from a numerical study of multi-time-scale bistable dynamics for CO oxidation on a catalytic surface in a flowing, well-mixed gas stream. The problem is posed in terms of surface and gas-phase submodels that dynamically interact in the presence of stochastic perturbations, reflecting the impact of molecular-scale fluctuations on the surface and turbulence in the gas. Wavelet-based methods are used to encode and characterize the temporal dynamics produced by each submodel and detect the onset of sudden state shifts (bifurcations) caused by nonlinear kinetics. When impending state shifts are detected, a more accurate but computationally expensive integration scheme can be used. This appears to make it possible, at least in some cases, to decrease the net computational burden associated with simulating multi-time-scale, nonlinear reacting systems by limiting the amount of time in which the more expensive integration schemes are required. Critical to achieving this is being able to detect unstable temporal transitions such as the bistable shifts in the example problem considered here. Lastly, our results indicate that a unique wavelet-based algorithm based on the Lipschitz exponent is capable of making such detections, even under noisy conditions, and may find applications in critical transition detection problems beyond catalysis.

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

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

Understanding Hydrogenation Chemistry at MgB 2 Reactive Edges from Ab Initio Molecular Dynamics

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB 2 without a priori assumption of reaction pathways. Focusing on highly reactive ($10\overline{1}0$) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH 4 – units and key chemical intermediates, involving H 2 dissociation, generation of functionalities and molecular complexes containing BH 2 and BH 3 motifs, and B–B bond breaking. Further, the genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B–H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg 2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB 2 -based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

08 HYDROGEN↗
DOE OSTI2017

A numerical analysis of a magnetocaloric refrigerator with a 16-layer regenerator

A numerical analysis was conducted to study a room temperature magnetocaloric refrigerator with a 16-layer parallel plates active magnetic regenerator (AMR). Sixteen layers of LaFeMnSiH having different Curie temperatures were employed as magnetocaloric material (MCM) in the regenerator. Measured properties data was used. A transient one dimensional (1D) model was employed, in which a unique numerical method was developed to significantly accelerate the simulation speed of the multi-layer AMR system. As a result, the computation speed of a multi-layer AMR case was very close to the single-layer configuration. The performance of the 16-layer AMR system in different frequencies and utilizations has been investigated using this model. To optimize the layer length distribution of the 16-layer MCMs in the regenerator, a set of 137 simulations with different MCM distributions based on the Design of Experiments (DoE) method was conducted and the results were analyzed. The results show that the 16-layer AMR system can operate up to 84% of Carnot cycle COP at a temperature span of 41 K, which cannot be obtained using an AMR with fewer layers. Here, the DoE results indicate that for a 16-layer AMR system, the uniform distribution is very close to the optimized design.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Investigation of ionization-induced electron injection in a wakefield driven by laser inside a gas cell

Ionization-induced electron injection was investigated experimentally by focusing a driving laser pulse with a maximum normalized potential of 1.2 at different positions along the plasma density profile inside a gas cell, filled with a gas mixture composed of 99%H 2 +1%N 2 . Changing the laser focus position relative to the gas cell entrance controls the accelerated electron bunch properties, such as the spectrum width, maximum energy, and accelerated charge. Simulations performed using the 3D particle-in-cell code WARP with a realistic density profile give results that are in good agreement with the experimental ones. Lastly, we discuss the interest of this regime for optimizing the bunch charge in a selected energy window.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

NONTHERMAL RADIATION FROM COSMIC-RAY MODIFIED SHOCKS

We calculate nonthermal radiation from cosmic-ray (CR) protons and electrons accelerated at CR modified plane and spherical shocks, using time-dependent, diffusive shock acceleration (DSA) simulations that include radiative losses of CR electrons. Strong non-relativistic shocks with physical parameters relevant for young supernova remnants (SNRs) are considered in both the plane-parallel and spherically symmetric geometries, and compared at times when their dynamical and CR properties are concordant. A thermal leakage injection model and a Bohm-like diffusion coefficient are adopted. After DSA energy gains balance radiative losses, the electron spectrum at the plane shock approaches a time-asymptotic spectrum with a super-exponential cutoff above the equilibrium momentum. The postshock electron spectrum cuts off at a progressively lower momentum downstream from the shock due to the energy losses. That results in the steepening of the volume integrated electron energy spectrum by one power of the particle energy. These features evolve toward lower energies in the spherical, SNR shocks. In a CR modified shock, pion decay gamma-ray emission reveals distinct signatures of nonlinear DSA due to the concave proton momentum spectrum. Although the electron momentum spectrum has a much weaker concavity, the synchrotron spectral slope at the shock may flatten by about 0.1-0.3 between radio and X-ray bands. The slope of the volume integrated emission spectrum behaves nonlinearly around the break frequency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Ion acceleration in electrostatic collisionless shock: on the optimal density profile for quasi-monoenergetic beams

Here, a numerical study on ion acceleration in electrostatic shock waves is presented, with the aim of determining the best plasma configuration to achieve quasi-monoenergetic ion beams in laser-driven systems. It was recently shown that tailored near-critical density plasmas characterized by a long-scale decreasing rear density profile lead to beams with low energy spread (Fiúza et al 2012 Phys. Rev. Lett. 109 215001). In this work, a detailed parameter scan investigating different plasma scale lengths is carried out. As result, the optimal plasma spatial scale length that allows for minimizing the energy spread while ensuring a significant reflection of ions by the shock is identified. Furthermore, a new configuration where the required profile has been obtained by coupling micro layers of different densities is proposed. Lastly, results show that this new engineered approach is a valid alternative, guaranteeing a low energy spread with a higher level of controllability.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Influence of the neutron-skin effect on nuclear isobar collisions at energies available at the BNL Relativistic Heavy Ion Collider

The unambiguous observation of a Chiral Magnetic Effect (CME)-driven charge separation is the core aim of the isobar program at RHIC consisting of $96\atop{40}$Zr+$96\atop{40}$Zr and $96\atop{44}$Ru+$96\atop{44}$Ru collisions at √s NN = 200 GeV. We quantify the role of the spatial distributions of the nucleons in the isobars on both eccentricity and magnetic field strength within a relativistic hadronic transport approach (SMASH, Simulating Many Accelerated Strongly-interacting Hadrons). In particular, we introduce isospin-dependent nucleon-nucleon spatial correlations in the geometric description of both nuclei, deformation for $96\atop{44}$Ru and the so-called neutron skin effect for the neutron-rich isobar i.e. $96\atop{40}$Zr. The main result of this study is a reduction of the magnetic field strength difference between $96\atop{44}$Ru+$96\atop{44}$Ru and $96\atop{40}$Zr+$96\atop{40}$Zr by a factor of 2, from 10% to 5% in peripheral collisions when the neutron-skin effect is included. Further, we find an increase of the eccentricity ratio between the isobars by up to 10% in ultra-central collisions as due to the deformation of $96\atop{44}$Ru while neither the neutron skin effect nor the nucleon-nucleon correlations result into a significant modification of this observable with respect to the traditional Woods-Saxon modeling. Our results suggest a significantly smaller CME signal to background ratio for the experimental charge separation measurement in peripheral collisions with the isobar systems than previously expected.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Strangeness production via resonances in heavy-ion collisions at energies available at the GSI Schwerionensynchrotron

Production of strange hadrons in elementary and heavy-ion reactions is studied with the hadronic transport approach called simulating many accelerated strongly interacting hadrons. The poorly known branching ratios of the relevant hadronic resonances are constrained from the known elementary hadronic cross sections and from invariant mass spectra of dileptons. The constrained model is employed as a baseline to compare to heavy-ion-collision experiments at low energies [E kin = (1-2) A GeV] and to predict some of the upcoming pion-beam results at the High-Acceptance Di-Electron Spectrometer, which are expected to be sensitive to the resonance properties. The employed vacuum-resonance approach proves to be viable for small systems at these energies, but for large systems additional in-medium effects might be required.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗