Orthogonal polynomials for the de Rham complex on the disk and cylinder

Titolo

Polinomi ortogonali per il complesso di de Rham su disco e cilindro.

1. Problema e Contesto

Il lavoro affronta la necessità di costruire basi polinomiali per la discretizzazione numerica delle equazioni alle derivate parziali (PDE) su domini geometrici specifici (disco e cilindro), rispettando la struttura topologica e geometrica del complesso di de Rham.
In particolare, l'obiettivo è sviluppare basi per gli spazi di Sobolev pertinenti ($\mathring{H}^1$, $\mathring{H}(\text{curl})$, $L^2$) che:

1. Rispettino le condizioni al contorno (funzioni nulle al bordo per $\mathring{H}^1$, vettori normali al bordo per $\mathring{H}(\text{curl})$).
2. Sfruttino la simmetria rotazionale del dominio, permettendo una decoupling (separazione) delle equazioni in modi di Fourier.
3. Consentano relazioni di ricorrenza semplici e sparse tra gli operatori differenziali (gradiente, rotore/curl, divergenza), essenziali per l'efficienza computazionale nel Calcolo Esterno agli Elementi Finiti (FEEC).

Le basi esistenti (come i polinomi di Zernike scalari) non sono direttamente adattabili ai campi vettoriali con condizioni al contorno di tipo "normale" o per la costruzione di sottocomplessi esatti nel contesto vettoriale/matriciale.

2. Metodologia

L'autore utilizza un approccio basato sulla teoria delle rappresentazioni dei gruppi di simmetria e sull'analisi degli operatori differenziali in coordinate polari.

• Basi Adattate alla Simmetria (Symmetry-Adapted Bases):
  Si definiscono funzioni scalari, vettoriali e matriciali che si trasformano secondo rappresentazioni irriducibili del gruppo di rotazione $SO(2)$. Per un modo $m$, una rotazione di angolo $\phi$ corrisponde alla moltiplicazione per $e^{im\phi}$ (scalare) o a una combinazione di rotazione vettoriale e fase (vettoriale). Questo garantisce l'ortogonalità automatica tra modi diversi ($m \neq n$).

• Costruzione di Polinomi Vettoriali e Matriciali:
  Partendo dai polinomi armonici generalizzati (analoghi ai polinomi di Zernike scalari $z_{mj}^{(\lambda)}$), l'autore costruisce:

  • Polinomi vettoriali armonici ($\mathbf{y}_m$).
  • Polinomi matriciali armonici ($\mathbf{Y}_m$).
  • Polinomi ortogonali vettoriali e matriciali pesati rispetto a pesi ultrasferici $(1-r^2)^\lambda$.

• Il Ruolo del Peso Matriciale $N(x,y)$:
  Il contributo centrale è l'introduzione di una matrice peso specifica:
  $$N(x, y) = \begin{pmatrix} 1-y^2 & xy \\ xy & 1-x^2 \end{pmatrix}$$
  Questa matrice ha la proprietà fondamentale di annullare la componente tangenziale di un vettore al bordo del disco ($r=1$). Moltiplicando polinomi vettoriali per $N(x,y)$, si ottiene una base per lo spazio $\mathcal{N}$ dei polinomi vettoriali normali al bordo.

• Riduzione a Polinomi Univariate:
  Per gestire i pesi matriciali complessi, il problema su disco viene mappato su un modulo di polinomi vettoriali univariati $\mathcal{M}$ su $[0,1]$. Utilizzando fattorizzazioni LU (o Cholesky) su operatori di tipo Jacobi, si costruiscono esplicitamente polinomi ortogonali per pesi specifici (come $V^{(b)}(t) = \frac{1}{1-t}t^b$) in termini di polinomi di Jacobi.

3. Risultati Chiave e Costruzioni

Il paper definisce tre famiglie di basi che formano il complesso di de Rham discreto:

1. Basi Scalari ($w_{mj}$):
   Polinomi di Zernike pesati che si annullano al bordo. Costituiscono una base per $\mathring{H}^1$.
   $$w_{mj}(x,y) = (1-r^2) z_{mj}^{(1)}(x,y)$$

2. Basi Vettoriali Normali ($n_{mj}^\nu$):
   Una nuova famiglia di polinomi vettoriali ortogonali rispetto al peso matriciale $N(x,y)$, che sono normali al bordo.
   $$n_{mj}^\nu(x,y) = N(x,y) v_{N,mj}^\nu(x,y)$$
   Dove $v_{N,mj}^\nu$ sono polinomi vettoriali costruiti tramite i polinomi univariati $q^{(b),\nu}_n$.

3. Basi Scalari $L^2$ ($z_{mj}$):
   I classici polinomi di Zernike standard, base per $L^2$.

Relazioni di Ricorrenza Semplici (Teoremi 6 e 7):
La scoperta più significativa è la semplicità delle relazioni tra gli operatori differenziali applicati a queste basi:

• Gradiente: Il gradiente di una funzione scalare pesata è proporzionale a una combinazione lineare delle basi vettoriali normali.
  $$\nabla w_{mj} = -(j+1) n_{m, j+1}^+$$
• Rotore (Curl): Il rotore delle basi vettoriali normali è proporzionale ai polinomi di Zernike scalari standard (con un fattore costante dipendente da $m$ e $j$).
  $$\text{curl } n_{mj}^- = \kappa_{mj} z_{mj}$$
  Inoltre, $\text{curl } n_{mj}^+ = 0$.

Queste relazioni permettono di decouplare il complesso di de Rham in piccoli sottocomplessi esatti per ogni modo $m$ e grado $j$:
$$0 \to \text{span}\{w_{mj}\} \xrightarrow{\nabla} \text{span}\{n_{m,j+1}^+, n_{m,j+1}^-\} \xrightarrow{\text{curl}} \text{span}\{z_{m,j+1}\} \to 0$$

4. Estensione ai Cilindri

I risultati sono generalizzati al caso tridimensionale su cilindri (sia periodici che finiti):

• Cilindro Periodico: Si utilizza la serie di Fourier nella direzione $z$. Le basi 3D sono prodotti tensoriali tra le basi sul disco e $e^{ikz}$. Le relazioni di ricorrenza includono termini di derivata rispetto a $z$ (che agiscono come moltiplicazione per $ik$).
• Cilindro Finito: Si sostituiscono le serie di Fourier con polinomi ultrasferici pesati (o funzioni di Legendre integrate) nella direzione $z$, mantenendo la struttura del complesso di de Rham e le condizioni al contorno.

5. Significato e Impatto

• FEEC e Risoluzione di PDE: Le basi costruite sono ideali per il Finite Element Exterior Calculus. La struttura di sottocomplessi esatti garantisce la stabilità numerica e la conservazione delle proprietà topologiche (come la divergenza nulla o il rotore nullo) nella discretizzazione.
• Efficienza Computazionale: Le relazioni di ricorrenza semplici e sparse permettono di costruire matrici di massa e operatori differenziali in modo efficiente, portando potenzialmente a risolutori a complessità ottimale per equazioni come quelle di Maxwell o di Stokes in geometrie cilindriche.
• Generalizzabilità: Il metodo basato su simmetria adattata e pesi matriciali offre un quadro teorico solido per estendere queste costruzioni a domini più complessi (anelli, sfere, capsule sferiche) e ad altri complessi (Koszul, elasticità).
• Nuova Classe di Polinomi: L'introduzione dei polinomi vettoriali $n_{mj}^\nu$ riempie un vuoto nella letteratura sui polinomi ortogonali multivariati, fornendo una base naturale per i campi vettoriali con condizioni al contorno di tipo Neumann/Robin (normali al bordo).

In sintesi, il paper fornisce un ponte teorico e computazionale fondamentale tra la teoria dei polinomi ortogonali, la teoria delle rappresentazioni e l'analisi numerica delle PDE su domini con simmetria rotazionale.