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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Roll-to-Roll SOFC Manufacturing for Low Cost Energy Generation

Oak Ridge National Laboratory (ORNL) and Redox Power Systems, LLC. collaborated on this project to develop an integrated, continuous roll-to-toll (R2R) manufacturing processes to produce solid oxide fuel cells (SOFCs). SOFCs are ceramic-based electrochemical devices that efficiently produce electricity from hydrogen-containing fuels, including natural gas. Standard processes for manufacturing SOFCs include tape casting of the anode and electrolyte, followed by stacking of the layers, and lamination using a uniaxial press or isostatic press. After high temperature firing, the cathode and electrical contact are screen printed and fired once again to complete the cell. This effort aimed to further develop R2R tape fabrication processes for green cell fabrication with the DOE Battery Manufacturing Facility (BMF) at ORNL as next generation improvements to current methods, focusing on R2R calendering and slot die coating for lower cost, high volume approaches for multilayer green cell fabrication. Two paths were investigated in parallel in this project. The first path was to optimize R2R processing using traditional tape cast materials from the Redox production chain on the BMF’s heated calender for compressing and laminating multilayer materials. The second path was to use BMF’s pilot-scale slot-die coater for R2R manufacturing as an alternative to tape casting of some or all of the layers in the green cell. Processes have been developed and and process conditions have been identified. Additionally, basic process economics were estimated for comparison to Redox current manufacturing approaches for the applicable process modules. The SBV program results will enable Redox to advance lower cost R2R manufacturing processes to produce millions of green cells per year on a single production line. As a result of this work, one invention disclosure has been submitted. In addition, one paper is under preparation for submission to a peer-reviewed journal.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

A Platform Technology for High-throughput Atomically Precise Manufacturing: Mechatronics at the Atomic Scale

The main objective of this project is to invent the necessary enabling technologies for high throughput atomically precise manufacturing (APM). APM is an emerging technology that refers to any manufacturing capability that enables fabrication of atomically precise structures, components, and devices under programmable control. APM will require positional assembly at the atomic and/or molecular scale, as well as at the nano and microscales using hierarchical assembly to create products ranging from nanoscale and quantum devices to macroscale systems and materials. This project builds on the expectation that commercial viability of APM will depend on a high level of parallelism to achieve the required throughput, a capability that currently does not exist. This project is the first attempt to address this key technological bottleneck.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Electron Acceleration at Pulsar Wind Termination Shocks

We study the acceleration of electrons and positrons at an electromagnetically modified, ultrarelativistic shock in the context of pulsar wind nebulae. We simulate the outflow produced by an obliquely rotating pulsar in proximity of its termination shock with a two-fluid code that uses a magnetic shear wave to mimic the properties of the wind. We integrate electron trajectories in the test-particle limit in the resulting background electromagnetic fields to analyze the injection mechanism. We find that the shock-precursor structure energizes and reflects a sizable fraction of particles, which becomes available for further acceleration. We investigate the subsequent first-order Fermi process sustained by small-scale magnetic fluctuations with a Monte Carlo code. We find that the acceleration proceeds in two distinct regimes: when the gyroradius r{sub g} exceeds the wavelength of the shear λ, the process is remarkably similar to first-order Fermi acceleration at relativistic, parallel shocks. This regime corresponds to a low-density wind that allows the propagation of superluminal waves. When r{sub g}<λ, which corresponds to the scenario of driven reconnection, the spectrum is softer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Modeling UV Radiation Feedback from Massive Stars. I. Implementation of Adaptive Ray-tracing Method and Tests [Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Algorithm]

Here, we present an implementation of an adaptive ray-tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses nonblocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems, including the propagation of radiation in vacuum and the expansions of various types of H II regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $\sim {10}^{3}$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is 12% when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with those from the M 1 closure relation. Although the ART and M 1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Myth of Fungible FTE: A Quantitative Assessment of Matrixed Resource Allocation

Matrix organizations allow scientific facilities to share specialized personnel across projects, operations, maintenance, and strategic initiatives. Nominal staffing allocations, however, may not capture the schedule consequences of fragmented individual commitments, limited access to specialist groups, and intermittent availability of key decision makers. We developed a stochas- tic, daily-time-step simulation of a hypothetical medium-sized accelerator-facility project com- prising sequential phases and parallel tasks. Each task requires role-specific work measured in FTE-days. Ordinary personnel may be unavailable because they contribute concurrently to other institutional activities, while designated key roles have independently specified daily un- availability probabilities. An organization-wide priority factor scales the number of people from each functional group who can effectively contribute to the project. It is interpreted as a composite proxy for project access and workforce fragmentation across competing commit- ments. We examined project completion time as a function of this factor and Project Lead unavailability using 100 Monte Carlo runs per condition. Increasing priority factor from 0.1 to 1.0 reduced median completion time from 1708.5 days (interquartile range 1681.5–1735.25) to 390 days (interquartile range 379–399). At priority factor = 0.1, increasing Project Lead unavailability from 0.5 to 0.9 increased median completion time from 1713.5 days (interquartile range 1691–1733.25) to 4,417 days (interquartile range 4271.75–4550.5). The model quantifies the commonly expected sensitivity of project schedules to fragmented resource commitments and limited coordination availability. Within this model, the results also indicate a possible threshold regime in which small increases in workforce availability yield only modest sched- ule improvements until sufficient capacity becomes accessible, after which project performance improves sharply. With further validation and calibration, this quantitative framework could support resource-allocation decisions during initial project planning and subsequent schedule rebaselining.

