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

by   Yubo Yang, et al.

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.



There are no comments yet.


page 1

page 2

page 3

page 4


Convergence analysis of a fully discrete energy-stable numerical scheme for the Q-tensor flow of liquid crystals

We present a fully discrete convergent finite difference scheme for the ...

Backward Euler method for the equations of motion arising in Oldroyd model of order one with nonsmooth initial data

In this paper, a backward Euler method combined with finite element disc...

Numerical stability of Grünwald-Letnikov method for time fractional delay differential equations

This paper is concerned with the numerical stability of time fractional ...

Travelling wave solutions for fully discrete FitzHugh-Nagumo type equations with infinite-range interactions

We investigate the impact of spatial-temporal discretisation schemes on ...

Discrete gradient structures of BDF methods up to fifth-order for the phase field crystal model

The well-known backward difference formulas (BDF) of the third, the four...

Stability properties of a projector-splitting scheme for the dynamical low rank approximation of random parabolic equations

We consider the Dynamical Low Rank (DLR) approximation of random parabol...

Towards Hybrid Two-Phase Modelling Using Linear Domain Decomposition

The viscous flow of two immiscible fluids in a porous medium on the Darc...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.