Skip to Main Content
Skip Nav Destination

A parallel finite-element time-domain method for transient electromagnetic simulation

Fu Haohuan, Wang Yingqiao, Evan Schankee Um, Fang Jiarui, Wei Tengpeng, Huang Xiaomeng and Yang Guangwen
A parallel finite-element time-domain method for transient electromagnetic simulation
Geophysics (July 2015) 80 (4): E213-E224


The finite-difference time-domain (FDTD) method has been widely used to perform transient electromagnetic simulation. Compared to the traditional FDTD method, the finite-element time-domain (FETD) method with unstructured meshes and an adaptive time-stepping scheme has the potential to dramatically cut down the number of unknowns and the number of time steps. However, with a need to solve a large-scale unstructured sparse matrix at every time step, the FETD method is generally difficult to scale with the available parallel computing resources. To make efficient use of the multicore resources in current computing platforms, we designed a parallel FETD method that used multiple threads to accelerate the key steps for solving a large sparse matrix, including the construction of the preconditioner and matrix-vector multiplication. Using a customized incomplete Cholesky preconditioner, we achieved faster convergence than many existing preconditioners and a scalable parallel performance over multiple cores. To further reduce the total calculation time, we implemented an overlapping scheme and a preconditioner-reusing scheme. The overlapping scheme executed the solver computation and the preconditioner computation concurrently using different parts of the parallel resources. The preconditioner reusing scheme reduced the need to recompute preconditioners when increasing the time-step size, thus reducing the frequency of the less scalable preconditioner computations. It brought performance benefits by reducing the frequency of less scalable preconditioner computation. Computed parallel FETD results showed the same level of accuracy as FDTD results. In terms of parallel performance, using eight cores, the method cut execution time by a factor of 4 (i.e., 50 parallel efficiency). For a realistic salt model, we simulated TEM responses from 0.01 s to 30 s in 5 min, using hundreds of CPU cores.

ISSN: 0016-8033
EISSN: 1942-2156
Coden: GPYSA7
Serial Title: Geophysics
Serial Volume: 80
Serial Issue: 4
Title: A parallel finite-element time-domain method for transient electromagnetic simulation
Affiliation: Tsinghua University, Laboratory for Earth System Modeling, Beijing, China
Pages: E213-E224
Published: 201507
Other Date: 20150602
Text Language: English
Publisher: Society of Exploration Geophysicists, Tulsa, OK, United States
References: 29
Accession Number: 2015-098582
Categories: Applied geophysics
Document Type: Serial
Bibliographic Level: Analytic
Annotation: Includes appendices
Illustration Description: illus. incl. 3 tables
Secondary Affiliation: Lawrence Berkeley National Laboratory, USA, United States
Country of Publication: United States
Secondary Affiliation: GeoRef, Copyright 2017, American Geosciences Institute. Reference includes data from GeoScienceWorld, Alexandria, VA, United States. Reference includes data supplied by Society of Exploration Geophysicists, Tulsa, OK, United States
Update Code: 201542
Close Modal

or Create an Account

Close Modal
Close Modal