SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations
Cody J. Balos (Lawrence Livermore National Laboratory)
Marcus Day (National Renewable Energy Laboratory)
Lucas Esclapez (National Renewable Energy Laboratory)
Anne M. Felden (TU Delft - Mechanical Engineering)
David J. Gardner (Lawrence Livermore National Laboratory)
Malik Hassanaly (National Renewable Energy Laboratory)
Daniel R. Reynolds (Southern Methodist University)
Jon S. Rood (National Renewable Energy Laboratory)
Jean M. Sexton (Lawrence Berkeley National Laboratory)
Nicholas T. Wimer (National Renewable Energy Laboratory)
Carol S. Woodward (Lawrence Livermore National Laboratory)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world’s largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger “outer-loop” time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.