GPU-Accelerated Energy-Conserving Methods for the Two-Dimensional Hyperbolized Serre-Green-Naghdi Equations
Abstract
We develop energy-conserving numerical methods for a two-dimensional hyperbolic approximation of the Serre-Green-Naghdi equations with variable bathymetry and either periodic or reflecting boundary conditions. The hyperbolic formulation avoids the costly inversion of an elliptic operator present in the classical model. Our schemes combine split forms with summation-by-parts (SBP) operators to construct semi-discretizations that conserve the total water mass and the total energy. We provide analytical proofs of these conservation properties and also verify them numerically. While the framework is general, our implementation focuses on second-order finite-difference SBP operators. The methods are implemented in Julia for CPU and GPU architectures (AMD and NVIDIA) and achieve substantial speedups on modern accelerators. We validate the approach through convergence studies based on solitary-wave and manufactured-solution tests, and by comparisons to analytical, experimental, and existing numerical results. All source code to reproduce our results is available online.
Cite
@article{arxiv.2601.02540,
title = {GPU-Accelerated Energy-Conserving Methods for the Two-Dimensional Hyperbolized Serre-Green-Naghdi Equations},
author = {Collin Wittenstein and Vincent Marks and Mario Ricchiuto and Hendrik Ranocha},
journal= {arXiv preprint arXiv:2601.02540},
year = {2026}
}
Comments
46 pages, 29 figures. Collin Wittenstein and Vincent Marks contributed equally to this work