Error analysis for a Finite Element Discretization of a corotational harmonic map heat flow problem

1. Il Problema Matematico

Il lavoro si concentra sul flusso di calore delle mappe armoniche (HMHF), un problema di evoluzione che descrive il gradiente $L^2$ dell'energia di Dirichlet per campi vettoriali vincolati a giacere sulla sfera unitaria $S^2$. L'equazione generale è:
$$\partial_t u = \Delta u + |\nabla u|^2 u$$
con vincolo di lunghezza unitaria $|u|=1$.

Gli autori restringono l'analisi a un caso specifico chiamato caso corotazionale. Assumendo che il dominio $\Omega$ sia il disco unitario in $\mathbb{R}^2$ e che la soluzione abbia una simmetria rotazionale specifica, il problema tridimensionale si riduce a un'equazione scalare in coordinate polari $(r, t)$ su $I=[0,1]$:
$$\partial_t u = \partial_{rr} u + \frac{1}{r} \partial_r u - \frac{\sin(2u)}{2r^2}$$
con condizioni al contorno omogenee $u(t,0)=u(t,1)=0$.
Questo caso è fondamentale perché:

• È collegato alla stabilità e alla formazione di singolarità (blow-up) nelle soluzioni dell'HMHF bidimensionale.
• Per dati iniziali con $|u_0(1)| \le \pi$, esiste una soluzione globale liscia, mentre per $|u_0(1)| > \pi$ si verifica un blow-up in tempo finito.

2. Metodologia e Discretizzazione

Gli autori propongono e analizzano un metodo numerico basato su:

• Spazio: Discretizzazione agli Elementi Finiti (FEM) conformi $H^1$ su una griglia uniforme. Viene utilizzato uno spazio di funzioni continue a tratti polinomiali di grado $k \ge 1$.
• Tempo: Schemi di Euler semi-implicito. La parte lineare dell'equazione è trattata implicitamente, mentre il termine non lineare (legato a $\sin(2u)$) è trattato in modo semi-implicito (valutato al passo temporale precedente).
  • Questo approccio trasforma il problema non lineare in un problema lineare risolvibile univocamente ad ogni passo temporale, riducendo il costo computazionale rispetto a metodi completamente impliciti che richiederebbero iterazioni di Newton.

3. Contributi Chiave e Innovazioni Teoriche

L'analisi di errore presentata è innovativa per diversi aspetti tecnici:

1. Proprietà di Convessità Nascosta:
   L'energia del sistema $E(u)$ non è convessa. Tuttavia, gli autori identificano una "convessità nascosta" attraverso una decomposizione dell'energia:
   $$E(u) = \frac{1}{2} \|u\|_{H^1_r}^2 - \int_0^1 F(u(r)) \frac{dr}{r}$$
   dove $F(v) = \frac{1}{2}(v^2 - \sin^2 v)$ è una funzione convessa. Questa proprietà è cruciale per dimostrare la stabilità discreta e controllare l'errore di linearizzazione.

2. Stima dell'Energia e Stabilità Discreta:
   Viene dimostrata una stima di dissipazione dell'energia discreta che mima quella della soluzione continua. Questo garantisce che l'energia numerica non cresca artificialmente, fornendo un limite uniforme sulla soluzione discreta ($\|u_h\|_{L^\infty} \le \pi$) sotto opportune condizioni iniziali.

3. Controllo dell'Errore di Linearizzazione:
   La parte più difficile dell'analisi è il controllo dell'errore introdotto dalla linearizzazione del termine non lineare. Gli autori utilizzano una combinazione di:

   • Una stima $L^\infty$ sulla proiezione di Galerkin della soluzione continua scalata per $1/r$ (Lemma 4.3).
   • L'argomento di convessità menzionato sopra.
     Questo permette di ottenere un limite di errore ottimale senza richiedere ipotesi di regolarità eccessive o condizioni di stabilità restrittive sul passo temporale.

4. Norme Pesate:
   A causa della trasformazione in coordinate polari, l'analisi utilizza norme pesate con il fattore $r$ (es. $\|v\|_{0,r}$, $\|v\|_{H^1_r}$) per gestire correttamente la singolarità in $r=0$ e le proprietà di embedding dello spazio.

4. Risultati Teorici e Numerici

• Stime di Errore: Il teorema principale (Teorema 4.4) stabilisce che, assumendo che la soluzione continua sia sufficientemente regolare, l'errore di discretizzazione soddisfa limiti ottimali:
  • Errore nella norma $H^1_r$: $O(\tau + h^k)$.
  • Errore nella norma $L^2_r$: $O(\tau + h^{k+1})$ (osservato numericamente, anche se non coperto dalla teoria principale).
    Qui $\tau$ è il passo temporale e $h$ è la dimensione della mesh.
• Validazione Numerica: Gli autori presentano esperimenti numerici (utilizzando il pacchetto Netgen/NGSolve) che confermano le previsioni teoriche:
  • Si osservano tassi di convergenza ottimali per elementi lineari ($k=1$) e quadratici ($k=2$).
  • Viene mostrata la corretta dissipazione dell'energia nel tempo, in accordo con il Teorema 3.1.
  • I risultati confermano che il metodo è stabile e preciso per soluzioni lisce.

5. Significato e Impatto

Questo lavoro è significativo perché:

• Fornisce la prima analisi di errore rigorosa per una discretizzazione agli elementi finiti del caso corotazionale dell'HMHF, colmando un vuoto nella letteratura esistente (che si era concentrata su differenze finite o su analisi di stabilità senza stime di errore ottimali).
• Dimostra che un metodo semi-implicito lineare è sufficiente per ottenere stabilità e convergenza ottimale, rendendo il metodo computazionalmente efficiente per simulazioni pratiche.
• Offre strumenti analitici (la decomposizione convessa dell'energia) che potrebbero essere applicati ad altre varianti di problemi di flusso di calore su varietà, sebbene gli autori notino che l'estensione al caso LLG (Landau-Lifshitz-Gilbert) completo o al caso equivariante generale sia ostacolata dalla mancanza di tale struttura convessa.

In sintesi, il paper stabilisce un fondamento teorico solido per l'uso di metodi agli elementi finiti nella simulazione di dinamiche di mappe armoniche con simmetria rotazionale, garantendo sia la stabilità numerica che la precisione asintotica.