A Conservative Fully Discrete Numerical Method for the Regularized Shallow Water Wave Equations

Dimitrios Mitsotakis, Hendrik Ranocha, David I. Ketcheson, Endre Süli

Research output: Contribution to journalArticlepeer-review

7 Scopus citations


This paper proposes a new, conservative fully discrete scheme for the numerical solution of the regularized shallow water Boussinesq system of equations in the cases of periodic and reflective boundary conditions. The particular system is one of a class of equations derived recently and can be used in practical simulations to describe the propagation of weakly nonlinear and weakly dispersive long water waves, such as tsunamis. Studies of small-amplitude long waves usually require long-time simulations in order to investigate scenarios such as the overtaking collision of two solitary waves or the propagation of transoceanic tsunamis. For long-time simulations of nondissipative waves such as solitary waves, the preservation of the total energy by the numerical method can be crucial in the quality of the approximation. The new conservative fully discrete method consists of a Galerkin finite element method for spatial semidiscretization and an explicit relaxation Runge--Kutta scheme for integration in time. The Galerkin method is expressed and implemented in the framework of mixed finite element methods. This paper provides an extended experimental study of the accuracy and convergence properties of the new numerical method. The experiments reveal a new convergence pattern compared to the standard, nonconservative Galerkin methods.
Original languageEnglish (US)
Pages (from-to)B508-B537
Number of pages1
JournalSIAM Journal on Scientific Computing
Issue number2
StatePublished - Apr 26 2021

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'A Conservative Fully Discrete Numerical Method for the Regularized Shallow Water Wave Equations'. Together they form a unique fingerprint.

Cite this