Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallelization”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Criteria for Realizing Room-Temperature Electrical Transport Applications of Topological Materials

The unusual electronic states found in topological materials can enable a new generation of devices and technologies, yet a long-standing challenge has been finding materials without deleterious parallel bulk conduction. This can arise either from defects or thermally activated carriers. In this study, the criteria that materials need to meet to realize transport properties dominated by the topological states, a necessity for a topological device, are clarified. This is demonstrated for 3D topological insulators, 3D Dirac materials, and 1D quantum anomalous Hall insulators, though this can be applied to similar systems. The key parameters are electronic bandgap, dielectric constant, and carrier effective mass, which dictate under what circumstances (defect density, temperature, etc.) the unwanted bulk state will conduct in parallel to the topological states. As these are fundamentally determined by the basic atomic properties, simple chemical arguments can be used to navigate the phase space to ultimately find improved materials. This will enable rapid identification of new systems with improved properties, which is crucial to designing new material systems and push a new generation of topological technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enabling Modular Autonomous Feedback‐Loops in Materials Science through Hierarchical Experimental Laboratory Automation and Orchestration

Abstract Materials acceleration platforms (MAPs) operate on the paradigm of integrating combinatorial synthesis, high‐throughput characterization, automatic analysis, and machine learning. Within a MAP, one or multiple autonomous feedback loops may aim to optimize materials for certain functional properties or to generate new insights. The scope of a given experiment campaign is defined by the range of experiment and analysis actions that are integrated into the experiment framework. Herein, the authors present a method for integrating many actions within a hierarchical experimental laboratory automation and orchestration (HELAO) framework. They demonstrate the capability of orchestrating distributed research instruments that can incorporate data from experiments, simulations, and databases. HELAO interfaces laboratory hardware and software distributed across several computers and operating systems for executing experiments, data analysis, provenance tracking, and autonomous planning. Parallelization is an effective approach for accelerating knowledge generation provided that multiple instruments can be effectively coordinated, which the authors demonstrate with parallel electrochemistry experiments orchestrated by HELAO. Efficient implementation of autonomous research strategies requires device sharing, asynchronous multithreading, and full integration of data management in experimental orchestration, which to the best of the authors’ knowledge, is demonstrated for the first time herein.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Giant Apparent Optical Circular Dichroism in Thin Films of Bismuth-Based Hybrid Organic-Inorganic Metal Halide Semiconductor Through Preferred Orientation

Introducing chirality into organic/inorganic hybrid materials can impart chiroptical properties such as circular dichroism. The ability to tune chiroptical properties in self-assembled materials can have important implications for spintronic and optoelectronic applications. Here, a chiral organic cation, (R/S)-4-methoxy-a-methylbenzylammonium, is incorporated to synthesize the bismuth-based hybrid organic-inorganic metal halide semiconductor, (R/S-MeOMePMA)BiI4. Thin films of this Bi-based compound demonstrate large chiroptical responses, with circular dichroism anisotropy (gCD) values up to ˜0.1, close to the highest value observed in another chiral metal-halide semiconductor, (R-MBA2CuCl4). Detailed investigation reveals that this large gCD in (R/S-MeOMePMA)BiI4 is caused by the apparent CD effect. Careful selection of deposition conditions and the concomitant thin-film orientation enables the control of gCD, with maximum value observed when its thin film has a well-crystallized preferred (001) orientation parallel to the substrate. The results support a growing body of evidence that low symmetry plays an important role in achieving unusually large gCD in these chiral metal-halide materials and provides design rules for achieving large chiroptical response via morphology control.

chiroptic response↗
DOE OSTI2020

Tighter Confinement Increases Selectivity of d-Glucose Isomerization Toward l-Sorbose in Titanium Zeolites

Aqueous-phase isomerization of D-glucose to D-fructose and L-sorbose is catalyzed in parallel by Lewis acidic Ti sites in siliceous frameworks. Glucose isomerization rates (per Ti, 373 K) are undetectable when Ti sites are confined within mesoporous voids (Ti-MCM-41, TiO 2 -SiO 2 ) and increase to detectable values when Ti sites are confined within the smaller 12-membered ring (12-MR) micropores of Ti-Beta. Furthermore, isomerization rates decrease to lower values (by ≈20×) with further decreases in micropore size as Ti sites are confined within 10-MR pores (Ti-MFI, Ti-CON), likely because of intrapore reactant diffusion restrictions, and reach undetectable values within the 8-MR pores of Ti-CHA as size exclusion prevents glucose from accessing active sites. Remarkably, the selectivity toward l -sorbose over d -fructose increases systematically as spatial constraints around Ti sites become tighter, and is >10 on Ti-MFI. These findings demonstrate the marked influence of confinement around Ti active sites on the selectivity between parallel stereoselective sugar isomerization pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A step towards the final frontier: Lessons learned from acceptance testing of the first HPE/Cray EX 3000 system at ORNL

