Finite-difference time-domain simulation of ground penetrating radar on dispersive, inhomogeneous, and conductive soils

TLDR

The study presents a 3D time‑domain numerical scheme for simulating ground‑penetrating radar in dispersive, inhomogeneous, conductive soils. The method uses a finite‑difference time‑domain (FDTD) discretization with piecewise‑linear recursive convolution (PLRC) to model Lorentz/Debye soil dispersion, fitting parameters to experimental data, and extends a perfectly matched layer (PML) for dispersive media as an absorbing boundary. The formulation is validated with example simulations, shows parallel scalability, and achieves near‑linear speedup on a 32‑processor system.

Abstract

A three-dimensional (3D) time-domain numerical scheme for simulation of ground penetrating radar (GPR) on dispersive and inhomogeneous soils with conductive loss is described. The finite-difference time-domain (FDTD) method is used to discretize the partial differential equations for time stepping of the electromagnetic fields. The soil dispersion is modeled by multiterm Lorentz and/or Debye models and incorporated into the FDTD scheme by using the piecewise-linear recursive convolution (PLRC) technique. The dispersive soil parameters are obtained by fitting the model to reported experimental data. The perfectly matched layer (PML) is extended to match dispersive media and used as an absorbing boundary condition to simulate an open space. Examples are given to verify the numerical solution and demonstrate its applications. The 3D PML-PLRC-FDTD formulation facilitates the parallelization of the code. A version of the code is written for a 32-processor system, and an almost linear speedup is observed.