Sur un ouvert $\Omega$ régulier, l’ensemble des fonctions lisses $C^{\infty}(\overline{\Omega})$ est dense dans les espaces de Sobolev $W^{1,p}(\Omega)$ (avec $1\leq p <\infty$). Pourtant, minimiser une fonctionnelle du calcul des variations sur $C^{\infty}(\overline{\Omega})$ ou sur $W^{1,p}(\Omega)$ peut conduire à des résultats différents: c’est le phénomène de Lavrentiev.

Il s’agit dans cet exposé d’identifier une large classe de fonctionnelles pour laquelle ce phénomène ne se produit pas. La preuve repose sur de nouvelles techniques d’approximation pour des versions paramétriques des problèmes variationnels considérés.

Dans cet exposé, nous considérons un opérateur non-auto adjoint sur un espace de Hilbert de la forme $H_0+V$ où $H_0$ est un opérateur auto-adjoint et $V$ est un opérateur borné à valeurs complexes. Nous supposons que la résolvante de $H_0$ satisfait un principe d’absorption limite et nous définissons les singularités spectrales de $H$ comme l’ensemble des points de son spectre essentiel tel que la résolvante de $H$ ne satisfait pas le principe d’absorption limite. Nous montrons alors que les singularités spectrales de $H$ sont en bijection avec des valeurs propres associées à des vecteurs propres spécifiques d’un prolongement de $H$ à un espace de Hilbert plus gros. Dans un deuxième temps, nous montrons que les états qui disparaissent à l’infini pour $H$ correspondent aux vecteurs propres généralisés de $H$ associés à des valeurs propres de partie imaginaire négative. Enfin nous définirons le sous-espace spectral absolument continu de $H$ et montrerons qu’il est égal à l’orthogonal de l’espace vectoriel engendré par tous les vecteurs propres généralisés de l’adjoint de $H$.

Des divergences apparaissent lorsque l’on cherche à effectuer des calculs en théorie des champs quantiques (comme par exemple pour la section efficace de diffusion). Il est donc essentiel de recourir à une procédure, nommée renormalisation, pour retirer ces divergences de nos modèles et obtenir des prédictions vérifiables expérimentalement.

En mathématiques, on peut décrire les modèles de théorie des champs par un Hamiltonien agissant sur un espace de Fock et dont le noyau n’est pas de carré intégrable. La procédure généralement suivie est d’abord d’introduire des coupures ultraviolettes permettant de régulariser le noyau. Ensuite, il faut démontrer que cet opérateur régularisé auquel on a soustrait un terme de compensation, converge vers un opérateur limite. Les termes de compensation généralement utilisés en mathématiques proviennent du premier ou des deux premiers ordres dans la théorie de la perturbation de grandeurs physiques telles que l’énergie de l’état fondamental. Dans cet exposé nous présenterons, sur un modèle jouet, une méthode permettant d’utiliser un nombre fini mais quelconque de termes de compensation, et donc d’ordre dans la théorie de la perturbation. Il s’agit d’un travail en cours en collaboration avec Jacob Møller.

Motivated by Tadmor’s work dedicated to multiscale image representation using hierarchical (BV,L^2) decompositions, we propose transposing their approach to the case of registration, task which consists in determining a smooth deformation aligning the salient constituents visible in an image into their counterpart in another. The underlying goal is to obtain a hierarchical decomposition of the deformation in the form of a composition of intermediate deformations: the coarser one, computed from versions of the two images capturing the essential features, encodes the main structural/geometrical deformation, while iterating the procedure and refining the versions of the two images yields more accurate deformations that map faithfully small-scale features. The proposed model falls within the framework of variational methods and hyperelasticity by viewing the shapes to be matched as Ogden materials. The material behaviour is described by means of a specifically tailored strain energy density function, complemented by L^∞ penalisations ensuring that the computed deformation is a bi-Lipschitz homeomorphism. Theoretical results emphasising the mathematical soundness of the model are provided, among which the existence of minimisers, a Γ-convergence result and an analysis of a suitable numerical algorithm, along with numerical simulations demonstrating the ability of the model to produce accurate hierarchical representations of deformations.

