Zero-viscosity limit of the chemotaxis-Navier-Stokes equations with the Navier-slip boundary condition

Résumé technique : Limite de viscosité nulle des équations de chemotaxie-Navier-Stokes avec la condition aux limites de Navier

Énoncé du problème
Cet article examine la justification mathématique rigoureuse de la limite de viscosité nulle pour un système couplé de chemotaxie-Navier-Stokes dans un demi-espace bidimensionnel ($\mathbb{R}^2_+$). Le système modélise l'interaction entre la densité bactérienne ($n_\varepsilon$), la concentration en oxygène ($c_\varepsilon$) et la vitesse du fluide ($u_\varepsilon$), où les coefficients de viscosité et de diffusion chimique sont mis à l'échelle par un petit paramètre $\varepsilon$. L'étude se concentre sur le régime où la viscosité du fluide et la diffusion chimotactique s'annulent simultanément, un scénario motivé par des observations expérimentales de panaches bioconvectifs et de micro-turbulence dans des suspensions bactériennes (par exemple, Dombrowski et al. [3], Tuval et al. [39]).

Un défi central dans cette limite est la formation de couches limites due à la discordance entre les conditions aux limites de type non-glissement (ou glissement de Navier) du système visqueux et les conditions de glissement libre du système limite non visqueux. Alors que la limite de viscosité nulle pour les équations de Navier-Stokes pures a fait l'objet d'études approfondies (par exemple, la théorie de Prandtl, les travaux d'Oleinik, Sammartino, Caflisch, Masmoudi), le système couplé de chemotaxie et de fluide présente des singularités supplémentaires dues au couplage non linéaire fort entre la dynamique des fluides et l'agrégation chimotactique.

Méthodologie
Les auteurs emploient la méthode du développement asymptotique apparié pour construire des solutions approchées et analyser rigoureusement le taux de convergence. Les composantes méthodologiques centrales incluent :

1. Développement asymptotique : La solution est décomposée en une couche externe (non visqueuse) et une correction de couche limite interne. Les auteurs introduisent des coordonnées de couche limite $z = y/\varepsilon$ et développent les variables ($n, c, u, p$) en puissances de $\varepsilon$.

   • Le développement inclut des termes jusqu'aux corrections de couche limite d'ordre deux ($nb,2, cb,2, ub,2$).
   • Une observation structurelle clé est que les équations pour les termes de couche limite d'ordre supérieur (par exemple, $nb,j$) dégénèrent souvent, perdant leurs dérivées temporelles pour devenir des relations linéaires impliquant uniquement des dérivées spatiales et des termes d'ordre inférieur connus. Cela permet aux auteurs de résoudre ces termes de manière séquentielle sans dépendance circulaire.

2. Espaces fonctionnels : Pour traiter la nature anisotrope de la couche limite (où les dérivées normales sont grandes), les auteurs utilisent des espaces de Sobolev conormaux anisotropes ($Y^{l,m}$). Ces espaces attribuent des poids différents aux dérivées tangentielles ($x$) et normales ($y$), en utilisant une fonction de poids $\psi(y)$ qui s'annule à la frontière. Ce cadre est adapté des techniques utilisées pour les équations de Navier-Stokes avec des conditions de glissement de Navier (Masmoudi et Rousset [26]) mais est appliqué ici de manière innovante au système de chemotaxie-fluide.

3. Analyse d'erreur : La différence entre la solution réelle et la solution approchée construite est définie comme les variables d'erreur $(n, c, u)$. Les auteurs dérivent un système d'équations d'erreur et établissent des estimations d'énergie pour ces variables.

   • Les estimations impliquent le contrôle des termes non linéaires, des commutateurs résultant des dérivées conormales et du couplage entre les équations de chemotaxie et de fluide.
   • Une attention particulière est portée au terme de pression et à l'équation de tourbillon pour retrouver les dérivées normales, en utilisant une variable de tourbillon modifiée $\tilde{\omega} = \omega - \varepsilon^2 \partial_x f e^{-y}$ pour gérer les effets de frontière.

4. Bien-posé : L'article établit le bien-posé local des équations de la couche externe (un système de chemotaxie-Euler) et des équations de la couche limite (un système d'équations paraboliques et de transport) dans des espaces de Sobolev appropriés, garantissant que la solution approchée est bien définie.

