Tailored meshing for parallel 3D electromagnetic modeling using high-order edge elements
Abstract
We present numerical experiments for geophysics electromagnetic (EM) modeling based upon high-order edge elements and supervised refinement approaches on massively parallel computers. Our high-order refinement strategy is based on and extends the PETGEM code. We focus on the performance study in terms of accuracy, convergence rate, and computational effort to solve real-life 3D setups based on synthetic and experimental data for energy reservoir characterization. These test cases show variable resolution discretization needs and realistic physical parameters. In general, our numerical results are consistent theoretically. The use of adapted meshes was efficient to achieve a certain accuracy level in the synthetic EM responses. Regarding global refinement, exhibits the best accuracy/performance trade-off. Selective -refinement might offer a better compromise between accuracy and computational cost. However, for refinement at different entities, the best refinement scheme consists of using at the volume level with at faces and edges. Thus, refinement can be competitive if applied hierarchically. Nevertheless, we acknowledge that the performance of our supervised refinement strategy depends on the input model (e.g., conductivity, frequency, domain decomposition strategy, among others). Whatever the chosen configuration, our numerical results provide an in-depth understanding of EM modeling's pros and cons when supervised refinement schemes are applied.
Cite
@article{arxiv.2208.00688,
title = {Tailored meshing for parallel 3D electromagnetic modeling using high-order edge elements},
author = {Octavio Castillo-Reyes and Adrian Amor-Martin and Arnaud Botella and Pierre Anquez and Luis Emilio García-Castillo},
journal= {arXiv preprint arXiv:2208.00688},
year = {2022}
}