A supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids (Q853177)

The authors propose a finite difference algorithm for the Poisson equation that yields second order accuracy for the solutions and their gradients on non-graded grids. They employ quadtree (in 2D) and octree (in 3D) data structures as an efficient means to represent the Cartesian grid, allowing for constraint-free grid generation. The discretization at one cell's node only uses nodes of two (2D) or three (3D) adjacent cells, producing schemes that are straightforward to implement. The linear systems obtained are nonsymmetric but are shown to be diagonally dominant. Numerical results in 2D and 3D demonstrate supra-convergence in the \( L^\infty \) norm.