3D magnetotelluric modeling using high-order tetrahedral Nédélec elements on massively parallel computing platforms
We 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 imple- mentation extends the PETGEM modeller capabilities, initially developed for active-sourc...
| Autores: | , , , , , , , |
|---|---|
| Tipo de recurso: | artículo |
| Fecha de publicación: | 2022 |
| País: | España |
| Institución: | Universitat Politècnica de Catalunya (UPC) |
| Repositorio: | UPCommons. Portal del coneixement obert de la UPC |
| Idioma: | inglés |
| OAI Identifier: | oai:upcommons.upc.edu:2117/363026 |
| Acceso en línea: | https://hdl.handle.net/2117/363026 https://dx.doi.org/10.1016/j.cageo.2021.105030 |
| Access Level: | acceso abierto |
| Palabra clave: | Imaging systems in geophysics High performance computing Magnetotelluric prospecting Magnetotelluric method Geophysical electromagnetics Numerical solutions High-order edge finite element Càlcul intensiu (Informàtica) Prospecció magnetotel·lúrica Àrees temàtiques de la UPC::Informàtica::Arquitectura de computadors Àrees temàtiques de la UPC::Enginyeria civil::Geologia |
| Sumario: | We 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 imple- mentation 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. |
|---|