Improved Shift-Operator FDTD Method for Anisotropic Magnetized Cold Plasmas With Arbitrary Magnetic Field Declination

Abstract

Abstract-In this paper, a recently improved SO-FDTD (shiftoperator finite difference time-domain) method is proposed and applied to the numerical analysis of the anisotropic magnetized plasma with arbitrary magnetic declination. By using the constitutive relation between polarized current density vector J and electric vector E and bringing the shift operators, the difference iteration equations of field components for Maxwell equations are derived in detail. Furthermore, the memory requirement is decreased significantly through incorporating a memory-minimized algorithm into the FDTD iterative cycles.The reflection and transmission coefficients of electromagnetic wave through a magnetized plasma layer are calculated by using this method. It is shown that the new method not only improves accuracy but also produces speed and memory advantages over the SO-FDTD method in kDB coordinates system proposed in the recent reference. In addition, the recursion formulae of the improved SO-FDTD method are deduced and programmed easily and they involve no complex variables, so the computations for the magnetized plasma become very simple.

Analysis of Propagation and Polarization Characteristics of Electromagnetic Waves Through Nonuniform Magnetized Plasma Slab Using Propagator Matrix Method

“…It is known that there are two kinds of eigen wave in an anisotropic media such as magnetized plasma [20], one is called type I wave and the other is type II wave. These two types of waves are coupled to each other at the interface -in other words, a type I wave may produce transmitted and reflected waves of both type I and type II.…”

“…In order to check the validity of the proposed propagator matrix method (PMM), a test is done by considering a uniform magnetized plasma plate which is similar to the case studied by Yin et al [20]. Here, the plate parameters are selected as: ω p = 2π × 28.7 × 10 9 rad/s, ω b = 8.8 × 10 10 rad/s, ν en = 2 × 10 10 rad/s and the plate thickness is 9 mm.…”

“…One is ordinary wave (O-wave), and the other is extraordinary wave (Xwave) [30]. Subsequently, three methods, i.e., the PMM proposed in this paper, the analytical method and the SO-FDTD method introduced in [20], are used to calculate EM wave propagating through this plate. When PMM is used, note that O-wave and X-wave are linearly polarized waves and not coupled with each other, the reflection and transmission coefficients of X-wave (O-wave) can be obtained by R 11 (R 22 ) and T 11 (T 22 ), respectively.…”

“…It is very important to find out the mechanism and phenomena about the interactions between EM wave and plasma in the study of this technology. In the past twenty years, a number of analytical and numerical methods have been proposed to investigate the wave propagation in plasma [2,3,[12][13][14][15][16][17][18][19][20][21]. Hu et al [15] used a scattering matrix method to calculate the reflection, transmission and absorption characteristic of EM wave propagation in nonuniform magnetized plasma, where the external magnetized field direction is perpendicular to the direction of EM wave propagation.…”

“…Liu et al [13] and Lei Shi et al [14] carried out the studies on the transmission solution of EM waves propagation through plasma sheath by adopting the finite difference time-domain (FDTD) method. Yin et al [20,21] have utilized the shift-operator finite difference timedomain (SO-FDTD) method to study the reflection and transmission characteristics for a plane wave incident on a magnetized plasma slab. The FDTD method is a more powerful numerical tool for dealing with complicated plasma model compared to analytical methods mentioned above, but it is a highly time consuming numerical technique due to the usage of massive storage resources required by the problem.…”