Notions of optimal transport theory and how to implement them on a computer

This article gives an introduction to optimal transport, a mathematical theory that makes it possible to measure distances between functions (or distances between more general objects), to interpolate between objects or to enforce mass/volume conservation in certain computational physics simulations. Optimal transport is a rich scientific domain, with active research communities, both on its theoretical aspects and on more applicative considerations, such as geometry processing and machine learning. This article aims at explaining the main principles behind the theory of optimal transport, introduce the different involved notions, and more importantly, how they relate, to let the reader grasp an intuition of the elegant theory that structures them. Then we will consider a specific setting, called semi-discrete, where a continuous function is transported to a discrete sum of Dirac masses. Studying this specific setting naturally leads to an efficient computational algorithm, that uses classical notions of computational geometry, such as a generalization of Voronoi diagrams called Laguerre diagrams.Comment: 32 pages, 17 figure

Finite element simulation of three-dimensional free-surface flow problems

An adaptive finite element algorithm is described for the stable solution of three-dimensional free-surface-flow problems based primarily on the use of node movement. The algorithm also includes a discrete remeshing procedure which enhances its accuracy and robustness. The spatial discretisation allows an isoparametric piecewise-quadratic approximation of the domain geometry for accurate resolution of the curved free surface. The technique is illustrated through an implementation for surface-tension-dominated viscous flows modelled in terms of the Stokes equations with suitable boundary conditions on the deforming free surface. Two three-dimensional test problems are used to demonstrate the performance of the method: a liquid bridge problem and the formation of a fluid droplet