Higher order difference numerical analyses of a 2D Poisson equation by the interpolation finite difference method and calculation error evaluation

In a previous paper, a calculation system for a high-accuracy, high-speed calculation of a one-dimensional (1D) Poisson equation based on the interpolation finite difference method was shown. Spatial high-order finite difference (FD) schemes, including a usual second-order accurate centered space FD scheme, are instantaneously derived on the equally spaced/unequally spaced grid points based on the definition of the Lagrange polynomial function. The upper limit of the higher order FD scheme is not theoretically limited but is studied up to the tenth order, following the previous paper. In the numerical analyses of the 1D Poisson equation published in the previous paper, the FD scheme setting method, SAPI (m), m = 2, 4, ..., 10, was defined. Due to specifying the value of m, the setting of FD schemes is uniquely defined. This concept is extended to the numerical analysis of two-dimensional Poisson equations. In this paper, we focus on Poiseuille flows passing through arbitrary cross sections as numerical calculation examples. Over regular and irregular domains, three types of FD methods-(i) forward time explicit method, (ii) time marching successive displacement method, and (iii) alternative direction implicit method-are formulated, and their characteristics of convergence and numerical calculation errors are investigated. The numerical calculation system of the 2D Poisson equation formulated in this paper enables high-accuracy and high-speed calculation by the high-order difference in an arbitrary domain. Especially in the alternative direction implicit method using the band diagonal matrix algorithm, convergence is remarkably accelerated, and high-speed calculation becomes possible.