Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallelization efficiency”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Kokkos OpenMPTarget Backend: Implementation and Lessons Learned

As the supercomputing landscape diversifies, solutions such as Kokkos to write vendor agnostic applications and libraries have risen in popularity. Kokkos provides a programming model designed for performance portability, which allows developers to write a single source implementation that can run efficiently on various architectures. At its heart, Kokkos maps parallel algorithms to architecture and vendor specific backends written in lower level programming models such as CUDA and HIP. Another approach to writing vendor agnostic parallel code is using OpenMP’s directives based approach, which lets developers annotate code to express parallelism. It is implemented at the compiler level and is supported by all major high performance computing vendors, as well as the primary Open Source toolchains GNU and LLVM. Since its inception, Kokkos has used OpenMP to parallelize on CPU architectures. In this paper, we explore leveraging OpenMP for a GPU backend and discuss the challenges we encountered when mapping the Kokkos APIs and semantics to OpenMP target constructs. As an exemplar workload we chose a simple conjugate gradient solver for sparse matrices. We find that performance on NVIDIA and AMD GPUs varies widely based on details of the implementation strategy and the chosen compiler. Furthermore, the performance of the OpenMP implementations decreases with increasing complexity of the investigated algorithms.

Gayatri, Rahulkumar↗
DOE OSTI2019

Investigation of resonant mode sensitivity in PhC based ring resonators

Quality factor and refractive index sensitivity are significant parameters in designing optical devices based on ring resonators such as filters, demultiplexers, switches in which the ring resonator is considered as an elementary componant. In this paper, we propose a photonic crystal ring resonator with super-ellipse shaped core. The proposed ring resonator placed between two parallels waveguides (bus waveguide and drop waveguide). The transmission efficiency of the proposed ring resonator at λ = 1547.3 nm is about 95% with bandwidth and quality factor values equal to 0.7 nm and 2210, respectively. Simulation results show that the resonance sensitivity upon refractive index variation of coupling, outer layer, adjacent and scattering rods are Δλ/Δn = 0.2 nm/0.01, Δλ/Δn = 0.225 nm/0.01, Δλ/Δn = 0.125 nm/0.01 and Δλ/Δn = 0.05 nm/0.01, respectively. The finite different time domain and plane wave expansion methods are used to calculate the outputs spectrum and band gap, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Experimental Observation of the Effects of Translational and Rotational Electrode Misalignment on a Planar Linear Ion Trap Mass Spectrometer

The performance of miniaturized ion trap mass analyzers is limited, in part, by the accuracy with which electrodes can be fabricated and positioned relative to each other. Alignment of plates in a two-plate planar LIT is ideal to characterize misalignment effects, as it represents the simplest possible case, having only six degrees of freedom (DOF) (three translational and three rotational). High-precision motorized actuators were used to vary the alignment between the two ion trap plates in five DOFs—x, y, z, pitch, and yaw. A comparison between the experiment and previous simulations shows reasonable agreement. Pitch, or the degree to which the plates are parallel along the axial direction, has the largest and sharpest impact to resolving power, with resolving power dropping noticeably with pitch misalignment of a fraction of a degree. Lateral displacement (x) and yaw (rotation of one plate, but plates remain parallel) both have a strong impact on ion ejection efficiency, but little effect on resolving power. The effects of plate spacing (y-displacement) on both resolving power and ion ejection efficiency are attributable to higher-order terms in the trapping field. Varying the DC (axial) trapping potential can elucidate the effects where more misalignments in more than one DOF affect performance. Implications of these results for miniaturized ion traps are discussed. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗
DOE OSTI2019

Closed-loop simulation with Grad-Shafranov equilibrium evolution for plasma control system development

Plasma control design increasingly depends on fast simulations able to connect to operational plasma control system (PCS) software for iterative development. GSevolve is a nonlinear free-boundary simulation that evolves the Grad-Shafranov equilibrium including current and pressure profiles, and can connect to all versions of the DIII-D PCS operational on devices around the world. Its ability to simulate the evolution of current continuously from open to closed field lines enables GSevolve to simulate full-discharge control from breakdown through plasma termination. GSevolve is compliant with the Plasma Control System Simulation Platform (PCSSP), a Matlab/Simulink®-based environment that serves as the primary framework for development of the ITER Plasma Control System (PCS). The code employs a novel response function-updating scheme to achieve fast simulations. The plasma evolves linearly until a significant shape change triggers a nonlinear update of the response while enforcing a prescribed accuracy. Old responses are recalled if the plasma changes back to an old shape. The equilibrium is a function of a state vector containing axisymmetric conductor currents and parameters related to parallel and perpendicular currents, which evolves according to an efficient state space circuit formulation. Control design and validation results are presented for simulations of DIII-D and EAST, connected to their corresponding versions of the DIII-D PCS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A remarkable CHF of 345W/cm 2 is achieved in a wicked-microchannel using HFE-7100

