Hexagonal grid approximation of the solution of the heat equation on special polygons

We consider the first type boundary value problem of the heat equation in two space dimensions on special polygons with interior angles α ⱼπ αjπ , j=1,2,…,M j=1,2,…,M , where α ⱼ∈ \ 1/2,1/3,2/3 \ αj∈{12,13,23} . To approximate the solution we develop two difference problems on hexagonal grids using two layers with 14 points. It is proved that the given implicit schemes in both difference problems are unconditionally stable. It is also shown that the solutions of the constructed Difference Problem 1 and Difference Problem 2 converge to the exact solution on the grids of order O ( h²+τ ² ) O(h2+τ2) and O ( h⁴+τ ) O(h4+τ) respectively, where h and √32h 32h are the step sizes in space variables x₁ x1 and x₂ x2 respectively and τ is the step size in time. Furthermore, theoretical results are justified by numerical examples on a rectangle, trapezoid and parallelogram.