IMIS | Flanders Marine Institute
 

Flanders Marine Institute

Platform for marine research

IMIS

Publications | Institutes | Persons | Datasets | Projects | Maps
[ report an error in this record ]basket (0): add | show Printer-friendly version

High-order discontinuous Galerkin schemes on general 2D manifolds applied to the shallow water equations
Bernard, P.-E.; Remacle, J.-F.; Comblen, R.; Legat, V.; Hillewaert, K. (2009). High-order discontinuous Galerkin schemes on general 2D manifolds applied to the shallow water equations. J. Comput. Physics 228(17): 6514-6535. hdl.handle.net/10.1016/j.jcp.2009.05.046
In: Journal of Computational Physics. Academic Press: Amsterdam etc.. ISSN 0021-9991, more
Peer reviewed article  

Available in Authors 
    VLIZ: Open Repository 279730 [ OMA ]

Author keywords
    Shallow water equations; High-order finite elements; Discontinuous Galerkin method; Spherical geometry

Authors  Top 
  • Bernard, P.-E., more
  • Remacle, J.-F., more
  • Comblen, R., more
  • Legat, V., more
  • Hillewaert, K.

Abstract
    An innovating approach is proposed to solve vectorial conservation laws on curved manifolds using the discontinuous Galerkin method. This new approach combines the advantages of the usual approaches described in the literature. The vectorial fields are expressed in a unit non-orthogonal local tangent basis derived from the polynomial mapping of curvilinear triangle elements, while the convective flux functions are written is the usual 3D Cartesian coordinate system. The number of vectorial components is therefore minimum and the tangency constraint is naturally ensured, while the method remains robust and general since not relying on a particular parametrization of the manifold. The discontinuous Galerkin method is particularly well suited for this approach since there is no continuity requirement between elements for the tangent basis definition. The possible discontinuities of this basis are then taken into account in the Riemann solver on inter-element interfaces.The approach is validated on the sphere, using the shallow water equations for computing standard atmospheric benchmarks. In particular, the Williamson test cases are used to analyze the impact of the geometry on the convergence rates for discretization error. The propagation of gravity waves is eventually computed on non-conventional irregular curved manifolds to illustrate the robustness and generality of the method.

All data in IMIS is subject to the VLIZ privacy policy Top | Authors