Second-generation wavelet finite element based on the lifting scheme for GPR simulation
- PDF / 1,697,818 Bytes
- 11 Pages / 595.22 x 842 pts (A4) Page_size
- 55 Downloads / 155 Views
Second-generation wavelet finite element based on the lifting scheme for GPR simulation* Feng De-Shan1,2, Zhang Hua1,2, and Wang Xun1,2 Abstract: Ground-penetrating radar (GPR) is a highly efficient, fast and non-destructive exploration method for shallow surfaces. High-precision numerical simulation method is employed to improve the interpretation precision of detection. Second-generation wavelet finite element is introduced into the forward modeling of the GPR. As the finite element basis function, the second-generation wavelet scaling function constructed by the scheme is characterized as having multiple scales and resolutions. The function can change the analytical scale arbitrarily according to actual needs. We can adopt a small analysis scale at a large gradient to improve the precision of analysis while adopting a large analytical scale at a small gradient to improve the efficiency of analysis. This approach is beneficial to capture the local mutation characteristics of the solution and improve the resolution without changing mesh subdivision to realize the efficient solution of the forward GPR problem. The algorithm is applied to the numerical simulation of line current radiation source and tunnel non-dense lining model with analytical solutions. Result show that the solution results of the secondgeneration wavelet finite element are in agreement with the analytical solutions and the conventional finite element solutions, thereby verifying the accuracy of the second-generation wavelet finite element algorithm. Furthermore, the second-generation wavelet finite element algorithm can change the analysis scale arbitrarily according to the actual problem without subdividing grids again. The adaptive algorithm is superior to traditional scheme in grid refinement and basis function order increase, which makes this algorithm suitable for solving complex GPR forward-modeling problems with large gradient and singularity. Keywords: Ground penetrating radar, wave equation, second-generation wavelet finite element method, lifting scheme, forward modeling
Introduction The objective of Ground Penetrating Radar (GPR) forward modeling is to solve Maxwell equations by using numerical algorithms. The current numerical simulation
algorithms used in solving the Maxwell equation are mainly as follows: finite difference method (FDM) and finite element method (FEM). FDM is simple in principle and easy to implement (Giannopoulos, 2005; Irving and Knight, 2006; Feng and Dai, 2011; Li et al., 2012), but it is unable to combine with unstructured grids.
Manuscript received by the Editor October 1,2018; revised manuscript received February 27,2019 * This work was supported by the National Natural Science Foundation of China (Nos. 41574116 and 41774132) and Hunan Provincial Innovation Foundation for Postgraduate (Grant Nos. CX2017B052) and the Fundamental Research Funds for the Central Universities of Central South University (Nos. 2018zzts693) 1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China. 2. Key Labo
Data Loading...
