High-order upwind schemes for the wave equation on overlapping grids: Maxwell's equations in second-order form

• Published in: Journal of computational physics 2018-01, Vol.352, p.534-567
• Main Authors: Angel, Jordan B., Banks, Jeffrey W., Henshaw, William D.
• Format: Article
• Language: English
• Subjects: Accuracy; Approximation; Boundary conditions; Computational physics; Conductors; Dielectrics; Dissipation; Electromagnetics; Geometry; Hyperbolic systems; Mathematical analysis; MATHEMATICS AND COMPUTING; Maxwell's equations; Overlapping grids; Scattering; Stability analysis; Upwind methods; Upwind schemes (mathematics); Wave equations
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2017.09.037

High-order accurate upwind approximations for the wave equation in second-order form on overlapping grids are developed. Although upwind schemes are well established for first-order hyperbolic systems, it was only recently shown by Banks and Henshaw [1] how upwinding could be incorporated into the second-order form of the wave equation. This new upwind approach is extended here to solve the time-domain Maxwell's equations in second-order form; schemes of arbitrary order of accuracy are formulated for general curvilinear grids. Taylor time-stepping is used to develop single-step space-time schemes, and the upwind dissipation is incorporated by embedding the exact solution of a local Riemann problem into the discretization. Second-order and fourth-order accurate schemes are implemented for problems in two and three space dimensions, and overlapping grids are used to treat complex geometry and problems with multiple materials. Stability analysis of the upwind-scheme on overlapping grids is performed using normal mode theory. The stability analysis and computations confirm that the upwind scheme remains stable on overlapping grids, including the difficult case of thin boundary grids when the traditional non-dissipative scheme becomes unstable. The accuracy properties of the scheme are carefully evaluated on a series of classical scattering problems for both perfect conductors and dielectric materials in two and three space dimensions. The upwind scheme is shown to be robust and provide high-order accuracy. •Novel upwind schemes for wave equations in second-order form are developed.•Arbitrary-order space-time schemes for overlapping grids are constructed.•Schemes are developed for Maxwell's equations in 2nd-order form.•Stability on overlapping grids is analyzed and confirmed in practice.•Comprehensive numerical results in 2D and 3D confirm stability and accuracy.