摘要
ADI-FDTD could break through the restriction of the CFL stability condition and is unconditionally stable. Through the discretization of the three two-dimensional Maxwell equations, this paper put forward two sub-time step iterative finite difference formulae of ADI- FDTD and UPML boundary condition of GPR wave, and gave the computing steps in detail separately. Based on this, the authors compile the corresponding program, use the program to simulate the two GPR models, get the wiggle maps, scan maps, and snapshots of the two models. The result indicated that the ADI- FDTD algorithm based on the UPML boundary condition eliminated the strong reflection of truncating boundaries and it could choose larger time steps. So it could simulate GPR model efficiently, reliably and accurately.