Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel data structures”

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 235 records · Page 13

Task-based Parallel Computation of the Density Matrix in Quantum-based Molecular Dynamics using Graph Partitioning

Quantum-based molecular dynamics (QMD) is a highly accurate and transferable method for material science simulations. However, the time scales and system sizes accessible to QMD are typically limited to picoseconds and a few hundred atoms. These constraints arise due to expensive self-consistent ground-state electronic structure calculations that can often scale cubically with the number of atoms. Linearly scaling methods depend on computing the density matrix P from the Hamiltonian matrix H by exploiting the sparsity in both matrices. The second-order spectral projection (SP2) algorithm is an O(N) algorithm that computes P with a sequence of 40 - 50 matrix-matrix multiplications. In this paper, we present task-based implementations of a recently developed data-parallel graph-based approach to the SP2 algorithm, G-SP2. We represent the density matrix P as an undirected graph and use graph partitioning techniques to divide the computation into smaller independent tasks. The partitions thus obtained are generally not of equal size and give rise to undesirable load imbalances in standard MPI-based implementations. This load-balancing challenge can be mitigated by dynamically scheduling parallel computations at runtime using task-based programming models. We develop task-based implementations of the data-parallel G-SP2 algorithm using both Intel's Concurrent Collections (CnC) as well as the Charm++ programming model and evaluate these implementations for future use. Scaling and performance results of our implementations are investigated for representative segments of QMD simulations for solvated protein systems containing more than 10,000 atoms.

97 MATHEMATICS AND COMPUTING↗

Crystal structure of halofuginone hydrobromide, C 16 H 18 BrClN 3 O 3 Br

