Discrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational Science2658-46702658-7149Peoples' Friendship University of Russia3032710.22363/2658-4670-2022-30-1-62-78Research ArticleFinite-difference methods for solving 1D Poisson problemNdayisengaSerge<p>Student of Department of Applied Probability and Informatics</p>1032195775@rudn.ruhttps://orcid.org/0000-0002-9297-9839SevastianovLeonid A.<p>Doctor of Physical and Mathematical Sciences, Professor of Department of Applied Probability and Informatics of Peoples’ Friendship University of Russia (RUDN University), Leading Researcher of Bogoliubov Laboratory of Theoretical Physics, JINR</p>sevastianov-la@rudn.ruhttps://orcid.org/0000-0002-1856-4643LovetskiyKonstantin P.<p>Candidate of Physical and Mathematical Sciences, Associate Professor of Department of Applied Probability and Informatics</p>lovetskiy-kp@rudn.ruhttps://orcid.org/0000-0002-3645-1060Peoples’ Friendship University of Russia (RUDN University)Bogoliubov Laboratory of Theoretical Physics Joint Institute for Nuclear Research01042022301627825022022Copyright © 2022, Ndayisenga S., Sevastianov L.A., Lovetskiy K.P.2022<p style="text-align: justify;">The paper discusses the formulation and analysis of methods for solving the one-dimensional Poisson equation based on finite-difference approximations - an important and very useful tool for the numerical study of differential equations. In fact, this is a classical approximation method based on the expansion of the solution in a Taylor series, based on which the recent progress of theoretical and practical studies allowed increasing the accuracy, stability, and convergence of methods for solving differential equations. Some of the features of this analysis include interesting extensions to classical numerical analysis of initial and boundary value problems. In the first part, a numerical method for solving the one-dimensional Poisson equation is presented, which reduces to solving a system of linear algebraic equations (SLAE) with a banded symmetric positive definite matrix. The well-known tridiagonal matrix algorithm, also known as the Thomas algorithm, is used to solve the SLAEs. The second part presents a solution method based on an analytical representation of the exact inverse matrix of a discretized version of the Poisson equation. Expressions for inverse matrices essentially depend on the types of boundary conditions in the original setting. Variants of inverse matrices for the Poisson equation with different boundary conditions at the ends of the interval under study are presented - the Dirichlet conditions at both ends of the interval, the Dirichlet conditions at one of the ends and Neumann conditions at the other. In all three cases, the coefficients of the inverse matrices are easily found and the algorithm for solving the problem is practically reduced to multiplying the matrix by the vector of the right-hand side.</p>1D Poisson equationfinite difference methodtridiagonal matrix inversionThomas algorithmGaussian elimination1D уравнение Пуассонаметод конечных разностейобращение трехдиагональной матрицыалгоритм Томасаисключение Гаусса<p>1. Introduction Applied mathematical models are mainly based on the use of partial differential equations [1]. The solution must satisfy a given equation of mathematical physics and some additional relations, which are, first, boundary and initial conditions. The most important for applications [2] are secondorder equations - elliptic, parabolic, and hyperbolic. Currently for equations of mathematical physics, methods of numerical solution and the appropriate software [3], [4], as well as computer algebra systems (CASs) such as Sage, Mathematica, Maxima and Maple are actively developed to implement these methods. Many features of stationary problems of mathematical physics described by elliptic equations of the second order can be illustrated by considering the simplest boundary value problems for an ordinary differential equation of the second order. Perhaps the simplest second-order elliptic equation is the Poisson equation. Let us consider some methods for the numerical solution of this equation and compare the investigated methods. The Poisson equation [1] is a special case of the heat conduction equation describing the dependence of the temperature of a medium on spatial coordinates and time, and the heat capacity and thermal conductivity of the medium (in the general case, inhomogeneous) are considered to be given. We will consider the problem of finding the steady-state distribution of density or temperature (e.g., when the distribution of sources does not depend on time). In this case, terms with time derivatives are eliminated from the non-stationary equation and a stationary heat equation is obtained, which belongs to the class of elliptic equations. A two-point boundary value problem is the problem of finding a solution to an ordinary differential equation or second-order systems in the interval</p>[A. N. Tikhonov and A. A. Samarskii, Equations of Mathematical Physics [Uravneniya matematicheskoy fiziki], 7th ed. Moscow: Moscow State University, Nauka, 2004, in Russian.][D. A. Yakovlev, V. G. Chigrinov, and H. S. Kwok, Modeling and optimization of LCD optical performance. New York: Wiley, 2015. DOI: 10.1002/9781118706749.][L. N. Trefethen, Approximation theory and approximation practice. Philadelphia: SIAM - Society for Industrial and Applied Mathematics, 2019.][M. Planitz et al., Numerical Recipes: The Art of Scientific Computing, 3rd ed. New York: Cambridge University Press, 2007.][N. N. Kalitkin and P. V. Koryakin, “Numerical methods [Chislennyye metody],” in Methods of Mathematical Physics, 1st ed. Moscow: Academia, 2013, vol. 2, in Russian.][A. A. Abramov and V. B. Andreyev, “On the application of the method of successive substitution to the determination of periodic solutions of differential and difference equations,” USSR Computational Mathematics and Mathematical Physics, vol. 3, no. 2, pp. 498-504, 1963. DOI: 10.1016/0041-5553(63)90034-x.][A. A. Samarskiy and A. V. Gulin, Numerical Methods of Mathematical Physics [Chislennyye metody matematicheskoy fiziki]. Moscow: Scientific world, 2003, in Russian.][L. H. Thomas, Elliptic problems in linear difference equations over a network. New York: Waston Sci. Comput. Lab. Rept., Columbia University, 1949.][S. B. Gueye, K. Talla, and C. Mbow, “Generalization of the exact solution of 1D Poisson equation with robin boundary conditions, using the finite difference method,” Journal of Electromagnetic Analysis and Applications, vol. 6, no. 12, pp. 372-381, 2014. DOI: 10.4236/jemaa.2014.612038.][S. B. Gueye, K. Talla, and C. Mbow, “Solution of 1D Poisson equation with Neumann-Dirichlet and Dirichlet-Neumann boundary conditions, using the finite difference method,” Journal of Electromagnetic Analysis and Applications, vol. 6, no. 10, pp. 309-318, 2014. DOI: 10.4236/jemaa.2014.610031.][S. B. Gueye, “The exact formulation of the inverse of the tridiagonal matrix for solving the 1D Poisson equation with the finite difference method,” Journal of Electromagnetic Analysis and Applications, vol. 6, no. 10, pp. 303-308, 2014. DOI: 10.4236/jemaa.2014.610031.][N. N. Kalitkin and E. A. Alshina, “Numerical Methods [Chislennyye metody],” in Numerical analysis. Moscow: Academia, 2013, vol. 1, in Russian.][A. Amosov, Y. Dubinsky, and N. Kopchenova, Computational Methods [Vychislitel’nyye metody], 4th ed. St. Petersburg: Lan’, 2021, in Russian.]