Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerated computing”

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,153 records · Page 64

Status of DUNE Offline Computing

We summarize the status of Deep Underground Neutrino Experiment (DUNE) Offline Software and Computing program. We describe plans for the computing infrastructure needed to acquire, catalog, reconstruct, simulate and analyze the data from the DUNE experiment and its prototypes in pursuit of the experiment's physics goals of precision measurements of neutrino oscillation parameters, detection of astrophysical neutrinos, measurement of neutrino interaction properties and searches for physics beyond the Standard Model. In contrast to traditional HEP computational problems, DUNE's Liquid Argon Time Projection Chamber data consist of simple but very large (many GB) data objects which share many characteristics with astrophysical images. We have successfully reconstructed and simulated data from 4% prototype detector runs at CERN. The data volume from the full DUNE detector, when it starts commissioning late in this decade will present memory management challenges in conventional processing but significant opportunities to use advances in machine learning and pattern recognition as a frontier user of High Performance Computing facilities capable of massively parallel processing. Our goal is to develop infrastructure resources that are flexible and accessible enough to support creative software solutions as HEP computing evolves.

43 PARTICLE ACCELERATORS↗

The Legendre Polynomial Axial Expansion Method

This work presents a new formulation of the axial expansion transport method explicitly using Legendre polynomials for arbitrarily high-order expansions. This new formulation also features an alternative method of axial leakage calculation to allow for nonextruded flat source region meshes. This alternative axial leakage is introduced alongside a balance equation requirement to ensure that neutron balance is preserved in the coarse mesh for a given axial leakage formulation, which allows for effective coarse mesh finite difference acceleration. A matrix exponential table method is derived to allow for fast computations of arbitrarily high-order matrix exponentials for this work and precludes the need for further research into matrix exponential calculations for this method. Numerical results are presented that demonstrate the stability of the axial expansion method in systems with voidlike regions, showcase the speedup from matrix exponential tables, and investigate the axial convergence of the method in terms of both expansion order and mesh size.

Herring, Nicholas↗

A GPU accelerated mixed-precision Smoothed Particle Hydrodynamics framework with cell-based relative coordinates

Smoothed Particle Hydrodynamics (SPH) is essential for modeling complex large-deformation problems across various applications, requiring significant computational power. A major portion of SPH computation time is dedicated to the Nearest Neighboring Particle Search (NNPS) process. While advanced NNPS algorithms have been developed to enhance SPH efficiency, the potential efficiency gains from modern computation hardware remain underexplored. Here, this study investigates the impact of GPU parallel architecture, low-precision computing on GPUs, and GPU memory management on NNPS efficiency. Our approach employs a GPU-accelerated mixed-precision SPH framework, utilizing low precision float-point 16 (FP16) for NNPS while maintaining high precision for other components. To ensure FP16 accuracy in NNPS, we introduce a Relative Coordinated-based Link List (RCLL) algorithm, storing FP16 relative coordinates of particles within background cells. Our testing results show three significant speedup rounds for CPU-based NNPS algorithms. The first comes from parallel GPU computations, with up to a 1000x efficiency gain. The second is achieved through low-precision GPU computing, where the proposed FP16-based RCLL algorithm offers a 1.5x efficiency improvement over the FP64-based approach on GPUs. By optimizing GPU memory bandwidth utilization, the efficiency of the FP16 RCLL algorithm can be further boosted by 2.7x, as demonstrated in an example with 1 million particles. Our code is released at https://github.com/pnnl/lpNNPS4SPH.

97 MATHEMATICS AND COMPUTING↗

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗

Computer simulation study of multiple germanium gamma-ray sensor arrays

Design studies of large high resolution gamma-ray spectrometer systems employing multiple sensors are greatly facilitated by the use of computer simulation techniques. A study has been made by interfacing a detector geometry code with the EGS-3 version of the Stanford Linear Accelerator Center Monte-Carlo Code which has been modified to run on our VAX computer. The detector responses have been simulated for input photons with energies between 0.1 and 10 MeV. Examples are given to illustrate the performance of complex sensor arrays with different configurations of the anticoincidence system. Of particular interest is a design composed of an array of 7 detector elements each consisting of a thin planar Ge (1 cm thick) detector mounted in front of and in tandem with a large coaxial Ge detector.

