Novel Scheme for GPU-Accelerated Finite-Difference Time-Domain Simulation of Electromagnetic Wave Interaction With Magnetic Plasma

Based on graphical processing unit acceleration, a new method of finite-difference time-domain scheme is proposed to simulate the interaction between electromagnetic waves and magnetized plasma in two-dimensional conditions. In this study, transversely electric and transversely magnetic are computed in time to avoid matrix operations involving Lorentz equations of motion. Compared to Young's method, the new method reduces addition and multiplication by about 63% and 66%, respectively. The simulation results of ionospheric wave propagation show that the new method agrees well with Young's method and the calculation speed is improved significantly. Finite-difference time-domain (FDTD) has been widely used in electromagnetic waves propagation in plasma. FDTD method computation can be greatly accelerated by graphical processing unit (GPU). However, in magnetized plasma, the media information is described by a tensor, and the computing includes the complex matrix operation. Moreover, in nonlinear physical process, the complex operation brings a larger amount of computation. This shortage is more significant for multicore GPU with weaker single-core computing power. Thus, it is essential to find new FDTD schemes to avoid the matrix operation and it is meaningful for both scientific calculation and engineering application. We present a new two-dimension finite-difference time-domain (scheme for wave propagation in and interaction with magnetized plasma The new method avoid complex matrix operation in solving couple of three velocities induced by Lorentz term The method performance shows less operation and faster speed in computation, especially for multicore parallel computation with graphical processing unit