A mathematical study of an elastic-viscous-plastic sea-ice model with the Kelvin-Voigt rheology

Di seguito è presentata una sintesi tecnica dettagliata del articolo "A mathematical study of an elastic-viscous-plastic sea-ice model with the Kelvin-Voigt rheology" di Boutros, Liu, Thomas e Titi.

1. Enunciato del Problema e Motivazione

Il lavoro affronta l'analisi matematica del modello di ghiaccio marino Elastico-Visco-Plastico (EVP), una regolarizzazione numerica ampiamente utilizzata del classico modello visco-plastico di Hibler. Sebbene il modello EVP sia computazionalmente efficiente per le simulazioni climatiche, la sua corretta impostazione matematica rigorosa (esistenza, unicità e dipendenza continua dai dati iniziali) rimane una sfida, in particolare in presenza di termini di advezione non lineari e senza l'uso di tagli artificiali sui coefficienti di viscosità.

Sfide Chiave:

• Reologia Non-Newtoniana: Il ghiaccio marino è modellato come un fluido non newtoniano con relazioni complesse tra sforzo e deformazione che coinvolgono la plasticità.
• Termine di Advezione: L'inclusione del termine di advezione non lineare ($u \cdot \nabla u$) nell'equazione della quantità di moto complica le stime energetiche e le dimostrazioni di esistenza globale.
• Taglio della Viscosità: I modelli standard richiedono spesso un "taglio" (parametro di regolarizzazione $\epsilon$) sul tensore di velocità di deformazione per evitare singolarità. Rimuovere questo taglio (prendendo il limite $\epsilon \to 0$) è un problema aperto maggiore per i modelli multidimensionali.
• Limitazioni Precedenti: In un lavoro precedente degli autori (arXiv:2505.03080), hanno dimostrato la corretta impostazione per un modello regolarizzato di Voigt in cui la regolarizzazione era applicata alla relazione costitutiva (evoluzione del tensore degli sforzi). Tuttavia, tale approccio richiedeva un'alta regolarità dei dati iniziali e non poteva gestire il termine di advezione quando si rimuoveva il taglio della viscosità.

Obiettivo:
Gli autori propongono una nuova formulazione in cui la regolarizzazione di Kelvin-Voigt è applicata direttamente all'equazione della quantità di moto (campo di velocità) piuttosto che all'evoluzione del tensore degli sforzi. Il loro obiettivo è dimostrare:

1. La corretta impostazione locale per il modello con il termine di advezione.
2. La corretta impostazione globale per il modello senza il termine di advezione.
3. La capacità di rimuovere il taglio della viscosità ($\epsilon \to 0$) in entrambi i casi.
4. L'esistenza di soluzioni per dati iniziali meno regolari rispetto agli studi precedenti.

2. Formulazione Matematica

Lo studio considera il modello EVP su un toro piatto bidimensionale $\mathbb{T}^2$. Il sistema è composto da un'equazione della quantità di moto e un'equazione di evoluzione per il tensore degli sforzi $\sigma$.

Il Sistema Regolarizzato (Kelvin-Voigt nella Quantità di Moto):
Gli autori introducono un parametro di regolarizzazione $\alpha > 0$ e modificano l'equazione della quantità di moto (1.6a/1.7a) per includere un termine $-\alpha^2 \partial_t \Delta u$:
$$\partial_t(u - \alpha^2 \Delta u) + (u \cdot \nabla)u = \nabla \cdot \sigma + \text{termini di forzante}$$
$$\frac{1}{E} \partial_t \sigma + \text{termini non lineari degli sforzi} = D(u)$$
dove $D(u)$ è il tensore della velocità di deformazione.

Regolarizzazione della Velocità di Deformazione:
Per gestire la singolarità nella relazione sforzo-deformazione, il modulo della velocità di deformazione $D = |D(u)|$ è regolarizzato come $D_\epsilon = \sqrt{|D(u)|^2 + \epsilon^2}$. Le dimostrazioni comportano la costruzione di soluzioni per $\epsilon > 0$ e il successivo passaggio al limite $\epsilon \to 0$.

3. Metodologia

La strategia di dimostrazione si basa su un'approssimazione multi-step e su un argomento di compattezza:

1. Doppia Regolarizzazione:

   • Velocità di Deformazione: Utilizzo di $D_\epsilon$ per ammorbidire la non linearità.
   • Quantità di Moto: Introduzione di un termine di regolarizzazione di ordine superiore $\beta^4 \partial_t \Delta^2 u$ (con $\beta > 0$) nell'equazione della quantità di moto per garantire una regolarità sufficiente per lo schema di approssimazione di Galerkin.
   • Dati Iniziali: Mollificazione dei dati iniziali $(u_0, \sigma_0)$ utilizzando un mollificatore standard $\phi_\delta$.

2. Approssimazione di Galerkin:

   • Costruzione di soluzioni approssimate a dimensione finita $(u_N, \sigma_N)$ utilizzando modi di Fourier.
   • Dimostrazione dell'esistenza locale per questi sistemi di equazioni differenziali ordinarie (ODE) mediante il teorema di Picard-Lindelöf.

