Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel, scaling”

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,477 records · Page 82

FY 2016 Fusion Energy Sciences (FES) Theory and Simulation Performance Target

Annual target: Predicting the magnitude and scaling of the divertor heat load width in magnetically confined burning plasmas is a high priority for the fusion program and ITER. One of the key unresolved physics issues is what sets the heat flux width at the entrance to the divertor region. Perform massively parallel simulations using 3D edge kinetic and fluid codes to determine the parameter dependence of the heat load width at the divertor entrance and compute the divertor plate heat flux applicable to moderate particle recycling conditions. Comparisons will be made with data from DIII-D, NSTX-U, and C-Mod.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗

A comment on the position dependent diffusion coefficient representation of structural heterogeneity

Experimental studies of the variation of the mean square displacement (MSD) of a particle in a confined colloid suspension that exhibits density variations on the scale length of the particle diameter are not in agreement with the prediction that the spatial variation in MSD should mimic the spatial variation in density. The predicted behavior is derived from the expectation that the MSD of a particle depends on the system density and the assumption that the force acting on a particle is a point function of position. The experimental data are obtained from studies of the MSDs of particles in narrow ribbon channels and between narrowly spaced parallel plates and from new data, reported herein, of the radial and azimuthal MSDs of a colloid particle in a dense colloid suspension confined to a small circular cavity. In each of these geometries, a dense colloid suspension exhibits pronounced density oscillations with spacing of a particle diameter. We remove the discrepancy between prediction and experiment using the Fisher-Methfessel interpretation of how local equilibrium in an inhomogeneous system is maintained to argue that the force acting on a particle is delocalized over a volume with radius equal to a particle diameter. Our interpretation has relevance to the relationship between the scale of inhomogeneity and the utility of translation of the particle MSD into a position dependent diffusion coefficient and to the use of a spatially dependent diffusion coefficient to describe mass transport in a heterogeneous system.

74 ATOMIC AND MOLECULAR PHYSICS↗

Online data analysis and reduction: An important co-design motif for extreme-scale computers

A growing disparity between supercomputer computation speeds and I/O rates means that it is rapidly becoming infeasible to analyze supercomputer application output only after that output has been written to a file system. Instead, data-generating applications must run concurrently with data reduction and/or analysis operations, with which they exchange information via high-speed methods such as interprocess communications. The resulting parallel computing motif, online data analysis and reduction (ODAR), has important implications for both application and HPC systems design. Here we introduce the ODAR motif and its co-design concerns, describe a co-design process for identifying and addressing those concerns, present tools that assist in the co-design process, and present case studies to illustrate the use of the process and tools in practical settings.

Data Analysis↗

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗

Optimizing Distributed Training on Frontier for Large Language Models

Large language models (LLMs) have demonstrated remarkable success as foundational models, benefiting various downstream applications through fine-tuning. Loss scaling studies have demonstrated the superior performance of larger LLMs compared to their smaller counterparts. Nevertheless, training LLMs with billions of parameters poses significant challenges and requires considerable computational resources. For example, training a one trillion parameter GPT-style model on 20 trillion tokens requires a staggering 120 million exaflops. This research explores efficient distributed training strategies to extract this computation from Frontier, the world's first exascale supercomputer. We enable and investigate various model and data parallel training techniques, such as tensor parallelism, pipeline parallelism, and sharded data parallelism, to facilitate training a trillion-parameter model on Frontier. We empirically assess these techniques and their associated parameters to determine their impact on memory footprint, communication latency, and GPU's computational efficiency. We analyze the complex interplay among these techniques and find a strategy to combine them to achieve high throughput through hyperparameter tuning. We have identified efficient strategies for training large LLMs of varying sizes through empirical analysis and hyperparameter tuning. For 22 Billion, 175 Billion, and 1 Trillion parameters, we achieved GPU throughputs of 38.38%, 36.14%, and 31.96%, respectively. For the training of the 175 Billion parameter model and the 1 Trillion parameter model, we achieved 100% weak scaling efficiency on 1024 and 3072 Mi250X GPUs, respectively. We also achieved strong scaling efficiencies of 89% and 87% for these two models. We trained these models only tens of iterations instead of training till completion.

Yin, Junqi↗

