Coupled Vlasov and non-Newtonian fluid dynamics: existence and large-time behavior

Résumé Technique : Dynamique Couplée de Vlasov et de Fluides Non-Newtoniens

Formulation du Problème
Cet article étudie l'analyse mathématique d'un système couplé cinétique-fluide décrivant l'interaction entre un système de particules diluées et un fluide incompressible non-newtonien sur un domaine périodique $\mathbb{T}^3$. Le système est régi par l'équation de Vlasov pour la fonction de distribution des particules $f(t, x, v)$ et par une équation de Navier-Stokes incompressible modifiée pour la vitesse du fluide $u(t, x)$. Le tenseur des contraintes du fluide $\tau(Du)$ suit une relation de loi de puissance, $\tau(Du) = \mu|Du|^{p-2}Du$, où $p > 1$ est l'exposant de la loi de puissance. Ce modèle capture les comportements des fluides de Newton généralisés, incluant les phénomènes de rhéofluidification ($p < 2$) et de rhéoépaississement ($p > 2$). Le couplage s'opère par un terme de force de traînée dans l'équation de quantité de mouvement du fluide et par une force correspondante dans l'équation de transport des particules.

Les principaux défis analytiques abordés sont :

1. Existence Globale : Établir l'existence de solutions faibles globales dans le temps pour le système couplé, particulièrement dans le régime où la contrainte du fluide est purement non linéaire (dépourvue de composante newtonienne linéaire) et où l'exposant $p$ se situe dans la plage inférieure du régime de rhéofluidification.
2. Comportement à Temps Long : Quantifier les taux de décroissance du système vers un alignement de vitesse (état monocinétique) et identifier comment la non-linéarité de la viscosité du fluide influence ces taux.

Méthodologie
Les auteurs emploient une stratégie d'approximation multi-étapes et de compacité pour prouver l'existence globale, suivie d'arguments basés sur l'énergie pour les asymptotiques à temps long.

• Schéma d'Approximation : Pour construire les solutions faibles, les auteurs introduisent une régularisation à deux niveaux. Premièrement, ils mollifient les termes non linéaires (spécifiquement la force de traînée et les termes convectifs) et tronquent les données initiales dans l'espace des vitesses à l'aide d'un paramètre $\varepsilon > 0$. Deuxièmement, ils appliquent une projection de Galerkin sur un sous-espace de dimension finie de champs de vecteurs sans divergence.
• Argument de Point Fixe : Pour le système régularisé, les auteurs découplent les équations cinétiques et de fluide en figeant les grandeurs macroscopiques (densité et quantité de mouvement). Ils définissent une application de solution sur un espace de Banach et utilisent le théorème du point fixe de Schauder pour établir l'existence de solutions au système régularisé sur un intervalle de temps local.
• Compacité et Limites : L'existence globale est obtenue en passant à la limite en deux étapes :
  1. Limite de Galerkin ($n \to \infty$) : En utilisant des estimations d'énergie uniformes et des lemmes de moyennage de vitesse, les auteurs établissent une compacité forte pour les grandeurs macroscopiques (densité et quantité de mouvement) et la vitesse du fluide. Cela permet d'identifier la limite du terme de traînée non linéaire et du tenseur des contraintes.
  2. Limite de Mollification ($\varepsilon \to 0$) : En supprimant la mollification, ils récupèrent une solution faible du système original. Un aspect critique de cette étape est la gestion de l'identification du tenseur de contrainte non linéaire $\tau(Du)$ sans l'aide d'une dissipation linéaire. Les auteurs adaptent la méthode de troncature $L^\infty$ (méthode de capacité) combinée à des arguments de monotonicité pour prouver la convergence presque partout des gradients de vitesse.
• Analyse à Temps Long : Pour le comportement asymptotique, les auteurs définissent une fonctionnelle d'énergie modulée $E(t)$ mesurant l'écart du système par rapport à un état monocinétique. Ils dérivent une identité différentielle reliant la dérivée temporelle de $E(t)$ à une fonctionnelle de dissipation totale $D(t)$, qui inclut la viscosité du fluide et la dissipation de traînée. En établissant une inégalité de la forme $E(t) \lesssim D(t)^\gamma$, ils dérivent des inégalités différentielles qui produisent des taux de décroissance explicites.

