Maximum Likelihood Filtering for Particle Tracking in Turbulent Flows
Lagrangian Particle Tracking (LPT) enables practitioners to study various concepts in turbulence by measuring particle positions in flows of interest. This data is subject to measurement errors, and filtering techniques are applied to mitigate these errors and improve the accuracy of analyses utilizing the data. We develop a new type of position filter through use of maximum likelihood estimation by considering both measurement errors and stochastic process physics. The maximum likelihood estimation scheme we develop is general and can accommodate many different stochastic process models, enabling it to be applied to many different turbulent flows. In this work, we propose a process model similar to existing, complimentary work in the development of B-splines. We compare our filtering scheme to existing schemes and find that our filter out performs the scheme proposed by Mordant et al. (2004) considerably, and produces similar performance to spline filters, proposed by Gesemann (2018). In comparing to the latter, we note that the maximum likelihood treatment provides a general framework which is capable of producing different filters based on the physics of interest, whereas the spline filters are built on less specific filtering theory and are therefore more difficult to adapt across diverse use cases in fluids. We quantify the performance of each of the filtering methods using error metrics which consider both noise reduction as well as signal degradation, and together these are used to define a concept of filter efficiency. The maximum likelihood filter developed in this work is shown to be the most efficient among all the methods examined when applied to simulated isotropic turbulence data from the Johns Hopkins Database.