Numerical modeling of round-penetrating radar (GPR) using finite-element method

Q: How can I download an article?

To download an article from SID, first log in to the site, search for the article title, and click on the 'Download Article' option.

Q: How can I download an ISI article?

To download an ISI article on SID, enter the keyword or article title in the search bar, view the relevant results, click on the desired article, and select the 'Download Article' option.

Q: How can I access the SID database?

To access the SID database, visit SID.ir, create an account, and log in to access scientific resources.

Q: Is downloading articles from SID free?

Some articles on SID are available for free, while others require payment. Details are specified on the article's page.

Zarei Sajjad; Rahimi Dalkhani Amin; AMINI NAVID

مرکز اطلاعات علمی Scientific Information Database (SID) - Trusted Source for Research and Academic Resources

Journal Paper

Paper Information

Journal: IRANIAN JOURNAL OF GEOPHYSICS Year:2017 | Volume:11 | Issue:3 Page(s): 144-156

Download Full-Text

Persian Verion

View:

761

Download:

Cites:

Information Journal Paper

Title

Numerical modeling of round-penetrating radar (GPR) using finite-element method

Author(s)

Zarei Sajjad | Rahimi Dalkhani Amin | AMINI NAVID | Issue Writer Certificate

Keywords

ground-penetrating radar (GPR)Q1

finite-element method (FEM)Q1

numerical modelingQ1

absorbing boundary conditions (ABC)Q1

Abstract

Ground-penetrating radar (GPR) is a popular geophysical method for high-resolution imaging of the shallow subsurface structures. numerical modeling of radar waves plays a significant role in interpretation, processing, and imaging of GPR data. A number of different approaches have been presented for the numerical modeling of GPR data. The most common approach for GPR modeling is the finite-difference method (FDM) because the FDM approach is conceptually simple and easy to program. The difficulties in applying boundary conditions at non-linear boundaries and the lack of sufficient accuracy in complex geometries are the most important drawbacks of FDM. This paper presents a finite-element method, for simulation of ground-penetrating radar (GPR) in two dimensions in the time-domain. FEM is a powerful and versatile numerical technique for handling problems involving complex geometries and inhomogeneous media. The technique is based on a weak formulation of Maxwell’ s equations. In the FEM method, the wavefield is discretized on the elements using Lagrange interpolation, and integration over an element is accomplished based upon the Gaussian-quad integration rule. The major difference between the various numerical methods is in the spatial discretization. In the elementalbased methods, the complex geometry of the problem is divided into several smaller and simpler elements, then the integrals are calculated for each element. These methods have no with any regular or irregular geometry. The responses of the model in the finite-element methods are approximated in nodal points, so nodal polynomials of Lagrange are used for interpolation of the model response. Besides, the systematic generality of the method makes it possible to construct general-purpose computer programs for solving a wide range of problems. In this paper, at first, Maxwell’ s equations are discretized, then the boundary condition is applied to minimize artificial reflections from the edges of the computation domain. Although the governing equations and mechanisms are completely different between radar and seismic waves, most of GPR data processing approaches are derived from seismic data processing. Due to similarities in these two techniques, accordingly, we implement the first-order Clayton and Engquist absorbing boundary conditions (first order CE-ABC) introduced in the numerical finitedifference modeling of seismic wave propagation. This boundary condition is simple to apply. The presented formulations are in matrix notation that simplifies the implementation of the relations in computer programs, especially in MATLAB application. After spatial discretization with FEM, time discretization is done by Finite-Central Difference (FCD). The time discretization is the most massive and time-consuming step in modeling, which spatial discretization has an important role in this process. The stiffness, mass and damping matrices are sparse and symmetrical in FEM; so if we use the optimized numerical algorithms and storages strategies, computational costs and processing-time can be reduced significantly. To investigate the efficiency of FEM, the computer program has been written in MATLAB and has been tested on two models. The results show that the radar wave simulation via FEM is an accurate and effective approach in complex models.

Cites

No record.

References

No record.

Cite

APA: Copy

Zarei, Sajjad, Rahimi Dalkhani, Amin, & AMINI, NAVID. (2017). Numerical modeling of round-penetrating radar (GPR) using finite-element method. IRANIAN JOURNAL OF GEOPHYSICS, 11(3 ), 144-156. SID. https://sid.ir/paper/133717/en

Vancouver: Copy

Zarei Sajjad, Rahimi Dalkhani Amin, AMINI NAVID. Numerical modeling of round-penetrating radar (GPR) using finite-element method. IRANIAN JOURNAL OF GEOPHYSICS[Internet]. 2017;11(3 ):144-156. Available from: https://sid.ir/paper/133717/en

IEEE: Copy

Sajjad Zarei, Amin Rahimi Dalkhani, and NAVID AMINI, “Numerical modeling of round-penetrating radar (GPR) using finite-element method,” IRANIAN JOURNAL OF GEOPHYSICS, vol. 11, no. 3 , pp. 144–156, 2017, [Online]. Available: https://sid.ir/paper/133717/en

Related Journal Papers