A consistent-splitting generalized scalar auxiliary variable scheme for the perturbed Boussinesq system

Sintesi Tecnica: Uno schema di tipo Consistent-Splitting Generalized Scalar Auxiliary Variable per il sistema di Boussinesq perturbato

Enunciato del Problema
Questo lavoro affronta la simulazione numerica del sistema di Boussinesq perturbato bidimensionale, che modella flussi guidati dalla galleggiabilità in ambienti stabilmente stratificati (ad esempio, dinamiche atmosferiche e oceaniche). Il sistema risultante (Equazione 6) è derivato sottraendo un equilibrio idrostatico stabile e linearmente stratificato dalle equazioni standard di Boussinesq. Questo sistema perturbato presenta un accoppiamento specifico in cui la galleggiabilità guida il moto verticale nell'equazione del momento, mentre il moto verticale adveta il gradiente di temperatura lineare di fondo nell'equazione della temperatura. Questo accoppiamento genera onde di gravità interne e guida il rilassamento delle perturbazioni verso l'equilibrio idrostatico. La principale sfida matematica risiede nel vincolo di incompressibilità e nella struttura a punto di sella del sistema velocità-pressione, complicata dalla necessità di un'accuratezza di ordine superiore e stabilità in presenza di convezione non lineare e termini di stratificazione lineare.

Metodologia
Gli autori propongono uno schema di splitting consistente del secondo ordine basato sull'approccio Generalized Scalar Auxiliary Variable (GSAV). La metodologia integra diverse componenti chiave:

1. Framework GSAV: Una variabile ausiliaria scalare $r(t)$ viene introdotta per emulare la legge dell'energia continua del sistema. Ciò consente la costruzione di un'energia discreta modificata che rimane incondizionatamente stabile.
2. Discretizzazione Temporale Consistent-Splitting: Lo schema estende il framework BDF2 (Backward Differentiation Formula) del secondo ordine con spostamento temporale sviluppato da Huang e Shen [17] per le equazioni di Navier-Stokes.
   • Trattamento Esplicito: Convezione non lineare ($u \cdot \nabla u$), advezione ($u \cdot \nabla \theta$) e accoppiamenti lineari (galleggiabilità $\theta e_2$ e stratificazione $\alpha v$) sono trattati esplicitamente tramite estrapolazione del secondo ordine al livello temporale traslato $t_{n+k}$.
   • Decoupling: Questo trattamento esplicito riduce ogni passo temporale alla risoluzione di un piccolo numero di sistemi lineari disaccoppiati a coefficienti costanti, evitando la necessità di solver non lineari o sistemi accoppiati monolitici.
3. Correzione della Pressione: Un'equazione di Poisson per la pressione viene risolta utilizzando un approccio di splitting consistente che mantiene un'accuratezza ottimale del secondo ordine, evitando le inconsistenze delle condizioni al contorno tipiche dei metodi di proiezione standard.
4. Analisi della Stabilità e dell'Errore: Gli autori dimostrano un teorema di stabilità debole incondizionata e derivano stime di errore ottimali del secondo ordine per velocità, pressione e temperatura. L'analisi traccia rigorosamente la dipendenza delle costanti di errore dalla viscosità cinematica ($\nu$) e dalla diffusività termica ($\gamma$).

Contributi Chiave
Il documento apporta tre contributi distinti rispetto alla letteratura esistente sugli schemi di tipo SAV per i sistemi di Boussinesq:

1. Analisi del Sistema Perturbato: A differenza dei lavori precedenti che analizzano il sistema di Boussinesq standard, questo articolo affronta specificamente il sistema perturbato (6) contenente il termine di accoppiamento della stratificazione lineare $\alpha v$ nell'equazione della temperatura.
2. Integrazione del Framework: Mentre Wagner et al. [36] hanno applicato l'approccio GSAV al sistema di Boussinesq standard, questo lavoro adotta il framework specifico di splitting consistente del secondo ordine con spostamento temporale di Huang e Shen [17], mantenendo però il meccanismo di aggiornamento GSAV standard, garantendo una piena accuratezza del secondo ordine sia nelle equazioni di velocità che di temperatura attraverso l'estrapolazione coerente di tutti i termini di accoppiamento.
3. Dipendenza Esplicita della Costante di Errore: Un contributo critico è il tracciamento esplicito della dipendenza della costante di errore da $\nu^{-1}$ e $\gamma^{-1}$. L'analisi rivela che la costante di errore contiene potenze negative di viscosità e diffusività, portando a una dipendenza esponenziale quadruplamente annidata rispetto a questi parametri.

Risultati

• Stabilità Teorica: Si dimostra che lo schema è incondizionatamente debolmente stabile. Le sequenze della variabile ausiliaria e del fattore di scala rimangono non negative, e le soluzioni numeriche soddisfano vincoli uniformi indipendenti dalla dimensione del passo temporale, a condizione che la costante di stabilizzazione $C$ superi una soglia determinata dai termini di forzante.
• Convergenza: Si stabiliscono tassi di convergenza ottimali del secondo ordine per la velocità, la pressione e la temperatura in norme appropriate.
• Esperimenti Numerici:
  • Test di Convergenza: Una soluzione fabbricata su un dominio limitato conferma la convergenza teorica del secondo ordine sia nelle norme $L^2$ che $H^1$.
  • Simulazione di Flusso Stratificato: Una simulazione a lungo termine ($T=10^4$) di un flusso stratificato con perturbazioni iniziali di temperatura gaussiane riproduce la dinamica fisica attesa:
    • Un regime transitorio di galleggiabilità in cui l'energia potenziale si converte in energia cinetica.
    • Un regime di oscillazione di onde interne con un periodo coerente con le previsioni teoriche ($T \approx \pi\sqrt{2/\alpha}$).
    • Un regime di decadimento esponenziale in cui le perturbazioni si rilassano verso l'equilibrio idrostatico, con tassi di decadimento coerenti con i limiti analitici derivati nel Lemma 2.

Significatività e Limitazioni
Il documento sostiene che lo schema proposto offre un metodo computazionalmente efficiente per simulare flussi stratificati, riducendo ogni passo temporale a risoluzioni lineari disaccoppiate pur mantenendo l'accuratezza del secondo ordine. Tuttavia, gli autori sono modesti riguardo alla robustezza dello schema nel limite di viscosità o diffusività che tendono a zero.

L'analisi dimostra che la costante di errore cresce come un esponenziale quadruplamente annidato di $\nu^{-1}$ e $\gamma^{-1}$. Di conseguenza, lo schema non è robusto quando $\nu, \gamma \to 0$. Gli esperimenti numerici confermano che per valori piccoli di viscosità e diffusività termica, lo schema richiede risoluzioni spaziali e temporali via via più fini per rimanere stabile, e può divergere a risoluzioni più grossolane. Gli autori concludono che, sebbene lo schema sia efficace per regimi pienamente dissipativi, sarebbe necessario un trattamento completamente implicito o stabilizzato dei termini di convezione per ottenere limiti di errore uniformi in $\nu$ e $\gamma$. Questa limitazione evidenzia un compromesso tra l'efficienza computazionale dell'approccio di splitting consistente esplicito e la sua stabilità in regimi a bassa diffusione.