In this paper, we present an h-adaptive discontinuous Galerkin formulation of the shallow water equations. For a discontinuous Galerkin scheme using polynomials up to order p, the spatial error of discretization of the method can be shown to be of the order of hp+1, where h is the mesh spacing. It can be shown by rigorous error analysis that the discontinuous Galerkin method discretization error can be related to the amplitude of the inter-element jumps. Therefore, we use the information contained in jumps to build error metrics and size field. Results are presented for ocean modelling problems. A first experiment shows that the theoretical convergence rate is reached with the discontinuous Galerkin high-order h-adaptive method applied to the Stommel wind-driven gyre. A second experiment shows the propagation of an anticyclonic eddy in the Gulf of Mexico.
All data in the Integrated Marine Information System (IMIS) is subject to the VLIZ privacy policy
Web site hosted and maintained by Flanders Marine Institute (VLIZ)
Webmaster info@vliz.be Number of visitors: 2709211 - Total hits: 10371149 (since 2006-01-17)