Contributions Clés et Résultats

1. Existence Globale de Solutions Faibles :
   Le papier prouve l'existence de solutions faibles globales dans le temps pour le système couplé Vlasov-non-newtonien pour tous les exposants de la loi de puissance $p > 8/5$. Ce résultat étend les théories précédentes qui nécessitaient souvent $p \ge 11/5$ ou reposaient sur la présence d'un terme visqueux linéaire pour assurer la coercivité. Les auteurs traitent explicitement la difficulté de contrôler la moyenne spatiale de la vitesse du fluide dans le cadre périodique, où la quantité de mouvement totale est conservée mais où les moyennes individuelles peuvent évoluer. Ils y parviennent en décomposant la vitesse du fluide en une composante de moyenne nulle et une moyenne spatiale, permettant ainsi des estimations a priori raffinées.

2. Estimations de Décroissance à Temps Long :
   Sous une hypothèse supplémentaire que la densité macroscopique des particules reste uniformément bornée dans le temps, les auteurs établissent des estimations de décroissance quantitative pour l'énergie modulée $E(t)$. Le taux de décroissance est montré comme dépendant de manière critique de l'exposant $p$ :

   • Décroissance Algébrique : Pour $p > 2$ (régime de rhéoépaississement), l'énergie modulée décroît algébriquement.
   • Décroissance Exponentielle : Pour $6/5 \le p \le 2$ (régimes de rhéofluidification et de type newtonien), le système présente une décroissance exponentielle.
     Cette dichotomie souligne que, dans le régime purement non linéaire sans dissipation linéaire, la dynamique à temps long est gouvernée par l'interaction entre la diffusion non linéaire du fluide et la traînée cinétique. La décroissance exponentielle pour $p \le 2$ est attribuée aux effets dissipatifs plus forts du fluide dans cette plage, tandis que la décroissance algébrique pour $p > 2$ reflète une diffusion plus faible.

3. Convergence vers l'Alignement :
   En conséquence des estimations de décroissance, le papier démontre que les vitesses moyennes du fluide et des particules convergent vers une limite commune $v_\infty$ quand $t \to \infty$. De plus, la fonction de distribution des particules converge au sens de la distance de Lipschitz bornée vers un profil monocinétique $\rho_f(t, x) \otimes \delta_{v_\infty}(v)$.

Signification et Revendications
Les auteurs positionnent ce travail comme la première théorie d'existence globale dans le temps pour les solutions faibles de ce système spécifique couplé cinétique-fluide dans le domaine périodique avec une viscosité de loi de puissance purement non linéaire, couvrant spécifiquement la plage $8/5 < p < 11/5$. Ils notent que les résultats précédents reposaient typiquement sur une composante linéaire (newtonienne) dans le tenseur des contraintes pour garantir la coercivité $L^2$ et faciliter les arguments de compacité. En supprimant cette partie linéaire, l'article aborde un cadre analytique plus complexe où les estimations de coercivité standards impliquant $\|\nabla u\|_{L^p}$ ne peuvent être directement appliquées en raison de la conservation de la quantité de mouvement totale.

Concernant le comportement à temps long, le papier affirme combler une lacune dans la littérature où les études quantitatives existantes de taux de décroissance reposaient fondamentalement sur des termes visqueux linéaires pour dériver une convergence exponentielle. Ce travail constitue la première étude quantitative de la dynamique à temps long dans le régime purement non linéaire, révélant que la décroissance exponentielle n'est pas attendue de manière générale lorsque $p > 2$, et que la décroissance algébrique est le résultat naturel de l'interaction entre la diffusion non linéaire et la traînée. Les résultats sont confinés au domaine périodique $\mathbb{T}^3$, distinguant l'analyse des travaux sur des domaines bornés avec conditions de Dirichlet ou sur l'espace entier.