Nakano, G. H.↗

Study of shock-induced combustion using an implicit TVD scheme

The supersonic combustion flowfields associated with various hypersonic propulsion systems, such as the ram accelerator, the oblique detonation wave engine, and the scramjet, are being investigated using a new computational fluid dynamics (CFD) code. The code solves the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. It employs an iterative method and a second order differencing scheme to improve computational efficiency. The code is currently being applied to study shock wave/boundary layer interactions in premixed combustible gases, and to investigate the ram accelerator concept. Results obtained for a ram accelerator configuration indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outward and downstream. The combustion process creates a high pressure region over the back of the projectile resulting in a net positive thrust forward.

Yungster, Shayne↗

Empowering Scientific Innovation Through An Integrated Research Infrastructure: The Role of the Advanced Computing Ecosystem

As the landscape of computational science evolves, the Department of Energy (DOE) is reimagining the roles of its large-scale computing facilities to meet emerging research challenges. The Integrated Research Infrastructure (IRI) program aims to transform how experiments are designed, conducted, and shared, with significant impacts on all stakeholders. In response, the Oak Ridge Leadership Computing Facility (OLCF) has established the Advanced Computing Ecosystem (ACE), a strategic framework to prepare its hardware, software, and experimental capabilities for the IRI era. ACE focuses on integrating novel compute environments, orchestrating advanced workflows, and developing foundational technologies, ensuring a seamless transition to IRI while accelerating scientific discovery. This paper outlines ACE's role in advancing OLCF's mission and its impact on the future of computational science.

Widener, Patrick↗

Celeritas R&D Report: Accelerating Geant4

Celeritas is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a Graphics Processing Unit (GPU)-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.4, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling straightforward integration of Celeritas into the high energy physics (HEP) experimental frameworks CMSSW and ATLAS FullSimLight. On the Perlmutter supercomputer, Celeritas performs EM physics between 3× and 18× faster using the machine’s Nvidia GPUs compared to using only CPUs, corresponding to an electrical power efficiency up to a factor of 5. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.7–2.2× on GPU and 1.2× on CPU. In a CMS test application using tt¯ events and the prototype Run 4 configuration, compared to Geant4 CPU, Celeritas with a Nvidia A100 improves overall throughput up to a factor of 2.7× but cannot be efficiently shared with more than 8 cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Engineering of niobium surfaces through accelerated neutral atom beam technology for quantum applications

A major roadblock to scalable quantum computing is phase decoherence and energy relaxation caused by qubits interacting with defect-related two-level systems (TLSs). Native oxides present on the surfaces of superconducting metals used in quantum devices are acknowledged to be a source of TLS that decrease qubit coherence times. Reducing microwave loss by “surface engineering” (i.e., replacing the uncontrolled native oxide of superconducting metals with a thin, stable surface with predictable characteristics) can be a key enabler for pushing performance forward with devices of higher intrinsic quality factor. In this work, we present a novel approach to replace the native oxide of niobium (typically formed in an uncontrolled fashion when its pristine surface is exposed to air) with an engineered oxide, using a room-temperature process that leverages accelerated neutral atom beam (ANAB) technology at 300 mm wafer scale. This ANAB is composed of a mixture of argon and oxygen, with tunable energy per atom, which is rastered across the wafer surface. The ANAB-engineered Nb-oxide thickness was found to vary from 2 to 6 nm depending on ANAB process parameters. The modeling of variable-energy x-ray photoelectron spectroscopy data confirms the thickness and compositional control of Nb surface oxide by the ANAB process. Furthermore, these results correlate well with those from transmission electron microscopy and x-ray reflectometry. Since ANAB is broadly applicable to material surfaces, the present study indicates its promise for modification of the surfaces of superconducting quantum circuits to achieve longer coherence times.

