应用数学青年讨论班(午餐会)—— p-Multigrid method for elliptic problem
Speaker(s):张东航(中国科学院数学与系统科学研究院基础软件研究中心)
Time:2024-10-09 11:45-13:00
Venue:智华楼四元厅
摘要:
In this talk, we propose the two-level and W-cycle algorithms of p-multigrid method designed to solve the linear systems of equations generated from p-version symmetric interior penalty discontinuous Galerkin (SIPDG) discretizations for elliptic problems. This SIPDG discretization employs hierarchical Legendre polynomial basis functions, where we can design restriction and prolongation operators between different discrete polynomial spaces naturally. Inspired by the uniform convergence theory of the W-cycle algorithm of hp-multigrid method in [P. F. Antonietti, et.al., SIAM J. Numer. Anal., 53 (2015)], we extend their work by providing a more refined matrix-based analysis. Specifically, we estimate the spectral radius of the stiffness matrix and its diagonal matrix, assess the approximation property of coarsest level correction, and analyze the smoothing properties of polynomial smoother based on fourth-kind Chebyshev polynomial iterative method. Building on these foundations, we provide a rigorous matrix-based convergence analysis for the proposed p-multigrid method, considering both inherited and non-inherited bilinear forms of SIPDG discretization. Our theoretical results show significant improvement over [P. F. Antonietti, et.al., SIAM J. Numer. Anal., 53 (2015)], reducing the required number of smoothing steps from O(p^2) to O(p), where p is the polynomial degree of the discrete broken polynomial space. Moreover, the convergence rate remains independent of the mesh size. Finally, several numerical experiments are presented to validate our theoretical findings.
报告人简介:
张东航,2015-2020 中国科学院数学与系统科学院,计算数学博士;2020-2023北京大学北京国际数学研究中心,博士后; 2023年至今,中国科学院数学与系统科学研究院基础软件研究中心。
报名问卷:
我们从11:45开始按照问卷情况提供午餐,请需要预定午餐的老师同学填写此报名问卷,10月8日(周二)下午3点截止。