Physics-Aware Adaptive Checkpointing with Shadow Systems for Nonlinear PDE Simulations
Authors:
- Qian GongQGQian Gong
- Mark AinsworthMAMark Ainsworth
- Jieyang ChenJCJieyang Chen
- Jaemoon LeeJLJaemoon Lee
- Scott KlaskySKScott Klasky
Preprint
- 5 February 2026
- preprint
- Published by Elsevier in SSRN Electronic Journal
Abstract
Large-scale simulations of nonlinear partial differential equations (PDEs) that exhibit strongly transient behavior and pattern-forming dynamics produce enormous amount of data, which, even with modern storage system, cannot be stored for later curation. Current I/O strategies either write dense time series of snapshots, which is often prohibitive in I/O and storage, or store a few checkpoints that enable restart but incur expensive recomputation cost and provide no control on post-restart error growth, especially when lossy compression is used. Moreover, most if not all existing strategies take no account of the actual physical state of the system. Here, we present a simple physics-aware I/O framework in which a lowcost shadow system adaptively triggers lossy checkpoints when the shadow system deviates from the fine-scale simulation. The shadow system can be a coarsened replica of the fine-scale simulation that evolves concurrently. This means that checkpoints are taken based on the physical state of the system: fewer checkpoints are triggered when the system is quiescent while more are taken when the system undergoes a rapid change. This type of behavior is observed in many systems such as Brusselator and FitzHugh-Nagumo. We illustrate that our framework maintains stable restarts, keeps fine-scale restart errors bounded by shadow errors, and reconstructs the time history with significantly lower error and storage than interpolating fixed-interval snapshots, at only modest additional runtime overhead.Keywords
This publication has 39 references indexed in Scilit:
- MGARD: A multigrid framework for high-performance, error-controlled data compression and refactoringSoftwareX, 2023
- An unstructured adaptive mesh refinement approach for computational fluid dynamics of reacting flowsJournal of Computational Physics, 2022
- A popular reaction-diffusion model fractional Fitzhugh-Nagumo equation: analytical and numerical treatmentApplied Mathematics-A Journal of Chinese Universities, 2021
- A fully adaptive nonintrusive reduced-order modelling approach for parametrized time-dependent problemsComputer Methods in Applied Mechanics and Engineering, 2021
- A computational modeling of two dimensional reaction–diffusion Brusselator system arising in chemical processesJournal of Mathematical Chemistry, 2014
- Adaptive mesh refinement for characteristic gridsGeneral Relativity and Gravitation, 2010
- Adaptive mesh refinement for coupled elliptic-hyperbolic systemsJournal of Computational Physics, 2006
- A second-order scheme for the “Brusselator” reaction–diffusion systemJournal of Mathematical Chemistry, 1999
- Analytical and Simulation Results for Stochastic Fitzhugh-Nagumo Neurons and Neural NetworksJournal of Computational Neuroscience, 1998
- A basic norm equivalence for the theory of multilevel methodsNumerische Mathematik, 1993