×

Von Neumann stability analysis of globally constraint-preserving DGTD and PNPM schemes for the Maxwell equations using multidimensional Riemann solvers. (English) Zbl 1416.78020

Summary: The time-dependent equations of computational electrodynamics (CED) are evolved consistent with the divergence constraints on the electric displacement and magnetic induction vector fields. Respecting these constraints has proved to be very useful in the classic finite-difference time-domain (FDTD) schemes. As a result, there has been a recent effort to design finite volume time domain (FVTD) and discontinuous Galerkin time domain (DGTD) schemes that satisfy the same constraints and, nevertheless, draw on recent advances in higher order Godunov methods. This paper catalogues the first step in the design of globally constraint-preserving DGTD schemes. The algorithms presented here are based on a novel DG-like method that is applied to a Yee-type staggering of the electromagnetic field variables in the faces of the mesh. The other two novel building blocks of the method include constraint-preserving reconstruction of the electromagnetic fields and multidimensional Riemann solvers; both of which have been developed in recent years by the first author. The resulting DGTD scheme is linear, at least when limiters are not applied to the DG scheme. As a result, it is possible to carry out a von Neumann stability analysis of the entire suite of DGTD schemes for CED at orders of accuracy ranging from second to fourth. The analysis requires some simplifications in order to make it analytically tractable, however, it proves to be extremely instructive. A von Neumann stability analysis is a necessary precursor to the design of a full DGTD scheme for CED. It gives us the maximal CFL numbers that can be sustained by the DGTD schemes presented here at all orders. It also enables us to understand the wave propagation characteristics of the schemes in various directions on a Cartesian mesh. We find that constraint-preserving DGTD schemes permit CFL numbers that are competitive with conventional DG schemes. However, like conventional DG schemes, the CFL of DGTD schemes decreases with increasing order. To counteract that, we also present constraint-preserving PNPM schemes for CED. We find that the third and fourth order constraint-preserving DGTD and P1PM schemes have some extremely attractive properties when it comes to low-dispersion, low-dissipation propagation of electromagnetic waves in multidimensions. Numerical accuracy tests are also provided to support the von Neumann stability analysis. We expect these methods to play a role in those problems of engineering CED where exceptional precision must be achieved at any cost.

MSC:

78M10 Finite element, Galerkin and related methods applied to problems in optics and electromagnetic theory
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs

Software:

