Leapfrog Time-Stepping for Hermite Methods

We introduce Hermite-leapfrog methods for first order linear wave systems. The new Hermite-leapfrog methods pair leapfrog time-stepping with the Hermite methods of Goodrich and co-authors et al. (Math Comput 75(254):595–630, 2006). The new schemes stagger field variables in both time and space and are high-order accurate for equations with smooth solutions and coefficients. We provide a detailed description of the method and demonstrate that the method conserves variable quantities. Higher dimensional versions of the method are constructed via tensor products. 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.

[1]  Amir Yefet,et al.  A staggered fourth-order accurate explicit finite difference scheme for the time-domain Maxwell's equations , 2001 .

[2]  Thomas Hagstrom,et al.  Hermite Methods for Aeroacoustics: Recent Progress , 2011 .

[3]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[4]  Xi Chen,et al.  A hybrid Hermite-discontinuous Galerkin method for hyperbolic systems with application to Maxwell's equations , 2014, J. Comput. Phys..

[5]  W. Symes.  Discontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients , 2015 .

[6]  Timothy C. Warburton,et al.  GPU Acceleration of Hermite Methods for the Simulation of Wave Propagation , 2016, ArXiv.

[7]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[8]  Thomas Hagstrom,et al.  Variations on Hermite methods for wave propagation , 2015, 1509.08012.

[9]  Randall J. LeVeque,et al.  Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems , 2007 .

[10]  A. Tarantola,et al.  Two‐dimensional nonlinear inversion of seismic waveforms: Numerical results , 1986 .

[11]  Thomas Hagstrom,et al.  Hermite Methods for the Scalar Wave Equation , 2018, SIAM J. Sci. Comput..

[12]  Maarten van Walstijn,et al.  A comparison of nonstaggered compact FDTD schemes for the 3D wave equation , 2010, 2010 IEEE International Conference on Acoustics, Speech and Signal Processing.

[13]  Daniel Appelö,et al.  Flux-Conservative Hermite Methods for Simulation of Nonlinear Conservation Laws , 2017, J. Sci. Comput..

[14]  David Medina,et al.  OKL: A Unified Language for Parallel Architectures , 2015 .

[15]  Bo Zhang,et al.  An explicit fourth-order staggered finite-difference time-domain method for Maxwell's equations , 2002 .

[16]  Thomas Hagstrom,et al.  Hermite methods for hyperbolic initial-boundary value problems , 2005, Math. Comput..

[17]  T. Hagstrom,et al.  P-adaptive Hermite methods for initial value problems , 2012 .