Analysis of a Backward Euler-type Scheme for Maxwell's Equations in a Havriliak-Negami Dispersive Medium

07/21/2020
by   Yubo Yang, et al.
0

For the Maxwell's equations in a Havriliak-Negami (H-N) dispersive medium, the associated energy dissipation law has not been settled at both continuous level and discrete level. In this paper, we rigorously show that the energy of the H-N model can be bounded by the initial energy so that the model is well-posed. We analyse a backward Euler-type semi-discrete scheme, and prove that the modified discrete energy decays monotonically in time. Such a strong stability ensures that the scheme is unconditionally stable. We also introduce a fast temporal convolution algorithm to alleviate the burden of the history dependence in the polarisation relation involving the singular kernel with the Mittag-Leffler function with three parameters. We provide ample numerical results to demonstrate the efficiency and accuracy of a full-discrete scheme via a spectra-Galerkin method in two spatial dimension. Finally, we consider an interesting application in the recovery of complex relative permittivity and some related physical quantities.

READ FULL TEXT

Please sign up or login with your details

Forgot password? Click here to reset