A not-so-strange term coming from somewhere

Résumé technique : Homogénéisation de l'équation de Laplace dans des domaines périodiquement perforés avec des conditions de Robin échelonnées

Énoncé du problème
L'article étudie le comportement asymptotique des solutions de l'équation de Laplace (avec un terme de potentiel) dans un domaine borné $\Omega \subset \mathbb{R}^n$ ($n \ge 3$, avec des extensions à $n=2$) périodiquement perforé par un ensemble de trous $P_\varepsilon$. Les trous sont définis par $P_\varepsilon = \varepsilon \mathbb{Z}^n + \varepsilon H_\varepsilon$, où $H_\varepsilon$ représente un trou de référence redimensionné à l'intérieur de la cellule unité.

L'équation gouvernante est :
$$-\Delta u_\varepsilon + \alpha u_\varepsilon = f \quad \text{dans } \Omega_\varepsilon := \Omega \setminus P_\varepsilon$$
soumise à une condition aux limites de Robin sur les frontières des trous $\Gamma_\varepsilon = \partial P_\varepsilon \cap \Omega$ :
$$\partial_\nu u_\varepsilon + \frac{\beta}{\mu_\varepsilon} u_\varepsilon = \frac{g}{\mu_\varepsilon} \quad \text{sur } \Gamma_\varepsilon$$
et à des conditions de Dirichlet sur la frontière extérieure $\partial \Omega$.

La nouveauté fondamentale réside dans l'échelle du coefficient de Robin. Le paramètre $\mu_\varepsilon$ est choisi proportionnel à l'aire totale des perforations, $\sigma(\Gamma_\varepsilon)$. Plus précisément, $\mu_\varepsilon$ satisfait $\lim_{\varepsilon \to 0} \frac{|\Omega|\mu_\varepsilon}{\sigma(\Gamma_\varepsilon)} = 1$. Cette échelle garantit que le terme d'énergie de surface dans le fonctionnel variationnel associé est du même ordre que l'énergie volumique, créant un régime critique où les effets de surface et de volume entrent en compétition.

Méthodologie
Les auteurs emploient un cadre d'homogénéisation basé sur l'analyse des minimiseurs d'énergie et des inégalités de trace. La méthodologie procède par plusieurs étapes clés :

1. Formulation variationnelle : Le problème est formulé comme la minimisation d'un fonctionnel $F_\varepsilon(u)$ comprenant un terme d'énergie volumique (termes de gradient et de potentiel) et un terme d'énergie de surface échelonné par $1/\mu_\varepsilon$.
2. Cadre abstrait : Les auteurs établissent une théorie générale pour une famille de trous satisfaisant des hypothèses géométriques et fonctionnelles abstraites (C1–C5). Celles-ci incluent la régularité Lipschitzienne, la compacité et l'existence d'opérateurs d'extension uniformes.
3. Analyse spectrale (problèmes de type Steklov) : Un outil technique central est l'analyse d'un problème de valeurs propres spécifique (noté $EV_{\varepsilon, \kappa}$) sur la cellule unité perforée. Il s'agit d'un problème de type Steklov où le paramètre spectral $\lambda(\varepsilon, \kappa)$ apparaît à la fois dans l'équation différentielle et dans la condition aux limites :
   $$-\Delta v_{\varepsilon, \kappa} + \kappa \lambda(\varepsilon, \kappa) v_{\varepsilon, \kappa} = 0 \quad \text{dans } T^n \setminus H_\varepsilon$$
   $$\partial_\nu v_{\varepsilon, \kappa} = \lambda(\varepsilon, \kappa) v_{\varepsilon, \kappa} \quad \text{sur } \partial H_\varepsilon$$
   Les auteurs étudient la limite asymptotique des valeurs propres échelonnées $\lambda^*(\kappa) = \lim_{\varepsilon \to 0} \varepsilon^{-2}\lambda(\varepsilon, \kappa)$.
4. Inégalités de trace : L'analyse de convergence repose sur l'établissement d'inégalités de trace qui bornent l'intégrale de surface des fonctions sur $\Gamma_\varepsilon$ par leur énergie volumique. Les constantes optimales dans ces inégalités sont liées aux valeurs propres du problème spectral mentionné ci-dessus.
5. Construction de correcteurs : Pour prouver la convergence forte en $L^2$, les auteurs construisent des correcteurs $w_\varepsilon$ basés sur les minimiseurs du problème spectral. Ces correcteurs sont utilisés pour tester la formulation faible du problème et identifier l'équation limite.