On the parameter dependence of the whistler anisotropy instability

The evolution of the whistler anisotropy instability relevant to whistler-mode chorus waves in the Earth’s inner magnetosphere is studied using kinetic simulations and is compared with satellite observations. The electron distribution is constrained by the whistler anisotropy instability to a marginal stability state and presents an upper bound of electron anisotropy, which agrees with satellite observations. The electron beta β ∥e separates whistler waves into two groups: (i) quasi-parallel whistler waves for β ∥e ≳ 0.02 and (ii) oblique whistler waves close to the resonance cone for β ∥e ≲ 0.02. Landau damping is important in the saturation and relaxation stage of the oblique whistler wave growth. The saturated magnetic field energy of whistler waves roughly scales with the electron beta β 2 ∥e , shown in both simulations and satellite observations. Lastly, these results suggest the critical role of electron beta β ∥e in determining the whistler wave properties in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of electronic transport on electrochemical performance of (Cu,Mn) 3 O 4 solid oxide fuel cell cathodes

Alkaline Earth free spinel oxides provide a potential benefit over Sr-doped perovskite-based materials commonly used as electrodes in high-temperature electrochemical energy conversion devices, e.g., solid oxide fuel cells (SOFCs). Sr-segregation is a known issue leading to performance degradation. In this study, Cu x Mn 3-x O 4 (x = 1, 1.2, and 1.5) porous electrodes were examined as SOFC cathodes using electrochemical impedance spectroscopy to investigate the oxygen reduction reaction (ORR) kinetics in relation to the material's intrinsic conductivity, the extrinsic electrode structure, and the cell test design. Similar to the electronic conducting (La,Sr)MnO 3 SOFC cathodes, the ORR kinetics of Cu x Mn 3-x O 4 spinel electrodes was governed by the oxygen adsorption and diffusion at the particle surface as well as the charge transfer at the triple phase boundaries. The overall electrode polarization resistance was highly dependent on contact density with the metallic current collector, active material particle connectivity, electrode thickness, and the intrinsic electronic materials conductivity. Here, we describe the importance of effective electronic charge transport parallel to the electrode surface in maximizing the electrochemically active electrode volume and enhancing electrode performance. We discuss an approach to optimize cell and electrode design with respect to active materials properties. This aspect is critical to ensure reliable evaluation of new materials, since laboratory-scale button-cells typically exhibit a high degree of electrode microstructure (e.g. porosity, thickness) and electrical contact density variation from sample to sample.

(Cu,Mn)3O4↗

Estimating the randomness of quantum circuit ensembles up to 50 qubits

Random quantum circuits have been utilized in the contexts of quantum supremacy demonstrations, variational quantum algorithms for chemistry and machine learning, and blackhole information. The ability of random circuits to approximate any random unitaries has consequences on their complexity, expressibility, and trainability. To study this property of random circuits, we develop numerical protocols for estimating the frame potential, the distance between a given ensemble and the exact randomness. Our tensor-network-based algorithm has polynomial complexity for shallow circuits and is high-performing using CPU and GPU parallelism. We study 1. local and parallel random circuits to verify the linear growth in complexity as stated by the Brown–Susskind conjecture, and; 2. hardware-efficient ansätze to shed light on its expressibility and the barren plateau problem in the context of variational algorithms. Our work shows that large-scale tensor network simulations could provide important hints toward open problems in quantum information science.

97 MATHEMATICS AND COMPUTING↗

Parallel electron force balance and the L-H transition

In one popular paradigm for the L-H transition, energy transfer to the mean flows directly depletes turbulence fluctuation energy, resulting in suppression of the turbulence and a corresponding transport bifurcation. To quantitatively evaluate this mechanism, one must remember that electron parallel force balance couples nonzonal velocity fluctuations with electron pressure fluctuations on rapid timescales, comparable with the electron transit time. For this reason, energy in the nonzonal velocity stays in a fairly fixed ratio to the free energy in electron density fluctuations, at least for frequency scales much slower than electron transit. Furthermore, in order for direct depletion of the energy in turbulent fluctuations to cause the L-H transition, energy transfer via Reynolds stress must therefore drain enough energy to significantly reduce the sum of the free energy in nonzonal velocities and electron pressure fluctuations. At low k⊥, the electron thermal free energy is much larger than the energy in nonzonal velocities, posing a stark challenge for this model of the L-H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flow-through compression cell for small-angle and ultra-small-angle neutron scattering measurements

