A New Family of Exponential-Based High-Order DGTD Methods for Modeling 3-D Transient Multiscale Electromagnetic Problems

A new family of exponential-based time integration methods are proposed for the time-domain Maxwell’s equations discretized by a high-order discontinuous Galerkin (DG) scheme formulated on locally refined unstructured meshes. These methods, which are developed from the Lawson method, remove the stiffness on the time explicit integration of the semidiscrete operator associated with the fine part of the mesh, and allow for the use of high-order time explicit scheme for the coarse part operator. They combine excellent stability properties with the ability to obtain very accurate solutions even for very large time step sizes. Here, the explicit time integration of the Lawson-transformed semidiscrete system relies on a low-storage Runge-Kutta (LSRK) scheme, leading to a combined Lawson-LSRK scheme. In addition, efficient techniques are presented to further improve the efficiency of this exponential-based time integration. For the efficient calculation of matrix exponential, we employ the Krylov subspace method. Numerical experiments are presented to assess the stability, verify the accuracy, and numerical convergence of the Lawson-LSRK scheme. They also demonstrate that the DG time-domain methods based on the proposed time integration scheme can be much faster than those based on classical fully explicit time stepping schemes, with the same accuracy and moderate memory usage increase on locally refined unstructured meshes, and are thus very promising for modeling 3-D multiscale electromagnetic problems.

[1]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[2]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[3]  Din-Kow Sun,et al.  Construction of Nearly Orthogonal Nedelec Bases for Rapid Convergence with Multilevel Preconditioned Solvers , 2001, SIAM J. Sci. Comput..

[4]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[5]  S. Lanteri,et al.  An Unconditionally Stable Discontinuous Galerkin Method for Solving the 2-D Time-Domain Maxwell Equations on Unstructured Triangular Meshes , 2008, IEEE Transactions on Magnetics.

[6]  L. Fezoui,et al.  Convergence and stability of a discontinuous galerkin time-domain method for the 3D heterogeneous maxwell equations on unstructured meshes , 2005 .

[7]  Jin-Fa Lee,et al.  Interconnect and lumped elements modeling in interior penalty discontinuous Galerkin time-domain methods , 2010, J. Comput. Phys..

[8]  Jian-Ming Jin,et al.  The Finite Element Method in Electromagnetics , 1993 .

[9]  Bin Li,et al.  An Inverse-Based Multifrontal Block Incomplete LU Preconditioner for the 3-D Finite-Element Eigenvalue Analysis of Lossy Slow-Wave Structures , 2015, IEEE Transactions on Microwave Theory and Techniques.

[10]  J. Williamson.  Low-storage Runge-Kutta schemes , 1980 .

[11]  M. Okoniewski,et al.  Three-dimensional subgridding algorithm for FDTD , 1997 .

[12]  Qing Huo Liu,et al.  Discontinuous Galerkin Time-Domain Methods for Multiscale Electromagnetic Simulations: A Review , 2013, Proceedings of the IEEE.

[13]  V. Lebedev,et al.  Explicit difference schemes for solving stiff systems of ODEs and PDEs with complex spectrum , 1998 .

[14]  Marcus J. Grote,et al.  Explicit local time-stepping methods for Maxwell's equations , 2010, J. Comput. Appl. Math..

[15]  Joseph W. H. Liu,et al.  The Multifrontal Method for Sparse Matrix Solution: Theory and Practice , 1992, SIAM Rev..

[16]  B. Minchev,et al.  A review of exponential integrators for first order semi-linear problems , 2005 .

[17]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[18]  Andreas Sturm,et al.  Error Analysis of a Second-Order Locally Implicit Method for Linear Maxwell's Equations , 2016, SIAM J. Numer. Anal..

[19]  S. Piperno.  Symplectic local time-stepping in non-dissipative DGTD methods applied to wave propagation problems , 2006 .

[20]  S.,et al.  Numerical Solution of Initial Boundary Value Problems Involving Maxwell’s Equations in Isotropic Media , 1966 .

[21]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[22]  J. Verwer,et al.  Unconditionally stable integration of Maxwell's equations , 2009 .

[23]  Jan Verwer.  Component splitting for semi-discrete Maxwell equations , 2011 .

[24]  E. Montseny,et al.  Dissipative terms and local time-stepping improvements in a spatial high order Discontinuous Galerkin scheme for the time-domain Maxwell's equations , 2008, J. Comput. Phys..

[25]  Y. Saad,et al.  Numerical Methods for Large Eigenvalue Problems , 2011 .

[26]  Marco Vianello,et al.  Efficient approximation of the exponential operator for discrete 2D advection–diffusion problems , 2003, Numer. Linear Algebra Appl..

[27]  J. D. Lawson.  Generalized Runge-Kutta Processes for Stable Systems with Large Lipschitz Constants , 1967 .

[28]  Bin Li,et al.  Windows simulator: An advanced three-dimensional finite-element S-parameter-simulation tool for microwave tubes , 2015, 2015 IEEE International Vacuum Electronics Conference (IVEC).

[29]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[30]  Julien Diaz,et al.  Multi-level explicit local time-stepping methods for second-order wave equations , 2015 .

[31]  Jian-Ming Jin,et al.  A comparative study of three finite element-based explicit numerical schemes for solving maxwell's equations , 2012, 2011 IEEE International Symposium on Antennas and Propagation (APSURSI).

[32]  M. Liou.  A novel method of evaluating transient response , 1966 .

[33]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[34]  Marcus J. Grote,et al.  High-order explicit local time-stepping methods for damped wave equations , 2011, J. Comput. Appl. Math..

[35]  Francesca Vipiana,et al.  EFIE Modeling of High-Definition Multiscale Structures , 2010, IEEE Transactions on Antennas and Propagation.

[36]  Geneviève Dujardin,et al.  High Order Exponential Integrators for Nonlinear Schrödinger Equations with Application to Rotating Bose-Einstein Condensates , 2015, SIAM J. Numer. Anal..

[37]  김덕영.  [신간안내] Computational Electrodynamics (the finite difference time - domain method) , 2001 .

[38]  Mike A. Botchev,et al.  Krylov subspace exponential time domain solution of Maxwell's equations in photonic crystal modeling , 2016, J. Comput. Appl. Math..

[39]  Stéphane Lanteri,et al.  Locally Implicit Time Integration Strategies in a Discontinuous Galerkin Method for Maxwell’s Equations , 2013, J. Sci. Comput..

[40]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[41]  R. Ward.  Numerical Computation of the Matrix Exponential with Accuracy Estimate , 1977 .

[42]  Stéphane Lanteri,et al.  Locally implicit discontinuous Galerkin method for time domain electromagnetics , 2010, J. Comput. Phys..

[43]  Jin-Fa Lee,et al.  Nonconformal Domain Decomposition Methods for Solving Large Multiscale Electromagnetic Scattering Problems , 2013, Proceedings of the IEEE.

[44]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[45]  Bin Li,et al.  A multifrontal block ILU preconditioner for the 3D finite-element eigenvalue analysis of lossy slow-wave structures of traveling-wave tubes , 2014, IEEE International Vacuum Electronics Conference.

[46]  M. Hochbruck,et al.  Exponential integrators , 2010, Acta Numerica.

[47]  Kurt Busch,et al.  Efficient multiple time-stepping algorithms of higher order , 2015, J. Comput. Phys..

[48]  Patrick R. Amestoy,et al.  Multifrontal parallel distributed symmetric and unsymmetric solvers , 2000 .