High-order triangle-based discontinuous Galerkin methods for hyperbolic equations on a rotating sphere
MetadataShow full item record
High-order triangle-based discontinuous Galerkin (DG) methods for hyperbolic equations on a rotating sphere are presented. The DG method can be characterized its the fusion of finite elements with finite volumes. This DG formulation uses high-order Lagrange polynomials on the triangle using nodal sets up to 15th order. The finite element-type area integrals are evaluated using order 2N Gauss cubature rules. This leads to a full mass matrix which, unlike for continuous Galerkin (CG) methods such as the spectral element (SE) method presented in Giraldo and Warburton [A nodal triangle-based spectral element method for the shallow water equations oil the sphere, J. COMPLIt. Phys. 207 (2005) 129-150], is small, local and efficient to invert. Two types of finite volume-type flux integrals are studied: a set based on Gauss-Lobatto quadrature points (order 2N-1) and a set based oil Gauss quadrature points (order 2N. Furthermore, we explore conservation and advection forms as well as strong and weak forms. Seven test cases are used to compare the different methods including some with scale contractions and shock waves. All three strong forms performed extremely well with the strong conservation form with 2N integration being the most accurate of the four DG methods studied. The strong advection form with 2N integration performed extremely well even for flows with shock waves. The strong conservation form with 2N-1 integration yielded results almost as good as those with 2N while being less expensive. All the DG methods performed better than the SE method for almost all the test cases, especially for those with strong discontinuities. Finally, the DG methods required less computing time than the SE method due to the local nature of the mass matrix. Published by Elsevier Inc.
The article of record as published may be located at http://dx.doi.org/10.1016/j.jcp.2005.09.029
Showing items related by title, author, creator and subject.
Giraldo, F.X. (2003);The Lagrange-Galerkin spectral element method for the two-dimensional shallow water equations is presented. The equations are written in conservation form and the domains are discretized using quadrilateral elements. ...
Simulation of Earthquake Rupture Dynamics in Complex Geometries Using Coupled Finite Difference and Finite Volume Methods O'Reilly, Ossian; Nordstrom, Jan; Kozdon, Jeremy E.; Dunham, Eric M. (2013-10);A numerical method suitable for wave propagation problems in complex geometries is developed for simulating dynamic earthquake ruptures with realistic friction laws. The numerical method couples an unstructured, node-centered ...
Williams, Roger Terry (Monterey, California. Naval Postgraduate School, 1978); NPS63-78-006In this report the Galerkin method is presented and the spectral and finite element methods are shown to be special cases of the general method. The spectral and finite methods are applied to a simple linear equation and ...