Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallel Applications”

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 685 records · Page 38

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗

Laser spectroscopy for the study of exotic nuclei

Investigation into the properties and structure of unstable nuclei far from stability is a key avenue of research in modern nuclear physics. These efforts are motivated by the continual observation of unexpected structure phenomena in nuclei with unusual proton-to-neutron ratios. In recent decades, laser spectroscopy techniques have made significant contributions in our understanding of exotic nuclei in different mass regions encompassing almost the entire nuclear chart. This is achieved through determining multiple fundamental properties of nuclear ground and isomeric states, such as nuclear spins, magnetic dipole and electric quadrupole moments and charge radii, via the measurement of hyperfine structures and isotope shifts in the atomic or ionic spectra of the nuclei of interest. These properties offer prominent tests of recently developed state-of-the-art nuclear theory and help to stimulate new developments in improving the many-body methods and nucleon–nucleon interactions at the core of these models. With the aim of exploring more exotic short-lived nuclei located ever closer to the proton and neutron driplines, laser spectroscopy techniques, with their continuous technological developments towards higher resolution and higher sensitivity, are extensively employed at current- and next-generation radioactive ion beam facilities worldwide. Ongoing efforts in parallel promise to improve the availability of these even more exotic species at next-generation facilities. Very recently, an innovative application of laser spectroscopy on molecules containing short-lived nuclei has been demonstrated offering additional opportunities for several fields of research, e.g. fundamental symmetry studies and astrophysics. Here, in this review, the basic nuclear properties measurable with laser spectroscopy will be introduced. How these observables are associated with nuclear structure and nucleon–nucleon interactions will be discussed. Following this, a general overview of different laser spectroscopy methods will be given with particular emphasis on technical advancements reported in recent years. The main focus of this article is to review the numerous highlights that have resulted from studying exotic nuclei in different mass regions with laser spectroscopy techniques since the last edition in this series. Finally, the challenges facing the field in addition to future opportunities will be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A two-phase mixing layer between parallel gas and liquid streams: multiphase turbulence statistics and influence of interfacial instability

The two-phase mixing layer formed between parallel gas and liquid streams is an important fundamental problem in turbulent multiphase flows. The problem is relevant to many industrial applications and natural phenomena, such as air-blast atomizers in fuel injection systems and breaking waves in the ocean. The velocity difference between the gas and liquid streams triggers an interfacial instability which can be convective or absolute depending on the stream properties and injection parameters. In the present study, a direct numerical simulation of a two-phase gas–liquid mixing layer that lie in the absolute instability regime is conducted. A dominant frequency is observed in the simulation and the numerical result agrees well with the prediction from viscous stability theory. As the interfacial wave plays a critical role in turbulence transition and development, the temporal evolution of turbulent fluctuations (such as the enstrophy) also exhibits a similar frequency. To investigate the statistical response of the multiphase turbulence flow, the simulation has been run for a long physical time so that time-averaging can be performed to yield the statistically converged results for Reynolds stresses and the turbulent kinetic energy (TKE) budget. An extensive mesh refinement study using from 8 million to about 4 billions cells has been performed. The turbulent dissipation is shown to be highly demanding on mesh resolution compared with other terms in TKE budget. The results obtained with the finest mesh are shown to be close to converged results of turbulent dissipation which allow us to obtain estimations of the Kolmogorov and Hinze scales. The estimated Kolmogorov scale is found to be similar to the cell size of the finest mesh used here. The computed Hinze scale is significantly larger than the size of droplets observed and does not seem to be a relevant length scale to describe the smallest size of droplets formed in atomization.

Mechanics↗

Width-bending characteristic of REBCO HTS tape and flat-tape Rutherford-type cabling

