Stochastic compressible Navier-Stokes equations under location uncertainty and their approximations for ocean modelling

Dit artikel presenteert een gezamenlijke theoretische en numerieke studie van stochastische comprimeerbare Navier-Stokes-vergelijkingen onder locatie-onzekerheid, waarbij hun toepassing in oceaanmodellen wordt aangetoond via Boussinesq-benaderingen die significante compressie-effecten in potentiële energie onthullen en een verbeterde energetische consistentie bieden voor subgrid-schaal verticale mengingsmodellen.

De kern

Stel je de oceaan voor als een gigantische, kolkende soeppot. Al lang proberen wetenschappers het "recept" te schrijven voor hoe deze soep beweegt, met behulp van perfecte, deterministische regels. Zij gingen ervan uit dat als je de exacte positie en snelheid van elk watermolecuul kende, je precies kon voorspellen waar de soep als volgende naartoe zou gaan.

De oceaan is echter te complex daarvoor. Hij zit vol met kleine, chaotische draaikolken (turbulentie) die te klein zijn om direct te zien of te meten. Proberen elke enkele draaikolk te simuleren, is als proberen elke zandkorrel op een strand te tellen om de getijden te voorspellen – het is onmogelijk en computertechisch te duur.

Dit artikel stelt een nieuwe manier voor om het recept te schrijven. In plaats van te proberen elke zandkorrel te volgen, stellen de auteurs voor een "correctiefactor" toe te voegen die gebaseerd is op onzekerheid. Zij noemen dit het Location Uncertainty (LU)-kader.

Hier is de kernidee opgesplitst in eenvoudige concepten:

1. De Analogie van de "Drunkards Walk" (Dronkenmansstap)

Stel je voor dat je door een drukke markt loopt. Je hebt een duidelijk bestemming (de "grootschalige" stroming), maar mensen die tegen je aanlopen duwen je lichtjes van koers in willekeurige richtingen.

• De Oude Manier: Je probeert het exacte pad van elke persoon die tegen je aanloopt te berekenen.
• De Nieuwe Manier (LU): Je accepteert dat je gestoten zal worden. Je modelleert je beweging als een gladde wandeling plus een willekeurige, trillende "Brownse beweging" (zoals een drunkards walk). Je weet niet precies waar de stoten je zullen duwen, maar je kent de statistieken van de stoten (hoe sterk ze zijn en hoe ze correleren).

2. De "Compressibele" Soep

De meeste oceaanmodellen gaan ervan uit dat water "oncompressibel" is – wat betekent dat het als een vast blok gelatine is dat niet kan worden samengedrukt. Maar in werkelijkheid kan water wel lichtjes worden samengedrukt, vooral wanneer de druk verandert of de temperatuur verschuift.

• De auteurs beginnen met de volledige, complexe fysica van compressibel water (water dat kan worden samengedrukt).
• Zij passen vervolgens hun wiskunde van "willekeurige stoten" toe op dit complexe systeem.
• Het Resultaat: Zij leiden een nieuwe set vergelijkingen af die eruitzien als de oude, maar met extra termen. Deze extra termen vertegenwoordigen het "werk" dat wordt verricht door de willekeurige stoten. Denk hierbij aan de energie die wordt overgedragen wanneer de menigte je duwt; het is niet zomaar een willekeurige duw, het verandert daadwerkelijk je snelheid en de warmte van je lichaam.

3. De "Verborgen Warmte" in de Mix

Het artikel richt zich sterk op temperatuur en convectie (warm water stijgt, koud water zakt).

