High order edge-based finite elements for 3D magnetotelluric modeling with unstructured meshes

We present a 3D magnetotelluric (MT) forward modeling program which computes the distribution of the electric field on an unstructured mesh. We employ high order edge elements and implement the Edge Finite Element method to accurately simulate magnetotelluric data in the form of impedance tensor. The program is open-source and optimized for both high-performance computing (HPC) architectures and nonHPC architectures. The paper details the theory, method and code work-flow. We evaluate the accuracy of our algorithm with various synthetic test models and compare our results with other Finite Element (FE) and Finite Difference (FD) algorithms. The experiments show that the implemented high order elements are effective in providing highly accurate solutions for multiple frequencies with few elements independently of the model's geometry. Hence, leading to reduced computing times.

Automated summary

The study presents a 3D magnetotelluric forward modeling program, which uses high-order edge elements and the Edge Finite Element method to accurately simulate magnetotelluric data. The program is open-source and optimized for both HPC and nonHPC architectures, providing highly accurate solutions for multiple frequencies with reduced computing times.