Abstract The width-bending behaviors for a rare earth barium copper oxide (REBCO) high-temperature superconducting (HTS) tape have been investigated. It has been found experimentally that the width-bending strain of a Yttrium Barium Copper Oxide (YBCO) tape does not degrade the critical current as much as is expected from the axial strain. The critical current is not directly affected by the width-bending strain but by the axial strain, which the width-bending strain generates by the Poisson effect. Since Poisson’s ratio is about 0.3, the axial strain effective on the critical current is about 30% of the width-bending strain. Therefore, the width-bending strain of even 1.5% degrades the critical current by only 15% if the REBCO layer side is bent inward. However, the critical current degraded by more than 70% when the REBCO layer is outward. The width-bending effects on the critical current have been further examined, considering Poisson’s ratio changes and the neutral plane shift of the REBCO tape substrate. Those changes would occur when the material yields due to severe width-bending. Based on the width-bending of REBCO tapes, flat-tape Rutherford-type cabling is discussed. A REBCO Rutherford-type cable can provide various advantages, especially for narrow-tape cabling. The flat-tape Rutherford-type cable has better characteristics against an electromagnetic transverse Lorentz force. Furthermore, the tape length of the cable can be approximately the same as the cable length allowing for excellent tape usage. REBCO Rutherford-type cabling will be a promising high-current, high-field cabling method using thin substrate REBCO tapes. It will be useful for AC ramp-field and pulse-field applications with low AC losses and low shield currents. It notes that the filaments in the cable are symmetrically distributed in parallel. Therefore, each filament’s inductance is uniform, and a uniform current distribution among the filaments can be obtained.

Physics↗

Enhancing quantum annealing accuracy through replication-based error mitigation *

Abstract Quantum annealers like those manufactured by D-Wave Systems are designed to find high quality solutions to optimization problems that are typically hard for classical computers. They utilize quantum effects like tunneling to evolve toward low-energy states representing solutions to optimization problems. However, their analog nature and limited control functionalities present challenges to correcting or mitigating hardware errors. As quantum computing advances towards applications, effective error suppression is an important research goal. We propose a new approach called replication based mitigation (RBM) based on parallel quantum annealing (QA). In RBM, physical qubits representing the same logical qubit are dispersed across different copies of the problem embedded in the hardware. This mitigates hardware biases, is compatible with limited qubit connectivity in current annealers, and is well-suited for currently available noisy intermediate-scale quantum annealers. Our experimental analysis shows that RBM provides solution quality on par with previous methods while being more flexible and compatible with a wider range of hardware connectivity patterns. In comparisons against standard QA without error mitigation on larger problem instances that could not be handled by previous methods, RBM consistently gets better energies and ground state probabilities across parameterized problem sets.

Djidjev, Hristo N. (ORCID:0000000192868824)↗

The Ongoing Evolution of OpenMP

Here, this paper presents an overview of the past, present and future of the OpenMP application programming interface (API). While the API originally specified a small set of directives that guided shared memory fork-join parallelization of loops and program sections, OpenMP now provides a richer set of directives that capture a wide range of parallelization strategies that are not strictly limited to shared memory. As we look toward the future of OpenMP, we immediately see further evolution of the support for that range of parallelization strategies and the addition of direct support for debugging and performance analysis tools. Looking beyond the next major release of the specification of the OpenMP API, we expect the specification eventually to include support for more parallelization strategies and to embrace closer integration into its Fortran, C and, in particular, C++ base languages, which will likely require the API to adopt additional programming abstractions.

97 MATHEMATICS AND COMPUTING↗

PETSc DMNetwork: A Library for Scalable Network PDE-Based Multiphysics Simulations

We present DMNetwork, a high-level package included in the PETSc library for the simulation of multiphysics phenomena over large-scale networked systems. The library aims at applications that have networked structures such as those in electrical, gas, and water distribution systems. DMNetwork provides data and topology management, parallelization for multiphysics systems over a network, and hierarchical and composable solvers to exploit the problem structure. DMNetwork eases the simulation development cycle by providing the necessary infrastructure through simple abstractions to define and query the network components. This article presents the design of DMNetwork, illustrates its user interface, and demonstrates its ability to solve multiphysics systems, such as an electric circuit, a network of power grid and water subnetworks, and transient hydraulic systems over large networks with more than 2 billion variables on extreme-scale computers using up to 30,000 processors.

extreme-scale computing↗

Radiographic Imaging and Tomography

Radiographic imaging and tomography (RadIT) come in many flavors such as X-ray imaging and tomography (IT), proton IT, neutron IT, muon IT, neutrino IT, etc. We identify five RadIT themes: Physics, Sources, Detectors, Methods, and Data Science, which are integral parts of image interpretation and three-dimensional (3D) tomographic reconstruction. Traditionally, RadIT have been driven by medicine, non-destructive testing, material sciences, and security applications. The latest thrusts of growth come from automation, machine vision, additive manufacturing and virtual reality (‘metaverse’). The five RadIT themes parallel their counterparts in optical IT. Synergies among different forms of RadIT and with optical IT motivate further advances towards multi-modal IT and quantum IT.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Correction for Metastability in the Quantification of PID in Thin-film Module Testing: Preprint

