The use of the line-by-line recurrent method for solving systems of difference elliptic equations with nine-diagonal matrices

A. A. Fomin ; L. N. Fomina

Vestnik Ûžno-Uralʹskogo gosudarstvennogo universiteta. Seriâ Vyčislitelʹnaâ matematika i informatika, Tome 8 (2019) no. 2, pp. 5-21

Résumé

The applying of the line-by-line recurrent method for solving systems of difference elliptic equations with nine-diagonal matrices is the subject of the article. Such matrices take place in the case of difference approximation of 2D differential problems of a higher order of accuracy on a regular grid covering the area under consideration. The technology of the so-called compensatory transform which allows replacing the initial nine-diagonal matrix of the system with the five-diagonal one is offered in the article, due to the fact that originally the line-by-line recurrent method was designed for solving systems of difference equations with a five-diagonal matrix. The efficiency of this technology is analyzed by comparing the solutions of the test boundary value problem in a unit square. The solutions are found both with the help of different implementations of the compensatory transform technology and by other modern highly efficient iterative methods for solving the systems of difference equations. The problem is solved on the sequence of grids from coarse (501$\times$501) to fine (4001$\times$4001) nodes. The accuracy of the solution convergence is determined by the relative norm of the residual, which is equal to $10^{-12}$ in the present work. It is shown that the line-by-line recurrent method retains its high efficiency over the entire range of the grids under consideration despite the use of the intermediate technology of the compensatory transform.

Détail

Keywords: grid method, system of difference elliptic equations, iterative method
Mots-clés : convergence of solution.

Bibliographie

[1] S. V. Patankar, Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing Corporation, New York, 1980, 197 pp.

[2] V. P. Il’in, Iterative Incomplete Factorization Methods, World Scientific Publishing Co., Singapore, 1992, 300 pp.

[3] B. P. Leonard, “A Stable and Accurate Convective Modelling Procedure Based on Quadratic Upstream Interpolation”, Computer Methods in Applied Mechanics and Engineering, 19:1 (1979), 59–98 | DOI

[4] P. H. Gaskell, A. K. C. Lau, “Curvature–Compensated Convective Transport: SMART, A New Boundedness – Preserving Transport Algorithm”, International Journal for Numerical Methods in Fluids, 8 (1988), 617–641 | DOI

[5] B. P. Leonard, “The ULTIMATE Conservative Difference Scheme Applied to Unsteady One-Dimensional Advection”, Computer Methods in Applied Mechanics and Engineering, 88:1 (1991), 17–74 | DOI

[6] M. S. Darwish, “A New High-Resolution Scheme Based on the Normalized Variable Formulation”, Numerical Heat Transfer, Part B: Fundamentals, 24 (1993), 353–371 | DOI

[7] M. S. Darwish, F. Moukalled, “The Normalized Weighting Factor Method: a Novel Technique for Accelerating the Convergence of High–Resolution Convective Schemes”, Numerical Heat Transfer, Part B: Fundamentals, 30 (1996), 217–237 | DOI

[8] D. V. Chirkov, S. G. Chernyi, “Comparison of Accuracy and Convergence of Some TVD-Schemes”, Computational Technologies, 5:5 (2000), 86–107

[9] A. I. Sukhinov, A. E. Chistakov, M. V. Yakobovskii, “Accuracy of the Numerical Solution of the Equations of Diffusion–Convection Using the Difference Schemes of Second and Fourth Order Approximation Error”, Bulletin of the South Ural State University. Series: Computational Mathematics and Software Engineering, 5:1 (2016), 47–62 | DOI

[10] L. A. Prokudina, N. M. Yaparova, M. P. Vikhirev, “Numerical Simulation of the Oscillations of the Elements of the Pipe with the Flow of an Incompressible Fluid”, Bulletin of the South Ural State University. Series: Computational Mathematics and Software Engineering, 7:3 (2018), 55–64 | DOI

[11] G. E. Schneider, M. Zedan, “A Modified Strongly Implicit Procedure for the Numerical Solution of Field Problems”, Numerical Heat Transfer, 4:1 (1981), 1–19 | DOI

[12] V. G. Zverev, “Modified Line-by-Line Method for Difference Elliptic Equations”, Computational Mathematics and Mathematical Physics, 38:9 (1998), 1490–1498

[13] D. F. Sikovskii, Computational Thermal Physics Methods, NSU, Novosibirsk, 2013, 98 pp.

[14] A. A. Fomin, L. N. Fomina, “On the Solution of Fluid Flow and Heat Transfer Problem in a 2D Channel with Backward–Facing Step”, Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta. Seriya: Fiziko-matematicheskie nauki, 21:2 (2017), 362–375 | DOI

[15] A. A. Fomin, L. N. Fomina, “On the Convergence of the Implicit Iterative Line–by–Line Recurrence Method for Solving Difference Elliptical Equations”, Computer Research and Modeling, 9:6 (2017), 857–880 | DOI

[16] M. S. Darwish, F. H. Moukalled, “Normalized Variable and Space Formulation Methodology for High-Resolution Schemes”, Numerical Heat Transfer, Part B: Fundamentals, 26 (1994), 79–96 | DOI

[17] H. A. Van der Vorst, “BI-CGSTAB: a Fast and Smoothly Converging Variant of BI-CG for the Solution of Nonsymmetric Linear Systems”, SIAM Journal on Scientific and Statistical Computing, 13 (1992), 631–644 | DOI

[18] A. V. Starchenko, “Comparative Analysis of Some Iterative Methods for the Numerical Solution of a Spatial Boundary Value Problem for Elliptic Equations”, Bulletin of the Tomsk State University. The Bulletin of Operational Scientific Information, 2003, no. 10, 70–80

[19] A. A. Fomin, L. N. Fomina, “Acceleration of the Line-by-Line Recurrent Method in Krylov Subspaces”, Tomsk State University Journal of Mathematics and Mechanics, 2011, no. 2, 45–54

[20] V. N. Faddeeva, Computational Methods of Linear Algebra, Dover Publications, N.Y., 1959, 252 pp.