Method for realizing perfectly matched layer by using current density convolution in plasma

Abstract

The invention discloses a method for realizing a perfectly matched layer by using current density convolution in a plasma. The method comprises the following steps of inputting a model file; initializing parameters and setting a PML (perfectly matched layer) coefficient and an absorbing boundary parameter; respectively updating and calculating the electric field component coefficients E<q>(y) and E<q>(x) of the whole calculation area in the y direction and the x direction, adding a field source to an electric field component coefficient, and updating and calculating the magnetic field component coefficient of the whole calculation area; updating and calculating the polarized current densities J<q>(x) and J<q>(y) of the whole calculation area, and updating and calculating the auxiliary variables of the electric and magnetic field component coefficients of the whole calculation area; updating and calculating electric and magnetic field components at an observing point; assigning q+1 assignment to q, judging whether step number q of Laguerre polynomials reaches a preset value or not, and if q does not reach the preset value, returning to step 3; if q reaches the preset value, ending. According to the method for realizing the perfectly matched layer by using the current density convolution in the plasma, the calculation speed is high, the memory consumption is low, and the method has very good absorbing effect on low frequencies and evanescent waves.

Description

technical field [0001] The invention belongs to the technical field of computational electromagnetics, and relates to a method for realizing a complete matching layer using current density convolution in plasma. Background technique [0002] The Finite-difference time-domain (FDTD) method is widely used in the simulation of electromagnetic wave propagation in dispersive media because of its simple implementation. However, its time step is limited by the Cauchy stability condition, which limits the application of the FDTD method to fine-structure models. In order to eliminate the limitation of the Cauchy stability condition, unconditionally stable finite-difference time-domain methods have been proposed, such as: Alternating-Direction-Implicit (ADI) finite-difference time-domain (ADI-FDTD) method and weighted-based Laguerre polynomials Finite-difference time-domain (WLP-FDTD) method. Among these methods, the ADI-FDTD method will produce a large dispersion error when using a...

Images