SUMMARYAllowing discontinuous or non‐matching mesh spacing across zonal interfaces within a computational domain offers many advantages, particularly in terms of easing the mesh generation process, reduction of required mesh densities, and relative motion between mesh zones. This paper presents a numerical study of a universal method for interpolating solution data across such interfaces. The method utilises radial basis functions (RBFs) for n‐dimensional volume interpolation, and treats the available solution data points simply as arbitrary clouds of points, eliminating all connectivity requirements and making it applicable to a wide range of computational problems. Properties of the developed meshless interface interpolation are investigated using analytic functions, and three issues are considered: the achievable order of spatial accuracy of the RBF interpolation alone and comparison with a variable order polynomial; the effect of a combined RBF and polynomial interpolation; and the ability of the method to recover frequency content. RBF interpolation alone is shown to achieve fourth‐order to sixth‐order spatial accuracy in one and two dimensions, and in three dimensions, using a small number of data points, third‐order and above is achievable even for a 3 : 1 discontinuous cell spacing ratio, that is a 27 : 1 volume ratio, across the interface. Hence, it is inefficient to include polynomial terms, since improving on the RBF spatial accuracy results in a significant increase in the system size and deterioration in conditioning. It is also shown that only five points per wavelength are required to capture both frequency and amplitude content of periodic solutions to less than 0.01% error.Copyright © 2013 John Wiley & Sons, Ltd.