Weak solutions and singular limits for a compressible fluid-structure interaction problem with slip boundary conditions

Résumé technique : Solutions faibles et limites singulières pour un problème d'interaction fluide-structure compressible avec conditions aux limites de glissement

Énoncé du problème
Cet article étudie un modèle mathématique décrivant l'interaction entre un fluide barotrope compressible tridimensionnel et une coque (plaque) solide viscoélastique bidimensionnelle. Le fluide est régi par les équations de Navier–Stokes compressibles, tandis que la structure élastique est modélisée par une équation linéaire d'ordre quatre (flexion) dans des coordonnées lagrangiennes. Une caractéristique déterminante de cette étude est l'imposition de conditions aux limites de glissement de Navier à l'interface fluide-structure, permettant un glissement tangentiel et une continuité de la vitesse normale, plutôt que la condition standard d'adhérence (no-slip). L'étude aborde deux défis principaux :

1. L'existence de solutions faibles globales pour le système couplé, en particulier en présence de frontières mobiles et de conditions de glissement.
2. La justification rigoureuse de la limite singulière où le fluide passe de compressible/visqueux à incompressible/non visqueux (régimes de faible nombre de Mach et de haut nombre de Reynolds), convergeant vers un système d'interaction Euler-plaque.

Les auteurs distinguent deux configurations géométriques : une géométrie de référence générale (frontière courbe) et une géométrie de référence plate (domaine de type plaque), notant que cette dernière permet des résultats de régularité plus forts en raison de la forme spécifique du couplage cinématique.

Méthodologie
L'analyse utilise une combinaison de techniques avancées issues de la théorie des fluides compressibles et de l'interaction fluide-structure (FSI) :

• Existence de solutions faibles :

  • Extension du domaine et régularisation : Les auteurs étendent le domaine fluide vers un grand domaine fixe $B$ et approximent les coefficients de viscosité et la pression (en utilisant un terme de pression artificielle $\delta \varrho^\beta$) pour gérer la frontière mobile et la non-linéarité.
  • Discrétisation temporelle et décomposition d'opérateurs : En suivant le schéma de [45], le problème est découplé en sous-problèmes fluide et structure, résolus itérativement, puis la limite est prise sur le pas de temps.
  • Équi-intégrabilité de la densité : Un obstacle technique critique consiste à établir l'équi-intégrabilité de la densité près de la frontière mobile sous des conditions de glissement. Les auteurs exploitent la structure spécifique de la condition de glissement de Navier et, dans le cas de la géométrie plate, la relation résultante $\partial_t w = (u \cdot n_w) \circ \Phi_w$ pour obtenir une régularité spatiale supplémentaire pour le déplacement $w$. Cela permet d'assouplir les exigences sur l'exposant adiabatique $\gamma$ par rapport aux cas d'adhérence sans amortissement.
  • Flux visqueux effectif : Des techniques standard impliquant l'identité du flux visqueux effectif sont adaptées pour identifier la limite faible du terme de pression.

• Limites singulières (Limite incompressible non visqueuse) :

  • Méthode d'entropie relative modifiée : Pour comparer la solution faible du système compressible (définie sur un domaine dépendant du temps $\Omega_{w^\varepsilon}$) avec la solution forte du système limite Euler-plaque (définie sur $\Omega_\eta$), les auteurs emploient une inégalité d'énergie relative modifiée.
  • Transformation de Piola : Une transformation de coordonnées naïve échoue car elle détruit la propriété de divergence nulle de la vitesse limite et modifie les vecteurs normaux aux frontières. Les auteurs introduisent une transformation de Piola basée sur la composition des applications de déformation. Cette transformation préserve la propriété solénoïdale de la vitesse et mappe correctement les vecteurs normaux, garantissant que les fonctions tests satisfont les conditions aux limites nécessaires.
  • Contrôle des erreurs : L'analyse nécessite de contrôler les termes d'erreur issus de la transformation (impliquant le jacobien et le gradient de déformation). Cela nécessite des hypothèses de régularité supplémentaires sur le déplacement de la solution faible, justifiées par la présence d'amortissement structurel ($\nu_s > 0$) et de configurations géométriques spécifiques.

Contributions et résultats clés

1. Existence de solutions faibles (Théorème 1.4) :

   • Les auteurs prouvent l'existence de solutions faibles globales pour le système FSI compressible avec conditions aux limites de glissement de Navier.
   • Conditions :
     • Avec amortissement structurel ($\nu_s > 0$), l'existence est démontrée pour $\gamma > 3/2$.
     • Sans amortissement structurel ($\nu_s = 0$), l'existence est démontrée pour la géométrie de référence plate ($\Omega \in \mathcal{G}_2$) à condition que $\gamma > 12/7$.
   • Signification : Cela étend les résultats précédents (par exemple [7], [45]) qui se concentraient principalement sur les conditions d'adhérence ou exigeaient un $\gamma$ plus élevé pour les conditions de glissement. La preuve met en évidence que dans les géométries générales sans amortissement, le manque de régularité de la vitesse normale empêche la fermeture des estimations, nécessitant soit un amortissement, soit une géométrie plate.

2. Limite incompressible non visqueuse (Théorème 1.7) :

   • L'article justifie rigoureusement la convergence du système fluide-structure compressible visqueux vers le système incompressible non visqueux (Euler) couplé à une plaque viscoélastique.
   • Régime : Faible nombre de Mach ($\varepsilon \to 0$) et haut nombre de Reynolds ($\nu \to 0$).
   • Exigences : La limite est établie pour des données initiales bien préparées, une géométrie de référence plate, un amortissement structurel ($\nu_s > 0$) et $\gamma > 3/2$. Une hypothèse de régularité supplémentaire sur le déplacement de la solution faible (spécifiquement concernant $\nabla w$) est requise pour gérer les termes d'erreur dans l'estimation d'entropie relative.
   • Vitesse de convergence : Les auteurs dérivent un taux de convergence dépendant de la perturbation des données initiales et des paramètres $\varepsilon$ et $\nu$.

3. Limite incompressible (Théorème 4.14) :

   • Comme sous-produit, l'article aborde également la limite où seul le nombre de Mach tend vers zéro (en maintenant la viscosité fixe). Dans ce cas, le système converge vers les équations de Navier–Stokes incompressibles couplées à une plaque viscoélastique. Notamment, ce résultat vaut pour $\gamma > 3$ sans les hypothèses de régularité supplémentaires requises pour la limite non visqueuse.

Signification et affirmations
Les auteurs affirment que ce travail présente le premier résultat concernant les solutions faibles de fluides compressibles interagissant avec des structures élastiques sous conditions aux limites de glissement de Navier. De plus, il fournit la première justification rigoureuse de la limite singulière à faible nombre de Mach et haut nombre de Reynolds pour de tels problèmes d'interaction fluide-structure.

L'article souligne que l'interplay entre la condition de glissement et la géométrie de la frontière mobile crée des difficultés analytiques significatives, en particulier pour retrouver la condition de couplage cinématique et établir l'équi-intégrabilité de la densité. L'utilisation de la transformation de Piola pour gérer le décalage de domaine dans la méthode d'entropie relative est présentée comme une innovation technique cruciale pour dériver des limites singulières dans des domaines FSI dépendant du temps.

Les auteurs restent modestes concernant la nécessité d'hypothèses de régularité supplémentaires pour la limite non visqueuse, notant que la régularité uniforme des solutions faibles dans le régime non visqueux est insuffisante pour justifier la limite sans ces contraintes supplémentaires, une limitation partagée par des problèmes similaires dans la littérature.