Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel Architectures”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Enabling Execution of a Legacy CFD Mini Application on Accelerators Using OpenMP

We describe the process and outcome of our efforts to port a legacy Fortran benchmark code to heterogeneous GPU-accelerated computing architectures using OpenMP. The benchmark code is one of the multi-zone NAS Parallel Benchmarks (NPB-MZ) called SP-MZ. This “mini-app” mimics the computation and data movement that is found in popular legacy and modern implicit computational fluid dynamics (CFD)solvers. Our objective was to examine how efficiently legacy Fortran codes can be ported to accelerators by leveraging OpenMP directives. We describe the development and optimization process and demonstrate the performance impact of various code modifications. We show select profiling results from the Nvidia nvvp profiler to help others diagnose and overcome performance issues in their own applications. We present results for two compute systems endowed with Nvidia V100 accelerators.

Ioannis Nompelis↗
DOE OSTI2020

ExaHDF5: Delivering Efficient Parallel I/O on Exascale Computing Systems

Scientific applications at exascale generate and analyze massive amounts of data. A critical requirement of these applications is the capability to access and manage this data efficiently on exascale systems. As such, parallel I/O, the key technology enables moving data between compute nodes and storage, faces monumental challenges from new applications, memory, and storage architectures considered in the designs of exascale systems. As the storage hierarchy is expanding to include node-local persistent memory, burst buffers, etc., as well as disk-based storage, data movement among these layers must be efficient. Parallel I/O libraries of the future should be capable of handling file sizes of many terabytes and beyond. In this paper, we describe new capabilities we have developed in Hierarchical Data Format version 5 (HDF5), the most popular parallel I/O library for scientific applications. HDF5 is one of the most used libraries at the leadership computing facilities for performing parallel I/O on existing HPC systems. The state-of-the-art features we describe include: Virtual Object Layer (VOL), Data Elevator, asynchronous I/O, full-featured single-writer and multiple-reader (Full SWMR), and parallel querying. In this paper, we introduce these features, their implementations, and the performance and feature benefits to applications and other libraries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The Kokkos Ecosystem [Brief]

In 2016/2017, the field of High-Performance Computing (HPC) entered a new era driven by fundamental physics challenges to produce ever more energy and cost-efficient processors. Since the convergence on the Message-Passing Interface (MPI) standard in the mid-1990s, application developers enjoyed a seemingly static view of the underlying machine — that of a distributed collection of homogeneous nodes executing in collaboration. However, after almost two decades of dominance, the sole use of MPI to derive parallelism acted as a limiter to improved future performance. While MPI is widely expected to continue to function as the basic mechanism for communication between compute nodes for the immediate future, additional parallelism is required on the computing node itself if high performance and efficiency goals are to be realized. When reviewing the architectures of the top HPC systems today, the change in paradigm is clear: the compute nodes of the leading machines in the world are either powered by many-core chips with a few dozen cores each, or use heterogeneous designs, where traditional CPUs marshal work to massively parallel compute accelerators which has as many as 200,000 processing threads in flight simultaneously. Complicating matters further for application developers, each processor vendor has its own preferred way of writing code for their architecture.The Kokkos EcoSystem was released by Sandia in 2017 to address this new era in HPC system design by providing a vendor independent performance portable programming system for scientific, engineering, and mathematical software applications written in the C++ programming language. Using Kokkos, application developers can be more productive because they will not have to create and maintain separate versions of their software for each architecture, nor will they have to be experts in each architecture's peculiar requirements. Instead, they will have a single method of programming for the diverse set of modern HPC architectures. While Kokkos started in 2011 as a programming model only, it soon became clear that complex applications needed more. It is also critical to have a portable mathematical functions and developers need tools to debug their applications, gain insight into the performance characteristics of their codes and tune algorithm performance parameters through automated processes. The Kokkos EcoSystem addresses those needs through its three main components: the Kokkos Core programming model, the Kokkos Kernels math library, and the Kokkos Tools project.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

A system for routing arbitrary directed graphs on SIMD architectures