Summary In this article, we summarize the deployment of the Air Force Weather (AFW) HPC11 system at Oak Ridge National Laboratory (ORNL) including the process followed to successfully complete acceptance testing of the system. HPC11 is the first HPE/Cray EX 3000 system that has been successfully released to its user community in a federal facility. HPC11 consists of two identical 800‐node supercomputers, Fawbush and Miller, with access to two independent and identical lustre parallel file systems. HPC11 is equipped with Slingshot 10 interconnect technology and relies on the HPE Performance Cluster Manager software for system configuration. ORNL has a clearly defined acceptance testing process used to ensure that every new system deployed can provide the necessary capabilities to support user workloads. We worked closely with HPE and AFW to develop a set of tests that used the United Kingdom's Meteorological Office's Unified Model and 4‐dimensional variational data assimilation. We also included benchmarks and applications from the Oak Ridge Leadership Computing Facility portfolio to fully exercise the HPE/Cray programming environment and evaluate the functionality and performance of the system. Acceptance testing of HPC11 required parallel execution of each element on Fawbush and Miller. In addition, careful coordination was needed to ensure successful acceptance of the newly deployed lustre file systems alongside the compute resources. In this work, we present test results from specific system components and provide an overview of the issues identified, challenges encountered, and the lessons learned along the way.

Melesse Vergara, Verónica G.↗
DOE OSTI2023

Scaled ILU Smoothers for Navier-Stokes Pressure Projection

Incomplete LU (ILU) smoothers are effective in the algebraic multigrid (AMG) V-cycle for reducing high-frequency components of the error. However, the requisite direct triangular solves are comparatively slow on GPUs. Previous work has demonstrated the advantages of Jacobi iteration as an alternative to direct solution of these systems. Depending on the threshold and fill-level parameters chosen, the factors can be highly nonnormal and Jacobi is unlikely to converge in a low number of iterations. We demonstrate that row scaling can reduce the departure from normality, allowing us to replace the inherently sequential solve with a rapidly converging Richardson iteration. There are several advantages beyond the lower compute time. Scaling is performed locally for a diagonal block of the global matrix because it is applied directly to the factor. Further, an ILUT Schur complement smoother maintains a constant GMRES iteration count as the number of MPI ranks increases, and thus parallel strong-scaling is improved. Our algorithms have been incorporated into hypre, and we demonstrate improved time to solution for linear systems arising in the Nalu-Wind and PeleLM pressure solvers. For large problem sizes, GMRES+AMG executes at least five times faster when using iterative triangular solves compared with direct solves on massively parallel GPUs.

algebraic multigrid↗
DOE OSTI2024

Angle-Resolved Polarized Raman Study of Layered Cr 2 Se 3

The polarization-resolved Raman spectra of two-dimensional Cr 2 Se 3 synthesized via chemical vapor deposition (CVD) and chemical vapor transport (CVT) techniques were investigated in detail. The samples were characterized using X-ray diffraction (XRD), transmission electron microscopy (TEM), and energy-dispersive X-ray spectroscopy (EDS). A distinct polarization dependence was observed in the Raman intensity of all the Cr-Cr, Cr-Se, and Se-Se modes in both samples. The observed angle-dependent Raman intensities of each peak could be related to the crystal structure-specific Raman tensor. XRD results of the bulk Cr 2 Se 3 sample synthesized via CVT confirm its trigonal crystal structure, and the Raman peaks can be fitted using the Raman tensors for the A g and E g modes for both the parallel and crossed polarizations. However, for the Cr 2 Se 3 samples directly grown on Si/SiO 2 substrates by CVD, it was necessary to assume the triclinic crystal structure in order to explain the polarized Raman dependence of all the peaks in both parallel and crossed polarization directions. Furthermore, this is the first experimental result suggesting the existence of triclinic Cr 2 Se 3 crystal structure, which has been theoretically predicted in the Materials Project database.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Scalable Interior‐Point Gauss–Newton Method for PDE‐Constrained Optimization With Bound Constraints

