Leapfrog time-stepping for Hermite methods

From MaRDI portal



Abstract: 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.


The authors consider initial-boundary value problems for a class of linear first-order hyperbolic systems. Finite difference time domain methods yield numerical solutions of these partial differential equations. \textit{J. Goodrich} et al., [Math. Comput. 75, No. 254, 595--630 (2006; Zbl 1103.35065)] used Hermite interpolation and a staggered grid to derive methods of high order. Now, the authors investigate a variant, which combines Hermite interpolation, leapfrog time-stepping, and staggered grids in both space and time. The Hermite interpolation operator applies derivative data up to some order \(m\). It follows that the method is consistent of order (at least) \(2m\), provided that the exact solution is sufficiently smooth. The authors show the stability of the method in one as well as in multiple space dimensions, where the Fourier transform of the solution is applied. The convergence is proven only in one space dimension. Results of numerical computations are presented for wave equations in one space dimension, which confirm the predicted convergence orders. In the case of even \(m\), the estimated convergence rates indicate an actually higher order of \(2m+2\). Moreover, wave equations are considered in two space dimensions, where the results indicate that the method is also convergent with rates as in the one dimensional case. Furthermore, the authors examine computations on graphics processing units in an example with three space dimensions. The performance demonstrates a potential for accelerated computations.











This page was built for publication: Leapfrog time-stepping for Hermite methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2316219)