36 MATERIALS SCIENCE↗

Probing Non-Equilibrium Pair-Breaking and Quasiparticle Dynamics in Nb Superconducting Resonators Under Magnetic Fields

We conducted a comprehensive study of the non-equilibrium dynamics of Cooper pair breaking, quasiparticle (QP) generation, and relaxation in niobium (Nb) cut from superconducting radio-frequency (SRF) cavities, as well as various Nb resonator films from transmon qubits. Using ultrafast pump–probe spectroscopy, we were able to isolate the superconducting coherence and pair-breaking responses. Our results reveal both similarities and notable differences in the temperature- and magnetic-field-dependent dynamics of the SRF cavity and thin-film resonator samples. Moreover, femtosecond-resolved QP generation and relaxation under an applied magnetic field reveals a clear correlation between non-equilibrium QPs and the quality factor of resonators fabricated by using different deposition methods, such as DC sputtering and high-power impulse magnetron sputtering. These findings highlight the pivotal influence of fabrication techniques on the coherence and performance of Nb-based quantum devices, which are vital for applications in superconducting qubits and high-energy superconducting radio-frequency applications.

43 PARTICLE ACCELERATORS↗

Molecular fingerprint and machine learning to accelerate design of high-performance homochiral metal–organic frameworks

In this report computational screening was employed to calculate the enantioseparation capabilities of 45 functionalized homochiral metal–organic frameworks (FHMOFs), and machine learning (ML) and molecular fingerprint (MF) techniques were used to find new FHMOFs with high performance. With increasing temperature, the enantioselectivities for (R,S)-1,3-dimethyl-1,2-propadiene are improved. The “glove effect” in the chiral pockets was proposed to explain the correlations between the steric effect of functional groups and performance of FHMOFs. Moreover, the neighborhood component analysis and RDKit/MACCS MFs show the highest predictive effect on enantioselectivities among the four ML classification algorithms with nine MFs that were tested. Based on the importance of MF, 85 new FHMOFs were designed, and a newly designed FHMOF, NO 2 -NHOH-FHMOF, with high similarity to the optimal MFs achieved improved chiral separation performance, with enantioselectivities of 85%. The design principles and new chiral pockets obtained by ML and MFs could facilitate the development of new materials for chiral separation.

molecular fingerprint↗

An autonomous laboratory for the accelerated synthesis of novel materials

To close the gap between the rates of computational screening and experimental realization of novel materials, we introduce the A-Lab, an autonomous laboratory for the solid-state synthesis of inorganic powders. This platform uses computations, historical data from the literature, machine learning (ML) and active learning to plan and interpret the outcomes of experiments performed using robotics. Over 17 days of continuous operation, the A-Lab realized 41 novel compounds from a set of 58 targets including a variety of oxides and phosphates that were identified using large-scale ab initio phase-stability data from the Materials Project and Google DeepMind. Synthesis recipes were proposed by natural-language models trained on the literature and optimized using an active-learning approach grounded in thermodynamics. Analysis of the failed syntheses provides direct and actionable suggestions to improve current techniques for materials screening and synthesis design. The high success rate demonstrates the effectiveness of artificial-intelligence-driven platforms for autonomous materials discovery and motivates further integration of computations, historical knowledge and robotics.

36 MATERIALS SCIENCE↗

Accelerating Combustion and Surface Chemistry Simulations

Design of modern combustion systems relies on computer models to predict how changes in design will affect performance. These models have largely displaced previous methods that rely on the designer’s intuition or costly and time-consuming physical testing. By using improved models, design cycles can be shortened, and cleaner and more efficient combustion devices can be created. This project aims to improve computer simulations of transportation fuels with the goal of making these simulations faster and more accurate for predicting combustion in vehicles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Hvac: Removing I/O Bottleneck for Large-Scale Deep Learning Applications

