High-order upwind schemes for the wave equation on overlapping grids: Maxwell’s equations in second-order form. (English) Zbl 1375.78039
Summary: 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 the second and the third author [ibid. 231, No. 17, 5854–5889 (2012; Zbl 1277.76058)] 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.
MSC:
| 78M20 | Finite difference methods applied to problems in optics and electromagnetic theory |
| 78A25 | Electromagnetic theory (general) |
| 35L05 | Wave equation |
Citations:
Zbl 1277.76058References:
| [1] | Banks, J. W.; Henshaw, W. D., Upwind schemes for the wave equation in second-order form, J. Comput. Phys., 231, 17, 5854-5889 (2012) · Zbl 1277.76058 |
| [2] | Henshaw, W. D., A high-order accurate parallel solver for Maxwell’s equations on overlapping grids, SIAM J. Sci. Comput., 28, 5, 1730-1765 (2006) · Zbl 1127.78011 |
| [3] | Yee, K. S., Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media, IEEE Trans. Antennas Propag., 14, 3, 302-307 (1966) · Zbl 1155.78304 |
| [4] | Nédélec, J., Mixed finite elements in \(R^3\), Numer. Math., 35, 3, 315-341 (1980) · Zbl 0419.65069 |
| [5] | Courant, R.; Isaacson, E.; Rees, M., On the solution of nonlinear hyperbolic differential equations by finite differences, Commun. Pure Appl. Math., 5, 243-255 (1952) · Zbl 0047.11704 |
| [6] | Godunov, S. K., Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics, Mat. Sb., 47, 271-306 (1959) · Zbl 0171.46204 |
| [7] | Boris, J. P.; Book, D. L., Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works, J. Comput. Phys., 11, 38-69 (1973) · Zbl 0251.76004 |
| [8] | Shoucri, M.; Knorr, G., Numerical integration of the Vlasov equation, J. Comput. Phys., 14, 1, 84-92 (1974) · Zbl 0275.65029 |
| [9] | Cheng, C. Z.; Knorr, G., The integration of the Vlasov equation in configuration space, J. Comput. Phys., 22, 3, 330-351 (1976) |
| [10] | Robert, A., A semi-Lagrangian and semi-implicit numerical integration scheme for the primitive meteorological equations, J. Meteorol. Soc. Jpn., Ser. II, 60, 1, 319-325 (1982) |
| [11] | Colella, P.; Woodward, P. R., The piecewise parabolic method (PPM) for gas-dynamical simulations, J. Comput. Phys., 54, 1, 174-201 (1984) · Zbl 0531.76082 |
| [12] | Harten, A.; Engquist, B.; Osher, S.; Chakravarthy, S., Uniformly high order accurate essentially non-oscillatory schemes, III, J. Comput. Phys., 71, 231-303 (1987) · Zbl 0652.65067 |
| [13] | Shu, C.-W.; Osher, S., Efficient implementation of essentially non-oscillatory shock-capturing schemes, J. Comput. Phys., 77, 2, 439-471 (1988) · Zbl 0653.65072 |
| [14] | Cockburn, B.; Shu, C. W., TVB Runge-Kutta local projection discontinuous Galerkin finite-element method for conservation-laws 2: general framework, Math. Comput., 52, 411-435 (1989) · Zbl 0662.65083 |
| [15] | Jiang, G.-S.; Shu, C.-W., Efficient implementation of weighted ENO schemes, J. Comput. Phys., 126, 1, 202-228 (1996) · Zbl 0877.65065 |
| [16] | Cioni, J.; Fezoui, L.; Steve, H., A parallel time-domain Maxwell solver using upwind schemes and triangular meshes, Impact Comput. Sci. Eng., 5, 3, 215-247 (1993) · Zbl 0788.65119 |
| [17] | Balsara, D. S.; Amano, T.; Garain, S.; Kim, J., A high-order relativistic two-fluid electrodynamic scheme with consistent reconstruction of electromagnetic fields and a multidimensional Riemann solver for electromagnetism, J. Comput. Phys., 318, 169-200 (2016) · Zbl 1349.76425 |
| [18] | Balsara, D. S.; Käppeli, R., Von Neumann stability analysis of globally divergence-free RKDG schemes for the induction equation using multidimensional Riemann solvers, J. Comput. Phys., 336, 104-127 (2017) · Zbl 1375.76212 |
| [19] | Christensen, R., Godunov methods on a staggered mesh - an improved artificial viscosity (1990), Lawrence Livermore National Laboratory, Research Report UCRL-JC-105269 |
| [20] | Richtmyer, R. D., Proposed numerical method for calculation of shocks (1948), Los Alamos National Laboratory, Tech. Rep. LA-671 |
| [21] | Richtmyer, R. D., Proposed numerical method for calculation of shocks, II (1948), Los Alamos National Laboratory, Tech. Rep. LA-699 |
| [22] | VonNeumann, J.; Richtmyer, R. D., A method for the numerical calculation of hydrodynamic shocks, J. Appl. Phys., 21, 232-243 (1950) · Zbl 0037.12002 |
| [23] | Banks, J. W.; Henshaw, W. D., High-order upwind methods for wave equations on curvilinear and overlapping grids, (Spectral and High Order Methods for Partial Differential Equations - ICOSAHOM 14. Spectral and High Order Methods for Partial Differential Equations - ICOSAHOM 14, Lect. Notes Comput. Sci. Eng., vol. 107 (2015), Springer), 137-145 · Zbl 1352.65226 |
| [24] | Appelö, D.; Banks, J. W.; Henshaw, W. D.; Schwendeman, D. W., Numerical methods for solid mechanics on overlapping grids: linear elasticity, J. Comput. Phys., 231, 6012-6050 (2012) · Zbl 1277.74005 |
| [25] | Appelö, D.; Hagstrom, T., A new discontinuous Galerkin formulation for wave equations in second-order form, SIAM J. Numer. Anal., 53, 6, 2705-2726 (2015) · Zbl 1330.65145 |
| [26] | Hesthaven, J. S., High-order accurate methods in time-domain computational electromagnetics: a review, Adv. Imaging Electron Phys., 127, 59-123 (2003) |
| [27] | Taflove, A.; Hagness, S. C., Computational Electrodynamics: The Finite-Difference Time-Domain Method (2000), Artech House · Zbl 0963.78001 |
| [28] | Cohen, G. C., Higher-Order Numerical Methods for Transient Wave Equations (2002), Springer: Springer New York · Zbl 0985.65096 |
| [29] | Driscoll, T. A.; Fornberg, B., Block pseudospectral methods for Maxwell’s equations II: two-dimensional, discontinuous coefficients case, SIAM J. Sci. Comput., 21, 3, 1146-1167 (1999) · Zbl 0949.65107 |
| [30] | Kim, I.; Hoefer, W. J.R., A local mesh refinement algorithm for the time domain-finite difference method using Maxwell’s curl equations, IEEE Trans. Microw. Theory Tech., 38, 6, 812-815 (1990) |
| [31] | Zakharian, A. R.; Brio, M.; Dineen, C.; Moloney, J. V., Second-order accurate FDTD space and time grid refinement method in three space dimensions, IEEE Photonics Technol. Lett., 18, 11, 1237-1239 (2006) |
| [32] | Balsara, D. S., Divergence-free adaptive mesh refinement for magnetohydrodynamics, J. Comput. Phys., 174, 2, 614-648 (2001) · Zbl 1157.76369 |
| [33] | Ismagilov, T. Z., Second order finite volume scheme for Maxwell’s equations with discontinuous electromagnetic properties on unstructured meshes, J. Comput. Phys., 282, 33-42 (2015) · Zbl 1352.65280 |
| [34] | Barbas, A.; Velarde, P., Development of a Godunov method for Maxwell’s equations with adaptive mesh refinement, J. Comput. Phys., 300, 186-201 (2015) · Zbl 1349.78079 |
| [35] | Kreiss, H.-O., Stability theory of difference approximations of mixed initial boundary value problems. I, Math. Comput., 22, 703-714 (1968) · Zbl 0197.13704 |
| [36] | Gustafsson, B.; Kreiss, H.-O.; Sundström, A., Stability theory of difference approximations for mixed initial boundary value problems. II, Math. Comput., 26, 119, 649-686 (1972) · Zbl 0293.65076 |
| [37] | Chesshire, G.; Henshaw, W., Composite overlapping meshes for the solution of partial differential equations, J. Comput. Phys., 90, 1-64 (1990) · Zbl 0709.65090 |
| [38] | Fornberg, B., A Practical Guide to Pseudospectral Methods (1996), Cambridge University Press: Cambridge University Press Cambridge · Zbl 0844.65084 |
| [39] | Gustafsson, B.; Kreiss, H.-O.; Oliger, J., Time Dependent Problems and Difference Methods (1995), John Wiley & Sons: John Wiley & Sons New York · Zbl 0843.65061 |
| [40] | Henshaw, W. D.; Schwendeman, D. W., An adaptive numerical scheme for high-speed reactive flow on overlapping grids, J. Comput. Phys., 191, 2, 420-447 (2003) · Zbl 1134.76427 |
| [41] | Strikwerda, J. C., Initial boundary value problems for the method of lines, J. Comput. Phys., 34, 1, 94-107 (1980) · Zbl 0441.65085 |
| [42] | Buchberger, B., Gröbner-Bases: An Algorithmic Method in Polynomial Ideal Theory (1985), Reidel Publishing Company: Reidel Publishing Company Dordrecht, Boston, Lancaster · Zbl 0587.13009 |
| [43] | Li, T. Y., Numerical solution of multivariate polynomial systems by homotopy continuation methods, Acta Numer., 6, 399-436 (1997) · Zbl 0886.65054 |
| [44] | J. Verschelde, Polynomial homotopy continuation with PHCpack, 2010.; J. Verschelde, Polynomial homotopy continuation with PHCpack, 2010. · Zbl 1308.68198 |
| [45] | Balanis, C. A., Advanced Engineering Electromagnetics (1989), John Wiley & Sons |
| [46] | Bowman, J. J.; Senior, T. B.A.; Uslenghi, P. L.E., Electromagnetic and Acoustic Scattering by Simple Shapes (1996), John Wiley and Sons, Hemisphere Publishing Corporation |
| [47] | van de Hulst, H. C., Light Scattering by Small Particles (1957), Dover Publications, Inc. |
| [48] | Hagstrom, T., Radiation boundary conditions for the numerical simulation of waves, Acta Numer., 8, 47-106 (1999) · Zbl 0940.65108 |
| [49] | Henshaw, W. D., CgMx User Guide: an Overture Solver for Maxwell’s Equations on Composite Grids (2015) |
| [50] | Engquist, B.; Majda, A., Absorbing boundary conditions for numerical simulation of waves, Proc. Natl. Acad. Sci., 74, 5, 1765-1766 (1977) · Zbl 0378.76018 |
| [51] | Henshaw, W. D.; Schwendeman, D. W., Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement, J. Comput. Phys., 227, 16, 7469-7502 (2008) · Zbl 1213.76138 |
| [52] | Banks, J. W.; Henshaw, W. D.; Shadid, J. N., An evaluation of the FCT method for high-speed flows on structured overlapping grids, J. Comput. Phys., 228, 15, 5349-5369 (2009) · Zbl 1280.76031 |
| [53] | Henshaw, W. D., Mappings for Overture, a description of the mapping class and documentation for many useful mappings (1998), Lawrence Livermore National Laboratory, Research Report UCRL-MA-132239 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
