Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow

Résumé Technique : Reconstruction compatible avec l'entropie pour les écoulements viscoélastiques à nombre de Weissenberg élevé

Énoncé du Problème
Le problème du nombre de Weissenberg élevé demeure un obstacle central dans la simulation des écoulements viscoélastiques. Dans les modèles de tenseur de conformation, l'état polymérique est représenté par un tenseur symétrique défini positif (SPD) $A$. Bien que la positivité soit essentielle pour la signification physique de l'énergie libre élastique ($\log \det A$) et du logarithme matriciel, l'article soutient que la positivité seule est insuffisante pour la stabilité numérique aux nombres de Weissenberg élevés.

Les méthodes existantes, telles que la conformation logarithmique ($A = \exp \Psi$) et les reconstructions par racine carrée ($A = BB^T$), imposent avec succès la contrainte SPD. Cependant, les expériences numériques montrent que ces méthodes peuvent encore produire des dépassements de contrainte (overshoots) non physiques et une croissance non physique de l'énergie élastique. Le problème central identifié est un décalage entre la contrainte de cône (positivité) et la contrainte de bilan couplé (compatibilité de l'énergie libre). Plus précisément, si le tenseur utilisé dans le travail de contrainte de l'équation de quantité diffère du tenseur utilisé dans le travail d'étirement ou les variables d'entropie de l'équation de conformation, un défaut de travail à signe indéfini apparaît. De plus, les reconstructions strictement positives peuvent introduire trois mécanismes non physiques spécifiques :

1. Biais d'entropie de type Jensen : Les applications de cartes non linéaires (comme l'exponentielle) à des perturbations de moyenne nulle dans l'espace logarithmique génèrent un biais de moyenne positif dans le tenseur physique et l'entropie.
2. Amplification exponentielle : Dans les états fortement étirés, de petites perturbations logarithmiques sont amplifiées par l'exponentielle de matrice, avec une mise à l'échelle proportionnelle à la plus grande valeur propre $\lambda_{\max}(A)$.
3. Défauts de travail à signe indéfini : L'utilisation de tenseurs reconstruits discordants pour les termes de contrainte et d'entropie empêche l'annulation exacte du travail polymérique, entraînant une injection ou un retrait incontrôlé d'énergie.

Méthodologie
L'article propose une reconstruction logarithmique corrigée et compatible avec l'entropie. Plutôt que de remplacer les méthodes de préservation de la positivité existantes, cette approche agit comme un filtre d'admissibilité appliqué à une reconstruction logarithmique d'ordre élevé brute avant qu'elle n'entre dans le bilan couplé contrainte-travail.

La méthodologie implique :

• Définition du chemin logarithmique : Étant donné un prédicteur physique $\hat{A}$ et une reconstruction d'ordre élevé brute $\tilde{A}$, un chemin est défini dans l'espace logarithmique : $A(\theta) = \exp(\hat{\Psi} + \theta(\tilde{\Psi} - \hat{\Psi}))$ pour $\theta \in [0, 1]$.
• Contrainte du budget d'entropie : La méthode cherche le plus grand $\theta$ tel que l'entropie (énergie libre) discrète totale ne dépasse pas l'entropie du prédicteur physique plus un budget local $\tau_h$. Le budget est mis à l'échelle pour correspondre à l'erreur de troncature attendue du schéma d'ordre élevé (par exemple, $O(h^{2k+2})$).
• Algorithme de bisection : L'optimal $\theta$ est trouvé via une recherche par bisection locale, cellule par cellule. Si la reconstruction brute satisfait le budget, $\theta=1$ (aucun changement). Sinon, la méthode revient le long du chemin vers le point admissible le plus élevé.
• Imposition de la cohérence : Le tenseur sélectionné est utilisé de manière cohérente à travers tous les termes de couplage : travail de contrainte de quantité, étirement de conformation, variables d'entropie et quadrature.

Contributions Clés et Résultats Théoriques
L'article fournit une analyse rigoureuse de ce principe de reconstruction, prouvant plusieurs propriétés clés :

