2-D electromagnetic modelling by finite element method on GPU

Abstract 2-D electromagnetic model with a line source uses Preconditioned Conjugate Gradient to solve the stiffness matrix and Herrmann's pseudo-delta function to simulate the neighborhood of the singular source. We employ this Finite Element Method (FEM) on Fermi-GPU, which is the third generation GPU for the general computing purpose. Our algorithm takes advantage of the GPU's particular characteristics of Single Instruction Multiple Thread (SIMT) architecture and runs entirely on the GPU to minimize the transfer time between CPU and GPU. Meanwhile, we design a reasonable pre-conditioner (Diagonal Pre-conditioner) on the GPU and inject a mixed scheme into our algorithm for obtaining the good performance. In the numerical experiments, the correctness of the algorithm is acknowledged, and it achieves high performance by outperforming CPU solvers over hundreds of times. Indeed, it allows finer mesh generation so as to obtain more accurate resolutions and advance the ability of data interpretation of geophysical prospecting.

[1]  Yuta Inoue,et al.  GPGPU-FDTD method for 2-dimensional electromagnetic field simulation and its estimation , 2009, 2009 IEEE 18th Conference on Electrical Performance of Electronic Packaging and Systems.

[2]  P. Wannamaker,et al.  Three-dimensional magnetotelluric modeling using difference equations­ Theory and comparisons to integral equation solutions , 1993 .

[3]  Markus Clemens,et al.  GPU Accelerated Adams–Bashforth Multirate Discontinuous Galerkin FEM Simulation of High-Frequency Electromagnetic Fields , 2010, IEEE Transactions on Magnetics.

[4]  Di Qing-yun.  2.5 D CSAMT modeling with finite element method , 2004 .

[5]  Efficient 3-D controlled-source electromagnetic modelling using an exponential finite-difference method , 2015 .

[6]  Xingming Sun,et al.  Achieving Efficient Cloud Search Services: Multi-Keyword Ranked Search over Encrypted Cloud Data Supporting Parallel Computing , 2015, IEICE Trans. Commun..

[7]  E. Haber,et al.  Fast Simulation of 3D Electromagnetic Problems Using Potentials , 2000 .

[8]  Tsili Wang,et al.  3‐D electromagnetic anisotropy modeling using finite differences , 2001 .

[9]  Hiroshi Takenaka,et al.  2.5-D modelling of elastic waves using the pseudospectral method , 1996 .

[10]  Robert Strzodka,et al.  Cyclic Reduction Tridiagonal Solvers on GPUs Applied to Mixed-Precision Multigrid , 2011, IEEE Transactions on Parallel and Distributed Systems.

[11]  Analysis of anomalous electrical conductivity and magnetic permeability effects using a frequency domain controlled-source electromagnetic method , 2016 .

[13]  Oliver Ritter,et al.  CSEM Monitoring of Reservoir Oil-saturation Using a Novel Borehole-to-surface Configuration , 2015 .

[14]  L. Knizhnerman,et al.  Spectral approach to solving three-dimensional Maxwell's diffusion equations in the time and frequency domains , 1994 .

[15]  Jack J. Dongarra,et al.  Using Mixed Precision for Sparse Matrix Computations to Enhance the Performance while Achieving 64-bit Accuracy , 2008, TOMS.

[17]  Robert B. Herrmann,et al.  SH-wave generation by dislocation sources—A numerical study , 1979, Bulletin of the Seismological Society of America.

[18]  Jian-Ming Jin,et al.  A spectral Lanczos decomposition method for solving 3-D low-frequency electromagnetic diffusion by the finite-element method , 1999 .

[19]  Michal Mrozowski,et al.  Jacobi and Gauss-Seidel preconditioned complex conjugate gradient method with GPU acceleration for finite element method , 2010, The 40th European Microwave Conference.

[20]  Gregory A. Newman,et al.  Electromagnetic induction in a fully 3‐D anisotropic earth , 2000 .

[21]  He Li-li.  Optimizing Sparse Matrix-vector Multiplication Based on GPU , 2010 .

[22]  C. Weiss,et al.  Full-physics 3-D heterogeneous simulations of electromagnetic induction fields on level and deformed sea ice , 2015, Annals of Glaciology.

[23]  A. Lichtenberg,et al.  Linear electromagnetic excitation of an asymmetric low pressure capacitive discharge with unequal sheath widths , 2016 .

[24]  Rajesh Bordawekar,et al.  Optimizing Sparse Matrix-Vector Multiplication on GPUs using Compile-time and Run-time Strategies , 2008 .

[25]  Michael Commer,et al.  A parallel finite-difference approach for 3D transient electromagnetic modeling with galvanic sources , 2004 .

[26]  Wolfgang Paul,et al.  GPU accelerated Monte Carlo simulation of the 2D and 3D Ising model , 2009, J. Comput. Phys..

[27]  Alan C. Tripp,et al.  Electromagnetic scattering of large structures in layered earths using integral equations , 1995 .

[28]  Oliver Ritter,et al.  Controlled‐source electromagnetic monitoring of reservoir oil saturation using a novel borehole‐to‐surface configuration , 2015 .

[29]  T. Takahashi,et al.  GPU‐accelerated boundary element method for Helmholtz' equation in three dimensions , 2009 .

[30]  Karol Miller,et al.  Real-Time Nonlinear Finite Element Computations on GPU - Application to Neurosurgical Simulation. , 2010, Computer methods in applied mechanics and engineering.

[31]  Y. Mitsuhata.  2-D electromagnetic modeling by finite‐element method with a dipole source and topography , 2000 .

[32]  Leonard J. Srnka,et al.  Special Section — Marine Controlled-Source Electromagnetic Methods An introduction to marine controlled-source electromagnetic methods for hydrocarbon exploration , 2007 .

[33]  Satoshi Matsuoka,et al.  Fast Conjugate Gradients with Multiple GPUs , 2009, ICCS.

[34]  Sharon C. Glotzer,et al.  Efficient neighbor list calculation for molecular simulation of colloidal systems using graphics processing units , 2016, Comput. Phys. Commun..

[35]  Xavier Andrade,et al.  Application of Graphics Processing Units to Accelerate Real‐Space Density Functional Theory and Time‐Dependent Density Functional Theory Calculations , 2016 .

[36]  Jie Zhang,et al.  A GPU approach to FDTD for radio coverage prediction , 2008, 2008 11th IEEE Singapore International Conference on Communication Systems.

[37]  David B. Davidson,et al.  GPU-based LU decomposition for large method of moments problems , 2010 .