• Het Probleem: In standaardmodellen stopt koud water dat zakt vaak abrupt aan de onderkant van het gemengde laagje (het bovenste deel van de oceaan). In werkelijkheid prikken deze "pluimen" van koud water vaak, als een speer, door tot in het diepere, warmere water. Dit wordt penetrerende convectie genoemd.
• De Ontdekking: Toen de auteurs hun nieuwe stochastische model draaiden, ontdekten zij dat de termen voor "willekeurige stoten" dit doordringende effect op natuurlijke wijze nabootsten.
• De Metafoor: Stel je een menigte mensen (de oceaan) voor die probeert een zware doos (een koudwaterpluim) te verplaatsen. Standaardmodellen gedragen zich als een stijve muur die de doos stopt. Het nieuwe model gedraagt zich als een chaotische menigte; de willekeurige duwtjes geven de doos genoeg extra momentum om door de muur te glijpen en dieper te gaan dan verwacht.

4. Twee Manieren om Energie te Meten

De auteurs vonden iets interessants over hoe zij de energie van het systeem maten:

• Interne Energie (De "Warmte"): Toen zij alleen naar de warmte keken, waren de effecten van "samendrukken" (compressie) miniem en niet erg belangrijk. Dit kwam overeen met de oude, eenvoudigere modellen.
• Potentiële Energie (De "Hoogte"): Maar toen zij keken naar de energie gerelateerd aan hoogte (hoe hoog het water staat in het zwaartekrachtsveld), werden de effecten van "samendrukken" zeer belangrijk.
• De Conclusie: Het is als het meten van een stuiterende bal. Als je alleen meet hoe heet de bal wordt wanneer hij op de vloer landt, lijkt de stuiter niet belangrijk. Maar als je meet hoe hoog hij stuitert, is de impact enorm. De auteurs ontdekten dat de willekeurige druktermen in hun model werken als een verborgen veer, die beïnvloedt hoe hoog het water "stuiter" in de energiebegroting.

5. De "Drift" en de "Diffusie"

De wiskunde produceert twee specifieke nieuwe termen die fungeren als de "correctiefactoren":

• De Drift (Itô-Stokes Drift): Dit is een systematische duw veroorzaakt door het feit dat de willekeurige stoten niet perfect uniform zijn. Het is als een rivierstroom die iets anders stroomt omdat de rotsen (turbulentie) in een specifiek patroon zijn gerangschikt.
• De Diffusie: Dit is het verspreidingseffect veroorzaakt door de willekeurige stoten.

Samenvatting van de Prestatie

De auteurs hebben succesvol een brug gebouwd tussen de rommelige, chaotische realiteit van de oceaan en de schone, wiskundige modellen die we gebruiken om deze te voorspellen.

• Zij begonnen met de meest complexe fysica die mogelijk is (compressibel, willekeurig, thermodynamisch).
• Zij toonden aan dat wanneer je het terugbrengt tot het "standaard" oceaanbeeld (Boussinesq-benadering), je nieuwe vergelijkingen nog steeds werken en eigenlijk de voorspelling van hoe diep koud water zakt, verbeteren.
• Zij bewezen dat je niet elke kleine draaikolk hoeft te simuleren om het juiste antwoord te krijgen; je hoeft alleen maar wiskundig rekening te houden met de onzekerheid van waar het water naartoe gaat.

Kortom, zij vervingen de onmogelijke taak van "elke zandkorrel tellen" door een slimmere strategie: "rekening houden met de neiging van het zand om te verspreiden", en ontdekten dat deze benadering de diepe, doordringende stromingen van de oceaan veel beter vastlegt dan eerdere methoden.

Technische samenvatting

Technische Samenvatting: Stochastische Compressibele Navier-Stokes-vergelijkingen en Hun Benaderingen voor Oceaanmodellering

