A New 3D Marine Controlled-Source Electromagnetic Modeling Algorithm Based on Two New Types of Quadratic Edge Elements

This study proposes a new algorithm for a higher-order vector finite element method based on two new types of second-order edge elements to solve the electromagnetic field diffusion problem in a 3D anisotropic medium. To avoid source singularity in the quasistatic variant of Maxwell’s function, a secondary field formulation was adopted. The modeling domain was discretized using two types of quadratic edge hexahedral elements, which were obtained using the edge unification method to reduce variables on each side of two conventional quadratic edge elements. Compared with the traditional quadratic element, the number of unknowns that needed to be solved was significantly reduced. The sparse linear equation of the finite element system was solved using an open-source direct solver called MUMPS. The numerical results demonstrated that the proposed algorithm has the same level of accuracy as the conventional vector finite element method and has a significant advantage over it in terms of computational cost.