From Linear to Nonlinear: A Resolvente criterion for Polynomial Stability of Semigroups Generated by Monotone Operators

1. Problématique et Contexte

La stabilité asymptotique des équations d'évolution dissipatives est un sujet central en analyse. Dans le cadre linéaire, le théorème de Borichev-Tomilov (2010) établit une correspondance précise entre la croissance du résolvant sur l'axe imaginaire et le taux de décroissance polynomiale de l'énergie. Plus précisément, pour un semi-groupe $C_0$ borné $e^{-tA}$, la décroissance $\|e^{-tA}A^{-1}\| \le C t^{-1/\alpha}$ est équivalente à la condition de résolvant $\|(isI + A)^{-1}\| \le C|s|^\alpha$ lorsque $|s| \to \infty$.

Cependant, dans le cadre non linéaire, cette approche échoue car les opérateurs non linéaires n'ont pas de spectre au sens classique, et l'analyse sur l'axe imaginaire (via la transformée de Fourier) n'est pas applicable. De plus, de nombreux systèmes physiques (dissipation non locale, amortissement dégénéré, équations de type Kirchhoff ou p-Laplacien) présentent des mécanismes de dissipation qui s'affaiblissent à l'équilibre, empêchant la stabilité exponentielle et conduisant à une décroissance polynomiale.

Le défi principal est donc de développer un critère de stabilité polynomial pour les opérateurs monotones maximaux non linéaires, sans recourir à la théorie spectrale linéaire, tout en évitant les difficultés techniques liées aux méthodes classiques de multiplicateurs (qui exigent souvent une régularité supérieure que les solutions faibles ne possèdent pas).

2. Méthodologie : Le Principe Taubérien Non Linéaire

Les auteurs proposent un changement de paradigme fondamental : remplacer l'analyse spectrale sur l'axe imaginaire ($i\mathbb{R}$) par une analyse asymptotique du résolvant réel près de l'origine ($\lambda \to 0^+$).

A. Le Cadre du Résolvant Réel

Pour un opérateur monotone maximal $A$ dans un espace de Hilbert $H$, on considère l'équation de résolvant stationnaire :
$$\lambda x_\lambda + A(x_\lambda) \ni y, \quad \lambda > 0$$
L'hypothèse centrale est que le comportement singulier (ou "explosion") de la norme $\|x_\lambda\|$ lorsque $\lambda \to 0^+$ révèle l'échelle de non-linéarité effective de l'opérateur.

B. Le Lien avec la Décroissance Polynomiale

L'article établit un principe Taubérien non linéaire :

1. Homogénéité : Si l'opérateur $A$ est $\alpha$-homogène (c'est-à-dire $A(\lambda u) = \lambda^\alpha A(u)$), la solution de l'équation de résolvant satisfait $\|x_\lambda\| \sim \lambda^{-\frac{1}{\alpha-1}}$.
2. Coercivité : Cette croissance du résolvant implique une inégalité de coercivité pour la dissipation : $\langle \xi, u \rangle \ge m \|u\|^{\alpha+1}$ pour tout $\xi \in A(u)$.
3. Résultat de décroissance : Cette coercivité conduit à une inégalité différentielle pour l'énergie $E(t)$ :
   $$E'(t) + c E(t)^{\frac{\alpha+1}{2}} \le 0$$
   Ce qui implique une décroissance polynomiale optimale :
   $$E(t) \le C(1+t)^{-\frac{2}{\alpha-1}}$$

C. Alignement Dissipatif

