Free-surface Stokes problem: stability estimates and time-step improvements

Sintesi Tecnica: Problema di Stokes a superficie libera: stime di stabilità e miglioramenti del passo temporale

Definizione del Problema
La simulazione accurata di flussi con superficie libera altamente viscosi, come la dinamica delle calotte glaciali, la convezione del mantello e il flusso di lava, richiede la risoluzione delle equazioni accoppiate Stokes/superficie libera. Questi problemi sono caratterizzati da un dominio mobile in cui la velocità del fluido e la pressione sono accoppiate con l'evoluzione dell'altezza della superficie libera. Sebbene le discretizzazioni standard che utilizzano elementi finiti nello spazio e metodi di Eulero nel tempo siano comuni, esse spesso soffrono di severi vincoli di stabilità. Nello specifico, gli schemi di Eulero espliciti e semi-impliciti richiedono frequentemente passi temporali proibitivi per rimanere stabili, portando a costi computazionali elevati. Le tecniche di stabilizzazione esistenti, come l'Algoritmo di Stabilizzazione della Superficie Libera (FSSA), migliorano le dimensioni del passo temporale ma mancano di limiti teorici a priori sulla stabilità e non garantiscono la stabilità incondizionata o la conservazione del volume.

Metodologia
Gli autori presentano un'analisi teorica rigorosa della stabilità e delle proprietà di conservazione del sistema accoppiato Stokes/superficie libera sia per fluidi Newtoniani che non Newtoniani (viscosità con legge di potenza). L'analisi è condotta a due livelli:

1. Livello Continuo: Gli autori derivano un bilancio di energia totale che accoppia l'energia delle equazioni del momento di Stokes con l'energia dell'equazione dell'altezza della superficie libera. Questa derivazione si basa sull'instaurare specifiche relazioni tra il termine di avvezione della superficie libera e gli integrali di superficie, collegandoli successivamente al termine di gravità nelle equazioni di Stokes.
2. Livello Discreto: Lo studio esamina schemi completamente discreti utilizzando elementi finiti nello spazio combinati con l'integrazione temporale di Eulero esplicita, semi-implicita e implicita.
   • Regolarizzazione: Per affrontare la mancanza di regolarità intrinseca della funzione di altezza della superficie libera (necessaria per la ben posto-età del problema di Stokes), gli autori introducono un termine di stabilizzazione sui bordi di ordine superiore, $J(s^{n+1}, w)$, trattato implicitamente nel tempo.
   • Proposta di Stabilizzazione: Motivati dall'identificazione di termini "spuri" nelle stime di energia del metodo di Eulero esplicito, gli autori propongono un nuovo termine di stabilizzazione aggiunto alla forma debole di Stokes. Questo termine penalizza la velocità normale alla superficie libera ed è progettato per cancellare precisamente i termini destabilizzanti derivanti dal trattamento esplicito dell'equazione della superficie libera.

Contributi Chiave

• Bilancio dell'Energia Totale: Il documento deriva un nuovo bilancio dell'energia totale per il sistema accoppiato Stokes/superficie libera, accoppiando l'energia del momento con l'energia dell'altezza superficiale. Questo è identificato come un elemento fondamentale per comprendere la stabilità, che era precedentemente assente nella letteratura.
• Dimostrazioni di Stabilità Teorica:
  • Eulero Implicito: Dimostrato essere incondizionatamente stabile rispetto alla dimensione del passo temporale (assumendo la continuità Lipschitz della superficie libera) e capace di conservare il volume del dominio.
  • Eulero Esplicito (non stabilizzato): Dimostrato essere instabile nella norma dell'energia totale per qualsiasi passo temporale a causa di termini spuri relativi alla velocità normale superficiale.
  • Eulero Semi-implicito (FSSA): Mostrato essere solo condizionatamente stabile e, criticamente, dimostrato non conservare il volume del dominio.
  • Eulero Esplicito Stabilizzato Proposto: Gli autori dimostrano che il loro nuovo termine di stabilizzazione rende il metodo di Eulero esplicito incondizionatamente stabile rispetto alla dimensione del passo temporale (sotto l'assunzione di continuità Lipschitz) preservando al contempo il volume del dominio.
• Conservazione del Volume: Il documento fornisce prove teoriche e verifiche numeriche che lo schema esplicito stabilizzato proposto conserva il volume del dominio, mentre lo schema semi-implicito stabilizzato da FSSA non lo fa.

Risultati
Esperimenti numerici sono stati condotti utilizzando la libreria open-source FEniCS sia per casi di fluidi Newtoniani che non Newtoniani, incluso un caso con geometria complessa basata sui dati della calotta glaciale della Groenlandia.

• Stabilità: Lo schema di Eulero esplicito non stabilizzato era instabile per tutti i passi temporali testati. Gli schemi stabilizzati da FSSA (sia semi-impliciti che espliciti) erano condizionatamente stabili, fallendo a passi temporali più grandi. Al contrario, lo schema EE-Stabilized proposto è rimasto stabile per passi temporali arbitrariamente grandi in entrambi i casi, Newtoniano e non Newtoniano.
• Conservazione del Volume: Lo schema EE-Stabilized ha conservato il volume del dominio fino alla precisione di macchina. Lo schema SIE-FSSA non ha conservato il volume, confermando i risultati teorici.
• Accuratezza: I test di convergenza hanno mostrato che i termini di stabilizzazione proposti sono coerenti con il problema discretizzato. Lo schema EE-Stabilized ha mantenuto un'accuratezza del primo ordine nell'altezza superficiale e nella velocità sotto il raffinamento simultaneo di mesh e passo temporale, con errori comparabili o inferiori ad altri schemi stabilizzati.
• Oscillazioni: Il termine di stabilizzazione dei bordi implicito ha efficacemente soppresso le oscillazioni ad alta frequenza nell'altezza della superficie libera senza introdurre ulteriori restrizioni sul passo temporale.

Significatività e Rivendicazioni
Il documento sostiene che lo schema di Eulero esplicito stabilizzato proposto offre un'alternativa robusta ed efficiente ai risolutori implicitamente costosi. Garantendo teoricamente la stabilità incondizionata rispetto alla dimensione del passo temporale e assicurando la conservazione del volume, il metodo consente simulazioni stabili di flussi a superficie libera anche con geometrie complesse e termini sorgente. Gli autori sottolineano che la loro stabilizzazione è la prima a garantire teoricamente passi temporali arbitrari per la discretizzazione di Eulero esplicito di questo problema accoppiato. L'approccio è notato per essere simmetrico (benefico per i risolutori iterativi), semplice da implementare e privo di parametri da regolare. Gli autori concludono che, sebbene la loro analisi assuma la continuità Lipschitz della superficie libera, le stime di stabilità in norma $L^2$ derivate sono una base necessaria per la ben posto-età una volta applicata l'appropriata regolarizzazione. Si suggerisce come lavoro futuro l'esplorazione di termini di regolarizzazione per la ben posto-età e l'estensione di questa analisi di stabilità ai problemi Navier-Stokes/superficie libera e a integratori temporali di ordine superiore.