A novel CN–ICCG–FDTD algorithm research of plasma reflection and transmission characteristics of electromagnetic wave

Abstract The advantage of the ICCG method for solving large sparse matrix is taken in the CN–FDTD equation solving. The CN–ICCG–FDTD can accelerate iteration for numerical calculation, and reduce memory overhead. Maxwell's equation of the electromagnetic wave in dispersive medium plasma is deduced from the CN format after being differentiated in the time domain, the FDTD method become unconditionally stable. In this process, a large sparse matrix is generated. As for solving such matrix, ICCG method has sufficient advantages. In ICCG method, convergence is so fast and stable that it is quite easy for computer programming. In addition, data required by ICCG method are memorized in a computer in the format of the one-dimensional compression, which helps the computer memory save a large amount of capacity, especially when the issues are rather complex to cope with, this algorithm will display such strength. With these advantages, ICCG method is able to calculate the reflection coefficient and the transmission coefficient of electromagnetic waves in plasma as well as the phase angle of them. The calculation indicates that the time step required by CN–ICCG–FDTD method has eliminated the constraints of the CFL, thus shortening the required time, making the calculated result stable and accurate, and improving the efficiency of programming.

[1]  X Wu.  THE CALCULATION OF THREE-DIMENSIONAL GEOELECTRIC FIELD OF POINT SOURCE BY INCOMPLETE CHOLESKY CONJUGATE GRADIENT METHOD , 1998 .

[2]  Rushan Chen,et al.  An improved formalism for the Crank‐Nicolson finite‐difference time‐domain by conformal mapping method , 2009 .

[3]  Edward K. N. Yung,et al.  The unconditionally stable Crank Nicolson FDTD method for three-dimensional Maxwell's equations , 2006 .

[4]  Hong-Wei Yang,et al.  Analysis of the temperature influence on the plasma reflectance by a new FDTD method , 2012 .

[5]  A. Ishimaru.  Electromagnetic Wave Propagation, Radiation, and Scattering , 1990 .

[6]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[7]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[8]  Yan Liu,et al.  Shift operator method and Runge-Kutta exponential time differencing method for plasma , 2011 .

[9]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[10]  Guilin Sun,et al.  Approximate Crank-Nicolson schemes for the 2-D finite-difference time-domain method for TE/sub z/ waves , 2004, IEEE Transactions on Antennas and Propagation.

[11]  Michael A. Jensen,et al.  A massively parallel computation strategy for FDTD: time and space parallelism applied to electromagnetics problems , 1995 .

[12]  D. Kershaw.  The incomplete Cholesky—conjugate gradient method for the iterative solution of systems of linear equations , 1978 .

[13]  Yuan Nai-chang,et al.  An auxiliary differential equation FDTD method for anisotropic magnetized plasmas , 2005 .

[14]  Hong-Wei Yang.  Research on reflection and transmission characteristics of electromagnetic wave to corn leaf , 2012 .

[15]  C. W. Trueman,et al.  Unconditionally stable Crank-Nicolson scheme for solving two-dimensional Maxwell's equations , 2003 .

[16]  Dennis M. Sullivan,et al.  Frequency-dependent FDTD methods using Z transforms , 1992 .

[17]  Chen Rushan,et al.  SO-FDTD method and its application to the calculation of electromagnetic wave reflection coefficients of plasma , 2006 .

[18]  Han Liu,et al.  Analysis on the Calculation of Plasma Medium with Parallel SO‐FDTD Method , 2009 .