Probleemstelling
Oceaanmodellering staat voor aanzienlijke uitdagingen bij het weergeven van processen op subgrid-schaal, met name verticale menging en penetrerende convectie, die cruciaal zijn voor klimaatsimulaties (bijvoorbeeld de Atlantische Meridionale Omwentelingscirculatie). Traditionele benaderingen vertrouwen vaak op de Boussinesq-benadering, die oncompressibiliteit veronderstelt en compressie-effecten verwaarloost. Hoewel dit in veel contexten redelijk is, kan deze benadering de energetische consistentie tussen temperatuur- en impulsfluxen verduisteren en faalt hij in het vastleggen van de voortplanting van akoestische golven of specifieke thermodynamische arbeidstermen geassocieerd met compressie en dilatatie. Bovendien zijn bestaande stochastische modellen voor oceaan-dynamica grotendeels afgeleid voor oncompressibele stromingen, waardoor er een gat bestaat in het theoretisch kader voor compressibele, thermodynamisch actieve vloeistoffen. De auteurs beogen een rigoureuze stochastische formulering te ontwikkelen die locatie-onzekerheid (LU) in compressibele stromingen in rekening brengt, door expliciet de nodige thermodynamische termen af te leiden om de energetische consistentie van verticale modellen op subgrid-schaal te verbeteren.

Methodologie
Het artikel maakt gebruik van het Location Uncertainty (LU)-kader, dat de verplaatsing van vloeistofdeeltjes decomposeert in een in de tijd differentieerbare opgeloste snelheid en een sterk fluctuerend, in de tijd gedecorrelleerd component dat wordt gemodelleerd als een Brownse beweging.

1. Theoretische Afleiding:

   • De auteurs vertrekken vanuit de fundamentele behoudswetten voor massa, soortmassa, impuls en totale energie.
   • Zij passen de Stochastische Reynolds-Transportstelling (SRTT) toe, uitgebreid met stochastische brontermen. Dit leidt tot een algemeen systeem van stochastische compressibele Navier-Stokes-vergelijkingen in zowel conservatieve als niet-conservatieve (primitieve variabele) vormen.
   • Cruciaal is dat de afleiding expliciet rekening houdt met het werk van krachten (druk en viskeuze spanningen) en warmtestromen die inwerken op de stochastische verplaatsing. Dit resulteert in het ontstaan van specifieke stochastische termen, waaronder door ruis geïnduceerde diffusie, Itô-Stokes-drift en kwadratische covariatietermen (voortkomend uit het product van stochastische differentiaalvergelijkingen).
   • Het systeem wordt vervolgens onderworpen aan de Boussinesq- en hydrostatische benaderingen. De auteurs volgen zorgvuldig de orde van grootte van termen om een stochastisch Boussinesq-systeem af te leiden dat thermodynamische effecten behoudt, specifiek het werk van compressie uitgedrukt via potentiële energie.

2. Numerieke Toepassing:

   • Om het theoretisch kader te valideren, voeren de auteurs een Large Eddy Simulation (LES) uit van een door temperatuur gedreven vrije convectiegebeurtenis in een gestratificeerde oceaan met behulp van het CROCO-model (niet-hydrostatische, niet-Boussinesq-mogelijkheden).
   • Zij passen een tijdsfiltertechniek toe op de LES-gegevens (met een decorrelatietijdschaal van 3 uur) om opgeloste schalen te scheiden van fluctuaties op subgrid-schaal.
   • De residuen op subgrid-schaal worden gebruikt om de stochastische ruis (Brownse beweging) en zijn variantietensor te construeren.
   • Een offline 1D-stochastisch model wordt vervolgens aangedreven door deze gediagnosticeerde statistieken om de evolutie van het horizontaal gemiddelde temperatuurprofiel te simuleren. Dit maakt een directe vergelijking mogelijk tussen de output van het stochastische model en de referentie-LES-resultaten.

Belangrijkste Bijdragen

