Penelitian ini mengembangkan pemrograman paralel berbasis GPU CUDA untuk
solusi numerik Persamaan Maxwell TM 2-dimensi. Metode Numerik yang
digunakan adalah metode Discontinuous Galerkin (DG) dengan fungsi basis
linier, program hasil pengembangan digunakan untuk pemodelan awal Ground
Penetrating Radar (GPR). Transmitter (pemancar) radar dimodelkan sebagai
sumber (source) yang berupa fungsi waktu Blackman-Harris dengan frekuensi
pusat 500 MHz, sedangkan material berupa material heterogen yang terdiri atas
lapisan udara, tanah dan batu. Mesh yang digunakan adalah mesh tak berstruktur
(unstructured) dengan elemen berbentuk segitiga, sehingga bentuk domain ruang
dengan material yang heterogen dapat dimodelkan dengan relatif mudah. Lapisan
penyerap yang disebut Perfectly Matched Layer (PML) dipasang pada sisi luar
domain ruang, lapisan ini menggunakan pendekatan PML Berenger (1994).
Perbandingan hasil perhitungan metode DG dengan jawaban eksak dan metode
FDTD menunjukkan bahwa metode DG mempunyai penumpukan galat (error)
yang lebih sedikit dibanding metode FDTD, sehingga metode DG lebih stabil
untuk simulasi dengan durasi yang panjang. Selain itu simulasi GPR
menunjukkan bahwa metode yang dikembangkan mampu men-simulasikan GPR
dengan baik, PML mampu menyerap gelombang yang keluar dengan baik.
Fenomena perambatan gelombang berupa transmisi, refleksi dan difraksi juga
dapat disimulasikan dengan baik sesuai teori. Selain itu, penggunaan
pemrograman paralel berbasis GPU CUDA mampu mempercepat proses
perhitungan 7 kali lebih cepat pemrograman serial CPU.