In this article, we propose a mathematical model for one-dimensional heat conduction in a three-layered solid considering that an interfacial condition is present for the temperature and heat flux conditions between the layers. The numerical approach is developed by constructing a finite difference scheme to solve the initial boundary–interface problem. The numerical scheme is designed by considering the accuracy of the model on the inner part of each layer, then extending to the interfaces and boundaries by incorporating the continuous interfacial conditions. The finite difference scheme is unconditionally stable, convergent, and easy to implement since it consists of the solution of two algebraic systems. We provide three numerical examples to confirm that our numerical approximation is consistent with the analytical solution and the physical phenomenon.

In this paper, we consider the finite difference approximation for a one-dimensional mathematical model of heat conduction in a three-layered solid with interfacial conditions for temperature and heat flux between the layers. The finite difference scheme is unconditionally stable, convergent, and equivalent to the solution of two linear algebraic systems. We evaluate various methods for solving the involved linear systems by analyzing direct and iterative solvers, including GPU-accelerated approaches using CuPy and PyCUDA. We evaluate performance and scalability and contribute to advancing computational techniques for modeling complex physical processes accurately and efficiently.