Computational electrodynamics in material media with constraint-preservation, multidimensional Riemann solvers and sub-cell resolution - Part I, second-order FVTD schemes

Abstract While classic finite-difference time-domain (FDTD) solutions of Maxwell's equations have served the computational electrodynamics (CED) community very well, formulations based on Godunov methodology have begun to show advantages. We argue that the formulations presented so far are such that FDTD schemes and Godunov-based schemes each have their own unique advantages. However, there is currently not a single formulation that systematically integrates the strengths of both these major strains of development. While an early glimpse of such a formulation was offered in Balsara et al. [16] , that paper focused on electrodynamics in plasma. Here, we present a synthesis that integrates the strengths of both FDTD and Godunov-based schemes into a robust single formulation for CED in material media. Three advances make this synthesis possible. First, from the FDTD method, we retain (but somewhat modify) a spatial staggering strategy for the primal variables. This provides a beneficial constraint preservation for the electric displacement and magnetic induction vector fields via reconstruction methods that were initially developed in some of the first author's papers for numerical magnetohydrodynamics (MHD). Second, from the Godunov method, we retain the idea of upwinding, except that this idea, too, has to be significantly modified to use the multi-dimensionally upwinded Riemann solvers developed by the first author. Third, we draw upon recent advances in arbitrary derivatives in space and time (ADER) time-stepping by the first author and his colleagues. We use the ADER predictor step to endow our method with sub-cell resolving capabilities so that the method can be stiffly stable and resolve significant sub-cell variation in the material properties within a zone. Overall, in this paper, we report a new scheme for numerically solving Maxwell's equations in material media, with special attention paid to a second-order-accurate formulation. Several numerical examples are presented to show that the proposed technique works. Because of its sub-cell resolving ability, the new method retains second-order accuracy even when material permeability and permittivity vary by an order-of-magnitude over just one or two zones. Furthermore, because the new method is also unconditionally stable in the presence of stiff source terms (i.e., in problems involving giant conductivity variations), it can handle several orders-of-magnitude variation in material conductivity over just one or two zones without any reduction of the time–step. Consequently, the CFL depends only on the propagation speed of light in the medium being studied.

[1]  J. Hawley,et al.  Simulation of magnetohydrodynamic flows: A Constrained transport method , 1988 .

[2]  Dinshaw S. Balsara.  Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics , 2009, J. Comput. Phys..

[3]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[4]  C. Richard DeVore,et al.  Flux-corrected transport techniques for multidimensional compressible magnetohydrodynamics , 1989 .

[5]  Dinshaw S. Balsara.  Multidimensional HLLE Riemann solver: Application to Euler and magnetohydrodynamic flows , 2010, J. Comput. Phys..

[6]  J.V. Moloney,et al.  Second-order accurate FDTD space and time grid refinement method in three space dimensions , 2006, IEEE Photonics Technology Letters.

[7]  Dinshaw S. Balsara,et al.  Multidimensional Riemann problem with self-similar internal structure - part III - a multidimensional analogue of the HLLI Riemann solver for conservative hyperbolic systems , 2017, J. Comput. Phys..

[8]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[9]  Dinshaw S. Balsara,et al.  A high-order relativistic two-fluid electrodynamic scheme with consistent reconstruction of electromagnetic fields and a multidimensional Riemann solver for electromagnetism , 2016, J. Comput. Phys..

[10]  Michael Dumbser,et al.  ADER-WENO finite volume schemes with space-time adaptive mesh refinement , 2012, J. Comput. Phys..

[11]  Eleuterio F. Toro,et al.  Towards Very High Order Godunov Schemes , 2001 .

[12]  Paul R. Woodward,et al.  On the Divergence-free Condition and Conservation Laws in Numerical Simulations for Supersonic Magnetohydrodynamical Flows , 1998 .

[13]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[14]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[15]  Pedro Velarde,et al.  Development of a Godunov method for Maxwell's equations with Adaptive Mesh Refinement , 2015, J. Comput. Phys..

[16]  D. Balsara,et al.  A Staggered Mesh Algorithm Using High Order Godunov Fluxes to Ensure Solenoidal Magnetic Fields in Magnetohydrodynamic Simulations , 1999 .

[17]  Francesco Miniati,et al.  A Divergence-free Upwind Code for Multidimensional Magnetohydrodynamic Flows , 1998 .

[18]  Dinshaw S. Balsara,et al.  Multidimensional Riemann problem with self-similar internal structure. Part I - Application to hyperbolic conservation laws on structured meshes , 2014, J. Comput. Phys..

[19]  Dinshaw S. Balsara,et al.  An efficient class of WENO schemes with adaptive order , 2016, J. Comput. Phys..

[20]  Michael Dumbser,et al.  Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems , 2007, J. Comput. Phys..

[21]  Steven G. Johnson,et al.  Advances in FDTD Computational Electrodynamics: Photonics and Nanotechnology , 2013 .

