Thomas O. Gallouët

A JKO Splitting Scheme for Kantorovich--Fisher--Rao Gradient Flows

In this article we set up a splitting variant of the JKO scheme in order to handle gradient flows with respect to the Kantorovich-Fisher-Rao metric, recently introduced and defined on the space of positive Radon measure with varying masses. We perform successively a time step for the quadratic Wasserstein/Monge-Kantorovich distance, and then for the Hellinger/Fisher-Rao distance. Exploiting some inf-convolution structure of the metric we show convergence of the whole process for the standard class of energy functionals under suitable compactness assumptions, and investigate in details the case of internal energies. The interest is double: On the one hand we prove existence of weak solutions for a certain class of reaction-advection-diffusion equations, and on the other hand this process is constructive and well adapted to available numerical solvers.

Incompressible immiscible multiphase flows in porous media: a variational approach

We describe the competitive motion of (N + 1) incompressible immiscible phases within a porous medium as the gradient flow of a singular energy in the space of non-negative measures with prescribed mass endowed with some tensorial Wasserstein distance. We show the convergence of the approximation obtained by a minimization scheme a la [R. Jordan, D. Kinder-lehrer \& F. Otto, SIAM J. Math. Anal, 29(1):1--17, 1998]. This allow to obtain a new existence result for a physically well-established system of PDEs consisting in the Darcy-Muskat law for each phase, N capillary pressure relations, and a constraint on the volume occupied by the fluid. Our study does not require the introduction of any global or complementary pressure.

BLOW-UP PHENOMENA FOR GRADIENT FLOWS OF DISCRETE HOMOGENEOUS FUNCTIONALS

We investigate gradient flows of some homogeneous functionals in

The gradient flow structure for incompressible immiscible two-phase flows in porous media

\mathbb R^N

Second order models for optimal transport and cubic splines on the Wasserstein space

RN, arising in the Lagrangian approximation of systems of self-interacting and diffusing particles. We focus on the case of negative homogeneity. In the case of strong self-interaction (super critical case), the functional possesses a cone of negative energy. It is immediate to see that solutions with negative energy at some time become singular in finite time, meaning that a subset of particles concentrate at a single point. Here, we establish that all solutions become singular in finite time, in the super critical case, for the class of functionals under consideration. The paper is completed with numerical simulations illustrating the striking non linear dynamics when initial data have positive energy.

An unbalanced Optimal Transport splitting scheme for general advection-reaction-diffusion problems

We approximate the regular solutions of the incompressible Euler equations by the solution of ODEs on finite-dimensional spaces. Our approach combines Arnold’s interpretation of the solution of the Euler equations for incompressible and inviscid fluids as geodesics in the space of measure-preserving diffeomorphisms, and an extrinsic approximation of the equations of geodesics due to Brenier. Using recently developed semi-discrete optimal transport solvers, this approach yields a numerical scheme which is able to handle problems of realistic size in 2D. Our purpose in this article is to establish the convergence of this scheme towards regular solutions of the incompressible Euler equations, and to provide numerical experiments on a few simple test cases in 2D.