Search NASA⌕ Search

SEARCH · Search NASA

Results for “code portability”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of Heat Pipe Failures in Microreactors

Microreactors provide new opportunities for nuclear reactor applications given their long-lived fuel source, compact size, portability, reliability as it relates to safety and their self-regulating nature. To take advantage of these opportunities there are technological challenges that need to be overcome. One advancing technology often used in conjunction with microreactor designs are heat pipes, which are used for heat removal from the reactor core to the power conversion system. Heat pipes are a substitute for the traditional system in which the same fluid used as coolant in the reactor is then used as the working fluid in the power conversion system. This document describes how heat pipes operate, what happens when a heat pipe fails, and what codes can be used to simulate the results. Our analyses show that a single/double heat pipe failure results in a temperature increase of 15-50°C respectively in surrounding heat pipes. Heat pipe microreactors could be easily designed such that this increase is within the maximum parameters allowable for safety and does not result in cascading heat pipe failure. Thus, such failures do not lead to severe degradation of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Porting HEP Parameterized Calorimeter Simulation Code to GPUs

The High Energy Physics (HEP) experiments, such as those at the Large Hadron Collider (LHC), traditionally consume large amounts of CPU cycles for detector simulations and data analysis, but rarely use compute accelerators such as GPUs. As the LHC is upgraded to allow for higher luminosity, resulting in much higher data rates, purely relying on CPUs may not provide enough computing power to support the simulation and data analysis needs. As a proof of concept, we investigate the feasibility of porting a HEP parameterized calorimeter simulation code to GPUs. We have chosen to use FastCaloSim, the ATLAS fast parametrized calorimeter simulation. While FastCaloSim is sufficiently fast such that it does not impose a bottleneck in detector simulations overall, significant speed-ups in the processing of large samples can be achieved from GPU parallelization at both the particle (intra-event) and event levels; this is especially beneficial in conditions expected at the high-luminosity LHC, where extremely high per-event particle multiplicities will result from the many simultaneous proton-proton collisions. We report our experience with porting FastCaloSim to NVIDIA GPUs using CUDA. A preliminary Kokkos implementation of FastCaloSim for portability to other parallel architectures is also described.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

To Exascale and Beyond—The Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM), a Performance Portable Global Atmosphere Model for Cloud-Resolving Scales

The new generation of heterogeneous CPU/GPU computer systems offer much greater computational performance but are not yet widely used for climate modeling. One reason for this is that traditional climate models were written before GPUs were available and would require an extensive overhaul to run on these new machines. In addition, even conventional “high–resolution” simulations don't currently provide enough parallel work to keep GPUs busy, so the benefits of such overhaul would be limited for the types of simulations climate scientists are accustomed to. The vision of the Simple Cloud-Resolving Energy Exascale Earth System (E3SM) Atmosphere Model (SCREAM) project is to create a global atmospheric model with the architecture to efficiently use GPUs and horizontal resolution sufficient to fully take advantage of GPU parallelism. After 5 years of model development, SCREAM is finally ready for use. In this paper, we describe the design of this new code, its performance on both CPU and heterogeneous machines, and its ability to simulate real-world climate via a set of four 40 day simulations covering all 4 seasons of the year.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗
DOE OSTI2024

Portable self-contained pressure testing manifold and associated methods

A portable, case-enclosed pressure test manifold, configured to support both pneumatic and hydrostatic pressure testing, and comprising a downstream fluid path (e.g., a system test gauge, a ball valve, a system test port, a system relief valve, and other components/piping) and an upstream fluid path (e.g., a low-pressure gauge, a high-pressure gauge, and other components/piping). A regulator and/or the ball valve selectively operate the downstream fluid path, defining a hydrostatic pressure testing path, to receive a liquid; and the upstream fluid path in combination with the downstream fluid path, defining a pneumatic pressure testing path, to receive a gas. Manifold components are pressure rated based on use and are ASME Code B31.3 compliant. A flexible hose connects with either of the pneumatic pressure testing path and the hydrostatic pressure testing path of the manifold. The case includes parking receptables and a carrying handle configured to stow the hose during transport.

