Implicit Adaptive Mesh Refinement for Dispersive Tsunami Propagation

Abstract: We present an algorithm to solve the dispersive depth-averaged Serre-Green-Naghdi (SGN) equations using patch-based adaptive mesh refinement. These equations require adding additional higher derivative terms to the nonlinear shallow water equations. This has been implemented as a new component of the open source GeoClaw software that is widely used for modeling tsunamis, storm surge, and related hazards, improving its accuracy on shorter wavelength phenomena. The equations require the solution of an elliptic system at each time step. The adaptive algorithm allows different time steps on different refinement levels, and solves the implicit equations level by level. Computational examples are presented to illustrate the stability and accuracy on a radially symmetric test case and two realistic tsunami modeling problems, including a hypothetical asteroid impact creating a short wavelength tsunami for which dispersive terms are necessary.

Has companion code repository: https://github.com/rjleveque/implicitamr-paper

Mathematics Subject Classification ID

Seismology (including tsunami modeling), earthquakes (86A15) Mesh generation, refinement, and adaptive methods for the numerical solution of initial value and initial-boundary value problems involving PDEs (65M50) Finite volume methods for initial value and initial-boundary value problems involving PDEs (65M08)

This page was built for publication: Implicit Adaptive Mesh Refinement for Dispersive Tsunami Propagation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6511027)