Interpolation based on radial basis functions (RBF) is a standard data map-
ping method used in multi-physics coupling. It works on scattered data without requiring
additional mesh topology or neighborhood information of support points. However, sys- tem matrices
of the equations for the coefficients tend to be ill-conditioned. In this work, we illustrate the
problem by a simple example and discuss possible remedies. Furthermore, we investigate the
numerical performance of this method on uniform and non-uniform meshes with a particular
focus on the coupling of black-box components where typically no information about the underlying
discretization can be extracted. Radial basis func- tion interpolation usually uses an
enhancement of the radial basis functions by a global polynomial in order to properly
capture constant components and linear trends in the given data. We present a method that
determines this polynomial independent from the radial basis function ansatz, which
substantially improves the condition number of the remaining RBF system. Furthermore, we
show that a rescaling approach can be used to either increase the accuracy or improve the
condition number even further by choosing radial basis functions with a smaller support
radius. The results represent an intermediate state with the aim to be integrated into the
multi-physics coupling library preCICE.
