Analytic First Derivatives of SIDER Interpolation

Sintesi Tecnica: Derivate Prime Analitiche dell'Interpolazione SIDER

Definizione del Problema
Nella analisi numerica e nella modellazione geometrica, esiste una crescente necessità di interpolare dati vincolati a varietà non lineari, nello specifico la sfera unitaria $S^2 \subset \mathbb{R}^3$. L'interpolazione polinomiale euclidea standard non riesce a preservare il vincolo sferico tra i punti dati, richiedendo spesso una normalizzazione puntuale che complica l'accuratezza e l'analisi delle derivate. Sebbene l'Interpolazione Sferica di ordine $n$ (SIDER-$n$) sia stata precedentemente introdotta come uno strumento di ricostruzione ad alto ordine per i metodi Spherical Essentially Non-Oscillatory (SENO), essa mancava di una teoria analitica per le derivate prime. L'assenza di informazioni sulle derivate limita l'applicazione di SIDER in schemi numerici che richiedono pendenze locali, cambiamenti direzionali o quantità differenziali (ad esempio, risolutori di Riemann generalizzati, procedure inverse di Lax–Wendroff, o quadratura ad alto ordine).

Metodologia
Il documento sviluppa formule analitiche per le derivate prime di curve SIDER di ordine arbitrario sfruttando la struttura ricorsiva dell'interpolazione. La metodologia principale si basa sui seguenti principi:

1. Composizione Ricorsiva: SIDER-$n$ è costruito interamente tramite operazioni di Interpolazione Lineare Sferica (SLERP), formando un albero binario di mappe SLERP.
2. Propagazione della Regola della Catena: Inveve di derivare un'espressione ad alto ordine completamente espansa, gli autori applicano la regola della catena direttamente attraverso la definizione ricorsiva. La derivata di una curva SIDER-$n$ è espressa come la derivata totale di un'operazione SLERP applicata a due rami SIDER-$(n-1)$ di ordine inferiore.
3. Derivata Totale di SLERP: Un elemento fondamentale viene derivato per la derivata totale di una SLERP con endpoint mobili e un parametro di interpolazione mobile. Questa formula tiene conto della variazione dell'angolo geodetico e degli endpoint.
4. Specializzazione: La formula ricorsiva generale viene specializzata per:
   • Nodi di Interpolazione Interi: Dove i due rami di ordine inferiore si incontrano nello stesso punto dato, permettendo una combinazione pesata semplificata delle derivate.
   • Endpoint: Dove un ramo viene valutato al di fuori dell'intervallo di interpolazione standard (estrapolazione), introducendo termini di mappa logaritmica.
   • Punti Medi: Posizioni a metà strada tra i campionamenti dei dati, che sono critiche per gli algoritmi basati sulla ricostruzione ad alto ordine.

Contributi Chiave

• Formule di Derivata Ricorsiva: Il documento fornisce un algoritmo ricorsivo compatto e pronto per l'implementazione per calcolare $D_n(t) = \frac{d}{dt}S_n(t)$ per qualsiasi ordine $n$. Ciò evita espressioni espanse ingombranti per curve ad alto ordine.
• Casi Base Espliciti: Il caso base SIDER2 è trattato esplicitamente, e SIDER3 e SIDER4 sono utilizzati per illustrare il meccanismo ricorsivo e la gestione specifica di endpoint e nodi interi.
• Preservazione della Tangenza: Viene fornita una prova rigorosa che dimostra che il vettore derivata analiticamente derivato è intrinsecamente tangente alla sfera in ogni punto ricostruito (sia nei nodi di campionamento che nei punti medi). Questa proprietà ($S_n(t) \cdot D_n(t) = 0$) deriva dalla struttura algebrica della derivata SLERP e non richiede proiezioni post-elaborazione.
• Derivate dei Punti Medi: Il documento deriva formule specifiche per le derivate nei coordinati di indice semi-interi ($\xi = j + 1/2$), essenziali per valutare le ricostruzioni tra i campioni di dati in schemi a volumi finiti e di tipo WENO.

Risultati

• Ricorsione Generale: La derivata $D_n$ è calcolata valutando la derivata totale SLERP sui due rami di ordine inferiore $A(t)$ e $B(t)$, scalati dal fattore di riparametrizzazione $n/(n-1)$.
• Semplificazione dei Nodi: Ai nodi interi $t_j = j/n$, la derivata si semplifica in una somma pesata delle derivate dei due rami di ordine inferiore che incontrano $p_j$. I pesi sommano a $n/(n-1)$, riflettendo la scala della parametrizzazione.
• Comportamento agli Endpoint: Agli endpoint, la derivata include termini che coinvolgono la mappa logaritmica ($\text{Log}_p$) per tenere conto della direzione dei valori del ramo estrapolato.
• Consistenza Geometrica: Visualizzazioni numeriche confermano che i vettori derivata calcolati giacciono esattamente nello spazio tangente della sfera nei punti interpolati, validando la prova teorica di tangenza.

Significato e Ambito
Il documento afferma che queste formule estendono il framework originale SIDER/SENO fornendo le necessarie informazioni differenziali per le ricostruzioni a valori sferici. La significatività risiede nel rendere possibile l'uso di SIDER in metodi ad alto ordine a volumi finiti, ENO/WENO e di tipo SENO per leggi di conservazione e problemi di evoluzione, dove il campo ricostruito deve essere differenziato pur preservando rigorosamente il vincolo della varietà. Il lavoro è presentato come una teoria analitica che interpreta il calcolo della derivata come una differenziazione in modalità forward dell'algoritmo di valutazione SIDER.

Gli autori osservano che l'implementazione rimane compatta e segue direttamente la definizione ricorsiva. Il documento conclude riconoscendo che, sebbene la profondità matematica dei risultati sia significativa, il lavoro è inteso come una preprint per la disseminazione e la futura verifica formale tra pari, con potenziali lavori futuri riguardanti la spaziatura dei parametri non uniforme e l'estensione ad altre varietà come i gruppi di rotazione.