Analytic First Derivatives of SIDER Interpolation

Résumé Technique : Dérivées Analytiques de Premier Ordre de l'Interpolation SIDER

Énoncé du Problème
En analyse numérique et en modélisation géométrique, il existe un besoin croissant d'interpoler des données contraintes à des variétés non linéaires, spécifiquement la sphère unité $S^2 \subset \mathbb{R}^3$. L'interpolation polynomiale euclidienne standard ne parvient pas à préserver la contrainte sphérique entre les points de données, nécessitant souvent une normalisation ponctuelle qui complique la précision et l'analyse des dérivées. Bien que l'Interpolation Sphérique d'ordre $n$ (SIDER-$n$) ait été introduite précédemment comme un outil de reconstruction de haut ordre pour les méthodes Spherical Essentially Non-Oscillatory (SENO), elle manquait d'une théorie analytique pour les dérivées de premier ordre. L'absence d'informations sur les dérivées limite l'application de SIDER dans les schémas numériques nécessitant des pentes locales, des changements directionnels ou des quantités différentielles (par exemple, les solveurs de Riemann généralisés, les procédures de Lax–Wendroff inverses, ou le quadrature de haut ordre).

Méthodologie
L'article développe des formules de dérivées analytiques de premier ordre pour les courbes SIDER d'ordre arbitraire en exploitant la structure récursive de l'interpolation. La méthodologie centrale repose sur les principes suivants :

1. Composition Récursive : SIDER-$n$ est entièrement construit à partir d'opérations d'Interpolation Linéaire Sphérique (SLERP), formant un arbre binaire de cartes SLERP.
2. Propagation par la Règle de la Chaîne : Au lieu de différencier une expression de haut ordre entièrement développée, les auteurs appliquent la règle de la chaîne directement à travers la définition récursive. La dérivée d'une courbe SIDER-$n$ est exprimée comme la dérivée totale d'une opération SLERP appliquée à deux branches de SIDER-$(n-1)$ d'ordre inférieur.
3. Dérivée Totale de SLERP : Un bloc de construction fondamental est dérivé pour la dérivée totale de SLERP avec des extrémités mobiles et un paramètre d'interpolation mobile. Cette formule prend en compte la variation de l'angle géodésique et des extrémités.
4. Spécialisation : La formule récursive générale est spécialisée pour :
   • Les Nœuds d'Interpolation Intérieurs : Où les deux branches d'ordre inférieur se rejoignent au même point de donnée, permettant une combinaison pondérée simplifiée des dérivées.
   • Les Extrémités : Où une branche est évaluée en dehors de l'intervalle d'interpolation standard (extrapolation), introduisant des termes de carte logarithmique.
   • Les Points Médians : Emplacements situés à mi-chemin entre les échantillons de données, qui sont critiques pour les algorithmes de reconstruction de haut ordre.

Contributions Clés

• Formules de Dérivées Récursives : L'article fournit un algorithme récursif compact et prêt pour l'implémentation pour calculer $D_n(t) = \frac{d}{dt}S_n(t)$ pour tout ordre $n$. Cela évite les expressions de haut ordre fastidieuses et développées.
• Cas de Base Explicites : Le cas de base SIDER2 est traité explicitement, et SIDER3 ainsi que SIDER4 sont utilisés pour illustrer le mécanisme récursif et la gestion spécifique des extrémités et des nœuds intérieurs.
• Préservation de la Tangence : Une preuve rigoureuse est fournie démontrant que le vecteur dérivé analytiquement est intrinsèquement tangent à la sphère en chaque point reconstruit (tant aux nœuds d'échantillonnage qu'aux points médians). Cette propriété ($S_n(t) \cdot D_n(t) = 0$) découle de la structure algébrique de la dérivée de SLERP et ne nécessite aucune projection de post-traitement.
• Dérivées aux Points Médians : L'article dérive des formules spécifiques pour les dérivées aux coordonnées d'indice semi-entier ($\xi = j + 1/2$), qui sont essentielles pour évaluer les reconstructions entre les échantillons de données dans les schémas de type volume fini et WENO.

Résultats

• Récursion Générale : La dérivée $D_n$ est calculée en évaluant la dérivée totale de SLERP sur les deux branches d'ordre inférieur $A(t)$ et $B(t)$, mises à l'échelle par le facteur de reparamétrage $n/(n-1)$.
• Simplification des Nœuds : Aux nœuds intérieurs $t_j = j/n$, la dérivée se simplifie en une somme pondérée des dérivées des deux branches d'ordre inférieur rejoignant $p_j$. Les poids somment à $n/(n-1)$, reflétant la mise à l'échelle de la paramétrisation.
• Comportement aux Extrémités : Aux extrémités, la dérivée inclut des termes impliquant la carte logarithmique ($\text{Log}_p$) pour rendre compte de la direction des valeurs de la branche extrapolée.
• Cohérence Géométrique : Des visualisations numériques confirment que les vecteurs dérivés calculés résident exactement dans l'espace tangent de la sphère aux points interpolés, validant la preuve théorique de tangence.

Signification et Portée
L'article affirme que ces formules étendent le cadre SIDER/SENO original en fournissant les informations différentielles nécessaires pour les reconstructions à valeurs sphériques. La signification réside dans la capacité de Sider à être utilisé dans les méthodes de volume fini de haut ordre, ENO/WENO, et de type SENO pour les lois de conservation et les problèmes d'évolution où le champ reconstruit doit être différencié tout en préservant strictement la contrainte de la variété. Le travail est présenté comme une théorie analytique qui interprète le calcul de la dérivée comme une différenciation en mode direct de l'algorithme d'évaluation SIDER.

Les auteurs notent que l'implémentation reste compacte et suit directement la définition récursive. L'article conclut en reconnaissant que bien que la profondeur mathématique des résultats soit significative, le travail est destiné à être un préprint pour la diffusion et une future vérification formelle par les pairs, avec des travaux futurs potentiels impliquant un espacement de paramètres non uniforme et des extensions à d'autres variétés comme les groupes de rotation.