Here, we present a scalable approach to solve a class of partial differential equation (PDE)‐constrained optimization problems with bound constraints. This approach utilizes a robust full‐space interior‐point (IP)‐Gauss–Newton optimization method. To cope with the poorly‐conditioned IP‐Gauss–Newton saddle‐point linear systems that need to be solved approximately, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE‐constrained optimization problems. A block Gauss–Seidel preconditioner is proposed for the GMRES‐based solution of the IP‐Gauss–Newton linear systems. It is shown, for a large‐class of PDE‐ and bound‐constrained optimization problems, that the spectrum of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix is asymptotically independent of discretization and is not impacted by the ill‐conditioning that notoriously plagues interior‐point methods. We exploit symmetry of the IP‐Gauss–Newton linear systems and propose a regularization and log‐barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)‐based solution of the equivalent IP‐Gauss–Newton–Schur complement linear systems. The eigenvalues of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix, that are not equal to one, are identical to the eigenvalues of the regularization and log‐barrier Hessian preconditioned Schur complement matrix. The scalability of the approach is demonstrated on two example problems. The numerical solution of these optimization problems is shown to require a discretization independent number of IP‐Gauss–Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill‐conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of the preconditioner, achieved via algebraic multigrid component solvers when applicable, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE‐ and bound‐constrained problems.

PDE-constrained optimization↗
DOE OSTI2021

Efficient quadrature rules for finite element discretizations of nonlocal equations

Here, in this paper, we design efficient quadrature rules for finite element (FE) discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work, we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first show that the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient FE assembly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Nonlinear multiobjective and dynamic real-time predictive optimization for optimal operation of baseload power plants under variable renewable energy

Considering the increase of disruptive variable renewable energy penetration into the power grid, this article focuses on the investigation of a multiobjective and dynamic real-time optimization framework to address the cycling of large-scale power plants under renewable penetration. In this framework, a parallelized particle swarm optimization step is first performed to generate feasible initial points. Then, a multiobjective and dynamic real-time optimization formulation generates optimal trajectories. Further, the benefit of predictive capability is investigated for the dynamic component, which introduces the novel nonlinear multiobjective and dynamic real-time predictive optimization approach. Two multiobjective formulations to obtain Pareto front optimal in real time are explored: the modified Tchebycheff-based weighted metric and ϵ-constraint methods. Economic and environmental objectives are considered in this study. A novel topical discussion on the intersection of dynamic real-time optimization with model predictive control is also presented. The developed framework is successfully applied to a baseload coal-fired power plant with postcombustion CO 2 capture. Results indicate that the approach can be deployed for a large-scale system if automatic differentiation, model reduction, and parallelization are adopted to improve computational tractability, with computational improvement up to 120-folds after performing these steps. Finally, market and carbon policies showed an impact on the optimal compromise between the objectives with an additional 63 ton of CO 2 captured under favorable market conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Preferred Void Orientation in Uniaxially Pressed PBX 9502

Ultra- and small-angle neutron scattering techniques were used to study the morphology of voids in directions parallel and perpendicular to the compaction plane of a uniaxially pressed cylinder of PBX 9502. Two-dimensional scattering patterns from samples cut parallel to the compaction plane had circular equal-intensity contours, implying that along this view the voids are circular or randomly oriented. Identical measurements of samples cut perpendicular to the compaction plane produced elliptical equal-intensity contours, indicating preferred orientation of voids. Analysis of the two dimensional scattering patterns confirmed that the voids are non-spherical with an average aspect ratio of 1.20, flattened so that the larger dimensions of the voids lie within the compaction plane. Model analysis of the combined ultra- and small-angle neutron scattering curves revealed that the PBX 9502 porosity is interconnected, in the form of aggregates comprised of non-spherical, primary voids. The volume dimension was found to be the same in both directions, indicating that the aggregate is randomly oriented throughout the sample, and that it is the primary voids that are preferentially oriented. Furthermore, void interfaces with the surrounding matrix were found to be rough or polydisperse, with a common roughness dimension or polydispersity in all directions.

