Barotropic-Baroclinic Splitting for Multilayer Shallow Water Models with Exchanges

Technische Samenvatting: Barotrope-barocline splitsing voor meerlaagse ondiepe watermodellen met uitwisselingen

Probleemstelling
Modellen voor de wereldwijde oceaancirculatie (OGCMs) worden geconfronteerd met een ernstige computationele beperking: de tijdstap wordt beperkt door de snelste golfsnelheid, specifiek de oppervlakte-zwaartekrachtgolven (ongeveer 200 m/s), terwijl de primaire fenomenen van belang (stratificatie en stromingen) zich op veel tragere tijdschalen ontwikkelen (ongeveer 3 m/s). Hoewel eerdere benaderingen zoals de rigide deksel-benadering (rigid lid) of impliciete oppervlaktehoogte-methoden zijn gebruikt, hebben deze beperkingen. De barotrope-barocline decompositie biedt een oplossing door snelle oppervlakte-zwaartekrachtgolven (barotroop) te scheiden van tragere verticale aanpassingen (baroclien). Echter, bestaande splitsingsmethoden hebben vaak moeite met het handhaven van consistentie in meerlaagse kaders met massa-uitwisselingen, het behoud van energie-eigenschappen en het waarborgen van stabiliteit zonder nauwkeurigheid te verliezen, met name in terreinvolgende coördinaten.

Methodologie
De auteurs stellen een numerieke analyse voor van een barotrope-barocline splitsing binnen een niet-lineair meerlaags kader met massa-uitwisselingen tussen lagen, geformuleerd als een exacte operator-splitsing. Het model is gebaseerd op onviskeuze Euler-vergelijkingen met een passieve tracer in terreinvolgende ($\sigma$) coördinaten, waarbij een constante dichtheid wordt aangenomen.

De kernmethodologie omvat het decomponeren van de laaggemiddelde Euler-vergelijkingen in twee afzonderlijke operatoren:

1. Barotrope stap: Deze stap behandelt de evolutie van de vrije oppervlaktehoogte ($h$) en de dieptegemiddelde snelheid ($\bar{u}$). Het is geformuleerd als een well-balanced eenlaags ondiep watermodel dat druktermen en topografie-brontermen bevat, maar barocline forcering uitsluit. Cruciaal is dat deze stap het transport van individuele laag-snelheden ($u_\alpha$) en tracers ($T_\alpha$) behandelt bij de gemiddelde snelheid $\bar{u}$.
2. Barocline stap: Deze stap beheert de verticale uitwisselingen tussen lagen en de aanpassing van snelheden ten opzichte van het gemiddelde. Het lost de afwijking van de gemiddelde snelheid op ($u_\alpha - \bar{u}$) en dwingt de beperking af dat laaghoogtes proportioneel blijven aan de totale hoogte ($h_\alpha = l_\alpha h$) via Lagrange-multiplicatoren (uitwisselingstermen $G_{\alpha \pm 1/2}$). Dit systeem is volledig 3D maar bevat geen druk- en bathymetrie-termen, waardoor een grotere tijdstap mogelijk is.

Het numerieke schema maakt gebruik van een Lie-splitsing met een subcycling-strategie:

• De barocline stap wordt opgelost op een grote tijdstap $\Delta t$.
• De barotrope stap wordt op hetzelfde interval opgelost met behulp van $K$ kleine sub-stappen $\delta t$, die beperkt worden door de zwaartekrachtgolfsnelheid.
• Een correctiemechanisme zorgt ervoor dat de door de barotrope stap voorspelde gemiddelde snelheid overeenkomt met het gemiddelde van de individuele snelheden die in de barocline stap zijn bijgewerkt.

Belangrijkste Bijdragen

• Exacte Operator-splitsing Formulering: Het artikel herformuleert de splitsing als een exacte operator-som ($M_L = BT + BC$), wat de wiskundige structuur van het meerlaagse systeem met uitwisselingen verheldert.
• Behoud van Fysische Eigenschappen: De auteurs bewijzen dat de splitsing het behoud van totale energie (of dissipatie) waarborgt en een discrete entropie-ongelijkheid bevredigt. Specifiek voldoet de barotrope stap aan een energie-gelijkheid, terwijl de barocline stap aan een energie-ongelijkheid voldoet vanwege de dissipatieve aard van de uitwisselingstermen.
• Discrete Maximumprincipes: Er wordt aangetoond dat het schema voldoet aan discrete maximumprincipes voor snelheid en tracerconcentraties, wat garandeert dat waarden binnen fysieke grenzen blijven.
• Well-Balancing: Door een well-balanced solver (hydrostatische reconstructie) te gebruiken voor de barotrope stap, behoudt het totale schema het "meer-in-rust" (lake-at-rest) evenwicht en behoudt het, wanneer het wordt uitgebreid met recente lineaire well-balancing strategieën, geostrofisch evenwicht op grove roosters.
• Hyperboliciteitsanalyse: Het artikel onderzoekt de hyperboliciteit van de gesplitste systemen, bewijst dat het barotrope systeem hyperbolisch is en analyseert de eigenwaarden van het barocline systeem, waarbij wordt opgemerkt dat dit strikt hyperbolisch blijft onder specifieke condities met betrekking tot snelheid-afwijkingen.

Resultaten
Numerieke experimenten demonstreren de effectiviteit van het voorgestelde schema:

• Computationele Efficiëntie: In simulaties met een laag Froude-getal (waarbij de ratio van de stromingssnelheid tot de zwaartekrachtgolfsnelheid klein is), vermindert de splitsingsmethode de computationele kosten drastisch vergeleken met ongesplitste schema's. De winst is proportioneel aan het Froude-getal.
• Nauwkeurigheid: Ondanks de splitsing en subcycling behoudt de methode een nauwkeurigheid die vergelijkbaar is met ongesplitste schema's. Convergentiestudies tegen analytische oplossingen van de hydrostatische Euler-vergelijkingen tonen een eerste-orde convergentie aan.
• Diffusiekarakteristieken: In bepaalde testgevallen (bijv. windgedreven holte) werd de splitsingsmethode minder diffusief gevonden dan het ongesplitste schema, hoewel de methode in andere gevallen (stationaire oplossingen) iets hogere diffusie vertoonde.
• Geostrofisch Evenwicht: De toepassing van lineaire well-balancing strategieën in de barotrope stap verbetert de resolutie van westelijke randstromingen in langdurige integraties aanzienlijk, waarbij nauwere en snellere stromingen worden vastgelegd vergeleken met standaard schema's op grove roosters.

Betekenis en Claims
Het artikel claimt dat dit werk een rigoureuze numerieke analyse biedt van een barotrope-barocline splitsing die zowel computationeel efficiënt als fysisch consistent is. De primaire betekenis ligt in het vermogen om langdurige oceaancirculatie te simuleren met grote tijdstappen (beperkt door trage barocline golven in plaats van snelle zwaartekrachtgolven) zonder het behoud van energie of de preservatie van kritieke evenwichten (meer-in-rust en geostrofisch) op te offeren. De auteurs benadrukken dat de well-balancing strategie die in de barotrope stap wordt gebruikt, wordt overgenomen door het totale schema, wat nauwkeurige simulaties van geostrofische stromingen op grove roosters mogelijk maakt. De methode wordt gepresenteerd als een robuust framework voor zowel globale als regionale toepassingen waarbij oppervlaktegolven niet het primaire focuspunt zijn.