There are many problems which can be described in terms of directed graphs that contain a large number of vertices where simple computations occur using data from connecting vertices. A method is given for parallelizing such problems on an SIMD machine model that is bit-serial and uses only nearest neighbor connections for communication. Each vertex of the graph will be assigned to a processor in the machine. Algorithms are given that will be used to implement movement of data along the arcs of the graph. This architecture and algorithms define a system that is relatively simple to build and can do graph processing. All arcs can be transversed in parallel in time O(T), where T is empirically proportional to the diameter of the interconnection network times the average degree of the graph. Modifying or adding a new arc takes the same time as parallel traversal.

Tomboulian, Sherryl↗
NASA NTRS1999

Integrated Space Transportation Plan: Defining Technology Requirements and Next Generation Launch Systems to Meet Commercial and Government Needs

NASA continues to focus on improving safety and reliability while reducing the annual cost of meeting human space flight and unique ISS and exploration needs. NASA's Space Transportation Architecture Study (STAS) Phase 2 in early 1998 focused on space transportation options. Subsequently, NASA directed parallel industry and government teams to conduct the Integrated Space Transportation Plan effort (STAS Phase 3). The objective of ISTP was to develop technology requirements, roadmaps, and risk reduction portfolio that considered expanded definition of "clean-sheet" and Shuttle-derived second generation ETO transportation systems in support of a 2005 RLV competition for NASA missions beginning 2010. NASA provided top-level requirements for improvements in safety, reliability, and cost and a set of design reference missions representing NASA ISS, human exploration, commercial, and other civil and government needs. This paper addresses the challenges of meeting NASA's objectives while servicing the varied market segments represented in the ISTP design reference missions and provides a summary of technology development needs and candidate system concepts. A comparison of driving requirements, architectures and technology needs is discussed and descriptions of viable Shuttle-derived and next generation systems to meet the market needs are presented.

Davidoff, Larry D.↗
NASA NTRS1983

Parallelism and pipelining in high-speed digital simulators

The attainment of high computing speed as measured by the computational throughput is seen as one of the most challenging requirements. It is noted that high speed is cardinal in several distinct classes of applications. These classes are then discussed; they comprise (1) the real-time simulation of dynamic systems , (2) distributed parameter systems, and (3) mixed lumped and distributed systems. From the 1950s on, the quest for high speed in digital simulators concentrated on overcoming the limitations imposed by the so-called von Neumann bottleneck. Two major architectural approaches have made ig possible to circumvent this bottleneck and attain high speeds. These are pipelining and parallelism. Supercomputers, peripheral array processors, and microcomputer networks are then discussed.

Karplus, W. J.↗
NASA NTRS2004

Paramedir: A Tool for Programmable Performance Analysis

Performance analysis of parallel scientific applications is time consuming and requires great expertise in areas such as programming paradigms, system software, and computer hardware architectures. In this paper we describe a tool that facilitates the programmability of performance metric calculations thereby allowing the automation of the analysis and reducing the application development time. We demonstrate how the system can be used to capture knowledge and intuition acquired by advanced parallel programmers in order to be transferred to novice users.

Jost, Gabriele↗
NASA NTRS2003

Matrix Cracking in 3D Orthogonal Melt-Infiltrated SiC/SiC Composites with Various Z-Fiber Types

The occurrence of matrix cracks in melt-infiltrated SiC/SiC composites with a 3D orthogonal architecture was determined at room temperature for specimens tested in tension oriented in the X-direction (parallel to Z-bundle weave direction) and Y-direction (perpendicular to Z-bundle weave direction) and Y-direction (perpendicular to Z-bundle weave direction). The fiber-types were Sylramic and Sylramic-IBN in the X and Y-directions and lower modulus ZMI, T300, and rayon in the Z-direction. Acoustic emission (AE) was used to monitor the matrix cracking activity. For Y-direction composites, the AE data was used to determine the exact (+/- 0.25 mm) location where matrix cracks occurred in the 3D orthogonal architecture. This enabled the determination of the stress-dependent matrix crack distributions for small but repeatable matrix rich 'unidirectional' and the matrix poor 'cross-ply' regions within the architecture. It was found that matrix cracking initiated at very low stresses (approx. 40 MPa) in the 'unidirectional' regions for the largest z-direction fiber tow composites. Decreasing the size of the z-fiber bundle, increased the stress for matrix cracking in the 'unidirectional' regions. Matrix cracking in the 'cross-ply' regions always occurred at higher stresses than in 'unidirectional' regions, and the stress-dependent matrix crack distribution of the 'cross-ply' regions was always over a wider stress-range than the 'unidirectional' regions. For composites tested in the X-direction, a lower elastic modulus and a narrower and lower stress-range for matrix cracking were observed compared to composites tested in the Y-direction.

