Transport

in Physics, Biology and Urban traffic

CIRM, Marseille, France
July 18 - August 26, 2022

Seminars

Seminars take place every morning at 9h in the auditorium. Tutorials take place at 14h in the auditorium.

Week 2 (July 25-29)

Monday. Mehdi Badsi (Université de Nantes), A kinetic model of plasma-probe interaction: theory and numerics. [Slides]
Tuesday. Aline Lefebvre-Lepot (CNRS, Ecole Polytechnique), Contacts: Numerical schemes based on convex optimization problems. [Abstract] [Slides]


In this presentation we focus on the numerical simulation of granular materials: collections of rigid, non-Brownian macroscopic grains (sand, grain, sugar, rubble, etc.). The contacts between grains lead to singular interactions for which adapted numerical schemes must be developed. We place ourselves in the framework of "Contact Dynamics" type models developed by J.J. Moreau, based on non-smooth convex analysis. The frictional models between grains lead to complex non-convex optimisation problems. The objective of the presentation is to show how algorithms based, at each instant, on convex optimization problems allow to obtain numerical long time simulations of large numbers of particles. In these schemes, the contact forces are obtained implicitly, as Lagrange multipliers associated with the constraints. The fundamental principle of dynamics is then obtained (in a discretised version) from the Euler (optimality) equations of the minimisation problem. In the case of frictionless grains, we come back to minimisation problems under affine constraints. In order to introduce friction (Coulomb model) between grains, we have to consider minimisation problems under conical constraint, for which the constraint is therefore non-derivable. We will illustrate this work with numerical simulations, showing that the resulting schemes can be used to study the macroscopic behaviour of granular materials. Part of the work was done in collaboration with Sylvain Faure, Philippe Gondret, Yvon Maday, Anne Mangeney, Hugo Martin, Bertrand Maury and Antoine Seguin.

Wednesday. Philippe Hoch (CEA DAM), High-order composite multi-dimensional finite volume schemes on unstructured meshes. [Abstract] [Slides]


In this talk, we consider the nodal numerical flux extension of classical Eulerian edge numerical flux schemes for linear and nonlinear systems of two dimensional hyperbolic conservation laws. In a first part, we present properties based both on mesh geometry and state interaction implying local conservativity and consistency. Such criteria apply to flux based schemes (VFFC, Roe,..) as well as viscosity based schemes (Rusanov, HLL, VF-Roe,..). This enables us to define composite edges/nodal fluxes in a natural way. In a second part, we deal with the stability (of discrete unknowns) and admissibility (in flux evaluation) of Euler system of gas dynamics. We focus on a nonlinear reconstruction for massic variables and propose a way to circumvent the direct velocity limitation which is not physically relevant. Hence, we show how to deduce from algebraic relation of total energy E = eps + 1/2 |U|^2 a limited velocity reconstruction induced by the direct limitation of internal massic energy eps.

Thursday. Vincent Calvez (CNRS, Institut Camille Jordan, Lyon), Biological waves at various scales. [Slides]
Friday. Nicolas Crouseilles (Inria Rennes Bretagne Atlantique), High order numerical methods for hybrid fluid-kinetic models. [Slides]

Tutorials
Wednesday (14h). Pierre Navaro (CNRS, Université de Rennes). Introduction to Git. [Slides]
Thursday (14h). Pierre Navaro (CNRS, Université de Rennes). Introduction to Julia. [Slides]

Week 3 (August 1-5)

Monday. Bastien Polizzi (Université de Franche-Comté), Mixture models for ecosystems. [Slides]
Tuesday. Charlotte Perrin (CNRS, Université de Marseille), Fluid flows under a congestion constraint.
Wednesday. Julien Olivier (Université de Marseille), On shear-banding in the Arrhenius model. [Abstract]


A story about rheology, and the stability of solutions of some system of PDEs. Featuring Fourier, Chapman-Enskog expansions and the Poincaré-Bendixon theorem.

Thursday. Yohan Penel (INRIA Paris), Further in the complexity of models? Application to the modelling of free-surface flows. [Slides]
Friday. Bertrand Maury (Université Paris Sud).

Week 4 (August 8-12)

Monday. Romain Ducasse (Université Paris Cité), Models in mathematical epidemiology and reaction-diffusion equations. [Abstract]