A fundamental change in the analysis for the accelerated stress testing of thin-film modules is proposed, whereby power changes due to metastability and other effects that may occur due to the thermal history are removed from the power measurement that we obtain as a function of the applied stress factor. The power of reference modules normalized to an initial state - undergoing the same thermal and light- exposure history but without the applied stress factor such as humidity or voltage bias - is subtracted from that of the stressed modules. For better understanding and appropriate application in standardized tests, the method is demonstrated and discussed for potential-induced degradation testing in view of the parallel-occurring but unrelated physical mechanisms that can lead to confounding power changes in the module.

14 SOLAR ENERGY↗

Parallel multiple-chain DRAM MCMC for large-scale geosteering inversion and uncertainty quantification

Geosteering is the proactive control of a wellbore placement based on the downhole measurements, aiming at maximizing economic production from the well. We report the recent development of azimuthal resistivity logging-whiledrilling (LWD) tools have a much larger depth of detection, thus sets higher demands for geosteering inversion as more unknown parameters need to be taken into account in the earth model to be inverted. For complicated nonlinear problems, traditional deterministic inversion methods are more likely to be trapped near local optima. In general, Markov chain Monte Carlo (MCMC) inversion methods are more capable of finding the global optimal solution and providing additional uncertainty analyses by sampling from the target distribution. However, MCMC methods usually have the problem of slow convergence. Even though optimization methods like delayed rejection (DR) and adaptive Metropolis (AM) were proposed to speed up the convergence, using MCMC methods to solve geosteering inversion problems may still incur unacceptable time cost. To reduce the sampling cost of MCMC methods, we use parallel multiple-chain DRAM MCMC methods to solve geosteering inverse problems and estimate the corresponding uncertainty. A clustering method, density-based spatial clustering of applications with noise (DBSCAN), is applied to select the best solution among many results. Lastly, the simulation results show that running many relatively short MCMC chains for one problem can obtain a similar result as running a long single chain. Besides, avoiding communications between multiple Markov chains during the sampling process yields almost linear scalability.

97 MATHEMATICS AND COMPUTING↗

Cabana: A Performance Portable Library for Particle-Based Simulations

Particle-based simulations are ubiquitous throughout many fields of computational science and engineering, spanning the atomistic level with molecular dynamics (MD), to mesoscale particle-in-cell (PIC) simulations for solid mechanics, device-scale modeling with PIC methods for plasma physics, and massive N-body cosmology simulations of galaxy structures, with many other methods in between (Hockney & Eastwood, 1989). While these methods use particles to represent significantly different entities with completely different physical models, many low-level details are shared including performant algorithms for short- and/or long-range particle interactions, multi-node particle communication patterns, and other data management tasks such as particle sorting and neighbor list construction. Cabana is a performance portable library for particle-based simulations, developed as part of the Co-Design Center for Particle Applications (CoPA) within the Exascale Computing Project (ECP) (Alexander et al., 2020). The CoPA project and its full development scope, including ECP partner applications, algorithm development, and similar software libraries for quantum MD, is described in (Mniszewski et al., 2021). Cabana uses the Kokkos library for on-node parallelism (Edwards et al., 2014; Trott et al., 2022), enabling simulation on multi-core CPU and GPU architectures, and MPI for GPU-aware, multi-node communication. Cabana provides particle simulation capabilities on almost all current Kokkos backends, including serial execution, OpenMP (including OpenMP-Target for GPUs), CUDA (NVIDIA GPUs), HIP (AMD GPUs), and SYCL (Intel GPUs), providing a clear path for the coming generation of accelerator-based exascale hardware. Cabana builds on Kokkos by providing new particle data structures and particle algorithms resulting in a similar execution policy-based, node-level programming model that is intended to be used in addition to the core Kokkos library within an application. Cabana is designed as an application and physics agnostic, but particle-specific toolkit which can either be used to generate a new application, or to be used as needed in existing applications at various levels of invasiveness including through interfaces that wrap user memory in existing data structures.

97 MATHEMATICS AND COMPUTING↗

ICED: An Integrated CGRA Framework Enabling DFVS-Aware Acceleration