The crystal structure of one form of halofuginone hydrobromide has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Halofuginone hydrobromide crystallizes in space group P2 1 (#4) with a = 8.87398(13), b = 14.25711(20), c = 15.0153(3) Å, β = 91.6867(15)°, V = 1898.87(4) Å 3 , and Z = 4. The crystal structure consists of alternating layers (parallel to the ab-plane) of planar and nonplanar portions of the cations. N–H∙∙∙Br and O–H∙∙∙Br hydrogen bonds link the protonated piperidine rings and bromide anions into a two-dimensional network parallel to the ab-plane. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

Incorporating Human Readiness Levels at Sandia National Laboratories

Since 2010, the concept of human readiness levels has been under development as a possible supplement to the existing technology readiness level (TRL) scale. The intent is to provide a mechanism to address safety and performance risks associated with the human component in a system that parallels the TRL structure already familiar to the systems engineering community. Sandia National Laboratories in Albuquerque, New Mexico, initiated a study in 2015 to evaluate options to incorporate human readiness planning for Sandia processes and products. The study team has collected the majority of baseline assessment data and has conducted interviews to understand staff perceptions of four different options for human readiness planning. Preliminary results suggest that all four options may have a vital role, depending on the type of work performed and the phase of product development. Upon completion of data collection, the utility of identified solutions will be assessed in one or more test cases.

99 GENERAL AND MISCELLANEOUS↗

Crystal structure of diclazuril, C 17 H 9 Cl 3 N 4 O 2

The crystal structure of diclazuril has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Diclazuril crystallizes in space group P2 1 /a (#14) with a = 27.02080(18), b = 11.42308(8), c = 5.36978(5) Å, β = 91.7912(7)°, V = 1656.629(15) Å 3 , and Z = 4. The crystal structure consists of layers of molecules parallel to the ac-plane. A strong N–H∙∙∙O hydrogen bond links the molecules into dimers along the a-axis with a graph set R2,2(8). The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

GridOPTICS/GridPACK

GridPACK is a software framework consisting of a set of modules designed to simplify the development of programs that model the power grid and run on parallel, high performance computing platforms. It also contains several fully developed applications, including powerflow, dynamic simulation, state estimation, Kalman filter analysis (dynamic state estimation), contingency analysis and real time path rating. These applications can be used either standalone or as components in more complicated workflows that combine several different types of application together. The framework modules are available as a combination of libraries and software templates and consist of components for setting up and distributing power grid networks, support for modeling the behavior of individual buses and branches in the network, converting the network models to the corresponding algebraic equations, and parallel routines for manipulating and solving large algebraic systems. The framework also contains a module for distributing tasks evenly amongst computing resources, even if individual tasks vary widely in their execution times. Additional modules support input and output, basic statistical analysis of contingency based calculations, distributed data structures, as well as basic profiling and error management.

Palmer, Bruce↗

Distributed-Memory Breadth-First Search on Massive Graphs. IN: Parallel Graph Algorithms

This chapter studies the problem of traversing large graphs using the breadth-first search order on distributed-memory supercomputers. We consider both the traditional level-synchronous top-down algorithm as well as the recently discovered direction optimizing algorithm. We analyze the performance and scalability trade-offs in using different local data structures such as CSR and DCSC, enabling in-node multithreading, and graph decompositions such as 1D and 2D decomposition.

97 MATHEMATICS AND COMPUTING↗

Crystal structure of vismodegib, C 19 H 14 Cl 2 N 2 O 3 S

The crystal structure of vismodegib has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Vismodegib crystallizes in space group P2 1 /a (#14) with a = 16.92070(20), b = 10.20235(4), c = 12.16161(10) Å, β = 108.6802(3)°, V = 1988.873(9) Å 3 , and Z = 4. The crystal structure consists of corrugated layers of molecules parallel to the bc-plane. There is only one classical hydrogen bond in the structure, between the amide nitrogen atom and the N atom of the pyridine ring. Pairs of these hydrogen bonds link the molecules into dimers, with a graph set R2,2(14) > a > a. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

Crystal structure of oxibendazole, C 12 H 15 N 3 O 3

The crystal structure of oxibendazole has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Oxibendazole crystallizes in space group C2/c (#15) with a = 23.18673(22), b = 5.35136(5), c = 19.88932(13) Å, β = 97.0876(9)°, V = 2449.018(17) Å 3 , and Z = 8. The structure consists of hydrogen-bonded layers of planar molecules parallel to the bc-plane. Strong N–H∙∙∙N hydrogen bonds link the molecules into dimers, with a graph set R2,2(8). N–H∙∙∙O hydrogen bonds further link these dimers into layers parallel to the bc-plane. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

Crystal structure of elvitegravir Form II, C 23 H 23 ClFNO 5

The crystal structure of elvitegravir Form II has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Elvitegravir Form II crystallizes in space group P2 1 (#4) with a = 11.54842(7), b = 14.04367(5), c = 13.33333(8) Å, β = 90.0330(6)°, V = 2162.427(14) Å 3 , and Z = 4. The crystal structure consists of alternating layers of parallel molecules perpendicular to the b-axis. The mean planes of the oxoquinoline ring systems in molecules 1 and 2 are 1(22)-1 and -1(22)1. Between the stacks are layers of the halogenated phenyl rings. These exhibit herringbone stacking. In each molecule, the carboxylic acid group forms a strong intramolecular O–H∙∙∙O hydrogen bond to the nearby carbonyl group. The hydroxyl group of each molecule forms a strong hydrogen bond to the carbonyl group of the carboxylic acid of the other molecule. These O–H∙∙∙O hydrogen bonds link the molecules into dimers, with a graph set R2,2(18) > a > c. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

Experimental setup combining in situ hard X-ray photoelectron spectroscopy and real-time surface X-ray diffraction for characterizing atomic and electronic structure evolution during complex oxide heterostructure growth

We describe the next-generation system for in situ characterization of a complex oxide thin film and heterostructure growth by pulsed laser deposition (PLD) using synchrotron hard X-rays. The system consists of a PLD chamber mounted on a diffractometer allowing both real-time surface X-ray diffraction (SXRD) and in situ hard X-ray photoelectron spectroscopy (HAXPES). HAXPES is performed in the incident X-ray energy range from 4 to 12 keV using a Scienta EW4000 electron energy analyzer mounted on the PLD chamber fixed parallel with the surface normal. In addition to the standard application mode of HAXPES for disentangling surface from bulk properties, the increased penetration depth of high energy photoelectrons is used for investigation of the electronic structure changes through thin films grown deliberately as variable thickness capping layers. Such heterostructures represent model systems for investigating a variety of critical thickness and dead layer phenomena observed at complex oxide interfaces. In this new mode of operation, in situ HAXPES is used to determine the electronic structure associated with unique structural features identified by real-time SXRD during thin film growth. The system is configured for using both laboratory excitation sources off-line and on-line operation at beamline 33-ID-D at the Advanced Photon Source. We illustrate the performance of the system by preliminary scattering and spectroscopic data on oxygen vacancy ordering induced perovskite-to-brownmillerite reversible phase transformation in La 2/3 Sr 1/3 MnO 3 films capped with oxygen deficient SrTiO 3–δ (100) layers of varying thickness.

47 OTHER INSTRUMENTATION↗

Porting hypre to heterogeneous computer architectures: Strategies and experiences

We report that linear systems are occurring in many applications, and solving them can take a large amount of the total simulation time. The high performance library hypre provides a variety of interfaces and linear solvers, including various multigrid methods, that have achieved good scalability on a variety of homogeneous parallel computer architectures. Heterogeneous architectures with nodes that have both CPUs and accelerators provide new challenges, since they require more fine-grained parallelism and reduced data movement between different memories on a single node as well as across nodes. We will discuss our experiences and strategies to port hypre to heterogeneous computers with accelerators, including the design of a new memory model, the use of abstractions, the BoxLoop macros in the structured and semi-structured interfaces, and the restructuring of algebraic multigrid (AMG) into modular components. We present numerical experiments comparing CPU and GPU performance for several test problems.

97 MATHEMATICS AND COMPUTING↗

Bit-GraphBLAS: Bit-Level Optimizations of Matrix-Centric Graph Processing on GPU

In the graph data structure like adjacency matrix, the connectivity of two nodes can be sufficiently represented using only 1 bit, but they are generally treated as 32-bit full-precision in state-of-the-art graph frameworks to adopt common sparse format such as CSR. Meanwhile, bit-level parallelism has recently be explored to have high-performance potential and low storage requirement on GPUs with dense bit-tiles. To fill the gap, our solution is a hierarchical storage format that contains the bit-indexing base and dense bit-tile units. Inherently, the granularity of the bit-tile is an essential factor in achieving both storage compression and GPU parallelism. How to find a sweet spot that trades off between avoiding sparsity and exploiting is comprehensively researched in this work. In the experiment, we evaluate the proposed storage format and algorithms on modern generation GPUs, including Pascal and Volta, to figure out critical software co-designs in conjunction with existing hardware-specific optimization.

Chen, Jou-An↗

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗

Crystal structure of brimonidine hydrogen tartrate, (C 11 H 11 BrN 5 )(HC 4 H 4 O 6 )

The crystal structure of brimonidine hydrogen tartrate has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional techniques. Brimonidine hydrogen tartrate crystallizes in space groupP2 1 (#4) witha= 7.56032(2),b= 7.35278(2),c= 30.10149(9) Å,β= 90.1992(2)°,V= 1673.312(10) Å 3 , andZ= 4 at 295 K. The crystal structure consists of alternating layers of cations and anions parallel to theab-plane. Each of the hydrogen tartrate anions is linked to itself by very strong charge-assisted O–H⋯O hydrogen bonds into chains along thea-axis. Each hydroxyl group of each tartrate acts as a donor in an O–H⋯O or O–H⋯N hydrogen bond. One of these is intramolecular, but the other three are intermolecular. These hydrogen bonds link the hydrogen tartrate anions into layers parallel to theab-plane and also link the anion–cation layers. The protonated N atoms act as donors in N–H⋯O or N–H⋯N hydrogen bonds to the carboxyl groups of the tartrates and to a ring nitrogen atom. These link the cations and anions, as well as providing cation–cation links. The amino N atoms of the cations form N–H⋯O hydrogen bonds to hydroxyl groups of the anions. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®)

Materials Science↗

Observation of Skewed Electromagnetic Wakefields in an Asymmetric Structure Driven by Flat Electron Bunches

Relativistic charged -particle beams that generate intense longitudinal fields in accelerating structures also inherently couple to transverse modes. The effects of this coupling may lead to beam breakup instability and thus must be countered to preserve beam quality in applications such as linear colliders. Beams with highly asymmetric transverse sizes (flat beams) have been shown to suppress the initial instability in slab -symmetric structures. However, as the coupling to transverse modes remains, this solution serves only to delay instability. In order to understand the hazards of transverse coupling in such a case, we describe here an experiment characterizing the transverse effects on a flat beam, traversing near a planar dielectric lined structure. Further, the measurements reveal the emergence of a previously unobserved skew-quadrupolelike interaction when the beam is canted transversely, which is not present when the flat beam travels parallel to the dielectric surface. We deploy a multipole field fitting algorithm to reconstruct the projected transverse wakefields from the data. We generate the effective kick vector map using a simple two -particle theoretical model, with particle -in -cell simulations used to provide further insight for realistic particle distributions.

43 PARTICLE ACCELERATORS↗

Understanding the Design Space of Sparse/Dense Multiphase Dataflows for Mapping Graph Neural Networks on Spatial Accelerators

Graph Neural Networks (GNNs) have garnered a lot of recent interest because of their success in learning representations from graph-structured data across several critical applications in cloud and HPC. Owing to their unique compute and memory characteristics that come from an interplay between dense and sparse phases of computations, the emergence of reconfigurable dataflow (aka spatial) accelerators offers promise for acceleration by mapping optimized dataflows (i.e., computation order and parallelism) for both phases. The goal of this work is to characterize and understand the design-space of dataflow choices for running GNNs on spatial accelerators in order for the compilers to optimize the dataflow based on the workload. Specifically, we propose a taxonomy to describe all possible choices for mapping the dense and sparse phases of GNNs spatially and temporally over a spatial accelerator, capturing both the intra-phase dataflow and the inter-phase (pipelined) dataflow. Using this taxonomy, we do deep-dives into the cost and benefits of several dataflows and perform case studies on implications of hardware parameters for dataflows and value of flexibility to support pipelined execution.

97 MATHEMATICS AND COMPUTING↗

Mesoflow: An Open-Source Reacting Flow Solver for Catalysis at Mesoscale

We present the capabilities and software performance metrics of our open-source continuum solver for catalysis, Mesoflow, developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex catalyst surface morphologies directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous catalyst interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPUs) compared to a single processor for representative problem sizes (2 million cell mesh). We will also present a brief introduction on how to build and use this software for application problems pertaining to catalytic upgrading and gas transport within porous catalyst particles.

adaptive meshing↗

Spotlight: efficient automated global optimization in rietveld analysis of diffraction data

Performing reliable Rietveld analysis on tens or hundreds of powder diffraction datasets from parametric or time-resolved experiments often poses a bottleneck in extracting meaningful results from the data. While automated analysis of data has recently been demonstrated, high temperature annealing studies, during which phase transformations occur and lattice parameters may change due to repartitioning of elements, are prime examples where automation by a simple phase identification from a database of room temperature structures or automation by sequential refinements is likely to fail. To enable reliable, efficient, automated Rietveld analysis, we present a Python package named Spotlight , building on established Rietveld packages such as MAUD, GSAS , or GSAS-II , which extends the refinement of best fit parameters to a global optimization using an ensemble of optimizers leveraging hierarchical parallel execution on high-performance computing clusters. Spotlight further enables the efficient design of refinement plans through the iterative automated machine-learning of a surrogate for the refinement on which the global optimizations are performed until results from the surrogate converge to the response surface data. We demonstrate Spotlight with the analysis of uranium molybdenum and Ti–6Al–4V datasets, as well as in two open-source tutorials analyzing aluminium oxide and lead sulphate.

36 MATERIALS SCIENCE↗