Une distinction cruciale est faite entre les systèmes où la dégénérescence de l'opérateur est alignée avec le canal de dissipation (ex: amortissement de Kelvin-Voigt non local dépendant de la vitesse) et ceux où elle ne l'est pas (ex: amortissement de Kirchhoff dépendant de l'énergie élastique).

• Aligné : Le critère du résolvant réel fonctionne directement.
• Non aligné : Le résolvant reste borné (pas d'explosion), mais la décroissance est lente. Dans ce cas, le théorème du résolvant seul est insuffisant et des arguments géométriques (multiplicateurs, GCC) dans le domaine temporel sont nécessaires.

3. Contributions Clés et Résultats

L'article propose dix applications concrètes démontrant la robustesse de cette approche :

1. Équation des ondes avec amortissement de Kelvin-Voigt non local :

   • Modèle : $u_{tt} - \Delta u - \|\nabla u_t\|^2 \Delta u_t = 0$.
   • Résultat : Le résolvant croît comme $\lambda^{-1}$, correspondant à une homogénéité $\alpha=3$. La méthode prouve une décroissance optimale en $O(t^{-1})$ pour des solutions faibles, là où les méthodes de multiplicateurs échoueraient par manque de régularité.

2. Généralisations et Dégénérescence :

   • Extension à des coefficients non locaux dépendant de puissances $\gamma$ de la norme du gradient. Le taux de décroissance est $O(t^{-1/\gamma})$.
   • Application à l'équation des ondes avec amortissement dégénéré sur la vitesse (Application III), confirmant le taux $O(t^{-1/\gamma})$.

3. Limites Structurelles (Contre-exemple Kirchhoff) :

   • Pour l'équation de Kirchhoff ($u_{tt} - \Delta u + \rho(\|\nabla u\|^2)u_t = 0$), la dégénérescence affecte la partie élastique, pas la vitesse. Le résolvant reste borné ($\gamma=0$), ce qui suggérerait une stabilité exponentielle si l'on suivait aveuglément la théorie linéaire.
   • L'article montre que la décroissance est en réalité polynomiale ($O(t^{-1/2\gamma})$) et que le critère du résolvant seul ne suffit pas ; il faut des arguments de contrôlabilité géométrique (GCC) dans le temps. Cela définit les limites précises de la méthode.

4. Équations Paraboliques et p-Laplacien :

   • Pour le p-Laplacien ($u_t - \Delta_p u = 0$), le résolvant est borné ($\gamma=0$), mais la décroissance est polynomiale $O(t^{-1/(p-2)})$ due uniquement à l'homogénéité de l'opérateur, et non à une obstruction dynamique.
   • Cas de dégénérescence intérieure (Application VIII) : Si le coefficient de diffusion s'annule sur un ensemble de mesure positive, le résolvant explose en $\lambda^{-1}$, indiquant l'absence totale de stabilité uniforme (sous-espaces invariants).

5. Mémoire et Couplages :

   • Application à des équations de viscoélasticité avec mémoire dépendante de l'état (Application X). La méthode sépare la partie linéaire (qui assure l'observabilité, $\gamma=0$) de la non-linéarité (qui dicte le taux polynomial).

6. Robustesse face aux termes sources :

   • L'article démontre que l'ajout de termes sources semi-linéaires (sous-critiques) ne modifie pas les exposants de décroissance polynomiale prédits par le critère du résolvant.

4. Signification et Impact

• Unification Théorique : L'article offre un cadre unifié ("Principe Taubérien Non Linéaire") qui généralise le théorème de Borichev-Tomilov au-delà du cadre linéaire. Il transforme un problème spectral (impossible en non linéaire) en un problème d'analyse asymptotique d'équations statiques.
• Avantage Technique Majeur : La méthode permet d'établir des taux de décroissance optimaux pour des solutions faibles (mild solutions) sans avoir besoin de régularité supérieure (comme $H^2$) requise par les méthodes classiques de multiplicateurs géométriques. Cela est crucial pour les problèmes dégénérés où la régularité n'est pas préservée.
• Diagnostic Structurel : Le critère du résolvant agit comme un outil de diagnostic. Il identifie non seulement les taux de décroissance, mais aussi les obstructions structurelles (comme les sous-espaces invariants dans les cas de dégénérescence totale) où la stabilisation est impossible.
• Perspective Physique : Il clarifie que la décroissance polynomiale n'est pas une propriété spectrale mystérieuse, mais la manifestation directe d'un mécanisme de coercivité non linéaire, codé dans le comportement du résolvant réel près de l'origine.

En conclusion, cet article fournit un outil mathématique rigoureux et prêt à l'emploi pour analyser la stabilité polynomiale d'une large classe de systèmes dissipatifs dégénérés, comblant le fossé entre la théorie spectrale linéaire mature et la complexité des équations aux dérivées partielles non linéaires.