On the fourth order accurate interpolation operator for the difference solution of the 3-dimensional Laplace equation

A. A. Dosiyev^a, E. Celiker<sup>b

a</sup> Department of Mechanics and Mathematics, Western Caspian University, Baku, Azerbaijan
^b University of Leicester, University Road, Leicester, UK

Abstract: A three-dimensional (3D) matching operator is proposed for a fourth-order accurate solution of a Dirichlet problem of Laplace's equation in a rectangular parallelepiped. The operator is constructed based on homogeneous, orthogonal-harmonic polynomials in three variables, and employs a cubic grid difference solution of the problem for the approximate solution inbetween the grid nodes. The difference solution on the nodes used by the interpolation operator is calculated by a novel formula, developed on the basis of the discrete Fourier transform. This formula can be applied on the required nodes directly, without requiring the solution of the whole system of difference equations. The fourth-order accuracy of the constructed numerical tools is demonstrated further through a numerical example.

Key words: 3D Laplace equation, cubic grids on parallelepiped, 15-point scheme, interpolation for harmonic functions, discrete Fourier transform.

MSC: 35A35, 35A40, 65D06, 65N15, 65N22

Received: 07.10.2023
Revised: 30.10.2023
Accepted: 19.11.2023

DOI: 10.15372/SJNM20240103

 English version: Numerical Analysis and Applications, 2024, 17:1, 28–42