Calculating 2D Helmholtz equation based on multigrid method and application in magnetotelluric modeling
Yang Zhenwei1,2, Feng Lei1,2, Zhao Ning3, Zhao Qiufang1,2, Yang Shuang'an1,2
1. Institute of Resources & Environment, Henan Polytechnic University, Jiaozuo, Henan 454000, China;
2. Collaborative Innovation Center of Coalbed Methane and Shale Gas for Central Plains Economic Region, Henan Province, Jiaozuo, Henan 454000, China;
3. College of Computer Sciences and Technology, Henan Polytechnic University, Jiaozuo, Henan 454000, China
Abstract:In order to enhance the computational efficiency and accuracy of Helmholtz equation,we discuss in the paper 2D Helmholtz equation calculation with multi-grid method and other methods such as conjugate gradient (CG),preconditioning conjugate gradient (PCG),and successive over relaxation (SOR).Delaunay triangular mesh algorithm with mesh refinement is automatically adopted in mesh subdivision.Research results show that the multi-grid method has great advantage in computation time consuming and iterative convergence efficiency,but its computational error is larger than the other algorithms.This might be caused by irregular grid subdivision,which is related with interpolation and restriction operator.Finally,a preliminary test on 2D magnetotelluric modeling with multi-grid algorithm is shown in the paper.
Wang Qingping,Bai Wuming,Wang Hongliang.Multigrid finite element analysis for 2D modeling of Poisson equation.Progress in Geophysics,2010,25(4):1467-1474.
Han Bo,Hu Xiangyun,He Zhanxiang et al.Mathematical classification of magnetotelluric inversion methods.OGP,2012,47(1):177-187.
[7]
Airaksinen T,Heikkola E,Pennanen A et al.An algebraic multigrid based shifted-Laplacian preconditioner for the Helmholtz equation.Journal of Computational Physics,2007,226(1):1196-1210.
[8]
Mulder W A.A multigrid solver for 3D electromagnetic diffusion.Geophysical Prospecting,2006,54(5):633-649.
[9]
Plaks A,Tsukerman I,Painchaud S et al.Multigrid methods for open boundary problems in geophysics.IEEE Transactions on Magnetics,2000,36(4):633-638.
[10]
Zaslavsky L Y,Schlick T.An adaptive multigrid technique for evaluating long-range forces in biomolecular simulations.Applied Mathematics and Computation,1998,97(2-3):237-250.
[11]
Tang J T,Wang F Y,Ren Z Y et al.3-D direct current resistivity forward modeling by adaptive multigrid finite element method.Journal of Central South University,2010,17(3):587-592.
Liu Jianxin,Guo Rongwen,Tong Xiaozhong.Boundary truncation of magnetotelluric modeling based on multigrid method.Journal of Central South University (Science and Technology),2011,42(11):3429-3437.
[13]
Mitsuhata Y J,Uchida T.3D magnetotelluric modeling using the T-Ω finite element method.Geophysics,2004,69(1):108-119.
[14]
Moucha R,Bailey R C.An accurate and robust multigrid algorithm for 2D forward resistivity modeling.Geophysical Prospecting,2004,52(3):197-212.
[15]
William L B.A Multigrid Tutorial.Society for Industrial and Applied Mathematics,2000.
[16]
Liu Yuanqing,Yuan Jiansheng.A finite element do-main decomposition combined with algebraic multigrid method for large-scale electromagnetic field computation.IEEE Transactions on Magnetics,2006,42(4):655-658.