oarse-grained reconfigurable arrays (CGRAs) are a promising solution to enable energy-efficient acceleration of applications from different domains. By leveraging reconfiguration at the functional level, they can adapt to significantly different computational patterns. Existing CGRA mapping approaches extract instruction-level parallelism, exploit loop-pipelining opportunities, guarantee the data dependency, and target high throughput of a given loop. However, the recurrence data-dependency in the DFG and the mismatch between required and available computing/communication resources complicate the mapping, and might lead to significant unbalances in the utilization of the CGRA's tiles. This results in wasted power for tiles with low utilization. Applying dynamic voltage and frequency scaling (DVFS) can potentially solve this challenge and improve energy efficiency by adjusting voltage and frequency of different tiles independently. CGRAs have also been successful in accelerating data-dependent streaming applications. However, in these applications, the execution time of each kernel in the pipeline might dynamically vary depending on the characteristics of the input. This also leads to under-utilization of resources for the dynamically changing kernels that do not limit the application throughput. DVFS can also improve energy efficiency for these applications by dynamically changing the voltage and frequency levels of tiles that host non performance-constraining kernels. This paper proposes ICEDTEA -- an integrated DVFS-aware framework to map applications on CGRAs that support power islands. ICEDTEA proposes a CGRA architecture supporting DVFS islands at varying granularity (from a single tile to a group of tiles) and the related DVFS-aware compilation and mapping toolchain. ICEDTEA is the first work that introduces DVFS support for spatio-temporal CGRAs at power-island levels. The experimental evaluation shows that ICEDTEA improves average utilization by 2.3$\times$ and energy-efficiency by 1.32$\times$ over a conventional CGRA. With streaming applications, ICEDTEA improves energy efficiency by 1.12$\times$ over a state-of-the-art CGRA that introduces partial dynamic reconfiguration to adapt to variations in kernels' throughput.

Tan, Cheng↗

Design for Recycling Principles Applicable to Selected Clean Energy Technologies: Crystalline-Silicon Photovoltaic Modules, Electric Vehicle Batteries, and Wind Turbine Blades

Abstract The global growth of clean energy technology deployment will be followed by parallel growth in end-of-life (EOL) products, bringing both challenges and opportunities. Cumulatively, by 2050, estimates project 78 million tonnes of raw materials embodied in the mass of EOL photovoltaic (PV) modules, 12 billion tonnes of wind turbine blades, and by 2030, 11 million tonnes of lithium-ion batteries. Owing partly to concern that the projected growth of these technologies could become constrained by raw material availability, processes for recycling them at EOL continue to be developed. However, none of these technologies are typically designed with recycling in mind, and all of them present challenges to efficient recycling. This article synthesizes and extends design for recycling (DfR) principles based on a review of published industrial and academic best practices as well as consultation with experts in the field. Specific principles developed herein apply to crystalline-silicon PV modules, batteries like those used in electric vehicles, and wind turbine blades, while a set of broader principles applies to all three of these technologies and potentially others. These principles are meant to be useful for stakeholders—such as research and development managers, analysts, and policymakers—in informing and promoting decisions that facilitate DfR and, ultimately, increase recycling rates as a way to enhance the circularity of the clean energy economy. The article also discusses some commercial implications of DfR. Graphical Abstract

36 MATERIALS SCIENCE↗

An analysis of various solution strategies and perturbations on inputs of the reactor shielding problem - 25584

Previous work focused on the use of automated variance reduction parameters for use in a fixed source reactor shielding problem to solve for ex-core particle transport results. The current work assesses that strategy in light of a fully-coupled Monte Carlo problem run as a critical system with parallel processing and the use of moderate variance reduction parameters set manually. The results from the two strategies compare well. Also, results are compared between the MCNP, Virtual Environment for Reactor Applications (VERA) and Shift computer codes. Shift and VERA deliver similar results and compare well to MCNP. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A New Integrated Analysis Suite for Fast-Ion Study in KSTAR

Here, an integrated workflow for fast-ion analysis was developed by adapting the One Modeling Framework for Integrated Task (OMFIT) workflow manager to support a standard and unified analysis platform for KSTAR users. The newly established analysis suite offers a graphical user interface–based workflow to enable users to readily access and handle experimental data archived in various data formats and servers. Further, users can analyze the data by importing modules designed for conducting certain tasks, such as profile fitting, equilibrium reconstruction, and postprocessing of tokamak data. The procedures for preparing the inputs for fast-ion simulations are streamlined by a common workflow manager, which enables the parallel processing of various tasks to efficiently analyze large fast-ion datasets. The OMFIT platform comprises a flexible Python-based application that enables users to freely manipulate the Python scripts for applications that are unavailable in the standard workflow. The framework also offers mapping tools to translate the output data into the Integrated Modeling and Analysis Suite format to maintain application compatibility for future ITER burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic simulation of a low-pressure helium discharge with comparison to experimental measurements

