Two Stage Implicit Method on Hexagonal Grids for Approximating the First Derivatives of the Solution to the Heat Equation

Abstract

The first type of boundary value problem for the heat equation on a rectangle is considered. We propose a two stage implicit method for the approximation of the first order derivatives of the solution with respect to the spatial variables. To approximate the solution at the first stage, the unconditionally stable two layer implicit method on hexagonal grids given by Buranay and Arshad in 2020 is used which converges with O(h(2) + tau(2)) of accuracy on the grids. Here, h and root 3/2 h are the step sizes in space variables x(1) and x(2), respectively and t is the step size in time. At the second stage, we propose special difference boundary value problems on hexagonal grids for the approximation of first derivatives with respect to spatial variables of which the boundary conditions are defined by using the obtained solution from the first stage. It is proved that the given schemes in the difference problems are unconditionally stable. Further, for r = omega tau/h(2) <= 3/7, uniform convergence of the solution of the constructed special difference boundary value problems to the corresponding exact derivatives on hexagonal grids with order O(h(2) + tau(2)) is shown. Finally, the method is applied on a test problem and the numerical results are presented through tables and figures.