A two-dimensional high-precision iterative implementation method in non-magnetized plasma

Abstract

The invention discloses a two-dimensional high-precision iterative non-magnetized plasma implementation method, comprising: inputting model files; initializing parameters and setting parameters; adding field sources to the electric field component coefficients in the y direction, and setting the electric field component coefficients It is recorded as the initial field value update calculation of the electric field component coefficient in the y direction of the entire calculation area update calculation of the electric field component coefficient in the x direction of the entire calculation area to determine whether the number of iterations k reaches the preset value; update calculation of the magnetic field component coefficient update calculation of the entire calculation area The polarization current density coefficient of the entire calculation area is updated to calculate the auxiliary variable of the electromagnetic field component coefficient of the entire calculation area; the electromagnetic field component at the observation point is updated to be calculated; and it is judged whether the order q of the Laguerre polynomial reaches the preset value. The invention is a two-dimensional high-precision iterative implementation method in non-magnetized plasma, which has high calculation precision and fast calculation speed, and has good absorption effect on low frequency and litter waves.

Description

technical field [0001] The invention belongs to the technical field of computational electromagnetics, and in particular relates to a two-dimensional high-precision iterative implementation method in non-magnetized plasma. Background technique [0002] The Finite-difference time-domain (FDTD) method is widely used in the simulation of electromagnetic wave propagation in dispersive media due to its advantages of simple calculation and easy implementation. However, its time step is limited by the Cauchy stability condition and cannot be selected larger. In multi-scale complex fine structure models, the FDTD method has slow calculation speed and low calculation efficiency. 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...

AI Technical Summary

Problems solved by technology

Someone proposed a perfectly matched layer (Perfectly matched layer, PML) absorption boundary. Later, PML was widely used to calculate the truncation of the area, and it was proved to be very effective. However, the research found that the traditional PML has an absorption effect on low frequencies and litter waves. Not ideal; using a PML (CFS-PML) absorption boundary with a complex frequency shift (CFS) factor can effectively improve the traditional PML's absorption effect for low frequencies, litter waves and glancing situations

Recently, someone proposed a WLP-FDTD method using an auxiliary differential equation that approximates a perfectly matched absorption boundary to solve the electromagnetic field problem in a dispersive medium. The absorption effect of this approximately perfectly matched absorption boundary is very poor, and there are errors in calculation. Moreover, this algorithm takes a long time to calculate and consumes a lot of memory.

Images