Milholland, Roger A.↗
DOE OSTI2020

Three dimensions, two microscopes, one code: Automatic differentiation for x-ray nanotomography beyond the depth of focus limit

Conventional tomographic reconstruction algorithms assume that one has obtained pure projection images, involving no within-specimen diffraction effects nor multiple scattering. Advances in x-ray nanotomography are leading toward the violation of these assumptions, by combining the high penetration power of x-rays, which enables thick specimens to be imaged, with improved spatial resolution that decreases the depth of focus of the imaging system. We describe a reconstruction method where multiple scattering and diffraction effects in thick samples are modeled by multislice propagation and the 3D object function is retrieved through iterative optimization. We show that the same proposed method works for both full-field microscopy and for coherent scanning techniques like ptychography. Our implementation uses the optimization toolbox and the automatic differentiation capability of the open-source deep learning package TensorFlow, demonstrating a straightforward way to solve optimization problems in computational imaging with flexibility and portability.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Western States Building Energy & Controls Apprenticeship (BECA) Program

How does this project help to understand the challenges of workforce development in the commercial building energy management industry? The purpose of this project was to create a replicable, scalable, and portable apprenticeship program for building energy management and controls. This paper will demonstrate how local processes can be more adaptable and inclusive than federal processes in achieving workforce development goals. Our work provides insights into the technical effectiveness of the program, enabling others to achieve greater success in their workforce development initiatives. The project required more time and financial resources than initially anticipated and spent a year in a no-cost extension working to accomplish the Statement of Project Objectives. The original goal was to create an Industry Related Apprenticeship Program (IRAP). The development of this apprenticeship presented challenges due to the lack of existing programs for reference and the absence of relevant industry classification codes by the Department of Labor (DOL). Notably, the role of a commercial building energy analyst is not recognized by the DOL. The SIC (industry codes) do not have a good description of this job. There are many that may fall into related categories, but none are the actual duties of an energy analyst. Discussions with the DOL indicated that substantial groundwork was necessary before a national apprenticeship program could be implemented, causing delays in the program’s commencement. As a result, students experienced longer wait times before starting their apprenticeship component. The apprenticeship program officially launched on September 14, 2021. The success of the State of Oregon’s apprenticeship program, the first of its kind in the state, underscored the flexibility and effectiveness of local initiatives compared to federal efforts. The COVID-19 pandemic also significantly impacted the project’s success. Beginning in March 2020, the pandemic led to widespread closures of schools and colleges by fall 2020. By September 2021, when the apprenticeship option became available, enrollment in colleges and universities nationwide had decreased, affecting student participation in the program. Moreover, as employers transitioned their employees to remote work, there was limited interaction with external personnel, influencing the willingness of training agents (employers) to integrate additional workers into their teams. This paper addresses several challenges encountered during the project, with the hope that future workforce development efforts will benefit from these experiences. We Final Technical Report 4 | Page encourage others to engage with state and federal agencies to enhance and update pathways for workforce development and apprenticeship programs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2018

PyFly: A fast, portable aerodynamics simulator

Here, we present a fast, user-friendly implementation of a potential flow solver based on the unsteady vortex lattice method (UVLM), namely PyFly. UVLM computes the aerodynamic loads applied on lifting surfaces while capturing the unsteady effects such as the added mass forces, the growth of bound circulation, and the wake while assuming that the flow separation location is known a priori. This method is based on discretizing the body surface into a lattice of vortex rings and relies on the Biot–Savart law to construct the velocity field at every point in the simulated domain. We introduce the pointwise approximation approach to simulate the interactions of the far-field vortices to overcome the computational burden associated with the classical implementation of UVLM. The computational framework uses the Python programming language to provide an easy to handle user interface while the computational kernels are written in Fortran. The mixed language approach enables high performance regarding solution time and great flexibility concerning easiness of code adaptation to different system configurations and applications. The computational tool predicts the unsteady aerodynamic behavior of multiple moving bodies (e.g., flapping wings, rotating blades, suspension bridges) subject to incoming air. The aerodynamic simulator can also deal with enclosure effects, multi-body interactions, and B-spline representation of body shapes. Finally, we simulate different aerodynamic problems to illustrate the usefulness and effectiveness of PyFly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Improving the Field Portability of Neutron and Gamma Counting Techniques

