-
A meshless method to compute the proper orthogonal decomposition and its variants from scattered data
Authors:
Iacopo Tirelli,
Miguel Alfonso Mendez,
Andrea Ianiro,
Stefano Discetti
Abstract:
Complex phenomena can be better understood when broken down into a limited number of simpler "components". Linear statistical methods such as the principal component analysis and its variants are widely used across various fields of applied science to identify and rank these components based on the variance they represent in the data. These methods can be seen as factorizations of the matrix colle…
▽ More
Complex phenomena can be better understood when broken down into a limited number of simpler "components". Linear statistical methods such as the principal component analysis and its variants are widely used across various fields of applied science to identify and rank these components based on the variance they represent in the data. These methods can be seen as factorizations of the matrix collecting all the data, which are assumed to be a collection of time series sampled from fixed points in space. However, when data sampling locations vary over time, as with mobile monitoring stations in meteorology and oceanography or with particle tracking velocimetry in experimental fluid dynamics, advanced interpolation techniques are required to project the data onto a fixed grid before carrying out the factorization. This interpolation is often expensive and inaccurate. This work proposes a method to decompose scattered data without interpolating. The approach is based on physics-constrained radial basis function regression to compute inner products in space and time. The method provides an analytical and mesh-independent decomposition in space and time, demonstrating higher accuracy than the traditional approach. Our results show that it is possible to distill the most relevant "components" even for measurements whose natural output is a distribution of data scattered in space and time, maintaining high accuracy and mesh independence.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
A simple trick to improve the accuracy of PIV/PTV data
Authors:
Iacopo Tirelli,
Andrea Ianiro,
Stefano Discetti
Abstract:
Particle Image Velocimetry (PIV) estimates velocities through correlations of particle images within interrogation windows, leading to a spatial modulation of the velocity field. Although in principle Particle Tracking Velocimetry (PTV) estimates locally a non-modulated particle displacement, to exploit the scattered data from PTV it is necessary to interpolate these data on a structured grid, whi…
▽ More
Particle Image Velocimetry (PIV) estimates velocities through correlations of particle images within interrogation windows, leading to a spatial modulation of the velocity field. Although in principle Particle Tracking Velocimetry (PTV) estimates locally a non-modulated particle displacement, to exploit the scattered data from PTV it is necessary to interpolate these data on a structured grid, which implies a spatial modulation effect that biases the resulting velocity field. This systematic error due to finite spatial resolution inevitably depends on the interrogation window size and on the interparticle spacing. It must be observed that all these operations (cross-correlation, direct interpolation or averaging in windows) induce modulation on both the mean and the fluctuating part. We introduce a simple trick to reduce this systematic error source of PIV/PTV measurements exploiting ensemble statistics. Ensemble Particle Tracking Velocimetry (EPTV) can be leveraged to obtain the high-resolution mean flow by merging the different instantaneous realisations. The mean flow can be estimated with EPTV, and the fluctuating part can be measured from PIV/PTV. The high-resolution mean can then be superposed to the instantaneous fluctuating part to obtain velocity fields with lower systematic error. The methodology is validated against datasets with a progressively increasing level of complexity: two virtual experiments based on direct numerical simulations (DNS) of the wake of a fluidic pinball and a channel flow and the experimental data of a turbulent boundary layer. For all the cases both PTV and PIV are analysed.
△ Less
Submitted 11 February, 2023;
originally announced February 2023.
-
An end-to-end KNN-based PTV approach for high-resolution measurements and uncertainty quantification
Authors:
Iacopo Tirelli,
Andrea Ianiro,
Stefano Discetti
Abstract:
We introduce a novel end-to-end approach to improving the resolution of PIV measurements. The method blends information from different snapshots without the need for time-resolved measurements on grounds of similarity of flow regions in different snapshots. The main hypothesis is that, with a sufficiently large ensemble of statistically-independent snapshots, the identification of flow structures…
▽ More
We introduce a novel end-to-end approach to improving the resolution of PIV measurements. The method blends information from different snapshots without the need for time-resolved measurements on grounds of similarity of flow regions in different snapshots. The main hypothesis is that, with a sufficiently large ensemble of statistically-independent snapshots, the identification of flow structures that are morphologically similar but occurring at different time instants is feasible. Measured individual vectors from different snapshots with similar flow organisation can thus be merged, resulting in an artificially increased particle concentration. This allows to refine the interrogation region and, consequently, increase the spatial resolution. The measurement domain is split in subdomains. The similarity is enforced only on a local scale, i.e. morphologically-similar regions are sought only among subdomains corresponding to the same flow region. The identification of locally-similar snapshots is based on unsupervised K-nearest neighbours search in a space of significant flow features. Such features are defined in terms of a Proper Orthogonal Decomposition, performed in subdomains on the original low-resolution data, obtained either with standard cross-correlation or with binning of Particle Tracking Velocimetry data with a relatively large bin size. A refined bin size is then selected according to the number of "sufficiently close" snapshots identified. The statistical dispersion of the velocity vectors within the bin is then used to estimate the uncertainty and to select the optimal K which minimises it. The method is tested and validated against datasets with a progressively increasing level of complexity: two virtual experiments based on direct simulations of the wake of a fluidic pinball and a channel flow and the experimental data collected in a turbulent boundary layer.
△ Less
Submitted 2 September, 2022; v1 submitted 5 May, 2022;
originally announced May 2022.