A CELL-CENTERED MULTIGRID SOLVER FOR THE FINITE VOLUME DISCRETIZATION OF ANISOTROPIC ELLIPTIC INTERFACE PROBLEMS ON IRREGULAR DOMAINS  

作　　者：Kejia Pan Xiaoxin Wu Hongling Hu Zhilin Li 

机构地区：[1]School of Mathematics and Statistics,HNP-LAMA,Central South University,Changsha 410083,China [2]Key Laboratory of Computing and Stochastic Mathematics(Ministry of Education),School of Mathematics and Statistics,Hunan Normal University,Changsha 410081,China [3]Department of Mathematics,North Carolina State University,Raleigh,NC 27695-8205,USA

出　　处：《Journal of Computational Mathematics》2025年第1期18-42,共25页计算数学(英文)

基　　金：supported by the National Natural Science Foundation of China(Grant No.42274101);X.X.Wu was supported by the Fundamental Research Funds for the Central Universities of Central South University(Grant No.2020zzts354);H.L.Hu was supported by the National Natural Science Foundation of China(Grant No.12071128);by the Natural Science Foundation of Hunan Province(Grant No.2021JJ30434);Z.L.Li was supported by a Simons Grant No.633724.

摘　　要：The aim of this paper is to develop a fast multigrid solver for interpolation-free finite volume (FV) discretization of anisotropic elliptic interface problems on general bounded domains that can be described as a union of blocks. We assume that the curved interface falls exactly on the boundaries of blocks. The transfinite interpolation technique is applied to generate block-wise distorted quadrilateral meshes, which can resolve the interface with fine geometric details. By an extensive study of the harmonic average point method, an interpolation-free nine-point FV scheme is then derived on such multi-block grids for anisotropic elliptic interface problems with non-homogeneous jump conditions. Moreover, for the resulting linear algebraic systems from cell-centered FV discretization, a high-order prolongation operator based fast cascadic multigrid solver is developed and shown to be robust with respect to both the problem size and the jump of the diffusion coefficients. Various non-trivial examples including four interface problems and an elliptic problem in complex domain without interface, all with tens of millions of unknowns, are provided to show that the proposed multigrid solver is dozens of times faster than the classical algebraic multigrid method as implemented in the code AMG1R5 by Stüben.

关 键 词：Elliptic interface problem Discontinuous coefficients Anisotropic coefficients Cascadic multigrid method Richardson extrapolation 

分 类 号：O24[理学—计算数学]