I will present in this talk some variations around a classical model in mathematical epidemiology, the (spatial) SIR model. This model is a rather simple system of parabolic equations. The goal of SIR-like models is to understand how a disease spreads in a population. I will start with a thorough analysis of the standard SIR system (long-time behavior of solutions, convergence, speed of propagation...). Then, we will present other models and the mathematical challenges they bear (models with delay, heterogeneous models.). Finally, I will show some recent results on a more complex model that takes into accounts variant (or strains, or mutations) of a disease.

Tuesday. Hervé Guillard (INRIA Sophia Antipolis), Singular limits of hyperbolic Partial Differential Equations in Fluid and Magnetohydro Dynamics
Wednesday. Boniface Nkonga (Université de Nice Sophia Antipolis), Hermite-Bézier Finite Element : Application to MHD in Tokamak plasmas.
Thursday. Mathieu Laurière (NYU Shanghai), Machine learning methods for mean field games and mean field control problems. [Abstract] [Slides]


In this talk, we will present stochastic numerical methods for mean field games and mean field control problems (also called optimal control of McKean-Vlasov dynamics). These problems arise as the limit of Nash equilibria or social optima in games when the number of players grows to infinity. We will mostly focus on methods based on neural networks to compute solutions when the model is fully known, motivated by applications in high dimension or with common noise. If time permits, we will also discuss model-free methods based on reinforcement learning.

Friday. Jean-Baptiste Saulnier (CNRS, LCB, Université Aix-Marseille). Collective behavior in Myxococcus xanthus bacteria.

Week 5 (August 15-19)

Tuesday. Erwan Hingant (Universidad del Bío-Bío, Chili), Equation de Becker-Döring I: modèle déterministe.
Wednesday. Benoît Gaudeul (Institut Camille Jordan, Université de Lyon), Finite volumes cookpots for the heat equation. [Slides]
Thursday. Romain Yvinec (INRAE, Université de Tours), Equation de Becker-Döring II: modèle stochastique.
Friday. Anaïs Crestetto (Université de Nantes), Numerical approximation of non classical solutions of Riemann problems. [Slides]

Week 6 (August 22-26)

Monday. Virginie Ehrlacher (Cermics, Ecole Nationale des Ponts et Chaussée), Hydrodynamic limits of multi-species mixtures: computation of the auto-diffusion matrix by tensor methods (joint work with Jad Dabaghi and Christoph Strössner). [Abstract] [Slides]


The modelisation of diffusion phenomena in solid crystalline alloys is very important in a wide variety of applications, one of which being the modelisation of the fabrication process of thin film solar cells. At the microscopic scale, atoms of various chemical species are located at the sites of the underlying crystalline lattice and may jump to another neighbouring site following a random Markov process, which is usually modeled as a symmetric multi-species exclusion process. In most realistic cases, the number of particles composing the system of interest is much too large for the movement of each atom to be fully simulated, and hydrodynamic limits (limits as the number of particles go to infinity) are thus considered instead. The latter read as cross-diffusion systems, i.e. deterministic systems modeling the evolution of the local concentrations or volumic fractions of each chemical species in the mixture. These systems read as nonlinear degenerate parabolic systems, which make their simulation very intricate. Another difficulty, which is the main focus of this talk, is that they involve the computation of some non-explicit terms, the so-called auto-diffusion matrix of the system, which can be expressed as the solution of a very high-dimensional optimization problem. Standard numerical approaches for the computation of this auto-diffusion matrix suffer from the so-called curse of dimensionality. As an alternative, we consider here a new numerical strategy, based on tensor decompositions, which proves to be highly effective in the present context. The aim of this talk is to present various theoretical and numerical results linked to this new approach.

Tuesday. Luca Nenna (Université Paris Saclay), On Optimal Transport, variational Mean Field Games and beyond. [Abstract] [Slides]


The minimization of a relative entropy (with respect to the Wiener measure) is a very old problem which dates back to Schrödinger. Strong connections and analogies between this problem and the Monge-Kantorovich problem with quadratic cost (namely the standard Optimal Transport problem) have been recently established. In particular, the entropic interpolation leads to a system of PDEs which presents strong analogies with the Mean Field Game (MFG) system with a quadratic Hamiltonian. In this talk, we will explain how such systems can indeed be obtained by minimization of a relative entropy at the level of measures on paths with an additional term involving the marginal in time. Moreover we will review all the links between Optimal transport and Variational MFG focusing both on theoretical and numerical aspects. If time permits we will also offer a glimpse on the case of multi-population MFG and its connection with some equations arising in quantum mechanics.