Barotropic-Baroclinic Splitting for Multilayer Shallow Water Models with Exchanges

Sintesi Tecnica: Scomposizione Barotropica-Baroclinica per Modelli di Shallow Water Multistrato con Scambi

Problematica
I Modelli di Circolazione Oceanica Globale (OGCM) affrontano un severo vincolo computazionale: il passo temporale è limitato dalla velocità dell'onda più rapida, specificamente le onde di gravità superficiali (circa 200 m/s), sebbene i fenomeni primari di interesse (stratificazione e correnti) si evolvano su scale temporali molto più lente (circa 3 m/s). Sebbene siano stati utilizzati approcci precedenti come l'approssimazione del pelo libero rigido o metodi per l'altezza superficiale implicita, questi presentano dei limiti. La scomposizione barotropica-baroclinica offre una soluzione separando le rapide onde di gravità superficiali (barotropiche) dai più lenti aggiustamenti verticali (baroclinici). Tuttavia, i metodi di scomposizione esistenti spesso faticano a mantenere la coerenza in framework multistrato con scambi di massa, a preservare le proprietà energetiche e a garantire la stabilità senza perdere accuratezza, in particolare in coordinate che seguono il terreno (terrain-following).

Metodologia
Gli autori propongono un'analisi numerica di una scomposizione barotropica-baroclinica all'interno di un framework multistrato non lineare con scambi di massa, formulata come una scomposizione esatta degli operatori. Il modello si basa sulle equazioni di Euler non viscose con un tracciante passivo in coordinate $\sigma$ (terrain-following), assumendo densità costante.

La metodologia centrale consiste nel decomporre le equazioni di Euler medie per strato in due distinti operatori:

1. Passaggio Barotropico: Questo passaggio gestisce l'evoluzione dell'altezza del pelo libero ($h$) e della velocità media di profondità ($\bar{u}$). È formulato come un modello di shallow water a uno strato ben bilanciato (well-balanced) che include termini di pressione e termini sorgente di topografia, ma esclude il forcing baroclinico. Fondamentalmente, questo passaggio tratta il trasporto delle singole velocità di strato ($u_\alpha$) e dei traccianti ($T_\alpha$) alla velocità media $\bar{u}$.
2. Passaggio Baroclinico: Questo passaggio gestisce gli scambi verticali tra gli strati e l'aggiustamento delle velocità rispetto alla media. Risolve la deviazione dalla velocità media ($u_\alpha - \bar{u}$) e impone il vincolo per cui le altezze degli strati rimangano proporzionali all'altezza totale ($h_\alpha = l_\alpha h$) tramite moltiplicatori di Lagrange (termini di scambio $G_{\alpha \pm 1/2}$). Questo sistema è completamente 3D ma privo di termini di pressione e batimetria, consentendo un passo temporale più ampio.

Lo schema numerico impiega una scomposione di Lie con una strategia di sub-ciclo (subcycling):

• Il passaggio baroclinico viene risolto su un grande passo temporale $\Delta t$.
• Il passaggio barotropico viene risolto sullo stesso intervallo utilizzando $K$ piccoli sotto-passi $\delta t$, vincolati dalla velocità dell'onda di gravità.
• Un meccanismo di correzione assicura che la velocità media predetta dal passaggio barotropico coincida con la media delle singole velocità aggiornate nel passaggio baroclinico.

Contributi Chiave

• Formulazione di Scomposizione Esatta degli Operatori: L'articolo riformula la scomposizione come una somma esatta di operatori ($M_L = BT + BC$), chiarificando la struttura matematica del sistema multistrato con scambi.
• Preservazione delle Proprietà Fisiche: Gli autori dimostrano che la scomposizione preserva la conservazione dell'energia totale (o la dissipazione) e soddisfa una disuguaglianza di entropia discreta. Nello specifico, il passaggio barotropico soddisfa un'uguaglianza energetica, mentre il passaggio baroclinico soddisfa una disuguaglianza energetica a causa della natura dissipativa dei termini di scambio.
• Principi di Massimo Discreto: Si dimostra che lo schema soddisfa principi di massimo discreto per la velocità e le concentrazioni dei traccianti, garantendo che i valori rimangano entro limiti fisici.
• Ben Bilanciamento (Well-Balancing): Utilizzando un risolutore well-balanced (ricostruzione idrostatica) per il passaggio barotropico, lo schema complessivo preserva l'equilibrio "lago a riposo" e, quando esteso con recenti strategie di well-balancing lineare, mantiene l'equilibrio geostrofico su mesh grossolane.
• Analisi dell'Iperbolicità: Il lavoro investiga l'iperbolicità dei sistemi scomposti, dimostrando che il sistema barotropico è iperbolico e analizzando gli autovalori del sistema baroclinico, notando che rimane strettamente iperbolico sotto specifiche condizioni riguardanti le deviazioni di velocità.

Risultati
Gli esperimenti numerici dimostrano l'efficacia dello schema proposto:

• Efficienza Computazionale: Nelle simulazioni a basso numero di Froude (dove il rapporto tra velocità del flusso e velocità dell'onda di gravità è piccolo), il metodo di scomposizione riduce drasticamente il costo computazionale rispetto agli schemi non scomposti. Il guadagno è proporzionale al numero di Froude.
• Accuratezza: Nonostante la scomposizione e il sub-ciclo, il metodo mantiene un'accuratezza comparabile agli schemi non scomposti. Gli studi di convergenza rispetto a soluzioni analitiche delle equazioni di Euler idrostatiche mostrano una convergenza del primo ordine.
• Caratteristiche di Diffusione: In determinati casi di test (ad esempio, una cavità guidata dal vento), il metodo di scomposizione è risultato meno diffusivo rispetto allo schema non scomposto, sebbene in altri casi (soluzioni stazionarie) abbia mostrato una diffusione leggermente superiore.
• Equilibrio Geostrofico: L'applicazione di strategie di well-balancing lineare nel passaggio barotropico migliora significativamente la risoluzione delle correnti di confine occidentali in integrazioni a lungo termine, catturando correnti più strette e veloci rispetto agli schemi standard su mesh grossolane.

Significato e Rivendicazioni
L'articolo sostiene che questo lavoro fornisce un'analisi numerica rigorosa di una scomposizione barotropica-baroclinica che è sia computazionalmente efficiente che fisicamente coerente. Il significato primario risiede nella capacità di simulare la circolazione oceanica a lungo termine con grandi passi temporali (limitati dalle onde barocline lente piuttosto che dalle rapide onde di gravità) senza sacrificare la conservazione dell'energia o la preservazione di equilibri critici (lago a riposo e geostrofico). Gli autori sottolineano che la strategia di well-balancing utilizzata nel passaggio barotropico viene ereditata dallo schema complessivo, permettendo simulazioni accurate di flussi geostrofici su griglie grossolane. Il metodo è presentato come un framework robusto sia per applicazioni globali che regionali dove le onde del pelo libero non sono il focus primario.