Research And Application Of Three-dimensional FDTD Parallel Algorithm

With the rapid development of the computer technology, the numerical methods of electromagnetic field have been developed rapidly. During recent years, Numerical Algorithms in The finite difference time domain (FDTD) have been widely used in a variety of practical problems. The wide-band frequency information can be gained one time through this method of calculation, and it is easy to handle complex shape and complex medium of electromagnetic problems.This dissertation lays s strong emphasis on the essential of a parallel algorithm for the three-dimensional Finite Difference Time Domain (FDTD) method based on MPI platform. The requirements for the memory and computation time in a electrically-large-size or complicated objects are dramatically decreased by using this method.During implementing the 3D parallel FDTD program, compared to traditional optimization algorithm (One-dimensional Parallel Algorithm of Three-dimensional FDTD Method), a 3D Cartesian topology is defined and the FDTD computation domain is divided into some subspaces using a spatial decomposition of the regular grid structure. Then the fields inside each sun-domain are computed on an individual processor with a few of data being communicated from neighborhood sub-domain. For avoided frequently packing and unpacking data, the communications are optimized by the use of user-defined data types, which group the data when their memory addresses are not continuous, then the cost of communication is reduced when the discontinuous data can be transmitted only once. The transmission is using block communication, shortening the delay and improving the efficiency of communication between sub-domains. Thus the optimized 3D parallel FDTD program can be result.Based on MPI development environment, the 3D FDTD parallel algorithm is implementation on a network of parallel systems with the PC. This dissertation discusses with emphasis on the parallel processing method to the absorbing boundary condition including the second-order Mur's absorbing condition, the uniaxial perfectly matched layer. The load balance is also analyzed.