An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials

doi:10.1093/gji/ggac152

	An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials
	Pan, Kejia 1; Wang, Jinxuan 1; Hu, Shuanggui 1; Ren, Zhengyong 2; Cui, Tao 3; Guo, Rongwen 2; Tang, Jingtian 2
刊名	GEOPHYSICAL JOURNAL INTERNATIONAL
	2022-05-20
卷号	230 期号:3 页码:1834-1851
关键词	Electrical properties Magnetotellurics Numerical modelling Numerical solutions
ISSN号	0956-540X
DOI	10.1093/gji/ggac152
英文摘要	The fast and accurate 3-D magnetotelluric (MT) forward modelling is core engine of the interpretation and inversion of MT data. In this study, we develop an improved extrapolation cascadic multigrid method (EXCMG) to solve the large sparse complex linear system arising from the finite-element (FE) discretization on non-uniform orthogonal grids of the Maxwell's equations using potentials. First, the vector Helmholtz equation and the scalar auxiliary equation are derived from the Maxwell's equations using Coulomb-gauged potentials. The weighted residual method is adopted to discretize the weak formulation and assemble the FE equation. Secondly, carefully choosing the preconditioned complex stable bi-conjugate gradient method (BiCGStab) as multigrid smoother, we develop an improved EXCMG method on non-uniform grids to solve the resulting large sparse complex non-Hermitian linear systems. Finally, several examples including three standard testing models (COMMEMI3D-1, COMMEMI3D-2 and DTM1.0) and a topographic model are used to validate the accuracy and efficiency of the proposed multigrid solver. Numerical results show that the proposed EXCMG algorithm greatly improves the efficiency of 3-D MT forward modelling, is more efficient than some existing solvers, such as Pardiso, incomplete LU factorization preconditioned biconjugate gradients stabilized method (ILU-BiCGStab) and flexible generalized minimum residual method with auxiliary space Maxwell preconditioner (FGMRES-AMS), and capable to simulate large-scale problems with more than 100 million unknowns.
资助项目	National Natural Science Foundation of China[41922027] ; Excellent Youth Foundation of Hunan Province of China[2019JJ20032] ; Fundamental Research Funds for the Central Universities of Central South University[1053320212467] ; Natural Science Foundation of China[42174171] ; High Performance Computing Center of Central South University
WOS研究方向	Geochemistry & Geophysics
语种	英语
出版者	OXFORD UNIV PRESS
WOS记录号	WOS:000798346100003
内容类型	期刊论文
源URL	[http://ir.amss.ac.cn/handle/2S8OKBNM/61381]
专题	中国科学院数学与系统科学研究院
通讯作者	Wang, Jinxuan
作者单位	1.Cent South Univ, Sch Math & Stat, Changsha 410083, Peoples R China 2.Cent South Univ, Sch Geosci & Infophys, Changsha 410083, Peoples R China 3.Chinese Acad Sci, Acad Math & Syst Sci, LSEC, NCMIS, Beijing 100190, Peoples R China
推荐引用方式 GB/T 7714	Pan, Kejia,Wang, Jinxuan,Hu, Shuanggui,et al. An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials[J]. GEOPHYSICAL JOURNAL INTERNATIONAL,2022,230(3):1834-1851.
APA	Pan, Kejia.,Wang, Jinxuan.,Hu, Shuanggui.,Ren, Zhengyong.,Cui, Tao.,...&Tang, Jingtian.(2022).An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials.GEOPHYSICAL JOURNAL INTERNATIONAL,230(3),1834-1851.
MLA	Pan, Kejia,et al."An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials".GEOPHYSICAL JOURNAL INTERNATIONAL 230.3(2022):1834-1851.