• Existence et Convexité : La fonctionnelle d'entropie le long du chemin logarithmique est convexe, garantissant l'existence d'un paramètre admissible maximal $\theta^*$.
• Amortissement Minimal : Parmi tous les états admissibles sur le chemin, le $\theta^*$ sélectionné minimise la distance par rapport à la reconstruction brute, rendant la correction « à amortissement minimal ».
• Contrôle Spectral : La correction préserve l'indépendance de cadre et maintient les bornes spectrales dérivées des logarithmes des extrémités.
• Inactivité Asymptotique : Si la reconstruction brute d'ordre élevé est déjà admissible au sens de l'entropie (c'est-à-dire que son défaut est de l'ordre de l'erreur de discrétisation), la correction devient inactive ($\theta^* = 1$), préservant ainsi la précision formelle d'ordre élevé.
• Avantage de Résolution en Régime de Fort Étirement : Dans les régimes de fort étirement, la méthode conserve une précision d'ordre élevé dans les métriques log-relatives et d'entropie, là où les reconstructions directes dans l'espace physique souffrent de problèmes d'échelle de dérivée proportionnels à $\lambda_{\max}(A)$.
• Estimation d'Énergie : Les auteurs prouvent une estimation de l'énergie libre discrète pour le système couplé, montrant que le budget de reconstruction $\tau_h$ agit comme une source d'entropie contrôlée et bornée plutôt que comme une viscosité artificielle.

Résultats Numériques
L'article valide la théorie via une série de diagnostics :

• Tests Scalaires : Il a été démontré que les perturbations logarithmiques de moyenne nulle créent un biais d'entropie positif (biais de Jensen) et que les reconstructions linéaires, bien que de faible biais, ne préservent pas la positivité.
• Tests d'Amplification : Ils ont confirmé que les défauts de force de contrainte dans les reconstructions logarithmiques brutes croissent linéairement avec l'étirement de fond $\lambda_{\max}(A)$.
• Vérification du Défaut de Travail : L'utilisation de tenseurs discordants pour la contrainte et l'entropie a montré des défauts d'énergie à signe indéfini.
• Boucles de Haut Nombre de Weissenberg Dynamiques : Dans les simulations avec $Wi=80$ et un étirement élevé, les reconstructions logarithmiques standards accumulaient de l'entropie et des erreurs de force de contrainte au fil du temps, tandis que la reconstruction corrigée contrôlait ces défauts, réduisant les erreurs de force de plusieurs ordres de grandeur.
• Benchmarks de Résolution : La méthode corrigée a maintenu la compatibilité de l'entropie et contrôlé les défauts de force à travers les raffinements de maillage et les variations de l'étirement, alors que les méthodes brutes montraient des erreurs croissantes avec l'augmentation de l'étirement.

Signification et Revendications
L'article positionne ce travail comme comblant une lacune critique entre la préservation de la positivité (traitée par les méthodes de conformation logarithmique/racine carrée) et la compatibilité thermodynamique (traitée par les méthodes dissipatives d'énergie libre). Les auteurs affirment que :

1. La positivité est une condition nécessaire mais insuffisante pour la robustesse aux nombres de Weissenberg élevés.
2. La correction proposée est une couche locale, calculable et à contrôle spectral qui peut être ajoutée aux schémas d'ordre élevé existants sans altérer leur structure fondamentale.
3. La méthode empêche l'injection d'énergie non physique en garantissant que le tenseur entrant dans l'équation de quantité est compatible avec les variables d'entropie.
4. Elle offre un avantage de résolution conditionnel dans les régimes de fort étirement en opérant dans l'espace de la variable logarithmique lisse tout en respectant le budget d'entropie.

Les auteurs concluent que pour les calculs à nombre de Weissenberg élevé, surveiller la valeur propre minimale ne suffit pas ; il faut également suivre l'excès d'entropie de reconstruction et les défauts de force de contrainte pour garantir la fidélité physique. La reconstruction logarithmique corrigée fournit un mécanisme pour imposer cette compatibilité tout en restant asymptotiquement inactive pour les écoulements suffisamment résolus.