Update-Magnitude State Redistribution (UM-SRD): A Shut-off Extension of Weighted SRD for Cut-Cell Methods

Résumé technique : Redistribution d'état par magnitude de mise à jour (UM-SRD)

Énoncé du problème
Les méthodes de cellules coupées permettent aux schémas aux volumes finis de traiter des géométries complexes sans maillages adaptés au corps, mais elles introduisent des défis de stabilité sévères dus aux cellules coupées arbitrairement petites. Les méthodes de redistribution d'état (SRD), et plus spécifiquement le cadre SRD pondéré développé par Berger & Giuliani (2024), résolvent ces problèmes de stabilité en fusionnant les petites cellules en voisinages, en faisant évoluer les moyennes sur le maillage fusionné, puis en redistribuant les résultats vers le maillage d'origine.

Cependant, une limitation critique de la formulation SRD pondérée existante est l'absence d'un mécanisme de « coupure ». Contrairement aux méthodes classiques de redistribution de flux, où le post-traitement s'annule naturellement lorsque les déséquilibres de flux sont nuls, le SRD pondéré applique un opérateur de redistribution non trivial à chaque pas de temps, indépendamment de la magnitude de la mise à jour aux volumes finis. Cette activation continue entraîne deux déficiences pratiques :

1. Échec à préserver les états stationnaires : Même lorsque l'équilibre des flux est exactement nul, le SRD remplace les valeurs ponctuelles par des moyennes pondérées par le volume, empêchant la convergence vers la précision machine et introduisant des erreurs de $O(1)$ dans les régions statiques.
2. Dissipation inutile : Dans les régions lisses avec des mises à jour petites mais non nulles, l'application continue de la redistribution introduit une diffusion numérique proportionnelle au gradient du voisinage, indépendante de la dynamique physique réelle.

Méthodologie
Cet article présente la Redistribution d'état par magnitude de mise à jour (UM-SRD), une extension de coupure du SRD pondéré conçue pour combiner l'opérateur SRD avec l'opérateur identité en fonction de la magnitude locale de la mise à jour aux volumes finis.

La méthode opère dans le cadre standard de trois étapes de pas de temps :

1. Mise à jour aux volumes finis : Calcul de l'état pré-redistribution $U^*$ en utilisant le schéma de base.
2. Agrégation de voisinage : Calcul des moyennes pondérées de voisinage $\hat{Q}^*_j$ sur les voisinages fusionnés $M_j$.
3. Redistribution combinée : Au lieu d'appliquer directement l'opérateur SRD $S_j$, l'UM-SRD applique une combinaison convexe :
   $$R_j(s_j) = (1 - s_j)I + s_j S_j$$
   où $I$ est l'opérateur identité et $s_j \in [0, 1]$ est un paramètre de mélange.

Le paramètre de mélange $s_j$ est déterminé par une fonction de coupure lisse et localement définie $h(\eta_j)$ :
$$s_j = \frac{\eta_j^p}{\eta_j^p + \tau^p}$$
Ici, $\eta_j$ est un indicateur normalisé de la magnitude absolue maximale de la mise à jour au sein du voisinage $M_j$ (spécifiquement $\eta_j = \frac{\max \|\Delta U_i\|}{\epsilon + \max \|\Delta U_k\|}$).

• Régime de mise à jour nulle : Si toutes les mises à jour dans le voisinage sont à la précision machine ($\Delta U_i \approx 0$), alors $\eta_j = 0$, $s_j = 0$, et la méthode se réduit exactement à l'identité ($U^{n+1} = U^*$).
• Régime actif : Si des mises à jour significatives se produisent (par exemple, près de discontinuités ou lorsque la rigidité des petites cellules est active), $\eta_j \approx 1$, $s_j \approx 1$, et la méthode se comporte comme un SRD pondéré complet.

Contributions analytiques clés
L'article fournit une analyse rigoureuse de l'UM-SRD pour un modèle d'advection linéaire unidimensionnel avec une seule petite cellule coupée :

• Conservation : L'opérateur combiné est prouvé conservateur par construction pour tout $s_j \in [0, 1]$.
• Diminution de la variation totale (TVD) : Sous la même condition CFL que le schéma de base d'upwind d'ordre un, l'UM-SRD est prouvé TVD sur le maillage de cellules coupées d'origine. Cela repose sur la combinaison convexe de l'identité et de l'opérateur SRD TVD.
• Préservation des états stationnaires : Il est prouvé que l'UM-SRD préserve exactement les états stationnaires du schéma de base. Si la mise à jour aux volumes finis est nulle, $s_j=0$, et aucune redistribution ne se produit, évitant ainsi la dérive spurious inhérente au SRD standard.
• Précision : Dans les régions lisses, l'UM-SRD est montré comme une perturbation du schéma de base. Si le paramètre de mélange évolue de manière appropriée avec la taille du maillage (ou via l'indicateur normalisé utilisé en pratique), la contribution additionnelle à l'erreur de troncature locale est d'ordre supérieur par rapport au schéma de base d'ordre un. Par conséquent, l'UM-SRD conserve le taux de convergence d'ordre un de la méthode d'upwind sous-jacente.

Résultats numériques
Des expériences numériques sur des problèmes modèles en une et deux dimensions confirment les prédictions théoriques :

• Convergence : L'UM-SRD présente une convergence d'ordre un dans les normes $L_1$ et $L_\infty$ pour l'advection lisse, avec des erreurs indiscernables du schéma de base et du SRD pondéré standard (différences < 0,01 %).
• Comportement de coupure : Dans les tests avec des impulsions compactes, le paramètre de mélange $s_j$ chute exactement à zéro lorsque l'impulsion est absente, récupérant le schéma de base, et monte à $\approx 0,99$ lorsque l'impulsion interagit avec la petite cellule.
• Stabilité : L'UM-SRD stabilise le schéma de base sous des pas de temps de maillage complet (où le CFL local à la petite cellule est 10 fois la limite de stabilité), empêchant la divergence exponentielle observée dans le schéma d'upwind de base.
• Propriétés non oscillatoires : Des tests avec des données discontinues montrent que l'UM-SRD reste non oscillatoire, correspondant au comportement du SRD standard près de la petite cellule.
• Préservation des états stationnaires : Les expériences confirment que tandis que le SRD standard introduit des erreurs de $O(1)$ dans les profils statiques, l'UM-SRD maintient le profil initial à la précision machine.
• Extension 2D : L'algorithme a été appliqué avec succès à une grille cartésienne 2D avec une frontière immergée, démontrant une convergence d'ordre un sans modification de la logique de mélange.

Importance
L'article affirme que l'UM-SRD offre une modification pratique et minimale de l'infrastructure existante du SRD pondéré. En réutilisant les structures de voisinage et de poids établies et en ajoutant uniquement un mécanisme de mélange local, il confère au SRD une propriété de mise à jour s'annulant. Cela combine les avantages de stabilité de la redistribution d'état pour les cellules coupées arbitrairement petites avec les propriétés de préservation de la précision des méthodes de redistribution de flux, permettant spécifiquement une préservation exacte des états stationnaires et réduisant la dissipation inutile dans les régions lisses. Les auteurs notent que bien que l'analyse actuelle soit restreinte à l'advection linéaire scalaire 1D, la structure algorithmique est indépendante de la dimension et prête à être intégrée dans des codes de cellules coupées de production.