Boiling heat transfer of dielectric fluids is a promising cooling technique for thermal management of microelectronic systems. However, the critical heat flux achievable is generally low, because of the poor thermophysical properties of these fluids. Here, to address this dilemma, we propose a new cooling concept to substantially enhance global liquid supply during phase-change process using enhanced capillary-driven force. Additionally, dedicated vapor pathways are designed among a bank of micro-pillars to facilitate vapor removal. In this work, new wicks comprised of silicon micro-pinfin arrays are explored to significantly enhance the flow boiling heat transfer performance. To examine the functionalities of this wick, experiments on HFE-7100 were carried out with mass velocities varying from 247 to 3,465 kg/m 2 s. To explore the enhancement mechanisms and to analyze the capillary-assisted flow boiling process, visualization studies were conducted. The results indicate that sustainable evaporation induced by wick microstructures and efficient liquid supply are the enhancement mechanisms compared to parallel microchannels with solid walls. It is found that the overall heat transfer coefficient is substantially increased up to 75%. Remarkably, a high critical heat flux (CHF) of approximately 345 W/cm 2 is recorded at G = 3,465 kg/m 2 s at coolant inlet temperature of ~20 °C. Equally importantly, this noticeable enhancement of CHF value is associated with drastically decreased pressure drops compared to microchannels decorated with μ-pinfin fences.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Chebyshev polynomial filtered subspace iteration in the discontinuous Galerkin method for large-scale electronic structure calculations

The Discontinuous Galerkin (DG) electronic structure method employs an adaptive local basis (ALB) set to solve the Kohn-Sham equations of density functional theory in a discontinuous Galerkin framework. The adaptive local basis is generated on-the-fly to capture the local material physics and can systematically attain chemical accuracy with only a few tens of degrees of freedom per atom. A central issue for large-scale calculations, however, is the computation of the electron density (and subsequently, ground state properties) from the discretized Hamiltonian in an efficient and scalable manner. We show in this work how Chebyshev polynomial filtered subspace iteration (CheFSI) can be used to address this issue and push the envelope in large-scale materials simulations in a discontinuous Galerkin framework. We describe how the subspace filtering steps can be performed in an efficient and scalable manner using a two-dimensional parallelization scheme, thanks to the orthogonality of the DG basis set and block-sparse structure of the DG Hamiltonian matrix. The on-the-fly nature of the ALB functions requires additional care in carrying out the subspace iterations. We demonstrate the parallel scalability of the DG-CheFSI approach in calculations of large-scale twodimensional graphene sheets and bulk three-dimensional lithium-ion electrolyte systems. In conclusion, employing 55 296 computational cores, the time per self-consistent field iteration for a sample of the bulk 3D electrolyte containing 8586 atoms is 90 s, and the time for a graphene sheet containing 11 520 atoms is 75 s.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

P SI4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree–Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4 's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. We announce a rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

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

A lightweight, user-configurable detector ASIC digital architecture with on-chip data compression for MHz X-ray coherent diffraction imaging

Today, most X-ray pixel detectors used at light sources transmit raw pixel data off the detector ASIC. With the availability of more advanced ASIC technology nodes for scientific application, more digital functionalities from the computing domains (e.g., compression) can be integrated directly into a detector ASIC to increase data velocity. In this paper, we describe a lightweight, user-configurable detector ASIC digital architecture with on-chip compression which can be implemented in 130 nm technologies in a reasonable area on the ASIC periphery. In addition, we present a design to efficiently handle the variable data from the stream of parallel compressors. The architecture includes user-selectable lossy and lossless compression blocks. The impact of lossy compression algorithms is evaluated on simulated and experimental X-ray ptychography datasets. This architecture is a practical approach to increase pixel detector frame rates towards the continuous 1 MHz regime for not only coherent imaging techniques such as ptychography, but also for other diffraction techniques at X-ray light sources.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Computable and Operationally Meaningful Multipartite Entanglement Measures

Multipartite entanglement is an essential resource for quantum communication, quantum computing, quantum sensing, and quantum networks. The utility of a quantum state |ψ$\rangle$ for these applications is often directly related to the degree or type of entanglement present in |ψ$\rangle$. Therefore, efficiently quantifying and characterizing multipartite entanglement is of paramount importance. Here, in this work, we introduce a family of multipartite entanglement measures, called concentratable entanglements. Several well-known entanglement measures are recovered as special cases of our family of measures, and hence we provide a general framework for quantifying multipartite entanglement. We prove that the entire family does not increase, on average, under local operations and classical communications. We also provide an operational meaning for these measures in terms of probabilistic concentration of entanglement into Bell pairs. Finally, we show that these quantities can be efficiently estimated on a quantum computer by implementing a parallelized SWAP test, opening up a research direction for measuring multipartite entanglement on quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

DIII-D Electron Cyclotron Heating and Current Drive System Status and Plans

