Bdx
Programme Venue

Q. Mérigot : Résolution numérique de l'équation d'Euler par transport optimal semi-discret


Résumé : Dans une première partie, je montrerai comment certains outils issus de la géométrie algorithmique permettent de résoudre numériquement le problème du transport optimal entre une mesure de probabilité absolument continue et une mesure à support fini (transport optimal semi-discret). Ensuite, je présenterai une méthode numérique développée avec Jean-Marie Mirebeau permettant de calculer des géodésiques minimisantes dans le groupe des applications préservant le volume. Ces géodésiques résolvent formellement l'équation d'Euler pour les fluides incompressibles (formulation d'Arnold). Notre méthode repose sur le théorème de décomposition polaire de Brenier, implémenté numériquement par transport optimal semi-discret. Cette méthode permet est suffisamment robuste pour extraire des solutions généralisées et multi-valuées de l'équation d'Euler, pour lesquelles la dimension du flot est supérieure à la dimension du domaine (ce qui est inévitable d'après Schnirelman). Nous montrons la convergence de cette méthode sous une hypothèse de régularité de la pression, et présentons des expériences numériques en 2D. (Travail commun avec Jean-Marie Mirebeau)