Recently an efficient pseudospectral time-domain (PSTD) algorithm has been developed to solve partial differential equations in computational electromagnetics and acoustics. It uses the fast Fourier transform (FFT) algorithm to approximate spatial derivatives, and the perfectly matched layer (PML) to eliminate the wraparound effect. Due to its high accuracy in the spatial derivatives, this method requires a significantly smaller number of unknowns than a conventional finite-difference timedomain (FDTD) method when solving large-scale problems. In this work, we further extend the PSTD algorithm to frequencydependent media and apply the algorithm to simulate groundpenetrating radar (GPR) measurements in a dispersive earth. The dispersion of the soil is treated by the recursive convolution approaches. The convergence property of the PSTD algorithm is investigated for the scattering of a dispersive cylinder. Multidimensional large-scale problems in GPR measurements are presented to demonstrate the efficiency of this frequency-dependent PSTD algorithm.Index Terms-Dispersive media, FDTD, ground-penetrating radar (GPR), numerical analysis, pseudospectral time-domain (PSTD), perfectly matched layer (PML).
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.