Assessment of edge-based finite element technique for geophysical electromagnetic problems: efficiency, accuracy and reliability
Document typeConference lecture
PublisherInternational Center for Numerical Methods in Engineering (CIMNE)
Rights accessOpen Access
European Commission's projectGEAGAM - Geophysical Exploration using Advanced GAlerkin Methods (EC-H2020-644202)
In Finite Element Methods for solving electromagnetic field problems, the use of Edge Elements has become very popular. In fact, Edge Elements are often said to be a cure to many difficulties that are encountered (particularly eliminating spurious solutions) and are claimed to yield accurate results [9, 16, 17]. We will shortly describe the mathematical formulation of linear edge elements and we go through the particular issues related to the implementation of these elements in order to solve geophysical electromagnetic problems. In particular, we describe a simple, flexible and parallel Fortran 90 implementation for Edge Elements. The code is based on an abstract data structure, which allows to use different kinds of solvers with little effort. The result is an implementation that allows users to specify Edge-based Finite Element variational forms of H(curl). Finally, we also show the performance of the code in terms of efficiency, accuracy and reliability, which will shape our future line of work in order to solve more complex problems.
CitationCastillo-Reyes, Octavio [et al.]. Assessment of edge-based finite element technique for geophysical electromagnetic problems: efficiency, accuracy and reliability. A: 1st Pan-American Congress on Computational Mechanics (PANACM 2015) and XI Argentine Congress on Computational Mechanics (MECOM 2015), Buenos Aires, Argentina 27-29 April 2015. "Proceedings of the 1st Pan-American Congress on Computational Mechanics". Barcelona: International Center for Numerical Methods in Engineering (CIMNE), 2015, p. 984-995.