Scientific communities are increasingly adopting deep learning (DL) models in their applications to accelerate scientific discovery processes. However, with rapid growth in the computing capabilities of HPC supercomputers, large-scale DL applications have to spend a significant portion of training time performing I/O to a parallel storage system. Previous research works have investigated optimization techniques such as prefetching and caching. Unfortunately, there exist non-trivial challenges to adopting the existing solutions on HPC supercomputers for large-scale DL training applications, which include non-performance and/or failures at extreme scale, lack of portability and generality in design, complex deployment methodology, and being limited to a specific application or dataset. To address these challenges, we propose High-Velocity AI Cache (HVAC), a distributed read-cache layer that targets and fully exploits the node-local storage or near node-local storage technology. HVAC seamlessly accelerates read I/O by aggregating node-local or near node-local storage, avoiding metadata lookups and file locking while preserving portability in the application code. We deploy and evaluate HVAC on 1,024 nodes (with over 6000 NVIDIA V100 GPUS) of the Summit supercomputer. In particular, we evaluate the scalability, efficiency, accuracy, and load distribution of HVAC compared to GPFS and XFS-on-NVMe. With four different DL applications, we observe an average 25 % performance improvement atop GPFS and 9% drop against XFS-on-NVMe, which scale linearly and are considered the performance upper bound. We envision HVAC as an important caching library for upcoming HPC supercomputers such as Frontier.

Khan, Awais↗

Three-dimensional structural analysis using interactive graphics

The application of computer interactive graphics to three-dimensional structural analysis was described, with emphasis on the following aspects: (1) structural analysis, and (2) generation and checking of input data and examination of the large volume of output data (stresses, displacements, velocities, accelerations). Handling of three-dimensional input processing with a special MESH3D computer program was explained. Similarly, a special code PLTZ may be used to perform all the needed tasks for output processing from a finite element code. Examples were illustrated.

Biffle, J.↗

Constraint Embedding Technique for Multibody System Dynamics

Multibody dynamics play a critical role in simulation testbeds for space missions. There has been a considerable interest in the development of efficient computational algorithms for solving the dynamics of multibody systems. Mass matrix factorization and inversion techniques and the O(N) class of forward dynamics algorithms developed using a spatial operator algebra stand out as important breakthrough on this front. Techniques such as these provide the efficient algorithms and methods for the application and implementation of such multibody dynamics models. However, these methods are limited only to tree-topology multibody systems. Closed-chain topology systems require different techniques that are not as efficient or as broad as those for tree-topology systems. The closed-chain forward dynamics approach consists of treating the closed-chain topology as a tree-topology system subject to additional closure constraints. The resulting forward dynamics solution consists of: (a) ignoring the closure constraints and using the O(N) algorithm to solve for the free unconstrained accelerations for the system; (b) using the tree-topology solution to compute a correction force to enforce the closure constraints; and (c) correcting the unconstrained accelerations with correction accelerations resulting from the correction forces. This constraint-embedding technique shows how to use direct embedding to eliminate local closure-loops in the system and effectively convert the system back to a tree-topology system. At this point, standard tree-topology techniques can be brought to bear on the problem. The approach uses a spatial operator algebra approach to formulating the equations of motion. The operators are block-partitioned around the local body subgroups to convert them into aggregate bodies. Mass matrix operator factorization and inversion techniques are applied to the reformulated tree-topology system. Thus in essence, the new technique allows conversion of a system with closure-constraints into an equivalent tree-topology system, and thus allows one to take advantage of the host of techniques available to the latter class of systems. This technology is highly suitable for the class of multibody systems where the closure-constraints are local, i.e., where they are confined to small groupings of bodies within the system. Important examples of such local closure-constraints are constraints associated with four-bar linkages, geared motors, differential suspensions, etc. One can eliminate these closure-constraints and convert the system into a tree-topology system by embedding the constraints directly into the system dynamics and effectively replacing the body groupings with virtual aggregate bodies. Once eliminated, one can apply the well-known results and algorithms for tree-topology systems to solve the dynamics of such closed-chain system.

Woo, Simon S.↗

Spring 2014 Internship Diffuser Data Analysis

