Pemodelan Awal Ground Penetrating Radar dengan Metode Discontinuous Galerkin dan PML Berenger

Pranowo Pranowo


This paper discusses the development of Discontinuous Galerkin method, which has both linear shape function and weight function, for modeling Ground Penetrating Radar (GPR) in heterogeneous media. The triangular meshes are used due to their flexibility to deal with complex geometries. The Berenger Perfectly Matched Layer (PML) is used as absorbing boundary condition at the truncation boundaries. The numerical results of the DG method are compared with the exact solutions and the numerical results of FDTD method and the comparisons show that DG method has better accuracy than FDTD method and more stable for long time simulation. The simulation results of GPR show that the PML works well. Propagating waves at the edge of absorbing boundaries can be suppressed without any significant reflection. The results also show that various waves e.g., transmission waves, reflection waves, and diffraction waves produced by heterogeneous material can be simulated well.

Full Text:



D. Daniels, Ground Penetrating Radar, 2nd Edition, The Institution of Electrical Engineers, London, 2004.

K. S. Yee, "Numerical Solution of Initial Boundary Value Problems Involving Maxwell's Equations in Isotropic Media", IEEE Transactions on Antennas and Propagation, May 1966 p 302-307.

J. M. Bourgeois, "A Fully Three-Dimensional Simulation of a Ground-Penetrating Radar: FDTD Theory Compared with Experiment", IEEE Transactions on Geoscience and Remote Sensing, Vol 34 no. 1 Jan 1996.

J. Irving and R. Knight, " Numerical modeling of ground-penetrating radar in 2-D using MATLAB", Computers & Geosciences 32 , p. 1247–1258, 2006.

N. Diamanti and A. Giannopoulos, "Implementation of ADI-FDTD subgrids in ground penetrating radar FDTD models", Journal of Applied Geophysics 67 , p. 309–317, 2009.

F.D. Shan and D. Q. Wei, "GPR numerical simulation of full wave field based on UPML boundary condition of ADI-FDTD", NDT&E International 44, p. 495–504, 2011.

Y. Shi and C.H. Liang, 2007, The finite-volume time-domain algorithm using least square method in solving Maxwell’s equations, Journal of Computational Physics 226 (2007) p. 1444–1457.

K. Sankaran et al., "Cell-Centered Finite-Volume-Based Perfectly Matched Layer for Time-Domain Maxwell System", IEEE Transactions on Microwave Theory and Techniques, Volume 54 issue 3 (2006) p.1269- 1276.

D. Baumann, et al., 2004, "Finite-volume time-domain (FVTD) modelling of a broadband double-ridged horn antenna", Int. J. Numer. Model, 17 (2004) , p.285–298.

D. Dosopoulos et al., "An MPI/GPU parallelization of an interior penalty discontinuous Galerkin time domain method for Maxwell’s equations", Radio Science, Vol. 46, RS0M05, 2011.

M. Min and P. Fischer, "An Efficient High-Order Time Integration Method for Spectral-Element Discontinuous Galerkin Simulations in Electromagnetics", J Sci Comput, p. 57:582–603,2013.

T. Lu et al., "Discontinuous Galerkin methods for dispersive and lossy Maxwell’s equations and PML boundary conditions", Journal of Computational Physics, 200 (2004) , p. 549–580.

T. Lu et al., "Discontinuous Galerkin Time-Domain Method for GPR Simulation in Dispersive Media", IEEE Transactions on Geoscience and Remote Sensing, Vol. 43, No. 1, January 2005.

J.P. Berenger, "A Perfectly Matched Layer for the Absorption of Electromagnetic Waves", Journal of Computational Physics 114, 1994.

R. W. Lewis, , P. Nithiarasu and K. N. Seetharamu, Fundamentals of the Finite Element Method for Heat Transfer and Fluid Flow, John Wiley & Sons, Chichester



  • There are currently no refbacks.

Copyright (c) 2016 Jurnal Nasional Teknik Elektro dan Teknologi Informasi (JNTETI)

JNTETI (Jurnal Nasional Teknik Elektro dan Teknologi Informasi)

Departemen Teknik Elektro dan Teknologi Informasi, Fakultas Teknik Universitas Gadjah Mada
Jl. Grafika No 2. Kampus UGM Yogyakarta 55281
+62 274 552305