Savannah River National Laboratory (SRNL) puts science to work by developing techniques that better meet customers needs for in-situ measurements. There is high demand for laboratory quality measurements in field situations, improving the cost and speed for first responders and environmental monitoring. This work focused on two improvements (1) reading a signal from a neutron detector using a Raspberry Pi (2) testing the Peak-to-Valley method for determining the depth of a source in soil contamination measurements. A Raspberry Pi is a cheap, compact computer. This affordable, compact computer is able to do simple tasks that a desktop computer is able to do, like spreadsheets, browse the internet, play games and even program a simple robot. Additionally, a Raspberry Pi can run off of a battery, so it could be deployed for long term measurements in adverse conditions like neutron monitors deployed in a shipping port. Using the Raspberry Pi and libraries easily accessible on the internet, a code was written in Python that allows the Pi to keep track of how many times a button is pushed. This work will be used with the future goal of creating low power-budget neutron detectors. 1) Make sure Python was correctly installed on the Raspberry Pi. 2) Write a code that will keep track of how many times a button connected to the Pi is pressed. A code was written to use the Raspberry Pi to count the number of times the button was pushed. Moving forward the Raspberry Pi will be tested with a neutron counter. Ideally, each signal from the neutron pulse will be equivalent of 'push the button' and be counted. In the characterization of a soil contamination area both the location and the depth of the source is important. Attenuation of the source by the soil can make determining its depth challenging, particularly for monoenergetic radionuclides like Cs-137. To address this issue, a method of calculating the 'peak-to-valley' ratio of the gamma count spectra is implemented to help calculate the known depth of the source. 1) Determine the 'peak-to-valley' ratio of a Cs-137 source at different depths. 2) Find an appropriate function to serve as a model for the peak-to-valley ratio as a function of depth. Several Petri dishes, filled with sand, were stacked with one containing a Cs-137 source. The individual Petri dishes allowed for the source to be placed at different depths. A bucket was placed over the sand and source, to keep the detector at a consistent height. A Petri dish containing a Eu-155 source, which acted as an internal calibration standard, was placed on top of the bucket. A NaI detector was place on top of the Eu-155 source and the assembly was surrounded by a tungsten shield. Count times were 20 minutes and the data was analyzed using PeakEasy. Data was recorded for Cs-137 at a series of depths in sand. Both exponential and logarithmic functions provided a reasonable fit to the data. The experimental data showed a decreasing 'peak-to-valley' ratio as the depth was increased. The next steps for this project involves determining the optimum fit function and combining the results from sand to different soil types to get a better representation of the soil found at SRS. Eventually, the peak-to-valley ratio will be used to estimate the depth of contamination in field measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

autoGEMM: Pushing the Limits of Irregular Matrix Multiplication on Arm Architectures

This paper presents an open-source library that pushes the limits of performance portability for irregular General Matrix Multiplication (GEMM) on the widely-used Arm architectures. Our library, autoGEMM, is designed to support a wide range of Arm processors: from edge devices to HPC-grade CPUs. autoGEMM generates optimized kernels for various hardware configurations by auto-combining fragments of autogenerated micro-kernels that employ hand-written optimizations to maximize computational efficiency. We optimize the kernel pipeline by tuning the register reuse and the data load/store overlapping. In addition, we use a dynamic tiling scheme to generate balanced tile shapes. Finally, we position autoGEMM on top of the TVM framework where our dynamic tiling scheme prunes the search space for TVM to identify the optimal combination of parameters for code optimization. Evaluations on five different classes of Arm chips demonstrate the advantages of autoGEMM. For small matrices, autoGEMM achieves 98% of peak and up to 2.0x speedup over state-of-the-art libraries such as LIBXSMM and LibShalom. For irregular matrices (i.e. tall skinny and long rectangles), autoGEMM is 1.3-2.0x faster than widely-used libraries such as OpenBLAS and Eigen. autoGEMM is publicly available at: https://github.com/wudu98/autoGEMM.

