A meshless finite difference method for elliptic interface problems based on pivoted QR decomposition

We propose to solve elliptic interface problems by a meshless finite difference method, where the second order elliptic operator and jump conditions are discretized with the help of the QR decomposition of an appropriately rescaled multivariate Vandermonde matrix with partial pivoting. A prescribed consistency order is achieved on irregular nodes with small influence sets, which allows to place the nodes directly on the unfitted interface and leads to sparse system matrices with the density of nonzero entries comparable to the density of the system matrices arising from the mesh-based finite difference or finite element methods. Numerical experiments on a number of standard test problems with known solutions demonstrate convergence orders up to O(h(6)) for both the approximate solution and its gradient, and a robust performance of the method in the case when the interface is known inaccurately. (c) 2020 IMACS. Published by Elsevier B.V. All rights reserved.

Automated summary

The study proposes a new method to solve elliptic interface problems using a meshless finite difference approach, achieving prescribed consistency order on irregular nodes and placing nodes directly on unfitted interfaces, resulting in sparse system matrices. Experimental results show convergence orders up to O(h(6) for both approximate solution and gradient, as well as robust performance when the interface is known inaccurately.