Contributions et résultats clés
Le résultat principal est la preuve rigoureuse de la limite de viscosité nulle avec un taux de convergence spécifique.

• Théorème principal (Théorème 2.2) : En supposant l'existence d'une solution approchée régulière satisfaisant des bornes spécifiques (Proposition 2.1), les auteurs prouvent que pour $\varepsilon$ suffisamment petit, le système de chemotaxie-Navier-Stokes admet une solution unique qui converge vers la solution approchée dans les espaces de Sobolev conormaux anisotropes.
• Taux de convergence : L'erreur entre la solution réelle et la solution approchée est bornée par $O(\varepsilon^{3/2})$ dans la norme $L^\infty(0, T; L^2 \cap L^\infty)$. Plus précisément :
  $$\|(n_\varepsilon - n_a, c_\varepsilon - c_a, u_\varepsilon - u_a)\|_{L^\infty_t Y^{1,2}} + \|\nabla(n_\varepsilon - n_a)\|_{L^2_t Y^{1,2}} + \dots \leq C \varepsilon^{3/2}$$
  Ce taux améliore les résultats précédents pour des systèmes similaires (par exemple, Hou [10] qui a obtenu $O(\varepsilon^{1/2})$ pour une configuration différente).
• Structure de la couche limite : L'article construit explicitement les profils de la couche limite. Il démontre que l'effet de couche limite dans ce système est induit par la dynamique des fluides plutôt que par le mécanisme chimotactique seul. En l'absence de couplage fluide ($u=0$), les corrections de couche limite pour les variables chimotactiques s'annulent identiquement.
• Gestion de la dégénérescence : Une contribution technique significative est la méthode utilisée pour résoudre les corrections de couche limite $nb,j$. En exploitant la dégénérescence des équations (où les dérivées temporelles s'annulent pour certains ordres), les auteurs dérivent des relations linéaires permettant de déterminer $nb,j$ à l'aide de termes d'ordre inférieur, contournant ainsi la difficulté de résoudre directement des EDP dépendantes du temps d'ordre élevé.

Portée et affirmations
L'article revendique fournir la première vérification mathématique rigoureuse de la structure asymptotique de la couche limite pour une classe d'équations de chemotaxie-fluide sous la limite de viscosité nulle.

• Nouveauté mathématique : Il étend la théorie des limites de viscosité nulle de la dynamique des fluides pure aux couplages complexes biologie-fluide. Le travail aborde la « région hautement singulière » créée par l'annulation simultanée des viscosités du fluide et de la chimotaxie.
• Pertinence physique : Les résultats offrent un fondement mathématique pour comprendre la formation des couches limites et de la micro-turbulence observées dans les suspensions bactériennes (par exemple, le phénomène de « ligne blanche » dans les boîtes de Pétri). Il confirme que la viscosité du fluide joue un rôle critique dans la détermination du comportement de la couche limite chimotactique.
• Avancement méthodologique : L'application réussie des espaces de Sobolev conormaux anisotropes à ce système couplé démontre un cadre robuste pour gérer l'interplay entre les termes non linéaires de chimotaxie et les couches limites des fluides.

Les auteurs restent modestes concernant le taux de convergence, notant que le taux $O(\varepsilon^{3/2})$ est sous-optimal par rapport au potentiel théorique du développement (qui va jusqu'à l'ordre deux) en raison de la perte inhérente de mise à l'échelle $\varepsilon$ lors de l'estimation des dérivées normales et de la troncature du développement à la correction d'ordre deux. Ils reconnaissent que l'assouplissement des hypothèses de régularité sur les données initiales reste un défi pour les travaux futurs.