Affiliations: [a] National Key Laboratory on Electromagnetic Environmental Effects and Electro-optical Engineering, PLA University of Science and Technology, Nanjing, Jiangsu, China | [b] School of Mathematics and Information Technology, Jiangsu Second Normal University, Nanjing, Jiangsu, China | [c] Jiangsu Regulatory Bureau of Nuclear and Radiation Safety, Nanjing, Jiangsu, China
Correspondence:
[*]
Corresponding author: Shu-Wen Chen, School of Mathematics and Information Technology, Jiangsu Second Normal University, Nanjing, Jiangsu, China. E-mail:[email protected]
Abstract:
This paper proposes the complex
frequency shifted (CFS) perfectly matched layer (PML) of dispersive media
for finite-difference time-domain (FDTD) method combined with weighted
Laguerre polynomials (WLP). According to the property of Laguerre basis
function, the relative dielectric constant ε _r (ω
) of dispersive media and complex frequency shifted (CFS) PML
parameters, auxiliary differential equation (ADE) technique is introduced.
Based on the ADE technique, the relationship between field components and
auxiliary differential variables is derived in Laguerre domain. Using ADE
scheme, the relationship between electric flux density D and electric field
E is derived in Laguerre domain. Substituting
auxiliary differential variables into CPML absorbing boundary conditions,
using auxiliary differential variables, electric flux density D of order q can be expressed
directly by magnetic field H in Laguerre
domain. Using the same procedure, magnetic field H of order q can be expressed
directly by electric field E in Laguerre
domain. Inserting H of order q$ into D , using the relationship between D and E , and using
central difference scheme, the formulations for dispersive media are
obtained. In order to validate the efficiency of the presented method, two
numerical examples are simulated. Numerical results show that, compared with
the Berenger PML (BPML) and nearly PML (NPML), the CFS-PML has about more
than 24 dB improvement in terms of the maximum relative error and much lower
reflection error for the late-time region.