J-2X engine testing on the A-2 test stand at the NASA John C. Stennis Space Center (SSC) has recently concluded. As part of that test campaign, the engine was operated at lower power levels in support of expanding the use of J-2X to other missions. However, the A-2 diffuser was not designed for engine testing at the proposed low power levels. To evaluate the risk of damage to the diffuser, computer simulations were created of the rocket engine exhaust plume inside the 50ft long, water-cooled, altitude-simulating diffuser. The simulations predicted that low power level testing would cause the plume to oscillate in the lower sections of the diffuser. This can possibly cause excessive vibrations, stress, and heat transfer from the plume to the diffuser walls. To understand and assess the performance of the diffuser during low power level engine testing, nine accelerometers and four strain gages were installed around the outer surface of the diffuser. The added instrumentation also allowed for the verification of the rocket exhaust plume computational model. Prior to engine hot-fire testing, a diffuser water-flow test was conducted to verify the proper operation of the newly installed instrumentation. Subsequently, two J-2X engine hot-fire tests were completed. Hot-Fire Test 1 was 11.5 seconds in duration, and accelerometer and strain data verified that the rocket engine plume oscillated in the lower sections of the diffuser. The accelerometers showed very different results dependent upon location. The diffuser consists of four sections, with Section 1 being closest to the engine nozzle and Section 4 being farthest from the engine nozzle. Section 1 accelerometers showed increased amplitudes at startup and shutdown, but low amplitudes while the diffuser was started. Section 3 accelerometers showed the opposite results with near zero G amplitudes prior to and after diffuser start and peak amplitudes to +/- 100G while the diffuser was started. Hot-Fire Test 1 strain gages showed different data dependent on section. Section 1 strains were small, and were in the range of 50 to 150 microstrain, which would result in stresses from 1.45 to 4.35 ksi. The yield stress of the material, A-285 Grade C Steel, is 29.7 ksi. Section 4 strain gages showed much higher values with strains peaking at 1600 microstrain. This strain corresponds to a stress of 46.41 ksi, which is in excess of the yield stress, but below the ultimate stress of 55 to 75 ksi. The decreased accelerations and strain in Section 1, and the increased accelerations and strain in Sections 3 and 4 verified the computer simulation prediction of increased plume oscillations in the lower sections of the diffuser. Hot-Fire Test 2 ran for a duration of 125 seconds. The engine operated at a slightly higher power level than Hot-Fire Test 1 for the initial 35 seconds of the test. After 35 seconds the power level was lowered to Hot-Fire Test 1 levels. The acceleration and strain data for Hot-Fire Test 2 was similar during the initial part of the test. However, just prior to the engine being lowered to the Hot-Fire Test 1 power level, the strain gage data in Section 4 showed a large decrease to strains near zero microstrain from their peak at 1500 microstrain. Future work includes further strain and acceleration data analysis and evaluation.

Laigaie, Robert T.↗

VTK-m: Visualization for the Exascale Era and Beyond

A recent trend in modern high-performance computing is the increasing use of hybrid architectures, where the vast majority of performance comes from accelerators. Modern accelerators are based on Graphics Processing Units (GPU) that contain many low power cores that in their aggregate provides an extremely high computation rate. Current and future CPU processors are requiring more explicit parallelism as each successive version of the hardware packs in more cores, and technologies like hyperthreading and vector operations require even more parallel processing to leverage each core’s full potential. As an example, the Frontier supercomputer installed at Oak Ridge National Laboratories recently hit a record breaking 1.1 exaflops1 on the LINPACK HPC benchmark [Shoemaker 2022]. The system contains 37632 AMD MI250x GPUs which requires more than half a billion threads to keep the system fully utilized [Khizeran 2022].VTK-m is a toolkit of scientific visualization algorithms for these emerging processor architectures. VTK-m supports the fine-grained concurrency for data analysis and visualization algorithms required to drive extreme scale computing by providing abstract models for data and execution that can be applied to a variety of algorithms across many different processor architectures.

Bolstad, Mark↗