42 ENGINEERING↗
DOE OSTI2022

Manipulation of the Magnetoresistance by Strain in Topological TaSe 3

1D TaSe 3 exhibits many unusual physical properties due to its distorted type-II chains. Ribbon-shaped single crystals can be easily bended along the $b-axis$, forming rings. This study investigates the magnetoresistance (MR) of TaSe 3 up to 60 T in both unbended (ribbon shape) and ring-shaped (bended ribbon) samples. Notable changes are found in the magnetotransport properties between the two different shaped samples. One is that the MR in ring-shaped samples is three orders lower than that in unbended samples under the same sample environment. In addition, linear MR is observed above ≈20 T in ring-shaped samples when the magnetic field is parallel or perpendicular to the rings. Quantum oscillations are also observed as a function of the magnetic field when the magnetic field is applied parallel to rings, possibly due to the Altshuler–Aronov–Spivak effect or the inversion of the lowest Landau level beyond the quantum limit. All these results are related to strain-induced electronic structure change in TaSe 3 , an effective way to tune physical properties in low-dimensional materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Rule of Thirds: Controlling Junction Chirality and Polarity in 3D DNA Tiles

The successful self-assembly of tensegrity triangle DNA crystals heralded the ability to programmably construct macroscopic crystalline nanomaterials from rationally-designed, nanoscale components. This 3D DNA tile owes its “tensegrity” nature to its three rotationally-stacked double helices locked together by the tensile winding of a center strand segmented into 7 base pair (bp) inter-junction regions, corresponding to two thirds of a helical turn of DNA. All reported tensegrity triangles to date have employed [Z + 2/3] turn inter-junction segments, yielding right-handed, antiparallel, “J1” junctions. Here, we report here a minimal DNA triangle motif consisting of 3-bp inter-junction segments, or one third of a helical turn. Here, we find that the minimal motif exhibits a reversed morphology with a left-handed tertiary structure mediated by a locally-parallel Holliday junction—the L1 junction. This parallel junction yields a predicted helical groove matching pattern that breaks the pseudosymmetry between tile faces, and the junction morphology further suggests a folding mechanism. We identify a Rule of Thirds by which supramolecular chirality can be programmed through inter-junction DNA segment length. These results underscore the role that global topological forces play in determining local DNA architecture and ultimately point to an under-explored class of self-assembling, chiral nanomaterials for topological processes in biological systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Surface Effects on Anisotropic Photoluminescence in One‐Dimensional Organic Metal Halide Hybrids

1D organic metal halide hybrids (OMHHs) exhibit strongly anisotropic optical properties, highly efficient light emission, and large Stokes shift, holding promise for novel photodetection and lighting applications. However, the fundamental mechanisms governing their unique optical properties and in particular the impacts of surface effects are not understood. Herein, 1D C 4 N 2 H 14 PbBr 4 by polarization‐dependent time‐averaged and time‐resolved photoluminescence (TRPL) spectroscopy, as a function of photoexcitation energy, is investigated. Surprisingly, it is found that the emission under photoexcitation polarized parallel to the 1D metal halide chains can be either stronger or weaker than that under perpendicular polarization, depending on the excitation energy. The excitation‐energy‐dependent anisotropic emission is attributed to fast surface recombination, supported by first‐principles calculations of optical absorption in this material. The fast surface recombination is directly confirmed by TRPL measurements, when the excitation is polarized parallel to the chains. The comprehensive studies provide a more complete picture for a deeper understanding of the optical anisotropy in 1D OMHHs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A localized ensemble of approximate Gaussian processes for fast sequential emulation

More attention has been given to the computational cost associated with the fitting of an emulator. Substantially less attention is given to the computational cost of using that emulator for prediction. This is primarily because the cost of fitting an emulator is usually far greater than that of obtaining a single prediction, and predictions can often be obtained in parallel. In many settings, especially those requiring Markov Chain Monte Carlo, predictions may arrive sequentially and parallelization is not possible. In this case, using an emulator procedure which can produce accurate predictions efficiently can lead to substantial time savings in practice. In this paper, we propose a global model approximate Gaussian process framework via extension of a popular local approximate Gaussian process (laGP) framework. Our proposed emulator can be viewed as a treed Gaussian process where the leaf nodes are laGP models, and the tree structure is learned greedily as a function of the prediction stream. The suggested method (called leapGP) has interpretable tuning parameters which control the time‐memory trade‐off. One reasonable choice of settings leads to an emulator with a training cost and makes predictions rapidly with an asymptotic amortized cost of .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fragme∩t: An Open‐Source Framework for Multiscale Quantum Chemistry Based on Fragmentation

