Robust Approaches to Handling Complex Geometries with Galerkin Difference Methods
arXiv:1806.06103 · doi:10.1016/j.jcp.2019.04.031
Abstract
The Galerkin difference (GD) basis is a set of continuous, piecewise polynomials defined using a finite difference like grid of degrees of freedom. The one dimensional GD basis functions are naturally extended to multiple dimensions using the tensor product constructions to quadrilateral elements for discretizing partial differential equations. Here we propose two approaches to handling complex geometries using the GD basis within a discontinuous Galerkin finite element setting: (1) using non-conforming, curvilinear GD elements and (2) coupling affine GD elements with curvilinear simplicial elements. In both cases the (semidiscrete) discontinuous Galerkin method is provably energy stable even when variational crimes are committed and in both cases a weight-adjusted mass matrix is used, which ensures that only the reference mass matrix must be inverted. Additionally, we give sufficient conditions on the treatment of metric terms for the curvilinear, nonconforming GD elements to ensure that the scheme is both constant preserving and conservative. Numerical experiments confirm the stability results and demonstrate the accuracy of the coupled schemes.
arXiv comment30 pages, 14 figures
References in corpus (3)
- GPU-accelerated discontinuous Galerkin methods on hybrid meshes
- Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
- An Energy Stable Approach for Discretizing Hyperbolic Equations with Nonconforming Discontinuous Galerkin Methods
Cited by in corpus (4)
- A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
- A weight-adjusted discontinuous Galerkin method for wave propagation in coupled elastic-acoustic media
- Efficient Tensor-Product Spectral-Element Operators with the Summation-by-Parts Property on Curved Triangles and Tetrahedra
- Energy-based discontinuous Galerkin difference methods for second-order wave equations