Journal
JOURNAL OF COMPUTATIONAL PHYSICS
Volume 228, Issue 17, Pages 6514-6535Publisher
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2009.05.046
Keywords
Shallow water equations; High-order finite elements; Discontinuous Galerkin method; Spherical geometry
Ask authors/readers for more resources
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. (C) 2009 Elsevier Inc. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available