A Numerical Study of Lossy Multipole Debye Dispersive Media Using a Recursive Integral-FDTD Method

摘要

A dispersive finite-difference time-domain (FDTD) method based on the recursive integration (RI) technique for the modeling of the lossy multipole Debye dispersive media is described in this article. The interaction between the electromagnetic (EM) field and the human tissues is simulated by means of the RI method, and the frequency-dependent formulations which possess good compatibility with the main FDTD algorithm are achieved. Next, the expression of the multipole Debye dispersion model is similar to the multipole complex frequency shift (CFS) stretching function which is utilized to build the multipole perfectly matched layer (MPML). Therefore, a significant feature of using the RI method is its overall modeling of both multipole Debye dispersive media and MPML boundary conditions. Furthermore, the stability analysis of the RI method in solving the lossy multipole Debye dispersive model indicates that the Courant–Friedrich–Levy (CFL) stability limit of the regular FDTD can be maintained. At last, the numerical cases performed in this work demonstrate the correctness of the RI method in modeling dispersive media.