In this talk we present the problem of designing infinite-dimensional observers for infinite-dimensional systems. We consider the situation where only boundary measurement is available. Infinite-dimensional Luenberger-like observers are elaborated to infinite-dimensional dynamical systems in the abstract framework of semigroup systems. Exponential convergence of the proposed observers is guaranteed in the abstract framework by using the Lyapunov techniques. As examples explicit observers are worked out for PDE systems such as Euler-Bernoulli elastic beam systems and water wave systems. Thanks to observability property exponential or strong convergence of the designed observers is established with the convergence rate estimated in function of parameters of the observed systems, which is a desirable property for practical applications.

1. Lien entre l’hydrodynamique quantique et les équations de Schrödinger non linéaires
2. Vitesse de propagation maximale pour les équations de Schrödinger

The standard approach for photoacoustic imaging with variable speed of sound is time reversal, which consists of solving a well-posed final-boundary value problem for the wave equation backwards in time. We present a gradient based approach which consists of the iterative Landweber regularization algorithm, where convergence is guaranteed by standard regularization theory, notably also in cases of trapping sound speed or for short measurement times.
We formulate and solve the direct and inverse problem on the whole Euclidean space, which is common in standard photoacoustic imaging, but not for time reversal algorithms, where the problems are considered on a domain enclosed by the measurement devices. We formulate both the direct and adjoint photoacoustic operator as the solution of an interior and an exterior differential equation which are coupled by transmission conditions. The former is solved numerically using a Galerkin scheme in space and finite difference discretization in time, while the latter consists of solving a boundary integral equation. We therefore use a boundary element method/finite element method approach for numerical solution of the forward operators.
We analyze this method, prove convergence, and provide numerical tests. Moreover, we compare the approach to time reversal.

We consider the Dirichlet Laplacian in a two-dimensional strip composed of segments translated along a straight line with respect to a rotation angle with velocity diverging at infinity. We show that this model exhibits a “raise of dimension” at infinity leading to an essential spectrum determined by an asymptotic three-dimensional tube of annular cross section. If the cross section of the asymptotic tube is a disc, we also prove the existence of discrete eigenvalues below the essential spectrum. Joint work with David Krejcirik (Prague).

Le but de ce travail est de montrer la contrôlabilité à zéro d’un système modélisant l’interaction entre un fluide incompressible et une structure solide en utilisant un contrôle distribué localement situé dans le domaine du fluide. Le fluide est modélisé par le système de Navier-Stokes avec les conditions de Navier sur le bord, tandis que le corps rigide est gouverné par les lois de Newton. Le résultat principal montre qu’on peut amener les vitesses du fluide et de la structure à zéro et qu’on peut contrôler exactement la position du solide qui est supposé être assez régulier et de forme géométrique quelconque. Le point clé consiste à l’élaboration d’une nouvelle inégalité de Carleman pour le système linéarisé associé à notre problème couplant les équations de Stokes et des équations différentielles ordinaires avec les conditions de Navier sur le bord.

Nous appliquons une méthode de commutateurs à l’opérateur de Schrödinger discret sur . Nous proposons un nouvel opérateur conjugué qui nous permet de prouver la présence de spectre purement absolument continue à certains niveaux d’énergie. Nous imposons aux potentiels une décroissance de type longue portée plus générale que dans le traitement classique. Typiquement, on suppose une décroissance sur où κ>1 et τ_j est une translation dans la direction j. Le cas classique correspond à κ=1.

Travail réalisé en collaboration avec Marc-Adrien Mandich

On s’intéresse aux équations paraboliques linéaires, donc aux équations de la forme

(avec une condition initiale).

Nous nous intéressons principalement à l’obtention d’approximations de sa solution fondamentale K_t(x, y). Les approximations de la solution fondamentale ont une longue histoire en mathématiques et dans d’autres domaines, sous le nom d’approximation WKB (en théorie quantique) ou approximations du noyau de la chaleur.

