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.
Lu Zhang, Thomas Hagstrom, Daniel Appelo
An energy-based discontinuous Galerkin method for the advective wave equation is proposed and analyzed. Energy-conserving or energy-dissipating methods follow from simple, mesh-independent choices of the inter-element fluxes, and both subsonic and supersonic advection is allowed. Error estimates in the energy norm are established, and numerical experiments on structured grids display optimal convergence in the $L^2$ norm for upwind fluxes. The method generalizes earlier work on energy-based discontinuous Galerkin methods for second order wave equations which was restricted to energy forms written as a simple sum of kinetic and potential energy.
Vladimir Druskin, Murthy Guddati, Thomas Hagstrom
We suggest a unified spectrally matched optimal grid approach for finite-difference and finite-element approximation of the PML. The new approach allows to combine optimal discrete absorption for both evanescent and propagative waves.
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.
Charles L. Epstein, Leslie Greengard, Thomas Hagstrom
We give a principled approach for the selection of a boundary integral, retarded potential representation for the solution of scattering problems for the wave equation in an exterior domain.