The DIII-D 110 GHz Electron Cyclotron Heating and Current Drive (ECH/ECCD) System supports a wide range of experiments at the DIII-D National Fusion Facility and is pioneering innovative new techniques for future fusion reactors. Here, the system currently includes 8 sockets with 4 installed gyrotrons, 31.75 mm circular corrugated transmission lines, four dual real time steerable launchers plus two Top Launchers. After the successful testing of the first Top Launcher, that injects the EC waves downward nearly parallel to the resonance layer with double current drive efficiency and a long absorption path, a second new Top Launcher was designed and installed with an upgraded design. Some of the operational limitations encountered for the 110 GHz X-mode second harmonic resonance scenario launched in over-dense plasmas were overcome with the alternative use of O-mode second harmonic launch. Two new gyrotrons in a new series manufactured by CPI were installed and tested, bringing up the total installed power to 3 MW. The installation of repaired and new gyrotrons will increase the number of installed gyrotrons to 7 in 2023. An addition to the DIII-D building will be completed in early 2022 to provide the space to expand the ECH/ECCD system to ten gyrotron systems, plus an extra gyrotron test socket.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A memory-driven mapping algorithm for heterogeneous systems

mpibind is a memory-driven algorithm to map parallel hybrid applications to the underlying hardware resources transparently, efficiently, and portably. There are two fundamental aspects of this algorithm. First, unlike existing mappings, its primary design point is the memory system. Compute elements are selected based on the identified memory components and not vice versa. Second, it embodies a global awareness of hybrid programming abstractions as well as heterogeneous devices.

Leon Borja, EdgarA↗
DOE OSTI2023

libEnsemble: A complete Python toolkit for dynamic ensembles of calculations

Almost all science and engineering applications eventually stop scaling: their runtime no longer decreases as available computational resources increase. Therefore, many applications will struggle to efficiently use emerging extreme-scale high-performance, parallel, and distributed systems. libEnsemble is a complete Python toolkit and workflow system for intelligently driving ensembles of experiments or simulations at massive scales. It enables and encourages multidisciplinary design, decision, and inference studies portably running on laptops, clusters, and supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2010

FY2010 Status Report on Advanced Neutronics Modeling and Validation

UNIC is the neutronics component of the massively parallel, multi-physics SHARP (Simulation for High-efficiency Advanced Reactor Prototyping) framework under development at Argonne National Laboratory. During this fiscal year, the SN2ND solver, MOCFE solver, and NODAL solver received significant development to meet the needs of the SHARP project. Additional follow-on analysis of the ZPR-6/6A from the previous year was performed in addition to new analysis of the ZPR-6/7 experiments where UNIC predictions of ZPR foil activation were made against experimentally measured values.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Scientific Data Management Center for Enabling Technologies

Managing scientific data has been identified by the scientific community as one of the most important emerging needs because of the sheer volume and increasing complexity of data being collected. Effectively generating, managing, and analyzing this information requires a comprehensive, end-to-end approach to data management that encompasses all of the stages from the initial data acquisition to the final analysis of the data. Fortunately, the data management problems encountered by most scientific domains are common enough to be addressed through shared technology solutions. Based on community input, we have identified three significant requirements. First, more efficient access to storage systems is needed. In particular, parallel file system and I/O system improvements are needed to write and read large volumes of data without slowing a simulation, analysis, or visualization engine.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2015

Benefits of Leapfrogging to Superefficiency and Low Global Warming Potential Refrigerants in Room Air Conditioning

Hydrofluorocarbons (HFCs) emitted from uses such as refrigerants and thermal insulating foam, are now the fastest growing greenhouse gases (GHGs), with global warming potentials (GWP) thousands of times higher than carbon dioxide (CO2). Because of the short lifetime of these molecules in the atmosphere,1 mitigating the amount of these short-lived climate pollutants (SLCPs) provides a faster path to climate change mitigation than control of CO2 alone. This has led to proposals from Africa, Europe, India, Island States, and North America to amend the Montreal Protocol on Substances that Deplete the Ozone Layer (Montreal Protocol) to phase-down high-GWP HFCs. Simultaneously, energy efficiency market transformation programs such as standards, labeling and incentive programs are endeavoring to improve the energy efficiency for refrigeration and air conditioning equipment to provide life cycle cost, energy, GHG, and peak load savings. In this paper we provide an estimate of the magnitude of such GHG and peak electric load savings potential, for room air conditioning, if the refrigerant transition and energy efficiency improvement policies are implemented either separately or in parallel.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

High Fidelity Simulations of Large-scale Wireless Networks (Part II - FY2017)

The ability to simulate wireless networks at large-scale for meaningful amount of time is considerably lacking in today's network simulators. For this reason, many published work in this area often limit their simulation studies to less than a 1,000 nodes and either over-simplify channel characteristics or perform studies over time scales much less than a day. In this report, we show that one can overcome these limitations and study problems of high practical consequence. This work presents two key contributions to high fidelity simulation of large-scale wireless networks: (a) wireless simulations can be sped up by more than 100X in runtime using ideas from spatial indexing algorithms and clipping of negligible signals and (b) clustering and task-oriented programming paradigm can be used to reduce inter- process communication in a parallel discrete event simulation resulting in a better scaling efficiency.

97 MATHEMATICS AND COMPUTING↗