An Algorithm of Linear Combinations: Thermal Conduction

This paper presents computational algorithms that make it possible to overcome some difficulties in the numerical solving boundary value problems of thermal conduction when the solution domain has a complex form or the boundary conditions differ from the standard ones. The boundary contours are assumed to be broken lines (the 2D case) or triangles (the 3D case). The boundary conditions and calculation results are presented as discrete functions whose values or averaged values are given at the geometric centers of the boundary elements. The boundary conditions can be imposed on the heat flows through the boundary elements as well as on the temperature, a linear combination of the temperature and the heat flow intensity both at the boundary of the solution domain and inside it. The solution to the boundary value problem is presented in the form of a linear combination of fundamental solutions of the Laplace equation and their partial derivatives, as well as any solutions of these equations that are regular in the solution domain, and the values of functions which can be calculated at the points of the boundary of the solution domain and at its internal points. If a solution included in the linear combination has a singularity at a boundary element, its average value over this boundary element is considered.