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

Blaise, Sébastien;St-Cyr, Amik
(2012) Monthly Weather Review — Vol. 104, n° 3, p. 978-996 (2012)

Files

blaise2012.pdf
  • Open Access
  • Adobe PDF
  • 2.72 MB

Details

Authors
  • Blaise, SébastienUCLouvain
    Author
  • St-Cyr, AmikRoyal Dutch Shell
    Author
Abstract
A discontinuous Galerkin model solving the shallow-water equations on the sphere is presented. 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. The proposed multiscale strategy is able to reduce the computational time by an order of magnitude on the tsunami simulation, clearly demonstrating its potential toward multiresolution three-dimensional oceanic and atmospheric applications.
Affiliations

Citations

Blaise, S., & St-Cyr, A. (2012). A dynamic hp-adaptive discontinuous Galerkin method for shallow water flows on the sphere with application to a global tsunami simulation. Monthly Weather Review, 104(3), 978-996. https://doi.org/10.1175/MWR-D-11-00038.1 (Original work published 2012)