$\Gamma$-convergence of convolution-type functionals for free discontinuity problems

1. Problématique et Contexte

L'article s'inscrit dans le domaine des problèmes de discontinuité libre, où l'on cherche à minimiser des fonctionnelles d'énergie qui dépendent à la fois du gradient d'une fonction et de la mesure de l'ensemble de ses sauts (discontinuités). L'exemple paradigmatique est la fonctionnelle de Mumford-Shah, utilisée en traitement d'images pour la segmentation :
$$\int_\Omega |\nabla u|^2 \, dx + \mathcal{H}^{d-1}(J_u)$$
où $J_u$ est l'ensemble des sauts de $u$.

Le défi principal réside dans la nature non-convexe et géométrique du terme de surface ($\mathcal{H}^{d-1}$), qui rend l'analyse numérique et théorique difficile. Une approche classique consiste à approximer ces fonctionnelles par des fonctionnelles non-locales plus régulières (souvent de type convolution), dont la minimisation est plus aisée.

L'objectif de ce papier est d'établir des résultats de compacité et de représentation intégrale pour une classe générale de fonctionnelles non-locales de type convolution, modélisées sur des travaux antérieurs (notamment ceux de Gobbino), mais avec des hypothèses de croissance et de noyaux plus générales.

2. Cadre Mathématique et Définitions

Les auteurs étudient une famille de fonctionnelles $\{F_\varepsilon\}_{\varepsilon > 0}$ définies sur $L^1_{loc}(\Omega; \mathbb{R}^m)$ par :
$$F_\varepsilon(u) := \frac{1}{\varepsilon^d} \int_\Omega \int_\Omega f_\varepsilon\left(x, \frac{y-x}{\varepsilon}, \frac{u(y)-u(x)}{\varepsilon}\right) \, dx \, dy$$
où les intégrandes $f_\varepsilon$ satisfont des conditions de croissance encadrées par des noyaux de convolution radiaux $\rho_1, \rho_2, \psi_\varepsilon$ et une fonction de perturbation $\eta_\varepsilon$.

Les conditions clés incluent :

• Bornes inférieures et supérieures : $f_\varepsilon$ est contrôlée par des termes quadratiques (liés à l'énergie de volume) et linéaires (liés à l'énergie de surface).
• Monotonie : $f_\varepsilon(x, \xi, z)$ est croissante par rapport à $|z|$.
• Comportement des noyaux : Les noyaux $\rho_i$ et $\psi_\varepsilon$ vérifient des conditions d'intégrabilité et de concentration (les interactions à courte distance dominent lorsque $\varepsilon \to 0$).

L'espace fonctionnel naturel pour la limite est l'espace des fonctions généralisées à variation bornée ($GSBV^2$), qui permet de traiter des fonctions dont le gradient est dans $L^2$ et dont l'ensemble de sauts a une mesure de Hausdorff finie.

3. Méthodologie

La preuve repose sur la théorie de la $\Gamma$-convergence dans la topologie $L^1_{loc}$. La stratégie se décompose en plusieurs étapes techniques majeures :

1. Réduction aux interactions à portée finie (Section 4) :
   Les auteurs montrent que le calcul de la $\Gamma$-limite de la fonctionnelle globale (avec des interactions à portée infinie) peut être obtenu en limitant d'abord l'analyse aux interactions à portée finie (troncature du noyau à une distance $T\varepsilon$), puis en faisant tendre $T \to +\infty$. Cela repose sur des estimations de type "extension" et des inégalités de contrôle pour les fonctions bornées.

2. Estimation fondamentale non-locale (Section 5) :
   Pour prouver la compacité, ils utilisent la méthode de localisation de $\Gamma$-convergence. Un obstacle majeur est l'absence de coercivité standard pour les fonctionnelles non-locales. Les auteurs établissent une estimation fondamentale (une version non-locale de l'estimation de De Giorgi-Letta) qui permet de "coller" deux fonctions minimisantes sur des domaines disjoints tout en contrôlant le coût énergétique de la transition. Cette étape utilise des arguments de troncature et de partition de l'unité.

