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 Oh2+τ2 of accuracy on the grids. Here, h and 32h are the step sizes in space variables x1 and x2, respectively and τ 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=ωτh2≤37, uniform convergence of the solution of the constructed special difference boundary value problems to the corresponding exact derivatives on hexagonal grids with order Oh2+τ2 is shown. Finally, the method is applied on a test problem and the numerical results are presented through tables and figures.

We give fourth-order accurate implicit methods for the computation of the first-order spatial derivatives and second-order mixed derivatives involving the time derivative of the solution of first type boundary value problem of two dimensional heat equation. The methods are constructed based on two stages: At the first stage of the methods, the solution and its derivative with respect to time variable are approximated by using the implicit scheme in Buranay and Arshad in 2020. Therefore, Oh4+τ of convergence on constructed hexagonal grids is obtained that the step sizes in the space variables x1, x2 and in time variable are indicated by h, 32h and τ, respectively. Special difference boundary value problems on hexagonal grids are constructed at the second stages to approximate the first order spatial derivatives and the second order mixed derivatives of the solution. Further, Oh4+τ order of uniform convergence of these schemes are shown for r=ωτh2≥116, ω>0. Additionally, the methods are applied on two sample problems.