Symplectic partitioned Runge–Kutta methods for two-dimensional numerical model of ground penetrating radar

Simulation of ground penetrating radar (GPR) wave propagation in two-dimensional (2-D) subsurface structure is developed using the symplectic partitioned Runge–Kutta (SPRK) method. A transmitting boundary is implemented to absorb waves at the edges of the modeling. For the 2-D case, the SPRK schemes require only two functions for the complete description of the electromagnetic field. To verify the performance of the proposed algorithms, results from comparisons between the SPRK schemes and the standard FDTD method are presented. In addition, a complicated subsurface structure model is considered. The wiggle trace profile of this model is obtained from forward simulation using a 2-order SPRK method.

► We use SPRK methods for simulation of GPR propagation in 2-D subsurface structure. ► The SPRK methods can preserve the global symplectic structure of the phase space. ► The SPRK methods are in a good agreement with the FDTD scheme. ► The SPRK schemes can save computer memory usage and CPU time significantly. ► We obtain the wiggle trace profile of the complicated subsurface structure model.