Application of Symplectic Partitioned Runge-Kutta Method and Total-Field Scattered-Field Formulation to Simulation of GPR Data

Abstract

The second-order Lobatto IIIA-IIIB symplectic partitioned RungeKutta (SPRK) method, combining with the first-order Mur absorbing boundary condition, is developed for the simulation of ground penetrating radar wave propagation in layered pavement structure. For 2-dimetional case, a significant advantage of this method is that only two functions need to be calculated at each time step. The total-field/scattered-field technique is used for plane wave excitation. Numerical examples are presented to verify the accuracy and efficiency of the proposed algorithm. The results illustrate that the reflected signal calculated by the SPRK method is in good agreement with that obtained using the finite difference time domain (FDTD) scheme, but the CPU time consumed by proposed algorithm is reduce about 20% of the FDTD scheme. In addition, an actual field test is conducted to evaluate the further performance of the SPRK method. It is found that the simulated waveform fits well with the measured signal in many aspects, especially in the peak amplitude and time delay.