Morscher, Gregory N.↗
DOE OSTI2025

Optical neural engine for solving scientific partial differential equations

Abstract Solving partial differential equations (PDEs) is the cornerstone of scientific research and development. Data-driven machine learning (ML) approaches are emerging to accelerate time-consuming and computation-intensive numerical simulations of PDEs. Although optical systems offer high-throughput and energy-efficient ML hardware, their demonstration for solving PDEs is limited. Here, we present an optical neural engine (ONE) architecture combining diffractive optical neural networks for Fourier space processing and optical crossbar structures for real space processing to solve time-dependent and time-independent PDEs in diverse disciplines, including Darcy flow equation, the magnetostatic Poisson’s equation in demagnetization, the Navier-Stokes equation in incompressible fluid, Maxwell’s equations in nanophotonic metasurfaces, and coupled PDEs in a multiphysics system. We numerically and experimentally demonstrate the capability of the ONE architecture, which not only leverages the advantages of high-performance dual-space processing for outperforming traditional PDE solvers and being comparable with state-of-the-art ML models but also can be implemented using optical computing hardware with unique features of low-energy and highly parallel constant-time processing irrespective of model scales and real-time reconfigurability for tackling multiple tasks with the same architecture. The demonstrated architecture offers a versatile and powerful platform for large-scale scientific and engineering computations.

Tang, Yingheng (ORCID:0009000153622546)↗
NASA NTRS1989

A language comparison for scientific computing on MIMD architectures

Choleski's method for solving banded symmetric, positive definite systems is implemented on a multiprocessor computer using three FORTRAN based parallel programming languages, the Force, PISCES and Concurrent FORTRAN. The capabilities of the language for expressing parallelism and their user friendliness are discussed, including readability of the code, debugging assistance offered, and expressiveness of the languages. The performance of the different implementations is compared. It is argued that PISCES, using the Force for medium-grained parallelism, is the appropriate choice for programming Choleski's method on the multiprocessor computer, Flex/32.

Jones, Mark T.↗
NASA NTRS1996

An Element-Based Concurrent Partitioner for Unstructured Finite Element Meshes

A concurrent partitioner for partitioning unstructured finite element meshes on distributed memory architectures is developed. The partitioner uses an element-based partitioning strategy. Its main advantage over the more conventional node-based partitioning strategy is its modular programming approach to the development of parallel applications. The partitioner first partitions element centroids using a recursive inertial bisection algorithm. Elements and nodes then migrate according to the partitioned centroids, using a data request communication template for unpredictable incoming messages. Our scalable implementation is contrasted to a non-scalable implementation which is a straightforward parallelization of a sequential partitioner.

concurrent partitioner finite element meshes distr↗
DOE OSTI2022

Molybdenum Carbide Electrocatalyst In Situ Embedded in Porous Nitrogen–Rich Carbon Nanotubes Promotes Rapid Kinetics in Sodium–Metal–Sulfur Batteries

