SUNDIALS Time Integrators for Exascale Applications with Many Independent ODE Systems
arXiv:2405.01713 · doi:10.1177/10943420241280060
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.
References in corpus (11)
- Nyx: A Massively Parallel AMR Code for Computational Cosmology
- Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers
- AMReX: Block-Structured Adaptive Mesh Refinement for Multiphysics Applications
- Self-Consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations
- Accelerating moderately stiff chemical kinetics in reactive-flow simulations using GPUs
- ARKODE: a flexible IVP solver infrastructure for one-step methods
- Improved Coupling of Hydrodynamics and Nuclear Reactions via Spectral Deferred Corrections
- Enabling GPU Accelerated Computing in the SUNDIALS Time Integration Library
- Modeling the Lyman- forest with Eulerian and SPH hydrodynamical methods
- The ACCEL project: simulating Lyman- forest in large-volume hydrodynamical simulations
- An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks