3D Magnetotelluric Modeling Using High-Order Tetrahedral Nédélec Elements on Massively Parallel Computing Platforms

dc.affiliation.dptoUC3M. Departamento de Teoría de la Señal y Comunicacioneses
dc.contributor.authorCastillo Reyes, Octavio
dc.contributor.authorModesto, David
dc.contributor.authorQueralt, Pilar
dc.contributor.authorMarcuello, Alex
dc.contributor.authorLedo, Juanjo
dc.contributor.authorAmor Martín, Adrián
dc.contributor.authorde la Puente, Josep
dc.contributor.authorGarcía Castillo, Luis Emilio
dc.contributor.funderEuropean Commissionen
dc.contributor.funderMinisterio de Ciencia e Innovación (España)es
dc.description.abstractWe present a routine for 3D magnetotelluric (MT) modeling based upon high-order edge finite element method (HEFEM), tailored and unstructured tetrahedral meshes, and high-performance computing (HPC). This implementation extends the PETGEM modeller capabilities, initially developed for active-source electromagnetic methods in frequency-domain. We assess the accuracy, robustness, and performance of the code using a set of reference models developed by the MT community in well-known reported workshops. The scale and geological properties of these 3D MT setups are challenging, making them ideal for addressing a rigorous validation. Our numerical assessment proves that this new algorithm can produce the expected solutions for arbitrarily 3D MT models. Also, our extensive experimental results reveal four main insights: (1) high-order discretizations in conjunction with tailored meshes can offer excellent accuracy; (2) a rigorous mesh design based on the skin-depth principle can be beneficial for the solution of the 3D MT problem in terms of numerical accuracy and run-time; (3) high-order polynomial basis functions achieve better speed-up and parallel efficiency ratios than low-order polynomial basis functions on cutting-edge HPC platforms; (4) a triple helix approach based on HEFEM, tailored meshes, and HPC can be extremely competitive for the solution of realistic and complex 3D MT models and geophysical electromagnetics in general.en
dc.description.sponsorshipThis project has been 65% cofinanced by the European Regional Development Fund (ERDF) through the Interreg V-A Spain–France– Andorra program (POCTEFA2014-2020). POCTEFA aims to reinforce the economic and social integration of the French–Spanish–Andorran border. Its support is focused on developing economic, social and environmental cross-border activities through joint strategies favoring sustainable territorial development. BSC authors received funding from the European Union’s Horizon 2020 programme, grant agreement N◦828947 and N◦777778, and from the Mexican Department of Energy, CONACYT-SENER Hidrocarburos grant agreement N◦B-S-69926.en
dc.identifier.bibliographicCitationCastillo-Reyes, O., Modesto, D., Queralt, P., Marcuello, A., Ledo, J., Amor-Martin, A., de la Puente, J., & García-Castillo, L. E. (2022). 3D magnetotelluric modeling using high-order tetrahedral Nédélec elements on massively parallel computing platforms. In Computers & Geosciences, 160, 105030-105042en
dc.identifier.publicationtitleCOMPUTERS & GEOSCIENCESen
dc.publisherELSEVIER BVen
dc.relation.projectIDGobierno de España. POCTEFA2014-2020es
dc.rights© 2022 The Author(s). Published by Elsevier Ltd.en
dc.rightsAtribución 3.0 España*
dc.rights.accessRightsopen accessen
dc.subject.otherMagnetotelluric methoden
dc.subject.otherGeophysical electromagneticsen
dc.subject.otherNumerical solutionsen
dc.subject.otherHigh-order edge finite elementen
dc.subject.otherHigh-performance computingen
dc.title3D Magnetotelluric Modeling Using High-Order Tetrahedral Nédélec Elements on Massively Parallel Computing Platformsen
dc.typeresearch article*
Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
1.77 MB
Adobe Portable Document Format