L’asymptotique du noyau de la chaleur en termes de la géométrie a conduit à une démonstration célèbre du théorème de l’indice local et, plus récemment, dans les travaux de Pierre-Henry Labordère, à des applications souvent utilisées en finance. Si les coefficients sont suffisamment bornés (y compris leurs dérivées), j’expliquerai une méthode différente pour obtenir des développements en série perturbative qui conduit à des calculs très explicites des coefficients du développement asymptotique.

Cette méthode combine la formule de Duhamel/séries perturbatives de Dyson habituelles avec des calculs explicites inspirés des calculs des groupes de Lie nilpotents. J’expliquerai comment ces calculs peuvent ensuite être utilisés pour obtenir des approximations de la solution u de l’équation parabolique et je montrerai quelques tests numériques. J’espère que cette méthode pourra être étendue à de nombreux autres types d’équations (quelques résultats récents sur Navier-Stokes sur la sphère nous donnent plus de raisons d’espérer).

Ces techniques peuvent être combinées avec d’autres techniques (dans les EDP, dans l’analyse numérique ou dans d’autres domaines) et, en principe, pourraient fonctionner pour d’autres équations d’évolution. Dans le premier exposé, Nabil Kazi-Tani expliquera également très brièvement comment des équations paraboliques de ce type découlent de processus stochastiques.

On s’intéresse aux équations paraboliques linéaires, donc aux équations de la forme

(avec une condition initiale).

Nous nous intéressons principalement à l’obtention d’approximations de sa solution fondamentale K_t(x, y). Les approximations de la solution fondamentale ont une longue histoire en mathématiques et dans d’autres domaines, sous le nom d’approximation WKB (en théorie quantique) ou approximations du noyau de la chaleur.

L’asymptotique du noyau de la chaleur en termes de la géométrie a conduit à une démonstration célèbre du théorème de l’indice local et, plus récemment, dans les travaux de Pierre-Henry Labordère, à des applications souvent utilisées en finance. Si les coefficients sont suffisamment bornés (y compris leurs dérivées), j’expliquerai une méthode différente pour obtenir des développements en série perturbative qui conduit à des calculs très explicites des coefficients du développement asymptotique.

Cette méthode combine la formule de Duhamel/séries perturbatives de Dyson habituelles avec des calculs explicites inspirés des calculs des groupes de Lie nilpotents. J’expliquerai comment ces calculs peuvent ensuite être utilisés pour obtenir des approximations de la solution u de l’équation parabolique et je montrerai quelques tests numériques. J’espère que cette méthode pourra être étendue à de nombreux autres types d’équations (quelques résultats récents sur Navier-Stokes sur la sphère nous donnent plus de raisons d’espérer).

Ces techniques peuvent être combinées avec d’autres techniques (dans les EDP, dans l’analyse numérique ou dans d’autres domaines) et, en principe, pourraient fonctionner pour d’autres équations d’évolution. Dans le premier exposé, Nabil Kazi-Tani expliquera également très brièvement comment des équations paraboliques de ce type découlent de processus stochastiques.

The aim of this talk is to study a stochastic Schrödinger equation with a quadratic nonlinearity and a space-time fractional perturbation, in space dimension $dleq 3$. When the Hurst index is large enough, we will prove local well-posedness of the problem using classical arguments. I will briefly talk about the case where we deal with a small Hurst index since even the interpretation of the equation needs some care. A renormalization procedure must come into the picture, leading to a Wick-type interpretation of the model. Our fixed-point argument then involves some specific regularization properties of the Schrödinger group, which allows us to cope with the strong irregularity of the solution. This is a joint work with Aurélien Deya and Laurent Thomann.

