Immersed Finite Element Method for Interface Problems with Algebraic Multigrid Solver  

作　　者：Wenqiang Feng Xiaoming He Yanping Lin Xu Zhang 

机构地区：[1]Department of Mathematics and Statistics,Missouri University of Science and Technology,Rolla,MO 65409,USA [2]Department of Applied Mathematics,Hong Kong Polytechnic University,Hung Hom,Hong Kong [3]Department of Mathematical and Statistics Science,University of Alberta,Edmonton,AB,T6G 2G1,Canada [4]Department of Mathematics,Virginia Tech,Blacksburg,VA 24061,USA

出　　处：《Communications in Computational Physics》2014年第4期1045-1067,共23页计算物理通讯（英文）

基　　金：supported by DOE grant DE-FE0009843;National Natural Science Foundation of China(11175052);GRF of HKSAR#501012 and NSERC(Canada).

摘　　要：This article is to discuss the bilinear and linear immersed finite element(IFE)solutions generated from the algebraic multigrid solver for both stationary and moving interface problems.For the numerical methods based on finite difference formulation and a structured mesh independent of the interface,the stiffness matrix of the linear system is usually not symmetric positive-definite,which demands extra efforts to design efficient multigrid methods.On the other hand,the stiffness matrix arising from the IFE methods are naturally symmetric positive-definite.Hence the IFE-AMG algorithm is proposed to solve the linear systems of the bilinear and linear IFE methods for both stationary and moving interface problems.The numerical examples demonstrate the features of the proposed algorithms,including the optimal convergence in both L 2 and semi-H1 norms of the IFE-AMG solutions,the high efficiency with proper choice of the components and parameters of AMG,the influence of the tolerance and the smoother type of AMG on the convergence of the IFE solutions for the interface problems,and the relationship between the cost and the moving interface location.

关 键 词：Interface problems immersed finite elements algebraic multigrid method 

分 类 号：O17[理学—数学]