Attention : le séminaire aura lieu en salle Döblin.

Abstract: Dispersive nonlinear partial differential equations can be used to describe a range of physical systems, from water waves to spin states in ferromagnetism. The numerical approximation of solutions with limited differentiability (low-regularity) is crucial for simulating fascinating phenomena arising in these systems including emerging structures in random wave fields and dynamics of domain wall states, but it poses a significant challenge to classical algorithms. Recent years have seen the development of tailored low-regularity integrators to address this challenge. Inherited from their description of physicals systems many such dispersive nonlinear equations possess a rich geometric structure, such as a Hamiltonian formulation and conservation laws. To ensure that numerical schemes lead to meaningful results, it is vital to preserve this structure in numerical approximations. This, however, results in an interesting dichotomy: the rich theory of existent structure-preserving algorithms is typically limited to classical integrators that cannot reliably treat low-regularity phenomena, while most prior designs of low-regularity integrators break geometric structure in the equation. In this talk, we will outline recent advances incorporating structure-preserving properties into low-regularity integrators. Starting from simple discussions on the nonlinear Schrödinger and the Korteweg–de Vries equation we will discuss the construction of such schemes for a general class of dispersive equations before demonstrating an application to the simulation of low-regularity vortex filaments. This is joint work with Yvonne Alama Bronsard, Valeria Banica, Yvain Bruned and Katharina Schratz.

On s’intéressera dans cet exposé à l’étude d’objets présentant un caractère singulier et des opérations mal posées, l’objet central étant l’opérateur d’Anderson, un opérateur de Schrödinger où le potentiel est un bruit blanc espace. Après avoir discuté de sa construction et des propriétés qui en découlent, on établira un contrôle de sa fonction de Green de l’opérateur afin de comprendre sa singularité.
Dans un deuxième temps, on s’intéressera à des dynamiques dirigées par cet opérateur et particulièrement à la construction de mesures invariantes dans ce cadre. On s’attardera sur la définition de la mesure gaussienne associée et son utilisation pour comprendre le modèle ${\mathrm{\Phi }}_2^4$ dans un environnement régi par un bruit blanc espace. Après avoir construit la mesure invariante par des méthodes variationnelles, on s’intéressera aux propriétés qui découlent du point de vue dynamique : solutions globales, propriété de Feller et propriété de Feller forte.
Les résultat présentés sont basés sur des travaux en collaboration avec Antoine Mouzard et Tristan Robert.

On va montrer comment on peut définir le carré du module de la solution de l’équation de Schrödinger fractionnaire sur le tore avec condition initiale dans un espace de Sobolev arbitrairement singulier. Ensuite on va montrer comment cela peut être utile dans des estimations d’observabilité.