Wu, Du↗
DOE OSTI2019

Parallel Computing for Module-Based Computational Experiment

Large-scale scientific code plays an important role in scientific researches. In order to facilitate module and element evaluation in scientific applications, we introduce a unit testing framework and describe the demand for module-based experiment customization. We then develop a parallel version of the unit testing framework to handle long-term simulations with a large amount of data. Specifically, we apply message passing based parallelization and I/O behavior optimization to improve the performance of the unit testing framework and use profiling result to guide the parallel process implementation. Finally, we present a case study on litter decomposition experiment using a standalone module from a large-scale Earth System Model. This case study is also a good demonstration on the scalability, portability, and high-efficiency of the framework.

Yao, Zhuo↗
DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Water Table Model (WTM) (v2.0.1): coupled groundwater and dynamic lake modelling

Abstract. Ice-free land comprises 26 % of the Earth's surface and holds liquid water that delineates ecosystems, affects global geochemical cycling, and modulates sea levels. However, we currently lack the capacity to simulate and predict these terrestrial water changes across the full range of relevant spatial (watershed to global) and temporal (monthly to millennial) scales. To address this knowledge gap, we present the Water Table Model (WTM), which integrates coupled components to compute dynamic lake and groundwater levels. The groundwater component solves the 2D horizontal groundwater flow equation using non-linear equation solvers from the C++ PETSc (Portable, Extensible Toolkit for Scientific Computation) library. The dynamic lake component makes use of the Fill–Spill–Merge (FSM) algorithm to move surface water into lakes, where it may evaporate or affect groundwater flow. In a proof-of-concept application, we demonstrate the continental-scale capabilities of the WTM by simulating the steady-state climate-driven water table for the present day and the Last Glacial Maximum (LGM; 21 000 calendar years before present) across the North American continent. During the LGM, North America stored an additional 14.98 cm of sea-level equivalent (SLE) in lakes and groundwater compared to the climate-driven present-day scenario. We compare the present-day result to other simulations and real-world data. Open-source code for the WTM is available on GitHub and Zenodo.

Callaghan, Kerry L. (ORCID:0000000226740838)↗
DOE OSTI2021

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

UEMtomaton: A Source-Available Platform to Aid in Start-up of Ultrafast Electron Microscopy Labs

