A parallel goal-oriented adaptive finite element method for 2.5D electromagnetic modeling
Kerry Key and Jeffrey Ovall
Contact the author: Please use for correspondence this email.
Submission date: 03. Aug. 2010 (revised version: June 2011)
published in: Geophysical journal international, 186 (2011) 1, p. 137-154
DOI number (of the published article): 10.1111/j.1365-246X.2011.05025.x
Keywords and phrases: marine electromagnetics, adaptive finite elements
Download full preprint: PDF (9359 kB)
We present a parallel goal-oriented adaptive finite element algorithm that can be used to rapidly compute highly accurate solutions for 2.5D controlled-source electromagnetic (CSEM) and 2D magnetotelluric (MT) modeling problems. We employ unstructured triangular grids to permit efficient discretization of complex modeling domains such as those containing topography, dipping layers and multiple scale structures. Iterative mesh refinement is guided by a goal-oriented error estimator based on a form of dual residual weighting, which is carried out using hierarchical basis computations. Our formulation of the error estimator considers the relative error in the strike aligned fields and their spatial gradients, and therefore results in a more efficient use of mesh vertices than previous error estimators based on absolute field errors. This algorithm is parallelized over frequencies, transmitters, receivers and wave-numbers, where adaptive refinement can be performed in parallel on subsets of these parameters while nearby parameters are able to share the refined grid, thus enabling our algorithm to achieve accurate solutions in run-times of seconds to tens of seconds for realistic models and data parameters when run on cluster computers containing about a thousand processors. Application of this new algorithm to a complex model that includes strong seafloor topography variations and multiple thin stacked reservoirs demonstrates the performance and scalability on a large cluster computer.