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 1,423 records · Page 79

A Methodology for the Rapid Development of Scalable HPC Data Services

There is growing evidence in the scientifc computing community that parallel file systems are not sufficient for all HPC storage workloads. This realization has motivated extensive research in new storage system designs. The question of which alternative we should turn towards implies that there could be a single answer satisfying a wide range of very diverse applications. We argue that such a generic solution does not exist. Instead, custom data services should be designed and tailored to the needs of specific applications on specific hardware. Furthermore, custom data services should be designed in close collaboration with users. In this paper, we present a methodology for the rapid development of such data services. This methodology promotes the design of reusable building blocks that can be composed together efficiently through a runtime based on high-performance threading, tasking, and remote procedure calls. We illustrate the success of our methodology by showcasing three examples of data services designed from the same building blocks, yet targeting entirely different applications.

Dorier, Mattieu↗

SimNet: Accurate and High-Performance Computer Architecture Simulation using Deep Learning

While cycle-accurate simulators are essential tools for architecture research, design, and development, their practicality is limited by an extremely long time-to-solution for realistic applications under investigation. This work describes a concerted effort, where machine learning (ML) is used to accelerate microarchitecture simulation. First, an ML-based instruction latency prediction framework that accounts for both static instruction properties and dynamic processor states is constructed. Then, a GPU-accelerated parallel simulator is implemented based on the proposed instruction latency predictor, and its simulation accuracy and throughput are validated and evaluated against a state-of-the-art simulator. Leveraging modern GPUs, the ML-based simulator outperforms traditional CPU-based simulators significantly.

97 MATHEMATICS AND COMPUTING↗

Near Real-Time Two and Three-Dimensional Electrical Geophysical Monitoring of Natural and Engineered Processes Associated with Subsurface Remediation - 16052

Legacy soil and groundwater contamination associated with nuclear weapons production is projected to be one of the largest cleanup liabilities to the U.S. Department of Energy for the foreseeable future. Remediation of subsurface contaminants is complicated by the combined effects of geochemical heterogeneity, geological heterogeneity, and the difficulty and expense of adequate borehole access. Each of these complications contributes to uncertainty in the distribution of contaminants and the performance of in situ treatments, which increases risk to human health and the environment and cost to closure in general. Electrical resistivity tomography (ERT) is a method of remotely imaging the electrical properties of the subsurface, which are governed by both geochemical and geological structure, thereby providing useful proxies for understanding contaminant distribution and the behavior in situ remediation processes. Recent advancements in both data collection and data processing capabilities are enabling ERT monitoring to be executed and provided to site operators in near real-time. To summarize the process, ERT surveys are rapidly and continuously collected during a time-sensitive operation such as an amendment injection. At the completion of each survey, data are transferred by wireless internet to offsite supercomputing resources for parallel tomographic inversion. Inversion results are then transferred back to onsite operators and/or other locations for visualization. In this paper, we demonstrate the utility of time-lapse ERT for monitoring fluid transport (both gas and liquid phase). We then demonstrate two real-time imaging applications for monitoring both natural and engineered processes at the Hanford Site within the vadose and the saturated zones. We expect this new capability for autonomous real-time imaging to enable ERT to be used as a cost-effective, short-term rapid feedback mechanism for guiding subsurface remediation injections, in addition to longer term monitoring of post-injection performance and environmental impact, all using the same electrode array and instrumentation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

PROTEUS-SN User Manual

PROTEUS-SN is a three-dimensional, highly scalable, high-fidelity neutron transport code developed at Argonne National Laboratory. The code is applicable to all spectrum reactor transport calculations, particularly those in which a high degree of fidelity is needed either to represent spatial detail or to resolve solution gradients. PROTEUS-SN solves the second order formulation of the transport equation using the continuous Galerkin finite element method in space, the discrete ordinates approximation in angle, and the multigroup approximation in energy. PROTEUS-SN’s parallel methodology permits the efficient decomposition of the problem by both space and angle, permitting large problems to run efficiently on hundreds of thousands of cores. PROTEUS-SN can also be used in serial or on smaller compute clusters (10’s to 100’s of cores) for smaller homogenized problems, although it is generally more computationally expensive than traditional homogenized methodology codes. PROTEUS-SN has been used to model partially homogenized systems, where regions of interest are represented explicitly and other regions are homogenized to reduce the problem size and required computational resources. PROTEUS-SN solves forward and adjoint eigenvalue problems and permits both neutron upscattering and downscattering. An adiabatic kinetics option has recently been included for performing simple time-dependent calculations in addition to standard steady state calculations. PROTEUS-SN handles void and reflective boundary conditions. Multigroup cross sections can be generated externally using the MC2-3 fast reactor multigroup cross section generation code or internally using the cross section application programming interface (API) which can treat the subgroup or resonance table libraries. PROTEUS-SN is written in Fortran 90 and also includes C preprocessor definitions. The code links against the PETSc, METIS, HDF5, and MPICH libraries. It optionally links against the MOAB library and is a part of the SHARP multi-physics suite for coupled multi-physics analysis of nuclear reactors. This user manual describes how to set up a neutron transport simulation with the PROTEUS-SN code. A companion methodology manual describes the theory and algorithms within PROTEUS-SN.

