Homogenization of an optimal control problem for nonlocal semilinear elasticity with soft inclusions

Résumé Technique : Homogénéisation d'un Problème de Contrôle Optimal pour l'Élasticité Sémilinéaire Non Locale avec Inclusions Molles

1. Formulation du Problème
Cet article étudie l'analyse asymptotique d'un Problème de Contrôle Optimal (PCO) défini sur un milieu élastique à haut contraste contenant des inclusions périodiques molles. Le domaine physique $\Omega \subset \mathbb{R}^N$ est décomposé en une phase de matrice rigide connectée ($\Omega^2_\varepsilon$) et une phase d'inclusion molle ($\Omega^1_\varepsilon$). Le système est régi par une équation d'élasticité sémilinéaire présentant un terme non local, modélisant des matériaux dont la réponse mécanique dépend d'une mesure globale de l'amplitude de la déformation (par exemple, en raison de la redistribution de solvants ou d'un rétroaction de champ moyen).

L'équation d'état directrice pour le déplacement $u_{\varepsilon\delta}$ est :
$$-\nabla \cdot \sigma_{\varepsilon\delta}(u_{\varepsilon\delta}) + \alpha \|u_{\varepsilon\delta}\|_{L^p(\Omega)}^p u_{\varepsilon\delta} = f(x) + \mathbb{1}_{\Omega^1_\varepsilon} \theta_\varepsilon \quad \text{dans } \Omega,$$
soumise à des conditions aux limites de Dirichlet et de Neumann. Ici :

• $\varepsilon > 0$ représente l'échelle de périodicité.
• $\delta > 0$ représente le paramètre de contraste, où la rigidité dans la phase molle évolue en $\delta^2$ par rapport à la phase rigide (d'ordre 1).
• $\theta_\varepsilon \in L^2(\Omega^1_\varepsilon)^N$ est le contrôle distribué agissant exclusivement au sein des inclusions molles.
• Le terme non local $\alpha \|u\|_{L^p}^p u$ (avec $p \ge 2$) couple tous les points spatiaux à travers l'amplitude globale de la $L^2$-norme du déplacement.

L'objectif est de minimiser une fonctionnelle de coût réduite $j_{\varepsilon\delta}(\theta_\varepsilon)$, qui équilibre l'erreur de suivi $L^2$ du déplacement par rapport à un profil désiré $u_d$ et un terme de régularisation pour l'énergie de contrôle :
$$j_{\varepsilon\delta}(\theta_\varepsilon) = \frac{1}{2} \int_\Omega |u_{\varepsilon\delta}(\theta_\varepsilon) - u_d|^2 dx + \frac{\gamma}{2} \int_{\Omega^1_\varepsilon} |\theta_\varepsilon|^2 dx.$$
L'analyse se concentre sur la limite quand $(\varepsilon, \delta) \to (0, 0)$ sous le régime d'échelle où le rapport $\kappa = \lim \delta/\varepsilon$ appartient à $(0, +\infty]$.

2. Méthodologie
Les auteurs emploient une combinaison de l'opérateur de dépliage périodique (periodic unfolding operator) et d'arguments de $\Gamma$-convergence pour dériver les systèmes limites.

• Opérateur de Dépliage : L'opérateur de dépliage périodique $T_\varepsilon$ est utilisé pour traiter les domaines et les coefficients oscillants. Cette technique permet de séparer les échelles macroscopiques et microscopiques, transformant les intégrales sur $\Omega$ en intégrales sur $\Omega \times Y$ (où $Y$ est la cellule périodique unitaire).
• Décomposition de l'Énergie et Estimations A Priori : En raison de la nature à haut contraste (dégénérescence de la rigidité dans la phase molle), les estimations $H^1$ standards échouent. Les auteurs introduisent une décomposition spécifique du champ de déplacement en une partie macroscopique et une partie correctrice microscopique (Lemmes 3.2–3.3). Cela permet de dériver les estimations uniformes nécessaires à la compacité.
• Convergence Forte via les Termes Non Locaux : Un défi méthodologique critique est le passage à la limite dans le terme sémilinéaire non local $\alpha \|u\|_{L^p}^p u$. La convergence faible est insuffisante pour ce terme. Les auteurs établissent la convergence forte $L^2$ des variables d'état dépliées en exploitant la convergence de l'énergie totale (via la $\Gamma$-convergence) et la convexité stricte de la fonctionnelle d'énergie.
• Analyse du Système Adjoint : Pour caractériser le contrôle optimal, les auteurs dérivent l'équation de l'état adjoint. Ils prouvent l'existence d'un état adjoint unique et établissent son comportement asymptotique, ce qui est essentiel pour dériver les conditions d'optimalité limites.

3. Contributions Clés et Résultats

• Dérivation du Système d'État Homogénéisé :
  L'article dérive le système d'état limite (homogénéisé) associé à l'échelle $\delta/\varepsilon \to \kappa$. Le système limite présente une structure à deux échelles :

  • Un déplacement macroscopique $u_0$.
  • Des correcteurs microscopiques $\tilde{W}$ (associés à la phase molle) et $\hat{U}$ (associés à la phase rigide).
  • La fonctionnelle d'énergie limite $T_\kappa$ dépend explicitement du paramètre $\kappa$.
  • Dans le cas spécifique $\kappa = +\infty$, le système se découple en un système homogénéisé totalement séparé des échelles.
  • Le problème limite est caractérisé par une inégalité variationnelle impliquant le tenseur de rigidité homogénéisé $A^{hom}$ et le terme non local agissant sur le champ limite combiné $u_0 + \frac{1}{\kappa}\tilde{W}$.

• Homogénéisation du Problème de Contrôle Optimal :
  Les auteurs prouvent que la séquence des contrôles optimaux microscopiques $\Theta_\varepsilon$ converge vers un contrôle optimal $\hat{\Theta}$ pour le problème de contrôle optimal (PCO) limite.

  • Existence et Unicité : L'existence de contrôles optimaux pour les problèmes microscopique et limite est établie. L'unicité est garantie pour le problème limite sous certaines conditions, bien que la non-convexité de la fonctionnelle de coût réduite (due à la non-linéarité) implique que l'unicité globale pour le problème microscopique n'est pas garantie de manière générale.
  • Convergence : En utilisant la $\Gamma$-convergence des fonctionnelles de coût réduites et la convergence forte des variables d'état et de l'adjoint, il est prouvé que tout point limite d'une séquence de contrôles optimaux microscopiques est un contrôle optimal pour le PCO limite.
  • Conditions d'Optimalité Limites : Le contrôle optimal limite est explicitement caractérisé via l'état adjoint limite : $\hat{\Theta} = -\frac{1}{\gamma} (v_0 + \frac{1}{\kappa}\hat{Z})$, où $v_0$ et $\hat{Z}$ sont les composantes de l'état adjoint limite.

• Gestion du Haut Contraste et de la Non-Localité :
  Le travail traite de la difficulté simultanée de l'élasticité à haut contraste (où l'énergie dégénère dans les phases molles) et de la sémilinéarité non locale. Les auteurs démontrent que le terme non local nécessite une convergence forte, qui est obtenue grâce à des arguments de convergence d'énergie, distinguant cette analyse de l'homogénéisation sémilinéaire locale standard où la convergence faible suffit souvent.

4. Signification et Nouveauté
L'article revendique sa signification principalement dans la structure de la sémilinéarité non locale et le cadre de haut contraste.

• Structure de l'Opérateur Non Local : Contrairement aux travaux précédents sur les PCO avec des équations sémilinéaires (par exemple, ceux utilisant des termes locaux $S(u)$), ce modèle utilise un opérateur d'espace de Hilbert $u \mapsto \alpha \|u\|_{L^p}^p u$. Cela couple tous les points spatiaux de manière globale. Les auteurs notent que cette structure nécessite une convergence forte $L^2$ (convergence des normes) plutôt qu'une simple convergence faible, une exigence non triviale dans les contextes à haut contraste où les bornes $H^1$ uniformes font généralement défaut.
• Élasticité à Haut Contraste : L'analyse étend la théorie de l'homogénéisation aux problèmes de contrôle optimal en élasticité avec des inclusions molles, un régime où l'énergie élastique dégénère. L'article fournit un cadre rigoureux pour dériver les systèmes limites dans ce régime, montrant comment le paramètre $\kappa$ influence le modèle effectif.
• Rigueur Méthodologique : Le travail combine la méthode indirecte (systèmes d'optimalité) avec des arguments directs de $\Gamma$-convergence pour gérer la convergence de la application contrôle-vers-état et de la fonctionnelle de coût, garantissant que la limite des contrôles optimaux est bien optimale pour le problème limite.

L'article conclut en organisant les résultats en preuves d'existence pour l'état et le contrôle, des estimations a priori, la dérivation des systèmes limites à deux échelles, et enfin la convergence du PCO, fournissant un cadre mathématique complet pour l'optimisation de forme/réponse dans ces métamatériaux complexes.