RIEMANN
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Balsara, D. S., Divergence-free adaptive mesh refinement for magnetohydrodynamics, J. Comput. Phys., 174, 614-648 (2001) · Zbl 1157.76369
[2] Balsara, D. S., Second-order-accurate schemes for magnetohydrodynamics with divergence-free reconstruction, Astrophys. J. Suppl. Ser., 151, 149-184 (2004)
[3] Balsara, D. S., Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics, J. Comput. Phys., 228, 5040-5056 (2009) · Zbl 1280.76030
[4] Balsara, D. S., Multidimensional HLLE Riemann solver; application to Euler and magnetohydrodynamic flows, J. Comput. Phys., 229, 1970-1993 (2010) · Zbl 1303.76140
[5] Balsara, D. S., A two-dimensional HLLC Riemann solver for conservation laws: application to Euler and magnetohydrodynamic flows, J. Comput. Phys., 231, 7476-7503 (2012) · Zbl 1284.76261
[6] Balsara, D. S.; Dumbser, M.; Abgrall, R., Multidimensional HLL and HLLC Riemann solvers for unstructured meshes - with application to Euler and MHD flows, J. Comput. Phys., 261, 172-208 (2014) · Zbl 1349.76426
[7] Balsara, D. S., Multidimensional Riemann problem with self-similar internal structure - part I - application to hyperbolic conservation laws on structured meshes, J. Comput. Phys., 277, 163-200 (2014) · Zbl 1349.76303
[8] Balsara, D. S.; Dumbser, M., Divergence-free MHD on unstructured meshes using high order finite volume schemes based on multidimensional Riemann solvers, J. Comput. Phys., 299, 687-715 (2015) · Zbl 1351.76092
[9] Balsara, D. S.; Dumbser, M., Multidimensional Riemann problem with self-similar internal structure - part II - application to hyperbolic conservation laws on unstructured meshes, J. Comput. Phys., 287, 269-292 (2015) · Zbl 1351.76091
[10] Balsara, D. S., Three dimensional HLL Riemann solver for structured meshes; application to Euler and MHD flow, J. Comput. Phys., 295, 1-23 (2015) · Zbl 1349.76584
[11] Balsara, D. S.; Vides, J.; Gurski, K.; Nkonga, B.; Dumbser, M.; Garain, S.; Audit, E., A two-dimensional Riemann solver with self-similar sub-structure - alternative formulation based on least squares projection, J. Comput. Phys., 304, 138-161 (2016) · Zbl 1349.76157
[12] Balsara, D. S.; Amano, T.; Garain, S.; Kim, J., High order accuracy divergence-free scheme for the electrodynamics of relativistic plasmas with multidimensional Riemann solvers, J. Comput. Phys., 318, 169-200 (2016) · Zbl 1349.76425
[13] Balsara, D. S.; Käppeli, R., von Neumann stability analysis of globally divergence-free RKDG and PNPM schemes for the induction equation using multidimensional Riemann solvers, J. Comput. Phys., 336, 104-127 (2017) · Zbl 1375.76212
[14] Balsara, D. S.; Taflove, A.; Garain, S.; Montecinos, G., Computational electrodynamics in material media with constraint-preservation, multidimensional Riemann solvers and sub-cell resolution - part I, second-order FVTD schemes, J. Comput. Phys., 349, 604-635 (2017) · Zbl 1380.78014
[15] Balsara, D. S.; Taflove, A.; Garain, S.; Montecinos, G., Computational electrodynamics in material media with constraint-preservation, multidimensional Riemann solvers and sub-cell resolution - part II, higher-order FVTD schemes, J. Comput. Phys., 354, 613-645 (2018) · Zbl 1380.78013
[16] Balsara, D. S.; Nkonga, B., Formulating multidimensional Riemann solvers in similarity variables - part III: a multidimensional analogue of the HLLI Riemann solver for conservative hyperbolic systems, J. Comput. Phys., 346, 25-48 (2017) · Zbl 1378.76056
[17] Barbas, A.; Velarde, P., Development of a Godunov method for Maxwell’s equations with adaptive mesh refinement, J. Comput. Phys., 300, 186-201 (2016) · Zbl 1349.78079
[18] Bokil, V.; Chang, Y.; Jiang, Y.; Li, F., Energy stable discontinuous Galerkin methods for Maxwell’s equations in nonlinear optical media, J. Comput. Phys., 350, 420-452 (2017) · Zbl 1380.78011
[19] Cockburn, B.; Shu, C.-W., TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws, II: general framework, Math. Comput., 52, 411-435 (1989) · Zbl 0662.65083
[20] Cockburn, B.; Hou, S.; Shu, C.-W., TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws, IV: the multidimensional case, J. Comput. Phys., 54, 545-581 (1990) · Zbl 0695.65066
[21] Cockburn, B.; Shu, C.-W., The Runge-Kutta discontinuous Galerkin method for conservation laws, V: multidimensional systems, J. Comput. Phys., 141, 199-224 (1998) · Zbl 0920.65059
[22] Cockburn, B.; Shu, C.-W., Runge-Kutta discontinuous Galerkin methods for convection dominated problems, J. Sci. Comput., 16, 173-261 (2001) · Zbl 1065.76135
[23] Cockburn, B.; Li, F.; Shu, C.-W., Locally discontinuous Galerkin methods for the Maxwell equations, J. Comput. Phys., 194, 588-610 (2004) · Zbl 1049.78019
[24] Dumbser, M.; Balsara, D. S.; Toro, E. F.; Munz, C.-D., A unified framework for the construction of quadrature-free one-step finite-volume and discontinuous Galerkin schemes, J. Comput. Phys., 227, 8209-8253 (2008) · Zbl 1147.65075
[25] Egger, H.; Kretzschmar, F.; Schnepp, S. M.; Tsukerman, I.; Weiland, T., Transparent boundary conditions for a discontinuous Galerkin Trefftz method, Appl. Math. Comput., 267, 42-55 (2015) · Zbl 1410.78027
[26] Gottlieb, S.; Shu, C.-W.; Tadmor, E., Strong stability-preserving higher order time discretization methods, SIAM Rev., 43, 1, 89-112 (2001) · Zbl 0967.65098
[27] Hesthaven, J.; Warburton, T., Nodal high-order methods on unstructured grids, I: time-domain solution of Maxwell’s equations, J. Comput. Phys., 181, 186-221 (2002) · Zbl 1014.78016
[28] 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
[29] Kretzschmar, F.; Schnepp, S. M.; Tsukerman, I.; Weiland, T., Discontinuous Galerkin methods with Trefftz approximations, Am. J. Comput. Appl. Math., 270, 211-222 (2014) · Zbl 1321.65155
[30] Liu, Y.; Shu, C.-W.; Tadmor, E.; Zhang, M., L2 stability analysis of the central discontinuous Galerkin method and comparison between the central and regular discontinuous Galerkin methods, Math. Model. Numer. Anal., 42, 593-607 (2008) · Zbl 1152.65095
[31] Munz, C.-D.; Omnes, P.; Schneider, R.; Sonnendrücker, E.; Voß, U., Divergence correction techniques for Maxwell solvers based on a hyperbolic model, J. Comput. Phys., 161, 484 (2000) · Zbl 0970.78010
[32] Reed, W. H.; Hill, T. R., Triangular Mesh Methods for the Neutron Transport Equation (1973), Los Alamos Scientific Laboratory: Los Alamos Scientific Laboratory Los Alamos, Technical Report LA-UR-73-479
[33] Spiteri, R. J.; Ruuth, S. J., A new class of optimal high-order strong-stability-preserving time-stepping schemes, SIAM J. Numer. Anal., 40, 469-491 (2002) · Zbl 1020.65064
[34] Spiteri, R. J.; Ruuth, S. J., Non-linear evolution using optimal fourth-order strong-stability-preserving Runge-Kutta methods, Math. Comput. Simul., 62, 125-135 (2003) · Zbl 1015.65031
[35] Shu, C.-W.; Osher, S. J., Efficient implementation of essentially non-oscillatory shock capturing schemes, J. Comput. Phys., 77, 439-471 (1988) · Zbl 0653.65072
[36] Shu, C.-W.; Osher, S. J., Efficient implementation of essentially non-oscillatory shock capturing schemes, II, J. Comput. Phys., 83, 32-78 (1989) · Zbl 0674.65061
[37] Taflove, A., Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures, Wave Motion, 10, 547-582 (1988) · Zbl 0672.73028
[38] Taflove, A.; Hagness, S., Computational Electrodynamics (2005), Artech House
[39] Taflove, A.; Oskooi, A.; Johnson, S., Advances in FDTD Computational Electrodynamics - Photonics and Nanotechnology (2013), Artech House
[40] Taflove, A.; Hagness, S., Finite difference time domain solution of Maxwell’s equations, WEEEE Rev., 8303 (2016)
[41] Yee, K. S., Numerical solution of initial boundary value problems involving Maxwell equation in an isotropic media, IEEE Trans. Antennas Propag., 14, 302 (1966) · Zbl 1155.78304
[42] Xu, Z.; Balsara, D. S.; Du, H., Divergence-free WENO reconstruction-based finite volume scheme for ideal MHD equations on triangular meshes, Commun. Comput. Phys., 19, 04, 841-880 (2016) · Zbl 1373.76149
[43] Yang, H.; Li, F., Stability analysis and error estimates of an exactly divergence-free method for the magnetic induction equations, Modél. Math. Anal. Numér., 50, 4, 965-993 (2016) · Zbl 1348.78028
[44] Zhang, M.; Shu, C.-W., An analysis of and a comparison between the discontinuous Galerkin and the spectral finite volume methods, Comput. Fluids, 34, 581-592 (2005) · Zbl 1138.76391
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.