• Algemene Stochastische Compressibele Formulering: Het artikel biedt een volledige afleiding van de stochastische compressibele Navier-Stokes-vergelijkingen, inclusief het transport van soorten en de volledige energiebalans. Dit breidt eerdere LU-werk uit, die voornamelijk beperkt was tot oncompressibele stromingen.
• Thermodynamische Consistentie: De afleiding omvat expliciet het werk van compressie en dilatatie. De auteurs tonen aan dat hoewel deze termen verwaarloosbaar zijn wanneer ze worden uitgedrukt in termen van inwendige energie (consistent met de Boussinesq-benadering), ze significant worden wanneer ze worden uitgedrukt in termen van potentiële energie. Dit onthult nieuwe stochastische druktermen binnen het gemengde laag die kinetische en potentiële energie koppelen.
• Stochastisch Boussinesq-model met Thermodynamica: De auteurs leiden een stochastisch Boussinesq-model af dat deze thermodynamische effecten integreert, waardoor een fysisch consistenter kader wordt geboden dan eerdere stochastische isochore modellen.
• Diagnose van Penetrerende Convectie: De studie identificeert dat de stochastische termen, met name de Itô-Stokes-drift en stochastische diffusie, in staat zijn de verticale temperatuurfluxen aan de basis van het gemengde laag te reproduceren. Dit is een gebied waar standaardparameteriseringen vaak falen in het vastleggen van effecten van penetrerende convectie (opwarming aan de basis van het gemengde laag).

Resultaten

• Energiebalans: In het numerieke experiment wordt de leidende orde balans voor temperatuurveranderingen beheerst door transport van inwendige energie, consistent met de Boussinesq-benadering. Echter, bij analyse van potentiële energie tonen de stochastische termen (specifiek de kwadratische covariatie $A_u$ en driftwerk $D_\sigma$ geassocieerd met stochastische druk) aanzienlijke bijdragen.
• Dynamica van het Gemengde Laag: Het stochastische model reproduceert succesvol de opwarming aan de basis van het gemengde laag geassocieerd met penetrerende convectie. Dit wordt toegeschreven aan het vermogen van het model om niet-lokale verticale dynamiek vast te leggen die wordt gedreven door horizontale variabiliteit en inertieprocessen, zonder te vertrouwen op expliciete diffusive sluitingsaannames.
• Beperkingen: Binnen het gemengde laag zelf toont het stochastische model beperkte vaardigheid in vergelijking met de referentie-LES. De auteurs merken op dat de relatief lage variabiliteit van het turbulente veld in dit gebied leidt tot verwaarloosbare bijdragen van de stochastische termen in hun huidige diagnostische opstelling. Het model worstelt ook nabij het oppervlak, waarschijnlijk vanwege het onvermogen van de stochastische representatie om de specifieke kleine schaal fluctuaties vast te leggen die door de LES worden opgelost.

Betekenis en Aanspraken
De auteurs beweren dat dit werk een gestructureerd modelleringkader biedt dat nieuwe perspectieven opent voor oceaan-dynamica binnen een stochastische setting.

• Fysische Rechtvaardiging: In tegenstelling tot fenomenologische benaderingen, worden de extra termen in de stochastische vergelijkingen rigoureus afgeleid uit behoudswetten en stochastische calculus, waardoor een duidelijke fysische rechtvaardiging wordt geboden voor hun vorm en rol (bijvoorbeeld analoog aan baropycnale arbeid in compressibele LES).
• Energetische Consistentie: Het kader biedt een weg om de energetische consistentie van verticale modellen op subgrid-schaal die worden gebruikt in Ocean General Circulation Models (OGCM's) te verbeteren door expliciet rekening te houden met fysieke onzekerheden en benaderingen, met name wat betreft de koppeling tussen dynamica en thermodynamica.
• Penetrerende Convectie: Het vermogen van de stochastische termen om kwalitatief ontbrekende Reynolds-spanningen aan de basis van het gemengde laag te reproduceren, suggereert een potentieel voor het weergeven van effecten van penetrerende convectie zonder de rekenkosten van volledige niet-hydrostatische oplosprogramma's.
• Toekomstige Richtingen: Het artikel suggereert bescheiden dat hoewel de huidige diagnostische benadering veelbelovend is, toekomstig werk zich moet richten op prognostische modellering van de ruiseigenschappen (bijvoorbeeld via TKE-gebaseerde modellen) en de opname van willekeurige oppervlakte-drukkrachten om de dynamica binnen het gemengde laag en nabij het oppervlak beter vast te leggen.