Fragment-based quantum chemistry offers a means to circumvent the nonlinear computational scaling of conventional electronic structure calculations, by partitioning a large calculation into smaller subsystems then considering the many-body interactions between them. Variants of this approach have been used to parameterize classical force fields and machine learning potentials, applications that benefit from interoperability between quantum chemistry codes. However, there is a dearth of software that provides interoperability yet is purpose-built to handle the combinatorial complexity of fragment-based calculations. To fill this void we introduce “Fragme∩t”, an open-source software application that provides a tool for community validation of fragment-based methods, a platform for developing new approximations, and a framework for analyzing many-body interactions. Fragme∩t includes algorithms for automatic fragment generation and structure modification, and for distance- and energy-based screening of the requisite subsystems. Checkpointing, database management, and parallelization are handled internally and results are archived in a portable database. Interfaces to various quantum chemistry engines are easy to write and exist already for Q-Chem, PySCF, xTB, Orca, CP2K, MRCC, Psi4, NWChem, GAMESS, and MOPAC. Applications reported here demonstrate parallel efficiencies around 96% on more than 1000 processors but also showcase that the code can handle large-scale protein fragmentation using only workstation hardware, all with a codebase that is designed to be usable by non-experts. Fragme∩t conforms to modern software engineering best practices and is built upon well established technologies including Python, SQLite, and Ray. The source code is available under the Apache 2.0 license.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Proximity Portability and in Transit , M-to-N Data Partitioning and Movement in SENSEI [Book Chapter]

In high-performance parallel in situ processing, the term in transit processing refers to those configurations where data must move from a producer to a consumer that runs on separate resources. In the context of parallel and distributed computing on an HPC platform one of the central challenges is to determine a mapping of data from producer ranks to consumer ranks. This problem is complicated by the heterogeneity that arises in producer-consumer pairs, such as when producer and consumer codes have different levels of concurrency, different scaling characteristics, or different data models. The resulting mapping and movement of data from M producer to N consumer ranks can have a significant impact on aggregate application performance, particularly when the data consumer requires only a subset of the overall data for its task. This chapter focuses on the design considerations that underlie SENSEI’s implementation to this challenging problem. These design considerations extend the core SENSEI architecture and include ideas like the need to accommodate flexibility in the choice of different partitioning methods, the ability for a data consumer to request and receive only the subset of data needed for its particular operation, and the ability to leverage any of several different data transport tools. The idea of proximity portability, being able to use different data transport methods as part of an in transit workflow, is illustrated through the use of three different transport layers where switching from one transport tool to another is accomplished with only a configuration file change. Here, the chapter also includes a performance analysis summary showing the performance gains that are possible in terms of multiple metrics, such as memory footprint, time to solution, and amount of data moved, when using optimized partitioners in an in transit setting, gains that are made possible by the implementation shaped by specific design considerations.

Bethel, E. Wes↗
DOE OSTI2024

CI/CD Efforts for Validation, Verification and Benchmarking OpenMP Implementations

Software developers must adapt to keep up with the changing capabilities of platforms so that they can utilize the power of High-Performance Computers (HPC), including exascale systems. OpenMP, a directive-based parallel programming model, allows developers to include directives to existing C, C++, or Fortran code to allow node level parallelism without compromising performance. This paper describes our CI/CD efforts to provide easy evaluation of the support of OpenMP across different compilers using existing testsuites and benchmark suites on HPC platforms. Our main contributions include (1) the set of a Continuous Integration (CI) and Continuous Development (CD) workflow that captures bugs and provides faster feedback to compiler developers, (2) an evaluation of OpenMP (offloading) implementations supported by AMD, HPE, GNU, LLVM, and Intel, and (3) evaluation of the quality of compilers across different heterogeneous HPC platforms. With the comprehensive testing through the CI/CD workflow, we aim to provide a comprehensive understanding of the current state of OpenMP (offloading) support in different compilers and heterogeneous platforms consisting of CPUs and GPUs from NVIDIA, AMD, and Intel.

Jarmusch, Aaron↗