Daniel Appelo, Thomas Hagstrom, Arturo Vargas
Arbitrary order dissipative and conservative Hermite methods for the scalar wave equation achieving $\mathcal{O}(2m)$ orders of accuracy using $\mathcal{O}(m^d)$ degrees of freedom per node in $d$ dimensions are presented. Stability and error analyses as well as implementation strategies for accelerators are also given.
Arturo Vargas, Thomas Hagstrom, Jesse Chan et al.
We introduce Hermite-leapfrog methods for first order wave systems. The new Hermite-leapfrog methods pair leapfrog time-stepping with the Hermite methods of Goodrich and co-authors. The new schemes stagger field variables in both time and space and are high-order accurate. We provide a detailed description of the method and demonstrate that the method conserves variable quantities in one-space dimension. Higher dimensional versions of the method are constructed via a tensor product construction. Numerical evidence and rigorous analysis in one space dimension establish stability and high-order convergence. Experiments demonstrating efficient implementations on a graphics processing unit are also presented.
Arturo Vargas, Jesse Chan, Thomas Hagstrom et al.
Hermite methods, as introduced by Goodrich et al., combine Hermite interpolation and staggered (dual) grids to produce stable high order accurate schemes for the solution of hyperbolic PDEs. We introduce three variations of this Hermite method which do not involve time evolution on dual grids. Computational evidence is presented regarding stability, high order convergence, and dispersion/dissipation properties for each new method. Hermite methods may also be coupled to discontinuous Galerkin (DG) methods for additional geometric flexibility. An example illustrates the simplification of this coupling of this coupling for the Hermite methods.