3. Représentation Intégrale (Section 6) :
   Une fois la compacité établie, le but est de montrer que la limite $\Gamma$ admet une représentation intégrale de la forme :
   $$F(u) = \int_\Omega f_{bulk}(x, \nabla u) \, dx + \int_{J_u} f_{surf}(x, [u], \nu_u) \, d\mathcal{H}^{d-1}$$
   Les auteurs ne peuvent pas appliquer directement les théorèmes classiques (comme ceux de Bouchitté, Fonseca, Leoni, Mascarenhas) car leurs conditions de croissance inférieure sont plus faibles (absence de terme linéaire en $|[u]|$ dans la borne inférieure).
   Solution innovante : Ils utilisent une méthode de perturbation. Ils ajoutent un terme $\delta \int |[u]|$ à la fonctionnelle limite pour satisfaire les hypothèses fortes des théorèmes existants, obtiennent la représentation pour $\delta > 0$, puis font tendre $\delta \to 0$ en utilisant des techniques de convergence dominée et de monotonie.

4. Caractérisation des densités d'énergie (Section 7) :
   Ils établissent des conditions nécessaires et suffisantes pour la convergence de la suite entière. Les densités d'énergie $f_{bulk}$ et $f_{surf}$ sont caractérisées comme des limites de problèmes de minimisation sur de petits cubes (cellules) pour les fonctionnelles approchantes.

4. Résultats Principaux

• Théorème de Compacité (Théorème 3.3) : Toute suite de fonctionnelles $\{F_{\varepsilon_n}\}$ admet une sous-suite $\Gamma$-convergente vers une fonctionnelle $F$. Cette limite est définie sur $GSBV^2(\Omega; \mathbb{R}^m)$ et possède une structure de fonctionnelle de discontinuité libre avec des densités de volume et de surface.
• Représentation Intégrale : La fonctionnelle limite s'écrit explicitement comme la somme d'une énergie volumique (dépendant du gradient) et d'une énergie surfacique (dépendant de l'amplitude du saut $[u]$ et de la normale $\nu_u$).
• Caractérisation des Densités (Théorème 3.4) : Les auteurs donnent des formules explicites pour les densités $f_{bulk}$ et $f_{surf}$ en termes de limites de problèmes de minimisation sur des cubes de taille $r$ tendant vers 0. Ces densités sont déterminées par les comportements asymptotiques des énergies minimales sur des configurations de référence (fonctions linéaires pour le volume, fonctions à saut pour la surface).
• Convergence des Minimiseurs : La convergence des valeurs minimales des problèmes de Dirichlet associés est garantie, ce qui est crucial pour les applications numériques et l'homogénéisation.

5. Signification et Contributions

• Généralisation des modèles existants : Ce travail étend les résultats de Gobbino (1998) et d'autres travaux antérieurs en autorisant des noyaux de convolution plus généraux (non nécessairement compacts, avec des comportements asymptotiques spécifiques) et des conditions de croissance plus larges, incluant des termes linéaires dans l'énergie de surface.
• Traitement du cas vectoriel : Les résultats sont établis pour des fonctions à valeurs dans $\mathbb{R}^m$ (cas vectoriel), ce qui est essentiel pour des applications en élasticité ou en mécanique des milieux continus.
• Nouvelle technique de perturbation : La méthode utilisée pour contourner les difficultés de coercivité faible via une perturbation $\delta$ et un passage à la limite est une contribution technique notable à la théorie de la représentation intégrale.
• Applications potentielles : Ces résultats fournissent une base théorique solide pour l'homogénéisation de fonctionnelles de discontinuité libre non-locales et pour la justification rigoureuse de schémas numériques basés sur des approximations non-locales (comme les méthodes de type "phase-field" ou les approximations par différences finies non-locales).

En résumé, cet article fournit un cadre robuste et général pour l'analyse asymptotique d'une large classe d'énergies non-locales, prouvant qu'elles convergent vers des modèles classiques de discontinuité libre et offrant des outils précis pour calculer les paramètres effectifs de ces modèles limites.