97 MATHEMATICS AND COMPUTING↗

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations.

97 MATHEMATICS AND COMPUTING↗

Wrap–Around Core–Shell Heterostructures of Layered Crystals

Engineered heterostructures create new functionality by integrating dissimilar materials. Combining different 2D crystals naturally produces two distinct classes of heterostructures, vertical van der Waals (vdW) stacks or 2D sheets bonded laterally by covalent line interfaces. Here, when joining thicker layered crystals, the arising structural and topological conflicts can result in more complex geometries. Phase separation during one–pot synthesis of layered tin chalcogenides spontaneously creates core–shell structures in which large orthorhombic SnS crystals are enclosed in a wrap–around shell of trigonal SnS 2 , forcing the coexistence of parallel vdW layering along with unconventional, orthogonally layered core–shell interfaces. Measurements of the optoelectronic properties establish anisotropic carrier separation near type II core–shell interfaces and extended long–wavelength light harvesting via spatially indirect interfacial absorption, making multifunctional layered core–shell structures attractive for energy–conversion applications.

36 MATERIALS SCIENCE↗

A novel approach to increase accuracy in remotely sensed evapotranspiration through basin water balance and flux tower constraints

Remote sensing-derived evapotranspiration (RSET) products capture the spatiotemporal variations of evapotranspiration (ET) from field to basin scales with unprecedented details. However, their accuracy varies across RSET estimation methods and diverse hydroclimate regions. While ET modeling efforts to account for biophysical processes and controlling parameters have made good progress in recent years, a parallel approach of integrating in-situ ET with RSET could reduce biases in RSET products. Basin water balance ET (WBET) and flux tower ET are widely applied to evaluate RSET accuracy, yet such ET measurements are rarely used for RSET bias corrections, especially for large area applications. To address this issue, we propose a novel approach: the water balance equivalence (WABE) method, which generates spatially continuous WBET for correcting biases in RSET products. The WABE method computes synthetic WBET by integrating observed WBET and flux tower-derived FLUXCOM ET, which fills the spatial gaps of observed WBET and generates a spatially continuous WBET dataset. Synthetic WBET (2002–2015 annual average) of eight-digit hydrologic unit code (HUC8) basins across the conterminous United States (CONUS), constituting 44 % (887 out of 2035 basins) of CONUS basins, was determined within 2.0 % (RMSE = 12 %) of observed WBET at CONUS and between 1–12 % (RMSE = 3–33 %) across 18 regions in CONUS. With WABE-based bias corrections, the overall annual bias of RSET decreased from 10 % (RMSE = 34 %) to 6 % (RMSE = 26 %) across 37 flux tower sites. The WABE method offers a new approach for RSET accuracy improvement and shows great promise for large area implementations with a potential to yield substantial benefits for building accurate basin water budgets and water management decisions.

Khand, Kul↗

OsCAldOMT1 is a bifunctional O-methyltransferase involved in the biosynthesis of tricin-lignins in rice cell walls

Abstract Lignin is a phenylpropanoid polymer produced in the secondary cell walls of vascular plants. Although most eudicot and gymnosperm species generate lignins solely via polymerization of p -hydroxycinnamyl alcohols (monolignols), grasses additionally use a flavone, tricin, as a natural lignin monomer to generate tricin-incorporated lignin polymers in cell walls. We previously found that disruption of a rice 5-HYDROXYCONIFERALDEHYDE O -METHYLTRANSFERASE ( OsCAldOMT1 ) reduced extractable tricin-type metabolites in rice vegetative tissues. This same enzyme has also been implicated in the biosynthesis of sinapyl alcohol, a monolignol that constitutes syringyl lignin polymer units. Here, we further demonstrate through in-depth cell wall structural analyses that OsCAldOMT1 -deficient rice plants produce altered lignins largely depleted in both syringyl and tricin units. We also show that recombinant OsCAldOMT1 displayed comparable substrate specificities towards both 5-hydroxyconiferaldehyde and selgin intermediates in the monolignol and tricin biosynthetic pathways, respectively. These data establish OsCAldOMT1 as a bifunctional O -methyltransferase predominantly involved in the two parallel metabolic pathways both dedicated to the biosynthesis of tricin-lignins in rice cell walls. Given that cell wall digestibility was greatly enhanced in the OsCAldOMT1 -deficient rice plants, genetic manipulation of CAldOMT s conserved in grasses may serve as a potent strategy to improve biorefinery applications of grass biomass.

