Optimal transport for particle image velocimetry
Martial Agueh, Boualem Khouider, Louis-Philippe Saumier · Communications in Mathematical Sciences · 2014
We present a new method for particle image velocimetry, a technique using successive laser images of particles immersed in a fluid to measure the velocity field of the fluid flow.The main idea is to recover this velocity field via the solution of the L 2 -optimal transport problem associated with each pair of successive distributions of tracers.We model the tracers by a network of Gaussian-like distributions and derive rigorous bounds on the approximation error in terms of the model's parameters.To obtain the numerical solution, we employ Newton's method, combined with an efficient spectral method, to solve the Monge-Ampère equation associated with the transport problem.We present numerical experiments based on two synthetic flow fields, a plane shear and an array of vortices.Although the theoretical results are derived for the case of a single particle in dimensions one and two, the results are valid in R d , d ≥ 1.Moreover, the numerical experiments demonstrate that these results hold for the case of multiple particles, provided the Monge-Ampère equation is solved on a fine enough grid.