We consider a nonlinear reaction–diffusion equation in a domain consisting of two bulk regions connected via small channels periodically distributed within a thin layer. The height and the thickness of the channels are of order $epsilon$, and the equation inside the layer depends on the parameter $epsilon$. We consider the critical scaling of the diffusion coefficients in the channels and nonlinear Neumann-boundary condition on the channels’ lateral boundaries. We derive effective models in the limit $epsilon to 0 $, when the channel-domain is replaced by an interface $Sigma$ between the two bulk-domains. Due to the critical size of the diffusion coefficients, we obtain jumps for the solution and its normal fluxes across $Sigma$, involving the solutions of local cell problems on the reference channel in every point of the interface $Sigma$. This is a joint work with Markus Gahn (University of Heidelberg)

Cet expose sera une introduction aux quelques résultats basiques pour les équations aux dérivées partielles sur domaine polyédrique. Par exemple, il est bien connu due aux travaux de Costabel, Dauge, Kondratiev, Mazya, Nicaise et d’autres que les solutions des équations elliptiques sur domaines non lisses souffrent d’une « perte de régularité. » Une question naturelle est qu’est-ce qu’on peut faire pour soulager ce problème, qui cause des ennuis pour les méthodes numériques. Les résultats que je vais exposer aident à  résoudre certains aspects de ce problème. L’exposé sera assez informel et basique.

Dans ce travail, on s’intéresse à  des configurations optimales de ressources (typiquement des denrées alimentaires) nécessaires à  la survie d’une espèce, dans un espace fermé. A cette fin, nous utilisons un modèle dit logistique pour décrire l’évolution de la densité d’individus constituant cette population. Cette équation fait intervenir une fonction représentant la répartition hétérogène (en espace) des ressources. La question principale traitée dans cet exposé peut se formuler ainsi : comment répartir de façon optimale des ressources dans un habitat ? Cette question peut se reformuler comme un problème de contrôle optimal ou d’optimisation de forme, dans lequel on cherche à  optimiser un critère mettant en jeu la valeur propre principale d’un opérateur par rapport au domaine occupé par les ressources ou encore une fonction de la densité de population dépendant implicitement d’un terme de ressources. Nous présenterons dans cet exposé de nouveaux résultats complétant une analyse ancienne de ces problèmes, ainsi que de nombreuses propriétés qualitatives.

We are used to self-adjoint Schrödinger operators with real potentials. In my talk I will try to convince you that the theory of 1 dimensional Schrödinger operators with complex potentials is very similar to the real case. For instance, the theory of their boundary value problem, formulas for their resolvents, etc. are essentially the same. The proofs are, however, often more difficult and one has to change the basic philosophy. For example, the concept of the self-adjointness should be replaced by the self-transposedness (called also the J-self-adjointness).

L’imagerie par résonance magnétique (IRM) est une technique d’imagerie médicale non-invasive et non-ionisante. L’examen IRM comporte plusieurs séquences d’acquisition, avec différents paramètres permettant de varier le contraste entre les tissus, et les données images sont acquises séquentiellement dans l’espace de Fourier. Lors du choix des paramètres, un compromis doit être choisi entre la résolution spatiale/temporelle des images, le rapport signal sur bruit, le type de contraste désiré et le temps d’acquisition. En imagerie thoracique et abdominale, les mouvements du patient (notamment cardio-respiratoires) imposent des contraintes additionnelles car ils peuvent dégrader la qualité des images (flous, artéfacts de type « fantômes ») et donc l’analyse, la quantification, et le diagnostic. Dans cet exposé nous présenterons quelques techniques de reconstruction des images permettant de lever certaines des limites rencontrées actuellement en imagerie clinique : la reconstruction conjointe de l’image et du mouvement du patient, la reconstruction super-résolution, ainsi que les algorithmes d’optimisation et techniques de régularisation associées.

Résumé

Nous étudions un système modélisant la dynamique d’une structure élastique immergée dans un fluide visqueux incompressible. Le mouvement du fluide est modélisé par les équations de Navier-Stokes et les déplacements élastiques suivent l’équation d’élasticité linéaire. Nous obtenons l’existence de solutions régulières locales en temps pour ce système.