Bai, Mei [SLAC National Accelerator Laboratory (SL↗
DOE OSTI2016

Parallel electric fields are inefficient drivers of energetic electrons in magnetic reconnection

We present two-dimensional kinetic simulations, with a broad range of initial guide fields, which isolate the role of parallel electric fields (E- ∥ ) in energetic electron production during collisionless magnetic reconnection. In the strong guide field regime, E- ∥ drives essentially all of the electron energy gains, yet fails to generate an energetic component. Here, we suggest that this is due to the weak energy scaling of particle acceleration from E- ∥ compared to that of a Fermi-type mechanism responsible for energetic electron production in the weak guide-field regime. Lastly, this result has important implications for energetic electron production in astrophysical systems and reconnection-driven dissipation in turbulence.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

An in situ shearing x-ray measurement system for exploring structures and dynamics at the solid–liquid interface

Revealing interfacial structure and dynamics has been one of the essential thematic topics in material science and condensed matter physics. Synchrotron-based x-ray scattering techniques can deliver unique and insightful probing of interfacial structures and dynamics, in particular, in reflection geometries with higher surface and interfacial sensitivity than transmission geometries. We demonstrate the design and implementation of an in situ shearing x-ray measurement system, equipped with both inline parallel-plate and cone-and-plate shearing setups and operated at the advanced photon source at Argonne National Laboratory, to investigate the structures and dynamics of end-tethered polymers at the solid-liquid interface. With a precise lifting motor, a micrometer-scale gap can be produced by aligning two surfaces of a rotating upper shaft and a lower sample substrate. A torsional shear flow forms in the gap and applies tangential shear forces on the sample surface. The technical combination with nanoscale rheology and the utilization of in situ x-ray scattering allow us to gain fundamental insights into the complex dynamics in soft interfaces under shearing. In this work, we demonstrate the technical scope and experimental capability of the in situ shearing x-ray system through the measurements of charged polymers at both flat and curved interfaces upon shearing. Furthermore, through the in situ shearing x-ray scattering experiments integrated with theoretical simulations, we aim to develop a detailed understanding of the short-range molecular structure and mesoscale ionic aggregate morphology, as well as ion transport and dynamics in soft interfaces, thereby providing fundamental insight into a long-standing challenge in ionic polymer brushes with a significant technological impact.

36 MATERIALS SCIENCE↗
DOE OSTI2017

AMPS: An Augmented Matrix Formulation for Principal Submatrix Updates with Application to Power Grids

We present an augmented matrix approach to update the solution to a linear system of equations when the coefficient matrix is modified by a few elements within a principal submatrix. This problem arises in the dynamic security analysis of a power grid, where operators need to perform $N-x$ contingency analysis, i.e., determine the state of the system when up to $x$ links from $N$ fail. Our algorithms augment the coefficient matrix to account for the changes in it, and then compute the solution to the augmented system without refactoring the modified matrix. We provide two algorithms, a direct method, and a hybrid direct-iterative method for solving the augmented system. We also exploit the sparsity of the matrices and vectors to accelerate the overall computation. Our algorithms are compared on three power grids with PARDISO, a parallel direct solver, and CHOLMOD, a direct solver with the ability to modify the Cholesky factors of the coefficient matrix. We show that our augmented algorithms outperform PARDISO (by two orders of magnitude), and CHOLMOD (by a factor of up to 5). Further, our algorithms scale better than CHOLMOD as the number of elements updated increases. The solutions are computed with high accuracy. Our algorithms are capable of computing $N-x$ contingency analysis on a $778K$ bus grid, updating a solution with $x=20$ elements in $1.6 \times 10^{-2}$ seconds on an Intel Xeon processor.

Yeung, Yu-Hong↗
DOE OSTI2016

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↗
DOE OSTI2021

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↗
DOE OSTI2018

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↗
DOE OSTI2021

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↗
DOE OSTI2021

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↗
DOE OSTI2024

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↗
DOE OSTI2017

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↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