Improving the accuracy of meshless methods via resolving power optimisation using multiple kernels

Meshless methods are commonly used to determine numerical solutions to partial differential equations (PDEs) for problems involving free surfaces and/or complex geometries, approximating spatial derivatives at collocation points via local kernels with a finite size. Despite their common use in turbulent flow simulations, the accuracy of meshless methods has typically been assessed using their convergence characteristics resulting from the polynomial consistency of approximations to operators, with little to no attention paid to the resolving power of the approximation. Here we provide a framework for the optimisation of resolving power by exploiting the non-uniqueness of kernels to provide improvements to numerical approximations of spatial derivatives. We first demonstrate that, unlike in finite-difference approximations, the resolving power of meshless methods is dependent not only on the magnitude of the wavenumber, but also its orientation, before using linear combinations of kernels to maximise resolving power over a range of wavenumbers. The new approach shows improved accuracy in convergence tests and has little impact on stability of time-dependent problems for a range of Eulerian meshless methods. Solutions to a variety of PDE systems are computed, with significant gains in accuracy for no extra computational cost per timestep in Eulerian frameworks. The improved resolution characteristics provided by the optimisation procedure presented herein enable accurate simulation of systems of PDEs whose solution contains short spatial scales such as flow fields with homogeneous isotropic turbulence.