二维边值问题有限元离散系统病态问题的病因抑制方法
A Method of Inhibiting Pathogeny for Ill-Condition Problems of the Finite Element Discrete Systems of 2D Boundary Value Problems
DOI: 10.12677/AAM.2021.1012442, PDF,    国家自然科学基金支持
作者: 张 衡:无损检测技术福建省高等学校重点实验室(福建技术师范学院),福建 福清;福建技术师范学院大数据与人工智能学院,福建 福清;湖南三一工业职业技术学院工程机械学院,湖南 长沙 ;张 宏:无损检测技术福建省高等学校重点实验室(福建技术师范学院),福建 福清;福建技术师范学院电子与机械工程学院,福建 福清;游文杰:福建技术师范学院大数据与人工智能学院,福建 福清;郑汉垣:龙岩学院传播与设计学院,福建 龙岩
关键词: 病因抑制方法病态结构病态因子去病因子最优预条件子Method of Inhibiting Pathogeny Ill-Condition Structure Ill-Condition Factor Ill Elimination Factor Optimal Preconditioner
摘要: 讨论二维Poisson方程边值问题离散系统的病态问题,针对三角剖分和和四角剖分,基于Lagrange形函数形成有限元离散系统的病态问题,提出病因抑制方法;给出该系统的病态结构、病态因子、去病因子;利用病态因子估计离散系统的条件数;利用去病因子为最优预条件子,精准抑制病态发作,该预条件子的使用,几乎不增加求解的计算量,预处理后离散系统保持正定对称,条件数关于离散系统规模一致有界。
Abstract: The ill conditioned problem of the finite element discrete system based on triangular subdivision, tetragonal subdivision and Lagrange shape function for 2D Poisson equation boundary value problem is discussed. The method of inhibiting pathogeny was proposed. The ill condition structure, ill condition factor and ill elimination factor of the system were given; and the condition number of the equation was estimated by the ill conditioned factor. The ill elimination factor was used as the optimal preconditioner to precisely suppress the ill condition. The use of the preconditioner hardly increases the calculation of iteration. After the pretreatment, the equations keep positive definite symmetry and the condition number is uniformly bounded with respect to the scale of equations.
文章引用:张衡, 张宏, 游文杰, 郑汉垣. 二维边值问题有限元离散系统病态问题的病因抑制方法[J]. 应用数学进展, 2021, 10(12): 4162-4171. https://doi.org/10.12677/AAM.2021.1012442

参考文献

[1] Zhong X.J. (1988) The Convergence of Krylov Subspace Methods for Large Nonsymmertric Linear Systems. Acta Mathematical Sinica, 14, 507-518. [Google Scholar] [CrossRef
[2] 吴勃英, 王德明, 丁效华, 李道华. 数值分析原理[M]. 北京: 科学出版社, 2003: 71-101.
[3] 张科. 两类线性方程组的预处理技术及数值求解方法[D]: [博士学位论文]. 上海: 上海大学, 2014: 7-9.
[4] 李荣华. 偏微分方程数值解法[M]. 第2版. 北京: 高等教育出版社, 2010: 139-151, 13-15.
[5] 于春肖, 苑润浩, 穆运峰. 新预处理ILUCG法求解稀疏病态线性方程组[J]. 数值计算与计算机应用, 2014, 35(1): 21-27.
[6] Xue, Q.F., Gao, X.B. and Liu, X.G. (2014) Comparison Theorems for a Class of Preconditioned AOR Iterative Methods. Journal of Mathematics, 34, 448-460.
[7] 潘春平, 马成荣, 曹文方, 王红玉. 一类预条件AOR迭代法的收敛性分析[J]. 数学杂志, 2013, 33(3): 479-484.
[8] 罗芳. L-矩阵的多参数预条件AOR迭代法[J]. 数学的实践与认识, 2013, 43(15): 277-282.
[9] 吴建平, 赵军, 马怀发, 宋君强, 张卫民, 李晓梅. 一般稀疏线性方程组的因子组合型并行预条件研究[J]. 计算机应用与软件, 2012, 29(5): 6-9, 108.
[10] 李继成, 蒋耀林. 预条件IMGS迭代方法的比较定理[J]. 数学物理学报, 2011, 31(4): 880-886.
[11] Pan, C.P. (2011) A New Effective Preconditioned Gauss-Seidel Iteration Method. Journal on Numerical Methods and Computer Applications, 32, 267-273.
[12] 张文生. 科学计算中的偏微分方程有限差分法[M]. 北京: 高等教育出版社, 2006: 205-209.
[13] Bai. Z.Z. and Li, G.Q. (2003) Restrictively Preconditioned Conjugate Gradient Methods for Systems of Linear Equations. IMA Journal of Numerical Analysis, 23, 561-580. [Google Scholar] [CrossRef
[14] Wang, Z.Q. (2015) Restrictively Preconditioned Chebyshev Method for Solving Systems of Linear Equations. Journal of Engineering Mathematics, 93, 61-76. [Google Scholar] [CrossRef
[15] Shen, H.L., Shao, X.H. and Zhang, T. (2012) Preconditioned Iterative Methods for Solving Weighted Linear Least Squares Problems. Applied Mathematics and Mechanics, 33, 375-384. [Google Scholar] [CrossRef
[16] 张衡, 郑汉垣. 二维泊松方程边值问题有限差分方程的病态结构和最优预条件子[J]. 福建师范大学福清分校学报, 2018(5): 1-6.
[17] 张衡. 一维问题m次lagrange形函数有限元方程的条件数与预处理[J]. 石河子大学学报(自然科学版), 2017, 35(4): 518-521.
[18] 张衡. 两点边值问题2次lagrange形函数有限元方程的条件数和预处理[J]. 福州大学学报(自然科学版), 2017, 45(5): 617-622.
[19] 张衡. 两点边值问题3次lagrange形函数有限元方程的条件数和预处理[J]. 计算力学学报, 2017, 34(5): 672-676.
[20] 张衡, 郑汉垣. 两点边值问题网格方程病态机理和预处理[J]. 福州大学学报(自然科学版), 2019, 47(3): 296-299, 306.
[21] 吕涛, 石济民, 林振宝. 区域分解算法——偏微分方程数值解新技术[M]. 北京: 科学出版社, 1997: 165-167.
[22] 林群. 微分方程数值解法基础教程[M]. 第2版. 北京: 科学出版社, 2003: 87-111.
[23] 李开泰, 黄艾香, 黄庆怀. 有限元方法及其应用[M]. 北京: 科学出版社, 2006: 20-27.
[24] 朱伯芳. 有限单元法原理与应用[M]. 第4版. 北京: 中国水利水电出版社, 2018: 120-121.