Global solutions of the 3D inhomogeneous incompressible viscoelastic system without structure assumptions

1. Problème Étudié

L'article s'intéresse à l'existence globale de solutions fortes pour le système de viscoélasticité inhomogène et incompressible en dimension 3 ($\mathbb{R}^3$). Le système modélise un fluide dont la densité $\tilde{\rho}$ n'est pas constante, avec une vitesse $u$, un gradient de déformation $F$ et une pression $p$.

Les équations gouvernantes sont :
$$\begin{cases} \tilde{\rho}_t + u \cdot \nabla \tilde{\rho} = 0, \\ \tilde{\rho} u_t + \tilde{\rho} u \cdot \nabla u + \nabla p = \mu \Delta u + c^2 \text{div}(\tilde{\rho} F F^\top), \\ F_t + u \cdot \nabla F = \nabla u F, \\ \text{div } u = 0, \end{cases}$$
avec des conditions initiales $(\tilde{\rho}_0, u_0, F_0)$ tendant vers $(1, 0, I)$ à l'infini.

Le défi principal : La plupart des résultats antérieurs sur l'existence globale nécessitaient des hypothèses structurelles supplémentaires sur le tenseur de déformation $F$, telles que les identités de type « div-curl » (ex: $\text{div}(\tilde{\rho} F^\top) = 0$ ou $\text{curl}(F^{-1}-I)=0$). L'objectif de ce travail est de prouver l'existence globale sans aucune hypothèse structurelle supplémentaire sur $F$, en ne supposant que la petitesse des données initiales dans des espaces de Sobolev appropriés.

2. Méthodologie

Les auteurs développent une approche combinant des techniques de transformation innovantes et une analyse spectrale fine.

A. Reformulation du Système

Pour contourner la difficulté liée au terme de pression $\nabla p$ (qui ne peut pas être éliminé directement par un projecteur de Helmholtz car le coefficient de viscosité est variable en raison de la densité non constante), les auteurs introduisent :

1. Un tenseur effectif $G$ : Défini par $G := \tilde{\rho} F F^\top - I$. Cela permet de transformer le système original en un système dissipatif couplé pour la vitesse $u$ et le tenseur $G$.
2. Une transformation de pression : Ils définissent une nouvelle pression $\tilde{p}$ via $\nabla \tilde{p} := \frac{1}{\tilde{\rho}} \nabla p$, ce qui permet de réécrire le système sous une forme plus maniable pour l'estimation d'énergie.

Le système réécrit (système (2.5) dans l'article) met en évidence une structure dissipative pour le couple $(u, G)$.

B. Analyse Spectrale et Décroissance Temporelle

La méthode clé pour obtenir des bornes uniformes est l'analyse spectrale de la partie linéaire du système réécrit :

• Les auteurs étudient la matrice du système linéarisé dans l'espace de Fourier.
• Ils divisent l'espace des fréquences $\mathbb{R}^3$ en deux sous-domaines ($S_1$ pour les basses fréquences et $S_2$ pour les hautes fréquences) basés sur les valeurs propres de la matrice.
• Ils établissent des estimations précises pour les noyaux de la solution (via le principe de Duhamel), montrant une décroissance temporelle améliorée pour le champ de vitesse $u$.

C. Remplacement des Conditions $H^{-1}$ par $L^1$

Contrairement aux travaux précédents (comme ceux de Zhu et al.) qui utilisaient des conditions de petitesse dans l'espace $H^{-1}$, les auteurs imposent des conditions de petitesse sur la norme $L^1$ des données initiales ($u_0$ et $G_0$).

• Cette condition $L^1$ est cruciale pour obtenir une décroissance temporelle plus rapide de la vitesse ($\sim (1+t)^{-3/4}$), ce qui est nécessaire pour contrôler les termes non linéaires dans le cas inhomogène.

3. Résultats Principaux

Théorème 1.1 (Existence Globale) :
Sous l'hypothèse que les données initiales $(\tilde{\rho}_0, u_0, F_0)$ satisfont :

1. Une condition de petitesse dans $H^2$ : $\|u_0\|_{H^2} + \|\tilde{\rho}_0 - 1\|_{H^2} + \|G_0\|_{H^2} \leq \delta$ (avec $\delta$ suffisamment petit).
2. Une condition d'intégrabilité $L^1$ : $\|u_0\|_{L^1} + \|G_0\|_{L^1} = M_0 < +\infty$ (la constante $M_0$ n'a pas besoin d'être petite).

Alors, le problème de Cauchy admet une solution forte globale unique $(\tilde{\rho}, u, F)$ telle que :

• Les normes de Sobolev de la densité et du tenseur de déformation restent uniformément bornées : $\sup_{t \geq 0} \|(\tilde{\rho}-1, F-I)(t)\|_{H^2} \lesssim \delta$.
• La vitesse $u$ décroit avec un taux optimisé pondéré par le temps :
  $$\sup_{t \geq 0} (1+t)^\sigma \|u(t)\|_{H^2}^2 + \int_0^\infty (1+t)^\sigma \|\nabla u(t)\|_{H^2}^2 dt < \infty$$
  où $\sigma \in (1, \sigma_0)$ et $1 < \sigma_0 < 3/2$.

4. Contributions Clés et Innovations

1. Suppression des hypothèses structurelles : C'est la première preuve d'existence globale pour le système viscoélastique inhomogène 3D sans utiliser les identités « div-curl » habituelles sur $F$. Cela élargit considérablement la classe de solutions admissibles.
2. Nouvelle technique de transformation : L'introduction du tenseur effectif $G$ et la transformation de la pression permettent de gérer la variable coefficient dans l'équation de la vitesse, un obstacle majeur dans le cas inhomogène.
3. Optimisation des conditions initiales : Le passage de la condition $H^{-1}$ à la condition $L^1$ permet d'obtenir un taux de décroissance temporelle plus fort pour la vitesse, essentiel pour fermer les estimations a priori dans le cas où la densité varie.
4. Optimalité du taux de décroissance : Les auteurs montrent que l'intervalle pour l'exposant de poids temporel $\sigma$ est optimal ($\sigma < 3/2$), ce qui est lié à la dimension de l'espace et à la nature du système.

5. Signification et Impact

Ce travail représente une avancée significative dans la théorie mathématique des fluides complexes.

• Théorique : Il résout un problème ouvert concernant la stabilité des solutions pour les fluides viscoélastiques inhomogènes sans contraintes géométriques artificielles sur la déformation.
• Méthodologique : Il démontre la puissance de l'analyse spectrale combinée à des techniques de transformation de variables pour traiter des systèmes couplés non linéaires avec des coefficients variables.
• Application : Bien que purement mathématique, ces résultats renforcent la compréhension fondamentale du comportement à long terme des fluides complexes (polymères, gels) dont la densité peut varier, un cas plus réaliste que les modèles à densité constante.

En résumé, l'article établit l'existence globale de solutions petites pour un système complexe et non linéaire en éliminant des hypothèses restrictives précédemment considérées comme nécessaires, grâce à une analyse fine de la décroissance temporelle via l'espace $L^1$.