This work is the first report of a molybdenum carbide-based electrocatalyst for sulfur-based sodium metal batteries (SMBs/NMBs). MoC/Mo 2 C is in-situ grown on nitrogen-doped carbon nanotubes in parallel with formation of extensive nanoporosity. Sulfur impregnation (50 wt% S) results in unique triphasic architecture termed MoC/Mo 2 C@PCNT-S. Quasi-solid-state phase transformation to Na 2 S is promoted in carbonate electrolyte, with in-situ time-resolved Raman, XPS and optical analysis demonstrating minimal soluble polysulfides. MoC/Mo 2 C@PCNT-S cathodes delivered among the most promising rate performance characteristics in literature, achieving 987 mAh g -1 at 1 Ag -1 , 818 mAh g -1 at 3 A g -1 , and 621 mAh g -1 at 5 A g -1 . The cells deliver superior cycling stability, retaining 650 mAh g -1 after 1000 cycles at 1.5 Ag -1 , corresponding to 0.028% capacity decay per cycle. High mass loading cathodes (64 wt% S, 12.7 mg cm -2 ) also show cycling stability, with anode degradation due to deep plating/stripping driving capacity decay. Density functional theory (DFT) demonstrates that formation energy of Na 2 S x (1 ≤ x ≤ 4) on surface of MoC/Mo 2 C is significantly lowered compared to analogous redox in liquid. Strong binding of Na 2 S x (1 ≤ x ≤ 4) on MoC/Mo 2 C surfaces results from charge transfer between the sulfur and Mo sites on carbides' surface.

25 ENERGY STORAGE↗
DOE OSTI2020

Scalable Data-Intensive Geocomputation: A Design for Real-Time Continental Flood Inundation Mapping

The convergence of data-intensive and extreme-scale computing enables an integrated software and data ecosystem for scientific discovery. Developments in this realm will fuel transformative research in data-driven interdisciplinary domains. Geocomputation provides computing paradigms in Geographic Information Systems (GIS) for interactive computing of geographic data, processes, models, and maps. Because GIS is data-driven, the computational scalability of a geocomputation workflow is directly related to the scale of the GIS data layers, their resolution and extent, as well as the velocity of the geo-located data streams to be processed. Unique in high user interactivity and low end-to-end latency requirements, geocomputation applications will dramatically benefit from the convergence of high-end data analytics (HDA) and high-performance computing (HPC). The application level challenge, however, is to identify and eliminate computational bottlenecks that arise along a geocomputation workflow. Indeed, poor scalability at any of the workflow components is detrimental to the entire end-to-end pipeline. Here, we study a large geocomputation use case in flood inundation mapping that handles multiple national-scale geospatial datasets and targets low end-to-end latency. We discuss benefits and challenges for harnessing both HDA and HPC for data-intensive geospatial data processing and intensive numerical modeling of geographic processes. We propose an HDA+HPC geocomputation architecture design that couples HDA (e.g., Spark)-based spatial data handling and HPC-based parallel data modeling. Key techniques for coupling HDA and HPC to bridge the two different software stacks are reviewed and discussed.

Liu, Yan↗
DOE OSTI2024

The data acquisition system of the LZ dark matter detector: FADR

The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals. This information is used to determine if the digitized waveforms should be preserved for offline analysis. The system is designed around the Kintex-7 FPGA. In addition to digitizing the PMT signals and providing basic event selection in real time, the flexibility provided by the use of FPGAs allows us to monitor the performance of the detector and the DAQ in parallel to normal data acquisition. Furthermore, the hardware and software/firmware of this FPGA-based Architecture for Data acquisition and Realtime monitoring (FADR) are discussed and performance measurements are described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

PyHydroGeophysX: An extensible open-source platform for integrating hydrological models with geophysical measurements

Hydrological models and geophysical measurements are widely used tools for understanding subsurface hydrological processes relevant to water resource management, yet they typically remain disconnected due to technical barriers. We present PyHydroGeophysX, an open-source Python platform bridging this gap by providing standardized interfaces between hydrological modeling software (MODFLOW, ParFlow) and geophysical simulation tools (PyGIMLi, SimPEG). The platform implements bidirectional workflows: translating hydrological outputs into simulated geophysical responses through petrophysical models, and extracting hydrological information from geophysical inversions. Key features include bidirectional workflow modules, configurable petrophysical models, time-lapse inversion with temporal regularization, parallel computing, and mesh utilities for property transfer between geophysical and hydrological grids. The modular architecture of PyHydroGeophysX enables researchers to incorporate additional models and methods, fostering broader adoption of integrated hydrogeophysical approaches. The software is freely available on GitHub and is intended for researchers and practitioners working at the intersection of hydrology and geophysics.

