Finite-difference methods for solving 1D Poisson problem

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.

В статье обсуждается постановка и анализ методов решения одномерного уравнения Пуассона на основе конечно-разностных аппроксимаций - важного и очень полезного инструмента численного исследования дифференциальных уравнений. По сути, это классический метод аппроксимации, основанный на разложении решения в ряд Тейлора. Развитие теоретических и практических результатов на базе этого метода в последние годы позволили повысить точность, стабильность и сходимость методов решения дифференциальных уравнений. Некоторые особенности этого анализа включают интересные расширения классического численного анализа начальных и граничных задач. В первой части излагается численный метод решения одномерного уравнения Пуассона, сводящийся к решению системы линейных алгебраических уравнений (СЛАУ) с ленточной симметричной положительно определённой матрицей. В качестве метода решения СЛАУ используется широко известный метод прогонки (метод Томаса). Во второй части представлен метод решения, основанный на аналитическом представлении точной обратной матрицы дискретизированного варианта уравнения Пуассона. Выражения для обратных матриц существенно зависят от типов граничных условий в исходной постановке. Представлены варианты обратных матриц для уравнения Пуассона с различными граничными условиями на концах исследуемого интервала - условиями Дирихле на обоих концах интервала, условиями Дирихле на одном из концов и Неймана на другом. Во всех трёх случаях коэффициенты обратных матриц легко вычисляются (выписываются) и алгоритм решения задачи практически сводится к умножению матрицы на вектор правой части.

Publisher
Федеральное государственное автономное образовательное учреждение высшего образования Российский университет дружбы народов (РУДН)
Number of issue
1
Language
English
Pages
62-78
Status
Published
Volume
30
Year
2022
Organizations
  • 1 Peoples’ Friendship University of Russia (RUDN University)
  • 2 Bogoliubov Laboratory of Theoretical Physics Joint Institute for Nuclear Research
Keywords
1D Poisson equation; finite difference method; Tridiagonal matrix inversion; Thomas algorithm; Gaussian elimination; 1D уравнение Пуассона; метод конечных разностей; обращение трехдиагональной матрицы; алгоритм Томаса; исключение Гаусса
Date of creation
06.07.2022
Date of change
06.07.2022
Short link
https://repository.rudn.ru/en/records/article/record/85021/
Share

Other records

Gambaryan M.M., Malykh M.D.
Discrete and Continuous Models and Applied Computational Science. Федеральное государственное автономное образовательное учреждение высшего образования Российский университет дружбы народов (РУДН). Vol. 30. 2022. P. 52-61
БАЛАН Е.Г.
Вестник Северного (Арктического) федерального университета. Серия: Гуманитарные и социальные науки. Федеральное государственное автономное образовательное учреждение высшего образования Северный (Арктический) федеральный университет имени М.В. Ломоносова. Vol. 22. 2022. P. 61-72