Application of the Method of Two-Sided Approximations for Finding Positive Axially Symmetric Solutions of Boundary Value Problems with Monotone Nonlinearities

• Datum: 2024
• Hauptverfasser: Пархоменко, Владислав, Сидоров, Максим
• Format: Artikel
• Sprache: Ukrainisch
• Veröffentlicht: Кам'янець-Подільський національний університет імені Івана Огієнка 2024
• Online Zugang: http://mcm-math.kpnu.edu.ua/article/view/313371
• Назва журналу: Mathematical and computer modelling. Series: Physical and mathematical sciences

This paper considers the application of the method of two-sided approximations for finding positive solutions to boundary value problems for nonlinear elliptic differential equations with axial symmetry. The case of Dirichlet boundary conditions (first kind) is considered, with a power-type monotone nonlinearity, where the exponent ranges from 0 to 1. By transitioning to polar coordinates in the boundary value problem for the elliptic equation, due to the axial symmetry of the solution, the problem is reduced to a boundary value problem for an ordinary differential equation on an interval with respect to a function that depends only on the polar radius, thus eliminating the dependence on the angle. The pole of the polar coordinate system becomes a singular point, where a boundedness condition must be imposed. For the boundary value problem, a Green's function is constructed to further reduce the problem to a Hammerstein integral equation. The integral equation is treated as a nonlinear operator equation in a Banach space of continuous functions on the interval, partially ordered by a cone of non-negative continuous functions on this interval. The operator is examined for properties such as monotonicity (isotonicity), positivity, boundedness, and concavity. Next, the initial approximation is found as the endpoints of the invariant conical segment for the isotone operator in such a way as to ensure the highest convergence rate of the iterative process. The following iterative sequences of two-sided approximations are constructed: the first sequence, which is non-decreasing with respect to the cone, and the second sequence, which is non-increasing with respect to the cone. At each iteration, the arithmetic mean of the upper and lower approximations is taken as the next approximation. The iterative process is terminated when the error estimate of the solution satisfies the specified accuracy. The theoretical results obtained in this work were validated through a computational experiment. The dependence of the solution on the parameters in the right-hand side was analyzed and illustrated with corresponding graphs.