Hydrogeophysics↗
DOE OSTI2022

GSoFa: Scalable Sparse Symbolic LU Factorization on GPUs

Decomposing a matrix $\mathbf {A}$ into a lower matrix $\mathbf {L}$ and an upper matrix $\mathbf {U}$, which is also known as LU decomposition, is an essential operation in numerical linear algebra. For a sparse matrix, LU decomposition often introduces more nonzero entries in the $\mathbf {L}$ and $\mathbf {U}$ factors than in the original matrix. A symbolic factorization step is needed to identify the nonzero structures of $\mathbf {L}$ and $\mathbf {U}$ matrices. Attracted by the enormous potentials of the Graphics Processing Units (GPUs), an array of efforts have surged to deploy various LU factorization steps except for the symbolic factorization, to the best of our knowledge, on GPUs. This article introduces gSoFa, the first GPU-based symbolic factorization design with the following three optimizations to enable scalable LU symbolic factorization for nonsymmetric pattern sparse matrices on GPUs. First, here we introduce a novel fine-grained parallel symbolic factorization algorithm that is well suited for the Single Instruction Multiple Thread (SIMT) architecture of GPUs. Second, we tailor supernode detection into a SIMT friendly process and strive to balance the workload, minimize the communication and saturate the GPU computing resources during supernode detection. Third, we introduce a three-pronged optimization to reduce the excessive space consumption problem faced by multi-source concurrent symbolic factorization. Taken together, gSoFa achieves up to 31× speedup from 1 to 44 Summit nodes (6 to 264 GPUs) and outperforms the state-of-the-art CPU project, on average, by 5×. Notably, gSoFa also achieves up to 47 percent of the peak memory throughput of a V100 GPU in the Summit Supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Creep and Fatigue Characterization of High Strength Alloy Thin Sections in Advanced CO2 Heat Exchangers

The objective of this work was to characterize and model elevated temperature creep and fatigue behavior for thin sheet and foil forms of gamma-prime strengthened alloys in wrought form and as-processed folded and brazed constructions. This work was motivated by the demanding temperature and pressure service conditions of the GEN3 Concentrated Solar Power (CSP) and supercritical CO2 (sCO2) power cycle working fluid. More specifically, the possibility of leveraging the superior creep strength of gamma-prime alloys in folded-fin and brazed-plate heat exchanger constructions. Gamma-prime alloys represent a step-change in raw-material strength over solid-solution strengthened alloys. And the folded-fin and brazed-plate heat exchanger architecture is lightweight and leverages cost-effective material stock forms. The investigation contained two parallel paths. (1) The first is referred to as a fundamental investigation where Oak Ridge National Laboratory conducts uniaxial creep testing on thin sheet and foil in wrought form. This effort aimed to serve as a benchmark against a relatively sparse existing database and a baseline comparison for path number 2. (2) The second path is referred to as the practical investigation where Brayton Energy manufactures plate-fin heat exchangers and performs pressurized creep and fatigue testing. This effort aimed to de-risk heat exchanger manufacturing process for service under sCO2 CSP conditions. A total of 14 uniaxial creep tests were completed using Haynes 282 thin sheet and foil. A variety of heat treatments were specified to coincide with path number 2. Baseline metallography of test samples and creep strength performance are contained. Benchmarks relatively to existing thick-form Haynes 282 are made, as well as to other thin-form Nickel-based superalloys. Description of a wrought-form modeling approach for thin gamma-prime alloys is also discussed. A total of 11 pressurized creep and fatigue tests were completed successfully with Haynes 282 heat exchanger prototypes. Manufacturing processing details, testing details, testing results, and failure analysis are discussed. Additionally, creep modeling techniques to predict failure are discussed, and modeling to support technological-to-market. In conclusion, Haynes 282 foils were demonstrated to yield rupture two-to-three orders of magnitude higher than Haynes 230 foils under similar conditions. And the manufactured heat exchanger prototypes demonstrated strength similar to the wrought constituents. Both of which contribute to elevated performance potential or cost savings in practice. Discussion is included.

14 SOLAR ENERGY↗