A dynamic hp-adaptive discontinuous Galerkin method for shallow-water flows on the sphere with application to a global tsunami simulation

It captures the dynamically varying key aspects of the flows by having the advantageous ability to locally modify the mesh as well as the order of interpolation within each element. The computational load is efficiently distributed among processors in parallel using a weighted recursive coordinate bisection strategy.

A simple error estimator, based on the discontinuity of the variables at the interfaces between elements, is used to select the elements to be refined or coarsened. The flows are expressed in three-dimensional Cartesian coordinates, but tangentially constrained to the sphere by adding a Lagrange multiplier to the system of equations. The model is validated on classic atmospheric test cases and on the simulation of the February 2010 Chilean tsunami propagation.