De nombreux travaux ont été consacrés à la dérivation d’approximations asymptotiques
des solutions d’une équation elliptique, lorsqu’on perturbe le milieu par des inhomogénéités
de petit volume. Les termes des développements asymptotiques des solutions contiennent
des informations sur la localisation, la forme et les propriétés physiques des inhomogénéités,
qui ont été exploitées avec bonheur dans le contexte des problèmes inverses d’identification
.
Dans cet exposé, nous présentons des résultats concernant le comportement des solutions
lorsque l’on perturbe la condition au bord sur un `petit’ ensemble ${\omega }_{\text{e}}$. Nous caractérisons
le terme de premier ordre du développement asymptotique en fonction de la mesure pertinente de
la taille de la perturbation ${\omega }_{\text{e}}$. Nous donnons des exemples explicites lorsque ${\omega }_{\text{e}}$
est une boule surfacique dans ${\text{R}}^d,d=2,3$.

Ce travail a été réalisé en collaboration avec Charles Dapogny et Michael Vogelius.

In this talk, I will briefly recall the historical Kalman filter and its ensemble form. Then I will show how the latter has been successfully implemented for data assimilation, in particular in numerical weather forecasting. More recently, the Ensemble Kalman Filter has been proposed as a methodology for solving very general inverse problems in high-dimensional contexts. I will present the theory, show some simple applications and point out the numerous open problems that remain.

Il est bien connu que la reconstruction d’une donnée initiale associée à une équation parabolique à partir de mesures internes de sa solution pendant un temps T, sur un domaine $\omega$ appelé domaine d’observation équivaut à la question de l’observabilité, ou plus précisément à la positivité de ce qu’on appelle la constante d’observabilité associée à $\omega$. Cette constante dépend du domaine d’observation $\omega$ mais aussi de façon cruciale de l’horizon temporel T.  

Dans cet exposé, nous nous intéressons au positionnement optimal de capteurs thermiques. Il est raisonnable de modéliser cette question apr la recherche des domaines extrémaux (lorsqu’ils existent) maximisant cette constante d’observabilité. Pour être physiquement pertinent, nous imposons une restriction sur la mesure du domaine observé. 

Après avoir introduit une relaxation convexe du problème d’optimisation de la forme, nous déterminons le comportement asymptotique des maximiseurs lorsque T tend vers $+\mathrm{\infty }$. En utilisant de façon cruciale un principe de la baignoire quantitatif, nous prouvons la forte convergence des maximiseurs vers la fonction caractéristique d’un ensemble mesurable que nous caractérisons précisément, et montrons en outre que cette convergence est exponentielle. 

Il s’agit d’un travail en collaboration avec Idriss Mazari (univ. Paris Dauphine) et Emmanuel Trélat (Sorbonne univ.)

! Attention ! Séminaire en visioconférence et à un horaire inhabituel.

Nous abordons l’approximation du flot de courbure moyenne des structures minces, pour lesquelles les méthodes classiques des champs de phase ne sont pas adaptées. Par structures minces, nous entendons soit des structures de codimension supérieure, typiquement des filaments, soit des surfaces non fermées et des surfaces non orientables.
Nous proposons une nouvelle approche qui consiste à introduire dans l’équation d’Allen-Cahn un terme de pénalisation localisé autour du squelette de l’ensemble en évolution. Cette approximation garantit une épaisseur minimale pendant l’évolution, prohibant ainsi les auto-intersections. L’efficacité numérique de notre approche est illustrée par des approximations du flot de courbure moyenne des filaments. Nous montrons son utilisation pour les approximations numériques aux problèmes de Steiner et de Plateau en dimension 3. Il s’agit d’un travail en collaboration avec Elie Bretin (INSA Lyon) et Simon Masnou (Lyon 1).

The mathematical modelling of cancer has been increasingly applying fluid-dynamics concepts to describe the mechanical properties of tissue growth. The biomechanical pressure plays a central role in these models, both as the driving force of cell movement and as an inhibitor of cell proliferation. In this talk, I will present how it is possible to build a bridge between models that have different pressure-velocity or pressure-density relations. In particular, I will focus on the inviscid limit from a Brinkman model to a porous medium-type model, and the incompressible limit that links the latter to a Hele-Shaw free boundary problem with density constraint.

Dans cet exposé, nous discutons de l’existence de solutions de norme L^2 prescrites pour des équations de Schrödinger non linéaires sur des graphes métriques. Une stratégie commune employée pour trouver une telle solution est de chercher un point critique sous contrainte de la fonctionnelle d’énergie associée. Certaines propriétés géométriques de la fonctionnelle varient en fonction de l’exposant du terme non linéaire de l’équation. Dans le cas dit de masse sous-critique, la fonctionnelle est bornée inférieurement et coercive sur la contrainte, de sorte que l’on peut rechercher un point critique en tant que minimum global. C’est pourquoi ce cas a été largement étudié ces dernières années.

Cependant, dans le cas complémentaire, connu sous le nom de masse sur-critique, la fonctionnelle d’énergie n’est plus bornée inférieurement sur la contrainte et présente un manque de d’estimation a priori sur les points critiques possibles. Par conséquent, on sait encore très peu de choses sur ce cas. A travers la présentation des quelques résultats existants, nous discuterons des principaux obstacles qui doivent être surmontés pour traiter ce cas sous des hypothèses générales. Nous présenterons également certains des outils qui ont déjà été développés à cette fin.

Cet exposé est basé sur des travaux communs avec J. Borthwick (Besançon puis Montréal), X. Chang (Changchun) et N. Soave (Turin).

We will discuss several models that might be regarded as migration models of populations moving from one part of a domain to the other and becoming part of the population living on the other side. Different situations assuming symmetry of movement between both sides of the domain, following a logistic model in their own environment and assuming spatial heterogeneities, are going to be discussed. Through such a common boundary both populations are coupled, acting as a permeable membrane on which their flow moves in and out. We will describe the precise interplay between the stationary solutions with respect to the parameters involved in the problem, in particular the growth rate of the populations and the coupling parameter involved on the boundary where the interchange of flux is taking place.

L’exposé se veut un résumé de mes travaux de thèse, qui portent une attention particulière aux schémas de Boltzmann sur réseau. Cette classe de schémas est utilisée depuis la fin des années ’80, en particulier en mécanique des fluides, et se caractérise par sa grande rapidité. Cependant, les méthodes de Boltzmann sur réseau sont très gourmandes en termes d’espace mémoire et conçues pour des maillages Cartésiens uniformes. De plus, nous manquons d’outils théoriques généraux qui permettent d’en analyser la consistance, la stabilité et enfin la convergence. Le travail s’articule autour de deux axes principaux. Le premier consiste à proposer une stratégie permettant d’appliquer les méthodes de Boltzmann sur réseau à des grilles de calcul non-uniformes adaptées dynamiquement en temps, afin de réduire le coût de calcul et de stockage. Le fait de pouvoir contrôler l’erreur commise et d’être en mesure d’employer la méthode quel que soit le schéma de Boltzmann sous-jacent sont des contraintes supplémentaires à prendre en compte. Pour cela, nous proposons d’adapter dynamiquement le réseau ainsi que d’ajuster toute méthode de Boltzmann à des maillages non-uniformes en nous appuyant sur la multirésolution. Cela a permis de proposer un cadre innovant pour des maillages mobiles en respectant les contraintes posées. Le second axe de recherche consiste à donner un cadre mathématiquement rigoureux aux méthodes de Boltzmann sur réseau, lié en particulier à leur consistance vis-à-vis des EDPs visées, leur stabilité et donc leur convergence. Pour cela, nous proposons une procédure, basée sur des résultats d’algèbre, pour éliminer les moments non-conservés de n’importe quel schéma de Boltzmann sur réseau, en le transformant en un schéma aux différences finies multi-pas sur les moments conservés. Les notions de consistance et stabilité pertinentes pour les méthodes de Boltzmann sur réseau sont donc celles des schémas aux différences finies. En particulier, tous les résultats concernant ces derniers, entre autres le théorème de Lax, se transpose naturellement aux schémas de Boltzmann sur réseau. Une étape ultérieure consiste à étudier la consistance et la stabilité directement sur le schéma de départ sans devoir calculer sa méthode aux différences finies « correspondante ». Cela permet d’en obtenir les équations modifiées et de montrer le bien-fondé des analyses de stabilité à la von Neumann couramment utilisées au sein de la communauté.

Depuis son introduction par Allen et Berkley en 1972, la méthode des
sources images est l’une des techniques les plus populaires pour la
modélisation des réponses impulsionnelles (RIR) en acoustique des
salles. Cette méthode modélise chaque réflexion d’une impulsion sonore
sur les murs d’une pièce rectangulaire (ou polyédrique) comme une source
impulsionnelle de type Dirac, obtenue à partir de critères géométriques
simples. Quelques travaux récents étudient l’estimation de la forme
d’une pièce tridimensionnelle en exploitant les temps d’arrivée des
échos dans l’enregistrement de la réponse impulsionnelle de salle.
Différentes limitations apparaissent dans ce type de méthode, notamment
la localisation temporelle des échos et leur labellisation. La méthode
présentée dans cet exposé permet la reconstruction des positions 3D des
sources images sans labellisation préalable des réflexions. Le problème
inverse est posé comme un problème convexe en dimension infinie de
reconstruction parcimonieuse en 3D des sources images, l’opérateur
linéaire d’observation à inverser faisant intervenir la solution de
l’équation des ondes avec un terme source mesure. Les dimensions d’une
pièce rectangulaire peuvent ensuite être estimées précisément à l’aide
du nuage de sources images ainsi reconstruites. L’exposé se conclura par
la présentation d’une approche alternative en cours de développement
basée sur l’optimisation de forme et la méthode des solutions
fondamentales, qui devrait permettre de dépasser le cas des pièces
rectangulaires.

The aim of this article is to analyze numerical schemes using two-layer neural networks with infinite width for the resolution of the high-dimensional Poisson partial differential equation (PDE) with Neumann boundary condition. Using Barron’s representation of the solution with a probability measure defined on the set of parameter values, the
energy is minimized thanks to a gradient curve dynamic on the 2-Wasserstein space of the set of parameter values defining the neural network. Inspired by the work from Bach and Chizat, we prove that if the gradient curve converges, then the represented function is the solution of the elliptic equation considered. In contrast to previous works, the activation function we use here is not assumed to be homogeneous to obtain global convergence of the flow. Numerical experiments are given to show the potential of the method.

On s’intéresse au problème d’optimisation de formes consistant à maximiser les valeurs propres du Laplacien avec conditions de Neumann homogènes. Ces valeurs propres interviennent notamment dans des problèmes acoustiques ou thermiques et sont en particulier liées à la « hot spot conjecture ». Contrairement aux valeurs propres de Dirichlet, celles associées au problème de Neumann sont de nature plutôt instable, ce qui rend le problème d’optimisation difficile. On verra comment certaines explorations numériques du problème pour des domaines du plan et de la sphère ont permis de mettre en évidence certaines propriétés des optima.

En fin de présentation, on fera une petite digression sur la capacité d’un réseau de neurones à apprendre les valeurs propres d’un opérateur.