Parallel finite element algorithm with domain decomposition for three-dimensional magnetotelluric modelling

We present a new finite element (FE) method for magnetotelluric modelling of three-dimensional conductivity structures. Maxwell's equations are treated as a system of first-order partial differential equations for the secondary fields. Absorbing boundary conditions are introduced, minimizing undesired boundary effects and allowing the use of small computational domains. The numerical algorithm presented here is an iterative, domain decomposition procedure employing a nonconforming FE space. It does not use global matrices, therefore allowing the modellization of large and complicated structures. The algorithm is naturally parallellizable, and we show results obtained in the IBM SP2 parallel supercomputer at Purdue University. The accuracy of the numerical method is verified by checking the computed solutions with the results of COMMEMI, the international project on the comparison of modelling methods for electromagnetic induction.