The problem of heat conduction in nonlinear media reduces to solving boundary value problems for nonlinear heat conduction equations, where the coefficients of the equation or the function of the power of heat sources depend on temperature according to some law. Among the numerical methods for solving problems for nonlinear equations of mathematical physics, one can distinguish finite difference methods, finite element methods, variational and projection methods, as well as iterative methods. Among the latter group, the two-sided approximation method is particularly attractive due to its ability to provide a convenient estimate for the error of the approximate solution and to prove the existence of a solution to the original problem. The theory of linear partially ordered spaces was developed by L. V. Kantorovich in the second half of the 1930s. Further development of this theory is associated with the works of M. A. Krasnosel’skii, H. Amann, V. I. Opojtsev, N. S. Kurpel, B. A. Shuvar, A. I. Kolosov etc. The aim of this article is to develop a two-sided approximation method based on the use of Green's functions for solving the first boundary value problem for a one-dimensional nonlinear heat conduction equation and to investigate its performance when solving test problems. To achieve this goal, the unknown function was replaced, and the boundary value problem was reduced to a Hammerstein integral equation, which was considered as a nonlinear operator equation in a partially ordered Banach space. Conditions for the existence of a unique positive solution to the problem and conditions for two-sided convergence of successive approximations to it were obtained. The developed method was implemented in software and tested on solving test problems. The results of the computational experiment are illustrated with graphical and tabular information
Read full abstract