Parallel 3-D marine controlled-source electromagnetic modelling using high-order tetrahedral Nédélec elements
Visualitza/Obre
Cita com:
hdl:2117/166189
Tipus de documentArticle
Data publicació2019-10
EditorOxford University Press
Condicions d'accésAccés obert
Tots els drets reservats. Aquesta obra està protegida pels drets de propietat intel·lectual i
industrial corresponents. Sense perjudici de les exempcions legals existents, queda prohibida la seva
reproducció, distribució, comunicació pública o transformació sense l'autorització del titular dels drets
ProjecteMATHROCKS - Multiscale Inversion of Porous Rock Physics using High-Performance Simulators: Bridging the Gap between Mathematics and Geophysics (EC-H2020-777778)
ChEESE - Centre of Excellence for Exascale in Solid Earth (EC-H2020-823844)
ChEESE - Centre of Excellence for Exascale in Solid Earth (EC-H2020-823844)
Abstract
We present a parallel and high-order Nédélec finite element solution for the marine controlled-source electromagnetic (CSEM) forward problem in 3-D media with isotropic conductivity. Our parallel Python code is implemented on unstructured tetrahedral meshes, which support multiple-scale structures and bathymetry for general marine 3-D CSEM modelling applications. Based on a primary/secondary field approach, we solve the diffusive form of Maxwell’s equations in the low-frequency domain. We investigate the accuracy and performance advantages of our new high-order algorithm against a low-order implementation proposed in our previous work. The numerical precision of our high-order method has been successfully verified by comparisons against previously published results that are relevant in terms of scale and geological properties. A convergence study confirms that high-order polynomials offer a better trade-off between accuracy and computation time. However, the optimum choice of the polynomial order depends on both the input model and the required accuracy as revealed by our tests. Also, we extend our adaptive-meshing strategy to high-order tetrahedral elements. Using adapted meshes to both physical parameters and high-order schemes, we are able to achieve a significant reduction in computational cost without sacrificing accuracy in the modelling. Furthermore, we demonstrate the excellent performance and quasi-linear scaling of our implementation in a state-of-the-art high-performance computing architecture.
CitacióCastillo-Reyes, O. [et al.]. Parallel 3-D marine controlled-source electromagnetic modelling using high-order tetrahedral Nédélec elements. "Geophysical Journal International", Octubre 2019, vol. 219, núm. 1, p. 39-65.
ISSN0956-540X
Col·leccions
Fitxers | Descripció | Mida | Format | Visualitza |
---|---|---|---|---|
Parallel 3D marine controlled-source.pdf | 5,896Mb | Visualitza/Obre |