[22]  Steven G. Johnson,et al.  Improving accuracy by subpixel smoothing in the finite-difference time domain. , 2006, Optics letters.

[23]  Dinshaw S. Balsara,et al.  Three dimensional HLL Riemann solver for conservation laws on structured meshes; Application to Euler and magnetohydrodynamic flows , 2015, J. Comput. Phys..

[24]  Dinshaw S. Balsara.  A two-dimensional HLLC Riemann solver for conservation laws: Application to Euler and magnetohydrodynamic flows , 2012, J. Comput. Phys..

[25]  John Lyon,et al.  A simulation study of east-west IMF effects on the magnetosphere , 1981 .

[26]  Michael Dumbser,et al.  A two-dimensional Riemann solver with self-similar sub-structure - Alternative formulation based on least squares projection , 2016, J. Comput. Phys..

[27]  Michael Dumbser,et al.  Efficient implementation of ADER schemes for Euler and magnetohydrodynamical flows on structured meshes - Speed comparisons with Runge-Kutta methods , 2013, J. Comput. Phys..

[28]  J. Bérenger.  Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[29]  A. Taflove,et al.  Finite-difference time-domain model of lasing action in a four-level two-electron atomic system. , 2004, Optics express.

[30]  Timur Z. Ismagilov.  Second order finite volume scheme for Maxwell's equations with discontinuous electromagnetic properties on unstructured meshes , 2015, J. Comput. Phys..

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

[32]  Steven J. Ruuth,et al.  Non-linear evolution using optimal fourth-order strong-stability-preserving Runge-Kutta methods , 2003, Math. Comput. Simul..

[33]  Willem Hundsdorfer,et al.  IMEX extensions of linear multistep methods with general monotonicity and boundedness properties , 2007, J. Comput. Phys..

[34]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[35]  Michael Dumbser,et al.  Multidimensional Riemann problem with self-similar internal structure. Part II - Application to hyperbolic conservation laws on unstructured meshes , 2015, J. Comput. Phys..

[36]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[37]  Rémi Abgrall,et al.  Multidimensional HLLC Riemann solver for unstructured meshes - With application to Euler and MHD flows , 2014, J. Comput. Phys..

[38]  Allen Taflove,et al.  Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures , 1988 .

[39]  Allen Taflove,et al.  Finite-Difference Time-Domain Solution of Maxwell's Equations , 2016 .

[40]  G. Russo,et al.  Implicit-explicit runge-kutta schemes and applications to hyperbolic systems with relaxation , 2005 .

[41]  E. Toro,et al.  Solution of the generalized Riemann problem for advection–reaction equations , 2002, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[42]  Steven J. Ruuth,et al.  A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods , 2002, SIAM J. Numer. Anal..

[43]  Zhiliang Xu,et al.  Divergence-Free WENO Reconstruction-Based Finite Volume Scheme for Solving Ideal MHD Equations on Triangular Meshes , 2011, 1110.0860.

[44]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[45]  Eleuterio F. Toro,et al.  ADER: Arbitrary High Order Godunov Approach , 2002, J. Sci. Comput..

[46]  Zach DeVito,et al.  Opt , 2017 .

[47]  Steven G. Johnson,et al.  Accurate finite-difference time-domain simulation of anisotropic media by subpixel smoothing. , 2009, Optics letters.

[48]  Dinshaw Balsara,et al.  Second-Order-accurate Schemes for Magnetohydrodynamics with Divergence-free Reconstruction , 2003, astro-ph/0308249.

[49]  Chi-Wang Shu.  Total-variation-diminishing time discretizations , 1988 .

[50]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[51]  Michael Dumbser,et al.  Divergence-free MHD on unstructured meshes using high order finite volume schemes based on multidimensional Riemann solvers , 2015, J. Comput. Phys..

[52]  Eleuterio F. Toro,et al.  ADER schemes for three-dimensional non-linear hyperbolic systems , 2005 .

[53]  D. Katz,et al.  Validation and extension to three dimensions of the Berenger PML absorbing boundary condition for FD-TD meshes , 1994, IEEE Microwave and Guided Wave Letters.

[54]  Inmaculada Higueras,et al.  Total-variation-diminishing implicit-explicit Runge-Kutta methods for the simulation of double-diffusive convection in astrophysics , 2011, J. Comput. Phys..

[55]  Claus-Dieter Munz,et al.  Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model , 2000 .

[56]  Dinshaw Balsara,et al.  Divergence-free adaptive mesh refinement for Magnetohydrodynamics , 2001 .

[57]  Michael Dumbser,et al.  Efficient, high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics , 2008, Journal of Computational Physics.

[58]  Dinshaw S. Balsara,et al.  Von Neumann stability analysis of globally divergence-free RKDG schemes for the induction equation using multidimensional Riemann solvers , 2017, J. Comput. Phys..