Contributions et résultats clés

• Identification d'un régime critique : L'article identifie une échelle spécifique pour le paramètre de Robin (l'inverse de l'aire totale de surface) qui comble le fossé entre la limite de Neumann classique (couplage nul) et la limite de Dirichlet (inclusions infiniment rigides).
• Terme « étrange » non linéaire : Le résultat principal (Théorème 2.8) démontre que l'équation homogénéisée contient un terme supplémentaire d'ordre zéro :
  $$-\Delta u_0 + \alpha u_0 + \beta \kappa^*(\beta) u_0 = f + \kappa^*(\beta) g \quad \text{dans } \Omega$$
  Le coefficient $\beta \kappa^*(\beta)$ est le « terme étrange ». Contrairement aux résultats classiques où le terme étrange est linéaire ou constant, ici le coefficient dépend de manière non linéaire du paramètre de Robin $\beta$.
• Caractérisation via un problème spectral : Le paramètre $\kappa^*(\beta) \in (0, 1)$ est déterminé de manière unique par l'équation $-\lambda^*(\kappa^*(\beta)) = \beta$, où $\lambda^*(\kappa)$ est la limite des valeurs propres échelonnées du problème de type Steklov.
• Interpolation des limites :
  • Lorsque $\beta \to 0$ (limite de Neumann), $\beta \kappa^*(\beta) \to 0$, retrouvant le Laplacien standard sans terme étrange.
  • Lorsque $\beta \to +\infty$ (limite de Dirichlet), $\beta \kappa^*(\beta)$ converge vers le terme étrange capacitaires classique $\Lambda^*_{Cap}$ (terme de Cioranescu-Murat), retrouvant le résultat connu pour les trous de Dirichlet.
• Formules explicites pour des trous sphériques : Dans le cas modèle spécifique où les trous sont sphériques ($D = B^n$) et échelonnés comme $r_\varepsilon \sim \varepsilon^{n/(n-2)}$, les auteurs dérivent une formule explicite pour le terme étrange :
  $$\kappa^*(\beta) = \frac{\sigma_n(n-2)}{\sigma_n(n-2) + \beta}$$
  conduisant à un coefficient de terme étrange de $\frac{\beta \sigma_n(n-2)}{\sigma_n(n-2) + \beta}$.
• Vitesses de convergence : Sous des hypothèses de régularité supplémentaires (domaine $\Omega \in C^2$, termes sources bornés), l'article fournit une estimation d'erreur pour la convergence de $u_\varepsilon$ vers la solution homogénéisée $u_0$ (Proposition 2.9).
• Robustesse vis-à-vis de l'aire de surface : L'article inclut un exemple (Section 6) où l'aire de surface des trous $\sigma(\Gamma_\varepsilon)$ tend vers l'infini lorsque $\varepsilon \to 0$. Malgré la disparition du coefficient de Robin $\beta/\mu_\varepsilon$, la grande aire de surface compense, et un terme étrange non trivial persiste dans la limite.

Portée et affirmations
Les auteurs affirment que ce travail fournit un cadre unifié pour les problèmes d'homogénéisation sur des domaines perforés. En identifiant l'échelle critique où les énergies de surface et de volume sont équilibrées, ils montrent que le comportement effectif n'est pas simplement une transition entre les extrêmes de Neumann et de Dirichlet, mais implique une interaction distincte et non linéaire caractérisée par un problème spectral de type Steklov.

L'article souligne que le « terme étrange » résultant interpole continûment entre les deux régimes classiques. Les auteurs notent que bien que l'existence d'un tel terme soit connue pour les conditions de Dirichlet, la dépendance non linéaire spécifique au paramètre de Robin et sa caractérisation via un problème spectral où le paramètre apparaît à la fois dans l'équation et dans la condition aux limites constituent une contribution nouvelle. Le travail est motivé par des applications en mécanique des milieux continus, telles que la modélisation de matériaux composites avec des effets d'énergie de surface (par exemple, des colloïdes de cristaux liquides), mais la contribution principale est mathématique, offrant une analyse asymptotique rigoureuse et des vitesses de convergence pour cette échelle critique spécifique.