59 BASIC BIOLOGICAL SCIENCES↗

Experience with the alpaka performance portability library in the CMS software

ion Library for Parallel Kernel Acceleration) is a header-only C++ library that provides performance portability across different back-ends, abstracting the underlying levels of parallelism. It supports serial and parallel execution on CPUs, and extremely parallel execution on NVIDIA, AMD and Intel GPUs.This contribution will show how alpaka is used in the CMS software to develop and maintain a single code base; to use different toolchains to build the code for each supported back-end, and link them into a single application; to seamlessly select the best backend at runtime, and implement portable reconstruction algorithms that run efficiently on CPUs and GPUs from different vendors. It will describe the validation and deployment of the alpaka-based implementation in the CMS High Level Trigger, and highlight how it achieves near-native performance.

Alawieh, Jaafar [CERN]↗

A poly-epoxy surface explored by Hartree-Fock ΔSCF simulations of C1s XPS spectra

Whereas poly-epoxy polymers represent a class of materials with a wide range of applications, the structural disorder makes them difficult to model. In the present work, we use good experimental model samples in the sense that they are pure, fully polymerized, flat and smooth, defect-free, and suitable for ultrahigh vacuum x-ray photoelectron spectroscopy, XPS, experiments. In parallel, we perform Hartree-Fock, HF, calculations of the binding energies, BEs, of the C1s electrons in a model molecule composed of the two constituents of the poly-epoxy sample. These C1s BEs were determined using the HF ΔSCF method, which is known to yield accurate values, especially for the shifts of the BEs, ΔBEs. We demonstrate the benefits of combining rigorous theory with careful XPS measurements in order to obtain correct assignments of the C1s XPS spectra of the polymer sample. Both the relative binding energies—by the ΔSCF method—and relative intensities—in the sudden approximation, SA, are calculated. It results in an excellent match with the experimental spectra. We are able to identify 9 different chemical environments under the C1s peak, where an exclusively experimental work would have found only 3 contributions. In addition, we observe that some contributions are localized at discrete binding energies, whereas others allow a much wider range because of the variation of their second neighbor bound polarization. Therefore, HF-ΔSCF simulations significantly increase the spectral resolution of XPS and thus offer a new avenue for the exploration of the surface of polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An extended magnetohydrodynamics model for relativistic weakly collisional plasmas

Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to GM/c{sup 2}, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel–Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

Rotation of polarization by a moving gravitational lens

We present a simple prescription for the rotation of polarization produced by a relativistically moving gravitational lens, applicable to arbitrary deflection angle and arbitrary velocity of the lens. When geometric optics is applicable, two independent components contribute to the total rotation of polarization: (i) in the frame of the lens the polarization vector experiences minimal rotation defined by the deflection angle (as measured by a set of remote observers, or no rotation if defined in terms of a parallel-propagated tetrad); (ii) the effect of the motion of the lens on the polarization can be taken into account exactly using special-relativistic Lorentz transformation of polarization. The effects of the gravitational lensing are thus parametrized by the deflection angle of the null geodesics (not necessarily small) and the motion of the lens (not necessarily with velocities much smaller than that of light).

79 ASTRONOMY AND ASTROPHYSICS↗

SU-G-TeP1-15: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation

Purpose: To develop a Graphic Processor Unit (GPU) accelerated deterministic solution to the Linear Boltzmann Transport Equation (LBTE) for accurate dose calculations in radiotherapy (RT). A deterministic solution yields the potential for major speed improvements due to the sparse matrix-vector and vector-vector multiplications and would thus be of benefit to RT. Methods: In order to leverage the massively parallel architecture of GPUs, the first order LBTE was reformulated as a second order self-adjoint equation using the Least Squares Finite Element Method (LSFEM). This produces a symmetric positive-definite matrix which is efficiently solved using a parallelized conjugate gradient (CG) solver. The LSFEM formalism is applied in space, discrete ordinates is applied in angle, and the Multigroup method is applied in energy. The final linear system of equations produced is tightly coupled in space and angle. Our code written in CUDA-C was benchmarked on an Nvidia GeForce TITAN-X GPU against an Intel i7-6700K CPU. A spatial mesh of 30,950 tetrahedral elements was used with an S4 angular approximation. Results: To avoid repeating a full computationally intensive finite element matrix assembly at each Multigroup energy, a novel mapping algorithm was developed which minimized the operations required at each energy. Additionally, a parallelized memory mapping for the kronecker product between the sparse spatial and angular matrices, including Dirichlet boundary conditions, was created. Atomicity is preserved by graph-coloring overlapping nodes into separate kernel launches. The one-time mapping calculations for matrix assembly, kronecker product, and boundary condition application took 452±1ms on GPU. Matrix assembly for 16 energy groups took 556±3s on CPU, and 358±2ms on GPU using the mappings developed. The CG solver took 93±1s on CPU, and 468±2ms on GPU. Conclusion: Three computationally intensive subroutines in deterministically solving the LBTE have been formulated on GPU, resulting in two orders of magnitude speedup. Funding support from Natural Sciences and Engineering Research Council and Alberta Innovates Health Solutions. Dr. Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗

