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

TL;DR

This paper introduces a kernel optimization framework for meshless methods that enhances resolving power and accuracy in PDE solutions, especially for turbulent flows, without increasing computational cost.

Contribution

It develops a novel approach to optimize kernel combinations in meshless methods, improving resolving power and accuracy for PDEs involving complex geometries and turbulence.

Findings

Enhanced accuracy in convergence tests

Little impact on stability of time-dependent problems

Significant gains in resolving short spatial scales

Abstract

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.