In situ measurements of geological materials under compression and with hydrostatic fluid pressure are important in understanding their behavior under field conditions, which in turn provides critical information for application-driven research. In particular, understanding the role of nano- to micro-scale porosity in the subsurface liquid and gas flow is critical for the high-fidelity characterization of the transport and more efficient extraction of the associated energy resources. In other applications, where parts are produced by the consolidation of powders by compression, the resulting porosity and crystallite orientation (texture) may affect its in-use characteristics. Small-angle neutron scattering (SANS) and ultra SANS are ideal probes for characterization of these porous structures over the nano to micro length scales. Here we show the design, realization, and performance of a novel neutron scattering sample environment, a specially designed compression cell, which provides compressive stress and hydrostatic pressures with effective stress up to 60 MPa, using the neutron beam to probe the effects of stress vectors parallel to the neutron beam. In conclusion, we demonstrate that the neutron optics is suitable for the experimental objectives and that the system is highly stable to the stress and pressure conditions of the measurements.

47 OTHER INSTRUMENTATION↗

A fast, dense Chebyshev solver for electronic structure on GPUs

Matrix diagonalization is almost always involved in computing the density matrix needed in quantum chemistry calculations. In the case of modest matrix sizes (≲4000), performance of traditional dense diagonalization algorithms on modern GPUs is underwhelming compared to the peak performance of these devices. This motivates the exploration of alternative algorithms better suited to these types of architectures. We newly derive, and present in detail, an existing Chebyshev expansion algorithm whose number of required matrix multiplications scales with the square root of the number of terms in the expansion. Focusing on dense matrices of modest size, our implementation on GPUs results in large speed ups when compared to diagonalization. Additionally, we improve upon this existing method by capitalizing on the inherent task parallelism and concurrency in the algorithm. Furthermore, this improvement is implemented on GPUs by using CUDA and HIP streams via the MAGMA library and leads to a significant speed up over the serial-only approach for smaller (≲1000) matrix sizes. Finally, we apply our technique to a model system with a high density of states around the Fermi level, which typically presents significant challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electric current paths in a Si:P delta-doped device imaged by nitrogen-vacancy diamond magnetic microscopy

The recently-developed ability to control phosphorous-doping of silicon at an atomic level using scanning tunneling microscopy, a technique known as atomic precision advanced manufacturing (APAM), has allowed us to tailor electronic devices with atomic precision, and thus has emerged as a way to explore new possibilities in Si electronics. In these applications, critical questions include where current flow is actually occurring in or near APAM structures as well as whether leakage currents are present. In general, detection and mapping of current flow in APAM structures are valuable diagnostic tools to obtain reliable devices in digital-enhanced applications. In this report, we used nitrogen-vacancy (NV) centers in diamond for wide-field magnetic imaging (with a few-mm field of view and micron-scale resolution) of magnetic fields from surface currents flowing in an APAM test device made of a P delta-doped layer on a Si substrate, a standard APAM witness material. We integrated a diamond having a surface NV ensemble with the device (patterned in two parallel mm-sized ribbons), then mapped the magnetic field from the DC current injected in the APAM device in a home-built NV wide-field microscope. The 2D magnetic field maps were used to reconstruct the surface current densities, allowing us to obtain information on current paths, device failures such as choke points where current flow is impeded, and current leakages outside the APAM-defined P-doped regions. Analysis on the current density reconstructed map showed a projected sensitivity of ~0.03 A m -1 , corresponding to a smallest-detectable current in the 200 μm wide APAM ribbon of ~6 μA. These results demonstrate the failure analysis capability of NV wide-field magnetometry for APAM materials, opening the possibility to investigate other cutting-edge microelectronic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Automated and highly parallelized Bayesian optimization scheme for direct drive fusion experiments on OMEGA

