Mario Ricchiuto

Collin Wittenstein, Vincent Marks, Mario Ricchiuto et al.

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.

Hendrik Ranocha, Mario Ricchiuto

We develop structure-preserving numerical methods for the Serre-Green-Naghdi equations, a model for weakly dispersive free-surface waves. We consider both the classical form, requiring the inversion of a non-linear elliptic operator, and a hyperbolic approximation of the equations, allowing fully explicit time stepping. Systems for both flat and variable topography are studied. Our novel numerical methods conserve both the total water mass and the total energy. In addition, the methods for the original Serre-Green-Naghdi equations conserve the total momentum for flat bathymetry. For variable topography, all the methods proposed are well-balanced for the lake-at-rest state. We provide a theoretical setting allowing us to construct schemes of any kind (finite difference, finite element, discontinuous Galerkin, spectral, etc.) as long as summation-by-parts operators are available in the chosen setting. Energy-stable variants are proposed by adding a consistent high-order artificial viscosity term. The proposed methods are validated through a large set of benchmarks to verify all the theoretical properties. Whenever possible, comparisons with exact, reference numerical, or experimental data are carried out. The impressive advantage of structure preservation, and in particular energy preservation, to resolve accurately dispersive wave propagation on very coarse meshes is demonstrated by several of the tests.

Moussa Ziggaf, Davide Torlo, Mario Ricchiuto

We develop arbitrarily high-order, stationarity-preserving stabilized finite element methods for multidimensional nonlinear hyperbolic balance laws on Cartesian grids. We aim at approximating all the steady states of the problem at hand, including non-trivial genuinely multidimensional equilibria, with a level of accuracy higher than the nominal one of the underlying scheme. We formalize more precisely the meaning of stationarity preservation, providing some technical conditions for its realizability. We then recast the multidimensional global-flux quadrature of Barsukow et al. (Num. Meth. PDEs, 2025) as a local preprocessing of the physical fluxes that maps continuous polynomial vector fields to a local space with Raviart--Thomas-type structure. Both the Galerkin and SUPG formulations are recast in this setting. The resulting methods extend the stationarity-preserving finite-volume approach of Barsukow et al. (J. Comput. Phys., 2026) to high-order continuous finite elements and Barsukow et al. (Num. Meth. PDEs, 2025) to nonlinear balance laws. We analyze key properties of the proposed schemes, including local conservation and nodal superconvergence of the discrete steady kernel, and we discuss their relation to low-Mach-compliant discretizations. We apply the framework to the compressible Euler equations with gravity. A simple source-term reformulation yields machine-precision preservation of isothermal hydrostatic equilibria. Extensive numerical benchmarks, including moving equilibrium, near-equilibrium, and instability-dominated regimes, demonstrate clear improvements in robustness and accuracy over standard SUPG and reference finite-volume methods.