Convergence of a two-parameter hyperbolic relaxation system toward the incompressible Navier-Stokes equations

Résumé technique : Convergence d'un système de relaxation hyperbolique à deux paramètres vers les équations de Navier-Stokes incompressibles

Énoncé du problème
L'article étudie la convergence rigoureuse d'un système de relaxation hyperbolique à deux paramètres (1.2) vers les équations de Navier-Stokes incompressibles (1.1) dans des domaines périodiques en deux dimensions ($T^2$) et en trois dimensions ($T^3$). Le système d'approximation (1.2) intègre deux paradigmes classiques : la méthode de compressibilité artificielle (AC), qui relaxe la contrainte de divergence nulle via une équation d'évolution temporelle pour la pression, et un modèle de relaxation de premier ordre pour un tenseur de type contrainte $U$. Le système est régi par de petits paramètres $\epsilon$ et $\delta$, représentant les temps de relaxation de la pression et du tenseur de contrainte, respectivement.

Le principal défi mathématique abordé est la limite simultanée $(\epsilon, \delta) \to 0$. Si les limites à un seul paramètre (soit $\epsilon \to 0$ pour l'AC, soit $\delta \to 0$ pour les modèles de relaxation) sont bien documentées, la limite simultanée présente des difficultés pour contrôler les erreurs de pression et gérer les perturbations des données initiales. Plus précisément, les auteurs visent à établir la convergence pour deux régimes distincts de perturbations de la vitesse initiale :

1. Perturbations vanishing (s'annulant) : Où la différence initiale $\|u_{\epsilon,\delta}^0 - u^0\|$ est de l'ordre de $O((\epsilon+\delta)^\alpha)$.
2. Perturbations d'ordre un ($O(1)$) : Où la différence initiale est de l'ordre de $O(1)$, un régime crucial pour la robustesse numérique mais traditionnellement difficile pour les méthodes basées sur l'énergie en raison de la dissipation limitée dans le système de relaxation.

Méthodologie
Les auteurs emploient des stratégies analytiques distinctes selon l'amplitude des perturbations initiales et la dimensionnalité du domaine.

• Système affine intermédiaire (Perturbations vanishing) : Pour dériver les taux de convergence pour la pression dans le régime de perturbation vanishing, les auteurs construisent un système affine intermédiaire. Ce système agit comme un pont entre le système de relaxation original (1.2) et le système cible de Navier-Stokes (1.1). En bornant les écarts entre le système intermédiaire et les systèmes limite et d'approximation, les auteurs obtiennent les estimations d'erreur nécessaires pour la pression, surmontant la limitation où les méthodes d'énergie directes ne produisent généralement que des bornes $O(1)$ pour les erreurs de pression.

• Énergie modulée et arguments de bootstrap (Perturbations $O(1)$) : Pour le régime de perturbation $O(1)$, les estimations d'énergie hyperboliques symétriques standards sont insuffisantes car la dissipation fournie par l'amortissement du tenseur $U_{\epsilon,\delta}$ n'est pas assez forte pour supprimer les effets non linéaires sur de longs intervalles de temps. À la place, les auteurs utilisent une énergie modulée (introduite dans les travaux antérieurs [4]) impliquant le tourbillon $\omega_{\epsilon,\delta} = \text{curl } u_{\epsilon,\delta}$.

  • Cette approche tire parti de la dissipation de $\nabla \omega_{\epsilon,\delta}$ pour contrôler le champ de vitesse.
  • Un argument de bootstrap est mis en œuvre pour fermer les estimations. Cela nécessite des relations d'échelle spécifiques entre les paramètres (spécifiquement $\epsilon = o(\delta)$) et des contraintes sur les données initiales où la divergence initiale de la vitesse est suffisamment petite par rapport au rotationnel initial.
  • Correction de couche initiale : Pour gérer le tenseur de type contrainte $U_{\epsilon,\delta}$ sans exiger qu'il soit $o(1)$ proche des données initiales limites, les auteurs introduisent des corrections de couche initiale (termes de décroissance exponentielle) dans les estimations d'énergie.

• Distinctions dimensionnelles :

  • En 3D ($T^3$), le tourbillon est un champ vectoriel. L'analyse identifie deux « cas extrêmes » distincts pour les taux de convergence basés sur l'ordre de la norme $L^2$ de la perturbation de vitesse initiale ($O(\delta^{3/4})$ vs. $O(\delta^{1/4})$), conduisant aux Théorèmes 2.6 et 2.7.
  • En 2D ($T^2$), le tourbillon est un scalaire. L'analyse révèle une seule construction viable pour l'ordre de convergence des perturbations initiales (Théorème 2.8).

Résultats clés
L'article établit les résultats de convergence rigoureux suivants :

1. Régime de perturbation vanishing (Théorèmes 2.1, 2.5, Corollaire 2.3) :

   • Pour des perturbations initiales suffisamment petites, la solution du système de relaxation existe sur le même intervalle de temps que la solution lisse de Navier-Stokes (ou un intervalle logarithmiquement étendu si la solution de Navier-Stokes est globale).
   • La convergence de la vitesse $u_{\epsilon,\delta} \to u$ est prouvée dans $L^\infty([0,T]; H^1)$ et $L^\infty([0,T]; L^2)$.
   • Crucialement, sous des conditions plus strictes (spécifiquement $\delta = o(\sqrt{\epsilon})$), la pression $p_{\epsilon,\delta}$ est montrée comme convergeant vers $p$ dans $L^\infty([0,T]; H^1)$.

2. Régime de perturbation d'ordre un (Théorèmes 2.6, 2.7, 2.8) :

   • Les auteurs prouvent la convergence globale de la vitesse (ou jusqu'à la durée de vie de la solution limite) même lorsque les perturbations de vitesse initiales sont de l'ordre de $O(1)$.
   • Résultats 3D : Deux théorèmes fournissent des bornes de convergence pour la vitesse en $L^2$ et le tourbillon en $H^1$ sous différentes hypothèses d'échelle pour les données initiales. Les bornes dépendent de $\delta$ et du rapport $\epsilon/\delta$.
   • Résultat 2D : Une convergence de vitesse globale similaire est établie pour le cas 2D avec une mise à l'échelle spécifique des perturbations initiales.
   • Dans tous les cas $O(1)$, la divergence du champ de vitesse $\text{div } u_{\epsilon,\delta}$ est montrée comme restant petite (contrôlée par $\epsilon$ et $\delta$), garantissant que la solution reste proche de la variété incompressible.

Signification et revendications
L'article prétend fournir une validation théorique fondamentale pour le modèle de relaxation hyperbolique à deux paramètres (1.2) dans les simulations numériques. Ses contributions spécifiques incluent :

• Limite simultanée à deux paramètres : Il comble la lacune dans la littérature existante concernant la limite simultanée $(\epsilon, \delta) \to 0$, un scénario moins exploré que les limites à un seul paramètre.
• Convergence de la pression : En introduisant le système affine intermédiaire, ce travail dérive rigoureusement les estimations d'erreur pour la pression dans le régime de perturbation vanishing, une étape souvent absente des analyses de modèles de relaxation précédentes.
• Robustesse face aux grandes perturbations : L'extension de la convergence de la vitesse au régime de perturbation $O(1)$ démontre la robustesse du modèle face à des incohérences initiales significatives, qui sont communes lors des démarrages numériques. Ceci est accompli en surmontant les limitations des estimations d'énergie standard grâce à la méthode d'énergie modulée et aux arguments de bootstrap.
• Aperçus structurels : L'analyse souligne la richesse structurelle du cas 3D par rapport au cas 2D, identifiant des régimes de taux de convergence distincts en trois dimensions.

Les auteurs notent que bien que la convergence de la vitesse soit établie pour les perturbations $O(1)$, la convergence de la pression dans ce régime spécifique demeure un défi ouvert pour la recherche future. Le travail est destiné à soutenir l'application de ce modèle dans des tâches en aval, telles que les statistiques des écoulements turbulents, en assurant la validité mathématique de l'approximation.