

















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.
| Comments: | 46 pages, 29 figures. Collin Wittenstein and Vincent Marks contributed equally to this work |
| Subjects: | Numerical Analysis (math.NA) |
| MSC classes: | 65M06, 65M12, 65M20 |
| Cite as: | arXiv:2601.02540 [math.NA] |
| (or arXiv:2601.02540v2 [math.NA] for this version) | |
| https://doi.org/10.48550/arXiv.2601.02540 arXiv-issued DOI via DataCite |
From: Collin Wittenstein [view email]
[v1]
Mon, 5 Jan 2026 20:26:56 UTC (5,807 KB)
[v2]
Tue, 26 May 2026 14:43:38 UTC (7,729 KB)
此内容由惯性聚合(RSS阅读器)自动聚合整理,仅供阅读参考。 原文来自 — 版权归原作者所有。