3-D direct current resistivity forward modeling by adaptive multigrid finite element method  

作　　者：汤井田 王飞燕 任政勇 郭荣文 

机构地区：[1]School of Info-physics and Geomatics Engineering, Central South University [2]Institute of Geophysics, ETH Zurich-Swiss Federal Institute of Technology [3]School of Earth and Ocean Sciences, University of Victoria

出　　处：《Journal of Central South University》2010年第3期587-592,共6页中南大学学报（英文版）

基　　金：Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of China;Projects(2007, 2008) supported by China Scholarship Council (CSC)

摘　　要：Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.

关 键 词：adaptive multigrid a-posteriori error estimator unstructured mesh V-CYCLE finite element method 

分 类 号：P631.322[天文地球—地质矿产勘探]