Abstract:In this paper, we consider Maxwell's equations in linear dispersive media described by a single-pole Lorentz model for electronic polarization. We study two classes of commonly used spatial discretizations: finite difference methods (FD) with arbitrary even order accuracy in space and high spatial order discontinuous Galerkin (DG) finite element methods. Both types of spatial discretizations are coupled with second order semi-implicit leap-frog and implicit trapezoidal temporal schemes studied in our previous research [5,6]. By performing detailed dispersion analysis for the semi-discrete and fully discrete schemes, we obtain rigorous quantification of the dispersion error for Lorentz dispersive dielectrics. In particular, comparisons of dispersion error can be made taking into account the model parameters, and mesh sizes in the design of the two types of schemes. The results for the numerical dispersion analysis can guide us in the optimal choice of discretization parameters for the more complicated and nonlinear models. The numerical dispersion analysis of the fully discrete FD and DG schemes, for the dispersive Maxwell model considered in this paper, clearly indicate the dependence of the numerical dispersion errors on spatial and temporal discretizations, their order of accuracy, mesh discretization parameters and model parameters. The results obtained here cannot be arrived at by considering discretizations of Maxwell's equations in free space. In particular, our results contrast the advantages and disadvantages of using high order FD or DG schemes and leap-frog or trapezoidal time integrators over different frequency ranges using a variety of measures of numerical dispersion errors. Finally, we highlight the limitations of the second order accurate temporal discretizations considered.

Two-Grid Methods for Maxwell’s Equations in a Cole-Cole Dispersive Medium

Discontinuous Galerkin Time Domain Electric Field Integral Equation Method For Solving Non-Conformal Structures

Two-grid finite element method with an H2N2 interpolation for two-dimensional nonlinear fractional multi-term mixed sub-diffusion and diffusion wave equation

Two‐grid Methods for Time‐harmonic Maxwell Equations

An Iterative Two-Grid Method of A Finite Element PML Approximation for the Two Dimensional Maxwell Problem

DEH Scheme DGTD-Based Transient Modeling Approach for the Cole–Cole Dispersive Media Using Tustin’s Method

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

A Weak Form Electromagnetic Transient Simulation Approach for Cole-Cole Medium

Discontinuous Galerkin Time-Domain Method for GPR Simulation in Dispersive Media

Superconvergence of a two‐grid method for fourth‐order dispersive‐dissipative wave equations

Dispersion Analysis of Finite Difference and Discontinuous Galerkin Schemes for Maxwell's Equations in Linear Lorentz Media

A Discontinuous Galerkin Time Domain Method for Solving Drude Dispersive Models

A new efficient 3D Discontinuous Galerkin Time Domain (DGTD) method for large and multiscale electromagnetic simulations

A postprocessing technique for a discontinuous Galerkin discretization of time-dependent Maxwell's equations

Recent Developments to the Discontinuous Galerkin Time-Domain Method for 3-D Multiscale Problems

Numerical Studies of Time Dependent Ginzburg-Landau Model by FEM with Moving Grid Deformation

Perfectly Matched Layer for Cole-Cole Dispersive Media in DGTD Method

Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions

An ADER Discontinuous Galerkin Method with Local Time-Stepping for Transient Electromagnetics.

Non-dissipative and spurious-free discontinuous Galerkin method for solving Maxwell equations:(Ⅰ) One-and two-dimensional cases

A two-grid discretization method for nonlinear Schrödinger equation by mixed finite element methods