Stability-preserving Lossy Compression for Large-scale Partial Differential Equations

Checkpoint/Restart (C/R) strategies are vital for fault tolerance in PDE-based scientific simulations, yet traditional checkpointing incurs significant I/O overhead. Lossy compression offers a scalable solution by reducing checkpoint data size, but conventional methods often lack control over physical invariants (e.g., energy), leading to instability such as oscillations or divergence in Partial Differential Equations (PDE) systems. This paper introduces a stability-preserving compression approach tailored for PDE simulations by explicitly controlling kinetic and potential energy perturbations to ensure stable restarts. Extensive experiments conducted across diverse PDE configurations demonstrate that our method maintains numerical stability with minimal error magnification—even across multiple checkpoint-restart cycles—outperforming state-of-the-art lossy compressors. Parallel evaluations on the Frontier supercomputer show up to 8.4× improvement in checkpoint write performance and 6.3× in read performance, while maintaining relative L2 errors ∼ 2e-6 throughout continued simulation. These results provide practical guidance for balancing compression accuracy, stability, and computational efficiency in large-scale PDE applications.

Gong, Qian [ORNL] (ORCID:0000000235704142)↗

Evolution of a minimal parallel programming model

Here, we take a historical approach to our presentation of self-scheduled task parallelism, a programming model with its origins in early irregular and nondeterministic computations encountered in automated theorem proving and logic programming. We show how an extremely simple task model has evolved into a system, asynchronous dynamic load balancing (ADLB), and a scalable implementation capable of supporting sophisticated applications on today’s (and tomorrow’s) largest supercomputers; and we illustrate the use of ADLB with a Green’s function Monte Carlo application, a modern, mature nuclear physics code in production use. Our lesson is that by surrendering a certain amount of generality and thus applicability, a minimal programming model (in terms of its basic concepts and the size of its application programmer interface) can achieve extreme scalability without introducing complexity.

97 MATHEMATICS AND COMPUTING↗

Exploring small-scale chemostats to scale up microbial processes: 3-hydroxypropionic acid production in S. cerevisiae

The physiological characterization of microorganisms provides valuable information for bioprocess development. Chemostat cultivations are a powerful tool for this purpose, as they allow defined changes to one single parameter at a time, which is most commonly the growth rate. The subsequent establishment of a steady state then permits constant variables enabling the acquisition of reproducible data sets for comparing microbial performance under different conditions. We performed physiological characterizations of a 3-hydroxypropionic acid (3-HP) producing Saccharomyces cerevisiae strain in a miniaturized and parallelized chemostat cultivation system. The physiological conditions under investigation were various growth rates controlled by different nutrient limitations (C, N, P). Based on the cultivation parameters obtained subsequent fed-batch cultivations were designed.We report technical advancements of a small-scale chemostat cultivation system and its applicability for reliable strain screening under different physiological conditions, i.e. varying dilution rates and different substrate limitations (C, N, P). Exploring the performance of an engineered 3-HP producing S. cerevisiae strain under carbon-limiting conditions revealed the highest 3-HP yields per substrate and biomass of 16.6 %C-mol and 0.43 g gCDW -1 , respectively, at the lowest set dilution rate of 0.04 h -1 . 3-HP production was further optimized by applying N- and P-limiting conditions, which resulted in a further increase in 3-HP yields revealing values of 21.1 %C-mol and 0.50 g gCDW -1 under phosphate-limiting conditions. The corresponding parameters favoring an increased 3-HP production, i.e. dilution rate as well as C- and P-limiting conditions, were transferred from the small-scale chemostat cultivation system to 1-L bench-top fermenters operating in fed-batch conditions, revealing 3-HP yields of 15.9 %C-mol and 0.45 g gCDW -1 under C-limiting, as well as 25.6 %C-mol and 0.50 g gCDW -1 under phosphate-limiting conditions. Small-scale chemostat cultures are well suited for the physiological characterization of microorganisms, particularly for investigating the effect of changing cultivation parameters on microbial performance. In our study, optimal conditions for 3-HP production comprised (i) a low dilution rate of 0.04 h -1 under carbon-limiting conditions and (ii) the use of phosphate-limiting conditions. Similar 3-HP yields were achieved in chemostat and fed-batch cultures under both C- and P-limiting conditions proving the growth rate as robust parameter for process transfer and thus the small-scale chemostat system as powerful tool for process optimization.

59 BASIC BIOLOGICAL SCIENCES↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