Modern computational validation efforts rely on comparison of known experimental quantities such as current, voltage, particle densities, and other plasma properties with the same values determined through simulation. A discrete photon approach for radiation transport was recently incorporated into a particle-in-cell/direct simulation Monte Carlo code. As a result, spatially and temporally resolved synthetic spectra may be generated even for non-equilibrium plasmas. The generation of this synthetic spectra lends itself to potentially new validation opportunities. In this work, initial comparisons of synthetic spectra are made with experimentally gathered optical emission spectroscopy. A custom test apparatus was constructed that contains a 0.5 cm gap distance parallel plane discharge in ultra high purity helium gas (99.9999%) at a pressure of 75 Torr. Plasma generation is initiated with the application of a fast rise-time, 100 ns full-width half maximum, 2.0 kV voltage pulse. Transient electrical diagnostics are captured along with time-resolved emission spectra. A one-dimensional simulation is run under the same conditions and compared against the experiment to determine if sufficient physics are included to model the discharge. To sync the current measurements from experiment and simulation, significant effort was undertaken to understand the kinetic scheme required to reproduce the observed features. Additionally, the role of the helium molecule excimer emission and atomic helium resonance emission on photocurrent from the cathode are studied to understand which effect dominates photo-feedback processes. Results indicate that during discharge development, atomic helium resonance emission dominates the photo-flux at the cathode even though it is strongly self-absorbed. A comparison between the experiment and simulation demonstrates that the simulation reproduces observed features in the experimental discharge current waveform. Here, the synthesized spectra from the kinetic method produces more favorable agreement with the experimental data than a simple local thermodynamic equilibrium calculation and is a first step towards using spectra generated from a kinetic method in validation procedures. The results of this study produced a detailed compilation of important helium plasma chemistry reactions for simulating transient helium plasma discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Engineering glycoside hydrolase stability by the introduction of zinc binding

The development of robust enzymes, in particular cellulases, is a key step in the success of biological routes to `second-generation' biofuels. The typical sources of the enzymes used to degrade biomass include mesophilic and thermophilic organisms. The endoglucanase J30 from glycoside hydrolase family 9 was originally identified through metagenomic analyses of compost-derived bacterial consortia. These studies, which were tailored to favor growth on targeted feedstocks, have already been shown to identify cellulases with considerable thermal tolerance. The amino-acid sequence of J30 shows comparably low identity to those of previously analyzed enzymes. As an enzyme that combines a well measurable activity with a relatively low optimal temperature (50°C) and a modest thermal tolerance, it offers the potential for structural optimization aimed at increased stability. Here, the crystal structure of wild-type J30 is presented along with that of a designed triple-mutant variant with improved characteristics for industrial applications. Through the introduction of a structural Zn 2+ site, the thermal tolerance was increased by more than 10°C and was paralleled by an increase in the catalytic optimum temperature by more than 5°C.

59 BASIC BIOLOGICAL SCIENCES↗

A High Performance Block Eigensolver for Nuclear Configuration Interaction Calculations

As on-node parallelism increases and the performance gap between the processor and the memory system widens, achieving high performance in large-scale scientific applications requires an architecture-aware design of algorithms and solvers. We focus on the eigenvalue problem arising in nuclear Configuration Interaction (CI) calculations, where a few extreme eigenpairs of a sparse symmetric matrix are needed. Here, we consider a block iterative eigensolver whose main computational kernels are the multiplication of a sparse matrix with multiple vectors (SpMM), and tall-skinny matrix operations. We then present techniques to significantly improve the SpMM and the transpose operation SpMM T by using the compressed sparse blocks (CSB) format. We achieve 3-4× speedup on the requisite operations over good implementations with the commonly used compressed sparse row (CSR) format. We develop a performance model that allows us to correctly estimate the performance of our SpMM kernel implementations, and we identify cache bandwidth as a potential performance bottleneck beyond DRAM. We also analyze and optimize the performance of LOBPCG kernels (inner product and linear combinations on multiple vectors) and show up to 15× speedup over using high performance BLAS libraries for these operations. The resulting high performance LOBPCG solver achieves 1.4× to 1.8× speedup over the existing Lanczos solver on a series of CI computations on high-end multicore architectures (Intel Xeons). We also analyze the performance of our techniques on an Intel Xeon Phi Knights Corner (KNC) processor.

97 MATHEMATICS AND COMPUTING↗