3. Stime A Priori:

   • Stima $L^\infty$ per gli Sforzi (Passaggio Cruciale): Gli autori derivano un nuovo limite $L^\infty$ per il tensore degli sforzi $\sigma$ (Lemma 2.3). Introducendo una variabile trasformata $\tau = \sigma + \frac{P}{2}I_2$ e sfruttando la struttura di smorzamento dell'equazione di evoluzione degli sforzi, dimostrano:
     $$\|\sigma(\cdot, t)\|_{L^\infty} \leq \|\sigma_0\|_{L^\infty} + 2P$$
     Questo limite è indipendente dai parametri di regolarizzazione ed è essenziale per chiudere le stime di ordine superiore.
   • Stime Energetiche: Derivazione di stime $H^1$ e $H^2$ per la velocità e stime $H^1$ per gli sforzi. Il limite $L^\infty$ su $\sigma$ permette il controllo dei termini non lineari che coinvolgono il tensore degli sforzi.

4. Passaggio al Limite:

   • Passo 1: Far tendere $\beta, \delta \to 0$ per recuperare la soluzione per il sistema con $D_\epsilon$ ma senza la regolarizzazione artificiale $\beta$.
   • Passo 2: Far tendere $\epsilon \to 0$ per rimuovere il taglio della velocità di deformazione.
   • Compattezza: Utilizzo del teorema di Banach-Alaoglu (per la convergenza debole-*) e del lemma di Aubin-Lions (per la convergenza forte) per estrarre sottosuccessioni convergenti.

5. Unicità e Dipendenza Continua:

   • Stabilire l'unicità analizzando la differenza tra due soluzioni e applicando la disuguaglianza di Grönwall, sfruttando il limite $L^\infty$ sul tensore degli sforzi.

4. Contributi e Risultati Chiave

Teorema 1.1 (Corretta Impostazione Locale con Advezione):

• Risultato: Dimostra l'esistenza di una soluzione forte locale nel tempo unica per il modello EVP di Kelvin-Voigt advettivo.
• Regolarità: $u \in C([0, T]; H^2)$ e $\sigma \in C([0, T]; H^1) \cap C([0, T]; L^\infty)$.
• Significato: Questo è il primo risultato rigoroso per il modello EVP multidimensionale con advezione in cui il taglio della viscosità è rimosso ($\epsilon \to 0$). Le stime precedenti non riuscivano a gestire il termine di advezione nel limite.

Teorema 1.2 (Corretta Impostazione Globale senza Advezione):

• Risultato: Dimostra l'esistenza di una soluzione forte globale nel tempo unica per il modello EVP di Kelvin-Voigt quando il termine di advezione è omesso.
• Significato: L'esistenza globale vale per dati iniziali con regolarità inferiore rispetto a quanto richiesto nel lavoro precedente degli autori (dove la regolarizzazione era sul tensore degli sforzi). Nello specifico, la velocità è regolarizzata, permettendo dati iniziali in $H^2$ invece di richiedere una maggiore regolarità.

Gestione del Taglio della Viscosità:

• Il lavoro giustifica con successo il limite $\epsilon \to 0$ (rimozione del taglio artificiale sui coefficienti di viscosità) in presenza di advezione (Teorema 1.1), un problema che rimane aperto per il modello standard di Hibler.

Interpretazione Fisica:

• Si sostiene che la regolarizzazione di Kelvin-Voigt nell'equazione della quantità di moto sia fisicamente più intuitiva rispetto alla regolarizzazione della relazione costitutiva, in quanto modella direttamente il comportamento elastico-viscoso del ghiaccio nel bilancio delle forze.

5. Significato e Impatto

• Rigor Matematico: Il lavoro fornisce la prima teoria rigorosa di corretta impostazione per il modello EVP con advezione e senza tagli di viscosità in più dimensioni.
• Superamento di Problemi Aperti: Risolve un importante problema aperto riguardante il limite $\epsilon \to 0$ in presenza di advezione, che in precedenza era intrattabile per i modelli di Hibler e Voigt-EVP.
• Requisiti di Regolarità Migliorati: Spostando la regolarizzazione sull'equazione della quantità di moto, gli autori abbassano i requisiti di regolarità dei dati iniziali, rendendo il modello applicabile a una classe più ampia di scenari fisici.
• Nuovi Strumenti Analitici: La derivazione della stima $L^\infty$ per il tensore degli sforzi (Lemma 2.3) è un contributo tecnico innovativo che permette la chiusura delle stime energetiche in presenza di forti non linearità.
• Modellazione Climatica: Questi risultati forniscono una solida base matematica per gli schemi numerici utilizzati nelle simulazioni climatiche su larga scala, garantendo che i modelli siano matematicamente fondati e stabili in condizioni realistiche.

In sintesi, questo lavoro avanza la teoria matematica della dinamica del ghiaccio marino stabilendo risultati robusti di corretta impostazione per una regolarizzazione motivata fisicamente del modello EVP, superando significativi ostacoli analitici legati alla non linearità, all'advezione e ai limiti di viscosità singolari.