Analysis of GPR Wave Propagation Using CUDA‐Implemented Conformal Symplectic Partitioned Runge‐Kutta Method

Abstract

Accurate forward modeling is of great significance for improving the accuracy and speed of inversion. For forward modeling of large sizes and fine structures, numerical accuracy and computational efficiency are not high, due to the stability conditions and the dense grid number. In this paper, the symplectic partitioned Runge-Kutta (SPRK) method, surface conformal technique, and graphics processor unit (GPU) acceleration technique are combined to establish a precise and efficient numerical model of electromagnetic wave propagation in complex geoelectric structures, with the goal of realizing a refined and efficient calculation of the electromagnetic response of an arbitrarily shaped underground target. The results show that the accuracy and efficiency of ground-penetrating radar (GPR) forward modeling are greatly improved when using our algorithm. This provides a theoretical basis for accurately interpreting GPR detection data and accurate and efficient forward modeling for the next step of inversion imaging.

Thermal Transport in Graphene Oxide Films: Theoretical Analysis and Molecular Dynamics Simulation

“…The annealing process [30] simulation is from 300 K up to 800 K. Then the system temperature down to 300 K is carried out under the NPT ensemble (constant number of atoms, constant pressure of 1.0 bar) to obtain a stable structure. The atomic coordinates of the pre-formed GO sheet and the periodic dimensions of the simulation cell are optimized using a gradient-based minimization method [43,44] implemented until pressure scaling at zero bar under periodic boundary to relax the structure and reduce the residue stress. The final overlap model of GO films with hydroxyl groups and vacancy defects is shown in Figure 3.…”

Vibration Analysis of Fluid Conveying Carbon Nanotubes Based on Nonlocal Timoshenko Beam Theory by Spectral Element Method

“…The solution of the ordinary differential Equations (9) and (10) with constant coefficients holds the following exponential form:η (ξ) = Ce ikξ ,φ(ξ) = De ikξ (11) where k is the dimensionless wave number, and C, D are undetermined coefficients. Substituting Equation (11) into Equations (9) and (10), the following matrix equation can be obtained:…”

“…Since the experiments at the nanoscale level are difficult to conduct [8] and have high computational cost, the molecular dynamics simulation method is more often used to deal with small-size atomic systems. Hence, the continuum-based solid mechanics models or methods [8][9][10][11] bring vitality for predicting the mechanical properties of a large CNT. As an important mechanical property, the vibration characteristics of a CNT have become another hot topic for research in recent years.…”

Forward Analysis of GPR for Underground Pipes Using CUDA-Implemented Conformal Symplectic Euler Algorithm

“…In order to accurately interpret the measured GPR images of underground pipes, it is necessary to study the GPR data numerical simulation of underground pipes, and evaluate what useful information can be reasonably extracted from the field data under various conditions [17]- [21]. In the previous study, we established an accurate and efficient forward model based on the symplectic Euler algorithm, surface conformal technology and GPU acceleration technology to simulate the circular cavity disease in pavement structure [22], [23]. However, the CPML boundary was not used in the previous simulation, and the absorption at the boundary was poor.…”