Finding the optimal implosion design on existing experimental facilities for inertial confinement fusion requires an exhaustive search of the vast design parameter space. This is infeasible both with experiments and with simulations. Consequently, a large fraction of the experimentally realizable design space remains unexplored, and new design schemes are challenging to optimize in a reasonable time frame. On the OMEGA laser facility, predictive machine learning models have been developed to accurately forecast the result of an experiment using only inexpensive simulations and the large dataset of prior experimental data. However, the full design space remains vast enough to be unassailable with simple optimization techniques. Here we develop an automated and optimally parallel Bayesian optimization algorithm that can entirely optimize the target and pulse shape of a direct-drive ICF implosion under a given design paradigm. We use this algorithm to find a markedly improved design for the performance implosions on OMEGA that is predicted to hydroequivalently scale to ignition at 2.15 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient Delaunay Tessellation through K-D Tree Decomposition

Delaunay tessellations are fundamental data structures in computational geometry. They are important in data analysis, where they can represent the geometry of a point set or approximate its density. The algorithms for computing these tessellations at scale perform poorly when the input data is unbalanced. We investigate the use of k-d trees to evenly distribute points among processes and compare two strategies for picking split points between domain regions. Because resulting point distributions no longer satisfy the assumptions of existing parallel Delaunay algorithms, we develop a new parallel algorithm that adapts to its input and prove its correctness. We evaluate the new algorithm using two late-stage cosmology datasets. The new running times are up to 50 times faster using k-d tree compared with regular grid decomposition. Moreover, in the unbalanced data sets, decomposing the domain into a k-d tree is up to five times faster than decomposing it into a regular grid.

Morozov, Dmitriy↗

Parallel electron force balance and the L-H transition

In one popular description of the L-H transition, energy transfer to the mean flows directly depletes turbulence fluctuation energy, resulting in suppression of the turbulence and a corresponding transport bifurcation. However, electron parallel force balance couples nonzonal velocity fluctuations with electron pressure fluctuations on rapid timescales, comparable with the electron transit time. For this reason, energy in the nonzonal velocity stays in a fairly fixed ratio to the free energy in electron density fluctuations, at least for frequency scales much slower than electron transit. In order for direct depletion of the energy in turbulent fluctuations to cause the L-H transition, energy transfer via Reynolds stress must therefore drain enough energy to significantly reduce the sum of the free energy in nonzonal velocities and electron pressure fluctuations. At low k, the electron thermal free energy is much larger than the energy in nonzonal velocities, posing a stark challenge for this model of the L-H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel electron force balance and the L-H transition

In one popular description of the L-H transition, energy transfer to the mean flows directly depletes turbulence fluctuation energy, resulting in suppression of the turbulence and a corresponding transport bifurcation. However, electron parallel force balance couples nonzonal velocity fluctuations with electron pressure fluctuations on rapid timescales, comparable with the electron transit time. For this reason, energy in the nonzonal velocity stays in a fairly fixed ratio to the free energy in electron density fluctuations, at least for frequency scales much slower than electron transit. In order for direct depletion of the energy in turbulent fluctuations to cause the L-H transition, energy transfer via Reynolds stress must therefore drain enough energy to significantly reduce the sum of the free energy in nonzonal velocities and electron pressure fluctuations. At low k, the electron thermal free energy is much larger than the energy in nonzonal velocities, posing a stark challenge for this model of the L-H transition.

L-H transition↗

MrHyDE v.1.0

SAND2024-01324O MrHyDE, which stands for Multi-resolution Hybridized Differential Equations, is a general-purpose C++ package for the solution of coupled multiphysics and multiscale systems on massively parallel computing systems. MrHyDE is designed to enable moving beyond forward simulation for multiscale applications which includes optimization, control, uncertainty quantification, and stochastic inversion. The framework provides interfaces to several packages within the Trilinos framework and leverages automatic differentiation to enable adjoint capabilities for large-scale, gradient-based optimization. MrHyDE provides automated multiscale capabilities through a subgrid model interface and multiscale Dirichlet-to-Neumann maps. For extreme-scale applications, MrHyDE provides in situ data-compression algorithms to reduce memory requirements while maintaining performance. MrHyDE is a general-purpose, computational framework for the solution of multiscale and multiphysics applications. It uses a combination of structure-preserving, physics-compatible discretizations, fully implicit methods, multi-resolution schemes, or fully explicit methods. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