The steady rise in the number of ultrafast electron microscopy (UEM) labs, in addition to the opacity and lack of detailed descriptions of current approaches that would enable point-by-point construction, has created an opportunity for sharing common methods and instrumentation for (for example) automating data acquisition to assist in efficient lab start-up and to learn about common and robust protocols. In the spirit of open sharing of methods, we provide here a description of an entry-level method and user interface (UI) for automating UEM experiments, and we provide access to the source code and scripts (source-available) for ease of implementation or as a starting reference point for those entering or seeking to enter the field (https://github.com/CEMSFlannigan/UEMtomaton/releases/tag/v1.0). Core instrumentation and physical connections in the UEM lab at Minnesota are described. Interface communication schemes consisting of duo server-client pairs between critical components – the optical delay stage and the UEM digital camera – are presented, with emphasis placed on describing the logic and communications sequence designed to conduct automated series acquisitions. An application designed and programmed with C++/CLI as Windows Forms in Microsoft Visual Studio – dubbed UEMtomaton – is also presented. Key to the UI layout is centralization of the automation tasks and establishment of communication within the software rather than by interfacing with each individual workstation. It is our hope that this note provides useful insight for current and future UEM researchers, particularly with respect to generalizability and portability of the approach to emerging labs. Here, we note that while this basic, entry-level approach is certainly not the most sophisticated or comprehensive of those currently in use, we feel there is nevertheless value in clearly communicating a proven straightforward method to hopefully lower the barrier to entry into the field.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

Dose Modeling for Harborview Research and Training Facility

The Harborview Research and Training (R&T) Facility at the University of Washington (UW) was accidently contaminated from a breached sealed source of Cesium-137. The release occurred during removal of the sealed source from its holder within a portable hot cell in the loading dock area near room 220. The loading dock and northeast corner of the building above the ground floor were impacted by the release, whereas the remaining portions of the R&T Facility were minimally affected. Detailed descriptions of the R&T Facility and the affected area/rooms are included in Section 2. The RESRAD Team in the Environmental Science Division of Argonne National Laboratory visited the site on August 12, 2019, and interviewed the R&T Facility building manager and staff working in the building at that time. Site- and building-specific data were collected including building floor plan CAD drawings, building occupants’ working hours, contaminated rooms photos and their dimensions, building air exchange rates, building materials and “curtain wall” details, and so forth. These data are being used in the dose modeling for potential occupants of the building after it is remediated. The dose modeling was performed using the RESRAD-BUILD computer code, which is designed for the evaluation of radiologically contaminated buildings (Yu et al. 2003). Two hypothetical building occupants of the Harborview R&T Facility were evaluated: a researcher occupying room 726 (a research laboratory on the seventh floor) and a custodian occupying room 220 (loading dock and trash compactor room with a rollup door). A maintenance worker that will be spending less time in impacted areas was not modeled because it will be bounded by other two hypothetical building occupants modeled. A summary description of the RESRAD-BUILD is provided in Appendix A, and the modeling of potential doses to a researcher and a custodian is described in detail in Section 3. A summary and a discussion of the dose modeling results are included in Section 4. The parameters used in the RESRAD-BUILD modeling are included in Appendix B. This report documents the dose modeling results; it does not describe measurements of the actual radioactivity levels at the Harborview R&T Facility, nor does it set any expectation for the actual final radioactivity levels. The radiation doses presented are modeling results based on the generic screening values of Cs-137 and several conservative assumptions and are meant only to assist in decision making regarding decontamination of the Harborview R&T Facility.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

Graphics processing unit accelerated phase field dislocation dynamics: Application to bi-metallic interfaces

We present the first high-performance computing implementation of the meso-scale phase field dislocation dynamics (PFDD) model on a graphics processing unit (GPU)-based platform. The implementation takes advantage of the portable OpenACC standard directive pragmas along with Nvidia's compute unified device architecture (CUDA) fast Fourier transform (FFT) library called CUFFT to execute the FFT computations within the PFDD formulation on the same GPU platform. The overall implementation is termed ACCPFDD-CUFFT. The package is entirely performance portable due to the use of OPENACC-CUDA inter-operability, in which calls to CUDA functions are replaced with the OPENACC data regions for a host central processing unit (CPU) and device (GPU). A comprehensive benchmark study has been conducted, which compares a number of FFT routines, the Numerical Recipes FFT (FOURN), Fastest Fourier Transform in the West (FFTW), and the CUFFT. The last one exploits the advantages of the GPU hardware for FFT calculations. The novel ACCPFDD-CUFFT implementation is verified using the analytical solutions for the stress field around an infinite edge dislocation and subsequently applied to simulate the interaction and motion of dislocations through a bi-phase copper-nickel (Cu–Ni) interface. It is demonstrated that the ACCPFDD-CUFFT implementation on a single TESLA K80 GPU offers a 27.6X speedup relative to the serial version and a 5X speedup relative to the 22-multicore Intel Xeon CPU E5-2699 v4 @ 2.20 GHz version of the code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