Construction of an iterative method for solving a nonlinear elliptic equation based on a mixed finite element method

Authors

• D.R. Baigereyev S. Amanzholov East Kazakhstan State University
• N.M. Temirbekov National Engineering Academy of the Republic of Kazakhstan
• D.A. Omariyeva D. Serikbayev East Kazakhstan State Technical University

80 54

Keywords:

mixed finite element method, nonlinear Poisson equation, a priori estimate, iterative method, Brezzi-Douglas-Marini elements

Abstract

This article is devoted to the construction and study of the finite element method for solving a two-dimensional nonlinear equation of elliptic type. Equations of this type arise in solving many applied problems, including problems of the theory of multiphase filtering, the theory of semiconductor devices, and many others. The relevance of the study of this problem is associated with the need to develop effective parallel methods for solving this problem. To discretize the equation, a mixed finite element method with Brezzi-Douglas-Marini elements is used. The issue of the convergence of the finite element method is investigated. To linearize the equation, the Picard iterative method is constructed. Two classes of basis functions of finite elements are used in the paper. A comparative analysis of the effectiveness of several direct and iterative methods for solving the resulting system of linear algebraic equations is carried out, including the method based on the Bunch-Kaufman LDLt factorization, the method of minimal residuals, the symmetric LQ method, the stabilized biconjugate gradient method, and a number of other iterative Krylov subspace algorithms with preconditioners based on incomplete LU decomposition. The method has been tested on several model problems by comparing an approximate solution with a known exact solution. The results of the analysis of the method error in various norms depending on the diameter of the mesh are presented.

References

[1] Gatica G., Baier R. and Tierra, G., "A mixed ﬁnite element method for Darcy’s equations with pressure dependent porosity" , Mathematics of Computation vol. 297 (2015): 1–33.

[2] Atkinson K. and Hansen O. "Solving the Nonlinear Poisson Equation on the Unit Disk" , Journal of Integral Equations
no. 3 (2005): 223-251.

[3] Auricchio F., Veiga L., Brezzi F. and Lovadina C. "Mixed ﬁnite element methods" , Encyclopedia of Computational Mechanics Second Edition no.1 (2017): 1–53.

[4] Puscas M. A., Enchery G. and Desroziers, S. "Application of the mixed multiscale ﬁnite element method to parallel simulations of two-phase ﬂows in porous media" , Oil and Gas Science and Technology vol. 73, no. 38 (2018): 1-14.

[5] Muzhinji K., Shateyi S. and Motsa S. "The Mixed Finite Element Multigrid Method for Stokes Equations" , The Scientific World Journal no. 460421 (2015): 1–12.

[6] Vorwerk J., Engwer C., Pursiainen S. and Wolters C. "A mixed ﬁnite element method to solve the EEG forward problem" ,
IEEE transations on medical imaging no. 4 (2016): 930–941.

[7] Spiridonov D., Huang J., Vasilyeva M., Huang Yu. and Chung, E. "Mixed generalized multiscale ﬁnite element method for Darcy-Forchheimer model" , Mathematics vol. 7, no. 1212 (2019): 1–13.

[8] Brown D. L. and Vasilyeva M. "Generalized multiscale ﬁnite element method for poroelasticity problems II: nonlinear coupling" , Journal of Computational and Applied Mathematics vol. 297 (2016): 132–146.

[9] Rebholz L., Viguerie A. and Xiao M. "Eﬃcient nonlinear iteration schemes based on algebraic splitting for the incom- pressible Navier-Stokes equations" , Mathematics of Computation vol. 88, no. 318 (2019): 1533–1557.

[10] Islam M., Hye A. and Mamun A. "Nonlinear Effects on the Convergence of Picard and Newton Iteration Methods in the Numerical Solution of One-Dimensional Variably Saturated–Unsaturated Flow Problems" , Hydrology vol. 4, no. 50 (2017): 1–18.

[11] Kuraz M., Mayer P. and Pech P. "Solving the nonlinear Richards equation model with adaptive domain decomposition" ,
Journal of Computational and Applied Mathematics vol. 270 (2014): 2-11.

[12] Nakshatrala K. and Turner D. "A mixed formulation for a modiﬁcation to Darcy equation based on Picard linearization and numerical solutions to large-scale realistic problems" , International Journal for Computational Methods in Engineering Science and Mechanics vol. 14, no. 6 (2013): 524–541.

[13] Madzvamuse A. and Chung A. "Fully implicit time-stepping schemes and non-linear solvers for systems of reaction– diﬀusion equations" , Applied Mathematics and Computation vol. 244 (2014): 361–374.

[14] Muccino J. and Luo H. "Picard iterations for a ﬁnite element shallow water equation model" , Ocean modeling vol. 10 (2005): 316–341.
[15] Zha Y., Yang J., Yin L., Zhang Y., Zheng, W. and Shi L. "A modiﬁed Picard iteration scheme for overcoming numerical diﬃculties of simulating inﬁltration into dry soil" , Journal of Hydrology vol. 551 (2017): 56–69.

[16] List F. and Radu F. "A study on iterative methods for solving Richards’ equation" , arXiv vol. 1507.07837v1 (2015): 1–16.

[17] Baboulin M., Dongarra J., Remy A., Tomov S. and Yamazaki I. "Solving dense symmetric indeﬁnite systems using GPU" ,
Concurrency and Computation vol. 29, no. 9 (2017): 1–17.

[18] Zhong-Zhi B. "Motivations and realizations of Krylov subspace methods for large sparse linear systems" , Journal of Computational and Applied Mathematics vol. 283, no. 1 (2015): 71–78.

[19] Tran H., Toh K. and Phoon K. "Preconditioned IDR(s) iterative solver for non-symmetric linear system associated with FEM analysis of shallow foundation" , International Journal for Numerical and Analytical Methods in Geomechanics vol. 37 (2013): 2972–2986.

[20] Choi S., Paige C. and Saunders M. "MINRES-QLP: A Krylov subspace method for indeﬁnite or singular symmetric systems" , SIAM Journal of Scientific Computing vol. 3, no. 4 (2011): 1810–1836.

[21] Chronopoulos A. T. and Kucherov A. "Block s-step Krylov iterative methods" , Numerical Linear Algebra with Applica- tions vol. 17, no. 1 (2010): 3–15.

[22] Cotter C. and Kirby R. "Mixed ﬁnite elements for global tide models" , Numerische Mathematik no. 133 (2015): 255–277.

[23] Chen Z. "Finite element methods and their applications"(Springer, 2007).

[24] Kozulik, V. and Gotovac, B. "Numerical solution of Poisson’s Equation in an arbitrary domain by using meshless R- function method" , Proceeding of the 27th DAAM International Symposium on Intellect Manufacturing and Automation (2016): 245–254.