The Landau--Lifshitz--Bloch equation on polytopal domains: Unique existence and finite element approximation

1. Il Problema e il Contesto

L'articolo si concentra sull'equazione di Landau-Lifshitz-Bloch (LLBE), un modello fondamentale per descrivere la dinamica della magnetizzazione nei ferromagneti a temperature elevate, in particolare sopra la temperatura di Curie ($T_c$). A differenza del modello standard di Landau-Lifshitz (valido solo a basse temperature), la LLBE interpola tra la dinamica a bassa temperatura e la teoria di Ginzburg-Landau delle transizioni di fase, permettendo di modellare fenomeni come la registrazione magnetica assistita dal calore (HAMR).

Il problema matematico studiato è la seguente equazione differenziale parziale non lineare su un dominio poliedrale limitato $D \subset \mathbb{R}^d$ ($d \in \{1, 2, 3\}$):
$$\begin{cases} \partial_t u = \kappa_1 \Delta u + \gamma u \times \Delta u - \kappa_2(1 + \mu|u|^2)u & \text{in } (0, T) \times D, \\ u(0, \cdot) = u_0 & \text{in } D, \\ \frac{\partial u}{\partial \nu} = 0 & \text{su } (0, T) \times \partial D. \end{cases}$$
Dove $u(t, x) \in \mathbb{R}^3$ è il campo di magnetizzazione macroscopica. Un aspetto cruciale è che, a differenza dell'equazione di Landau-Lifshitz standard, il modulo di $u$ non è conservato (non è vincolato alla sfera unitaria), il che rende inapplicabili molti metodi numerici esistenti progettati per il vincolo sferico.

Il dominio $D$ è considerato poliedrale (o poligonale in 2D), permettendo angoli rientranti (non convessi) nel caso bidimensionale, il che introduce sfide di regolarità per le soluzioni.

2. Metodologia

Gli autori adottano un approccio che combina analisi funzionale avanzata e metodi agli elementi finiti (FEM):

1. Analisi di Esistenza e Unicità: Utilizzano il metodo di Faedo-Galerkin per dimostrare l'esistenza e l'unicità di soluzioni forti (con regolarità superiore) per la LLBE su domini poliedrali. Vengono stabilite stime di decadimento temporale per la soluzione.
2. Schemi FEM Lineari: Propongono uno schema agli elementi finiti lineare e completamente discreto per l'approssimazione diretta della LLBE. Questo schema richiede la risoluzione di un sistema lineare ad ogni passo temporale, riducendo significativamente il costo computazionale rispetto agli schemi non lineari esistenti.
3. Regolarizzazione Viscosa ($\epsilon$-LLBE): Poiché lo schema diretto per la LLBE mostra tassi di convergenza sub-ottimali e richiede restrizioni severe sul rapporto tra passo temporale e dimensione della mesh, gli autori introducono una regolarizzazione viscosa (pseudo-parabolica). Aggiungendo il termine $-\epsilon \Delta \partial_t u$, ottengono l'equazione $\epsilon$-LLBE:
   $$\partial_t u_\epsilon - \epsilon \Delta \partial_t u_\epsilon = \kappa_1 \Delta u_\epsilon + \gamma u_\epsilon \times \Delta u_\epsilon - \kappa_2(1 + \mu|u_\epsilon|^2)u_\epsilon$$
4. Analisi di Convergenza:
   • Dimostrano che la soluzione $u_\epsilon$ della problema regolarizzato converge alla soluzione $u$ della LLBE originale quando $\epsilon \to 0^+$.
   • Sviluppano uno schema FEM lineare per l'$\epsilon$-LLBE e ne analizzano la stabilità e l'errore.
   • Derivano stime di errore uniformi nel tempo, superando il problema della crescita esponenziale dell'errore tipica delle stime a priori su intervalli temporali lunghi.

3. Contributi Chiave

I principali contributi teorici e numerici del lavoro sono:

• Esistenza e Unicità su Domini Non Lisci: Estendono i risultati di esistenza e unicità delle soluzioni forti della LLBE a domini poliedrali (inclusi quelli non convessi in 2D con angoli rientranti), un ambito precedentemente non esplorato in letteratura.
• Comportamento a Lungo Termine: Dimostrano rigorosamente che la magnetizzazione decade esponenzialmente verso lo stato zero ($u \to 0$) per $t \to \infty$, coerentemente con la fisica della transizione di fase sopra $T_c$.
• Schema FEM Lineare Efficiente: Propongono uno schema lineare per la LLBE che evita la risoluzione di sistemi non lineari ad ogni passo, offrendo un vantaggio computazionale significativo.
• Regolarizzazione e Ottimalità: Introducono la regolarizzazione $\epsilon$-LLBE per ottenere:
  • Stime di errore ottimali per lo schema numerico.
  • Stime di errore uniformi nel tempo, cruciali per simulazioni a lungo termine.
  • Un tasso di convergenza della soluzione regolarizzata $u_\epsilon$ verso $u$ dell'ordine di $O(\epsilon)$.
• Analisi su Domini con Singolarità: Gestiscono le singolarità associate agli angoli rientranti (corner singularities) nei domini non convessi, dimostrando che la regolarità della soluzione è limitata dall'angolo del dominio, influenzando il tasso di convergenza spaziale.

4. Risultati Teorici e Numerici

• Stime di Errore:
  • Per la LLBE diretta (senza regolarizzazione), lo schema FEM converge con un tasso sub-ottimale $O(\ell_h(h^\alpha + k))$, dove $\ell_h$ è un fattore logaritmico o di potenza dipendente dalla dimensione $d$, e $\alpha$ dipende dalla regolarità del dominio.
  • Per l'$\epsilon$-LLBE, lo schema FEM lineare ottiene stime di errore ottimali $O(h^\alpha + k)$ (o $O(h+k)$ in domini convessi con regolarità sufficiente).
  • Viene dimostrata la convergenza uniforme nel tempo dell'errore per l'approssimazione dell'$\epsilon$-LLBE.
• Simulazioni Numeriche:
  • Vengono presentate quattro simulazioni su domini diversi: quadrato (2D), cubo (3D), dominio a "L" (non convesso 2D) e dominio con angolo di Fichera (non convesso 3D).
  • I risultati confermano i tassi di convergenza teorici per $h$ (spaziale), $k$ (temporale) e $\epsilon$.
  • Le simulazioni mostrano la formazione di pareti di Bloch e il successivo decadimento della magnetizzazione, validando il modello fisico.
  • Viene osservato che la regolarizzazione $\epsilon$ non degrada la qualità della soluzione rispetto al caso originale, ma ne migliora le proprietà numeriche.

5. Significato e Impatto

Questo lavoro colma un divario significativo tra la teoria fisica della magnetizzazione ad alta temperatura e l'analisi numerica rigorosa.

• Rigorosità Matematica: Fornisce la prima prova di esistenza e unicità di soluzioni forti su domini poliedrali non lisci, essenziale per applicazioni pratiche dove i domini geometrici sono complessi.
• Efficienza Computazionale: Lo sviluppo di uno schema lineare stabile e convergente rende possibile la simulazione efficiente di fenomeni magnetici complessi su larga scala, superando le limitazioni degli schemi non lineari precedenti.
• Robustezza: L'uso della regolarizzazione viscosa offre un percorso teorico e numerico robusto per trattare problemi non lineari su domini con singolarità, garantendo stabilità e convergenza uniforme nel tempo, proprietà fondamentali per la simulazione di sistemi fisici reali.

In sintesi, l'articolo fornisce un quadro completo, dalla teoria all'implementazione numerica, per la risoluzione dell'equazione LLBE in contesti fisici realistici e geometricamente complessi.