A meshless data-tailored approach to compute statistics from scattered data with adaptive radial basis functions

TL;DR

This paper presents an advanced meshless RBF regression method that adaptively incorporates gradient information to improve velocity field reconstruction from scattered data, especially in regions with sharp gradients or anisotropic flows.

Contribution

It introduces an anisotropic, gradient-informed, and adaptive RBF framework that enhances accuracy and efficiency in reconstructing flow fields from scattered measurements.

Findings

Outperforms isotropic RBF in accuracy and smoothness

Reduces basis number by an order of magnitude

Improves reconstruction near shear layers and boundaries

Abstract

Constrained radial basis function (RBF) regression has recently emerged as a powerful meshless tool for reconstructing continuous velocity fields from scattered flow measurements, particularly in image-based velocimetry. However, existing formulations based on isotropic kernels often suffer from spurious oscillations in regions with sharp gradients or strong flow anisotropy. This work introduces an anisotropic, gradient-informed, and adaptively sampled extension of the constrained RBF framework for regression of scattered data. Gradient information is estimated via local polynomial regression at collocation points, smoothed, and used to (1) re-sample data, maximizing sampling density near steep gradients while downsampling in smooth regions, and (2) construct a local anisotropic metric that shapes each basis function according to the flow directionality. In addition, a gradient-informed regularization is introduced by embedding observed gradients into the least-squares system as weighted soft constraints. The resulting formulation is fully meshless, linear, and computationally efficient, while significantly improving reconstruction quality in challenging regions. The method is evaluated on both synthetic and experimental datasets, including direct numerical simulation (DNS) data of a turbulent channel and time-resolved particle tracking velocimetry of a turbulent jet. Results show that the proposed approach outperforms isotropic and gradient-free RBF formulations in accuracy, smoothness, and physical consistency -- particularly near shear layers and boundaries -- while reducing the number of bases by an order of magnitude. To support the application, we have created a repository (https://github.com/mendezVKI/SPICY_VKI) that provides access to the investigated datasets.