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

Questo lavoro presenta uno studio teorico e numerico congiunto delle equazioni di Navier-Stokes comprimibili stocastiche in condizioni di incertezza di posizione, dimostrandone l'applicazione alla modellazione oceanica attraverso approssimazioni di Boussinesq che rivelano effetti significativi di compressione nell'energia potenziale e offrono una migliore coerenza energetica per i modelli di mescolamento verticale su scala subgriglia.

L'essenziale

Immagina l'oceano come una gigantesca pentola di zuppa in ebollizione. Da lungo tempo, gli scienziati hanno cercato di scrivere la "ricetta" per il movimento di questa zuppa utilizzando regole perfette e deterministiche. Hanno assunto che, se si conoscesse la posizione esatta e la velocità di ogni singola molecola d'acqua, si potesse prevedere esattamente dove la zuppa si sarebbe mossa successivamente.

Tuttavia, l'oceano è troppo complesso per questo. È pieno di piccoli vortici caotici (turbolenza) troppo piccoli per essere visti o misurati direttamente. Tentare di simulare ogni singolo vortice è come cercare di contare ogni granello di sabbia su una spiaggia per prevedere la marea: è impossibile e computazionalmente troppo costoso.

Questo articolo propone un nuovo modo per scrivere la ricetta. Invece di cercare di tracciare ogni granello di sabbia, gli autori suggeriscono di aggiungere un "fattore di aggiustamento" basato sull'incertezza. Chiamano questo il framework Location Uncertainty (LU) (Incertezza di Posizione).

Ecco l'idea centrale scomposta in concetti semplici:

1. L'analogia del "Camminatore Ubriaco"

Immagina di camminare attraverso un mercato affollato. Hai una destinazione chiara (il flusso "su larga scala"), ma le persone che ti urtano ti spingono leggermente fuori rotta in direzioni casuali.

• Il vecchio modo: Cerchi di calcolare il percorso esatto di ogni persona che ti urta.
• Il nuovo modo (LU): Accetti di essere urtato. Modelli il tuo movimento come una camminata fluida più un movimento casuale e tremolante "moto browniano" (come il camminare di un ubriaco). Non sai esattamente dove gli urti ti spingeranno, ma conosci le statistiche degli urti (quanto sono forti e come sono correlati).

2. La zuppa "Comprimibile"

La maggior parte dei modelli oceanici assume che l'acqua sia "incomprimibile", il che significa che è come un blocco solido di gelatina che non può essere schiacciato. Ma in realtà, l'acqua può essere leggermente schiacciata, specialmente quando cambiano la pressione o la temperatura.

• Gli autori partono dalla fisica completa e complessa dell'acqua comprimibile (acqua che può essere schiacciata).
• Applicano poi la loro matematica degli "urti casuali" a questo sistema complesso.
• Il Risultato: Derivano un nuovo insieme di equazioni che assomigliano a quelle vecchie ma includono termini aggiuntivi. Questi termini extra rappresentano il "lavoro" compiuto dagli urti casuali. Pensa all'energia trasferita quando la folla ti spinge; non è solo una spinta casuale, ma cambia effettivamente la tua velocità e il calore del tuo corpo.

3. Il "Calore Nascosto" nel Mix

L'articolo si concentra pesantemente sulla temperatura e sulla convezione (acqua calda che sale, acqua fredda che scende).

• Il Problema: Nei modelli standard, quando l'acqua fredda scende, spesso si ferma bruscamente alla base dello strato misto (la parte superiore dell'oceano). In realtà, questi "plumi" di acqua fredda spesso penetrano, come una lancia, nell'acqua più calda e profonda. Questo è chiamato convezione penetrativa.
• La Scoperta: Quando gli autori hanno eseguito il loro nuovo modello stocastico, hanno scoperto che i termini degli "urti casuali" ricreavano naturalmente questo effetto di penetrazione.
• La Metafora: Immagina una folla di persone (l'oceano) che cerca di spostare una scatola pesante (un plume di acqua fredda). I modelli standard agiscono come un muro rigido che ferma la scatola. Il nuovo modello agisce come una folla caotica; lo spintonamento casuale dà alla scatola abbastanza momento extra per scivolare attraverso il muro e andare più in profondità del previsto.

4. Due Modi per Misurare l'Energia

Gli autori hanno scoperto qualcosa di interessante su come misuravano l'energia del sistema:

• Energia Interna (Il "Calore"): Quando hanno guardato solo il calore, gli effetti dello "schiacciamento" (compressione) erano minuscoli e non contavano molto. Questo corrispondeva ai vecchi modelli più semplici.
• Energia Potenziale (L'"Altezza"): Ma quando hanno guardato l'energia relativa all'altezza (quanto è alta l'acqua nel campo gravitazionale), gli effetti dello "schiacciamento" sono diventati molto importanti.
• La Conclusione: È come misurare una palla che rimbalza. Se misuri solo quanto la palla si scalda quando colpisce il pavimento, il rimbalzo non sembra importante. Ma se misuri quanto in alto rimbalza, l'impatto è enorme. Gli autori hanno scoperto che i termini di pressione casuali nel loro modello agiscono come una molla nascosta, influenzando quanto in alto l'acqua "rimbalza" nel bilancio energetico.

5. La "Deriva" e la "Diffusione"

La matematica produce due nuovi termini specifici che agiscono come i "fattori di aggiustamento":

• La Deriva (Deriva di Itô-Stokes): Questa è una spinta sistematica causata dal fatto che gli urti casuali non sono perfettamente uniformi. È come una corrente fluviale che scorre leggermente diversamente perché le rocce (turbolenza) sono disposte in uno schema specifico.
• La Diffusione: Questo è l'effetto di dispersione causato dagli urti casuali.

Sintesi del Raggiungimento

Gli autori hanno costruito con successo un ponte tra la realtà disordinata e caotica dell'oceano e i modelli matematici puliti che usiamo per prevederlo.

• Hanno iniziato con la fisica più complessa possibile (comprimibile, casuale, termodinamica).
• Hanno dimostrato che quando la si semplifica fino alla "visione standard" dell'oceano (approssimazione di Boussinesq), le nuove equazioni funzionano ancora e in realtà migliorano la previsione di quanto in profondità scende l'acqua fredda.
• Hanno dimostrato che non è necessario simulare ogni piccolo vortice per ottenere la risposta giusta; basta tenere conto matematicamente dell'incertezza di dove sta andando l'acqua.

In breve, hanno sostituito il compito impossibile di "contare ogni granello di sabbia" con una strategia più intelligente: "tenere conto della tendenza della sabbia a disperdersi", e hanno scoperto che questo approccio cattura le correnti profonde e penetranti dell'oceano molto meglio dei metodi precedenti.

Sommario tecnico

Sintesi Tecnica: Equazioni di Navier-Stokes Compressibili Stocastiche e Loro Approssimazioni per la Modellazione Oceanica

Enunciato del Problema
La modellazione oceanica affronta sfide significative nella rappresentazione dei processi a scala subgriglia, in particolare il rimescolamento verticale e la convezione penetrativa, che sono fondamentali per le simulazioni climatiche (ad esempio, la Circolazione Meridionale di Inversione Atlantica). Gli approcci tradizionali si basano spesso sull'approssimazione di Boussinesq, che assume l'incomprimibilità e trascura gli effetti di comprimibilità. Sebbene ragionevole in molti contesti, tale approssimazione può oscurare la coerenza energetica tra i flussi di temperatura e di quantità di moto e non riesce a catturare la propagazione delle onde acustiche o specifici termini di lavoro termodinamico associati alla compressione e alla dilatazione. Inoltre, i modelli stocastici esistenti per la dinamica oceanica sono stati derivati prevalentemente per flussi incomprimibili, lasciando un vuoto nel quadro teorico per i fluidi comprimibili e termodinamicamente attivi. Gli autori mirano a sviluppare una formulazione stocastica rigorosa che tenga conto dell'incertezza di posizione (LU) nei flussi comprimibili, derivando esplicitamente i termini termodinamici necessari per migliorare la coerenza energetica dei modelli verticali a scala subgriglia.

Metodologia
Il documento utilizza il quadro dell'Incertezza di Posizione (LU), che scompone lo spostamento della particella fluida in una velocità risolta differenziabile nel tempo e in una componente altamente fluttuante, decorrelata nel tempo, modellata come un moto browniano.

1. Derivazione Teorica:

   • Gli autori partono dalle leggi fondamentali di conservazione della massa, della massa delle specie, della quantità di moto e dell'energia totale.
   • Applicano il Teorema di Trasporto di Reynolds Stocastico (SRTT) esteso per includere termini sorgente stocastici. Ciò porta a un sistema generale di equazioni di Navier-Stokes compressibili stocastiche sia in forma conservativa che non conservativa (variabili primitive).
   • Crucialmente, la derivazione tiene esplicitamente conto del lavoro delle forze (pressione e sforzi viscosi) e dei flussi di calore agenti sullo spostamento stocastico. Ciò comporta l'emergere di termini stocastici specifici, inclusa la diffusione indotta dal rumore, la deriva di Itô-Stokes e i termini di covariazione quadratica (derivanti dal prodotto di differenziali stocastici).
   • Il sistema è quindi sottoposto alle approssimazioni di Boussinesq e idrostatica. Gli autori tracciano attentamente l'ordine di grandezza dei termini per derivare un sistema stocastico di Boussinesq che mantiene gli effetti termodinamici, in particolare il lavoro di compressione espresso attraverso l'energia potenziale.

2. Applicazione Numerica:

   • Per validare il quadro teorico, gli autori eseguono una Simulazione a Grande Vortice (LES) di un evento di convezione libera guidata dalla temperatura in un oceano stratificato utilizzando il modello CROCO (capacità non idrostatiche e non di Boussinesq).
   • Applicano una tecnica di filtraggio temporale ai dati LES (utilizzando una scala temporale di decorrelazione di 3 ore) per separare le scale risolte dalle fluttuazioni a scala subgriglia.
   • I residui a scala subgriglia sono utilizzati per costruire il rumore stocastico (moto browniano) e il suo tensore di varianza.
   • Un modello stocastico 1D offline è quindi forzato da queste statistiche diagnosticate per simulare l'evoluzione del profilo di temperatura mediato orizzontalmente. Ciò consente un confronto diretto tra l'output del modello stocastico e i risultati di riferimento LES.

Contributi Chiave

• Formulazione Stocastica Compressibile Generale: Il documento fornisce una derivazione completa delle equazioni di Navier-Stokes compressibili stocastiche, inclusa il trasporto delle specie e il bilancio energetico completo. Ciò estende il lavoro precedente sulla LU, che era principalmente limitato ai flussi incomprimibili.
• Coerenza Termodinamica: La derivazione include esplicitamente il lavoro di compressione e dilatazione. Gli autori dimostrano che, sebbene questi termini siano trascurabili quando espressi in termini di energia interna (coerentemente con l'approssimazione di Boussinesq), diventano significativi quando espressi in termini di energia potenziale. Ciò rivela nuovi termini di pressione stocastica all'interno dello strato misto che accoppiano l'energia cinetica e l'energia potenziale.
• Modello Stocastico di Boussinesq con Termodinamica: Gli autori derivano un modello stocastico di Boussinesq che incorpora questi effetti termodinamici, offrendo un quadro più coerente dal punto di vista fisico rispetto ai precedenti modelli stocastici isocori.
• Diagnostica della Convezione Penetrativa: Lo studio identifica che i termini stocastici, in particolare la deriva di Itô-Stokes e la diffusione stocastica, sono in grado di riprodurre i flussi verticali di temperatura alla base dello strato misto. Questa è una regione dove le parametrizzazioni standard spesso falliscono nel catturare gli effetti della convezione penetrativa (riscaldamento alla base dello strato misto).

Risultati

• Bilancio Energetico: Nell'esperimento numerico, il bilancio di ordine dominante per i cambiamenti di temperatura è governato dal trasporto di energia interna, coerentemente con l'approssimazione di Boussinesq. Tuttavia, analizzando l'energia potenziale, i termini stocastici (in particolare la covariazione quadratica $A_u$ e il lavoro di deriva $D_\sigma$ associati alla pressione stocastica) mostrano contributi significativi.
• Dinamica dello Strato Misto: Il modello stocastico riproduce con successo il riscaldamento alla base dello strato misto associato alla convezione penetrativa. Ciò è attribuito alla capacità del modello di catturare la dinamica verticale non locale guidata dalla variabilità orizzontale e dai processi inerziali senza fare affidamento su assunzioni esplicite di chiusura diffusiva.
• Limiti: All'interno dello strato misto stesso, il modello stocastico mostra abilità limitata rispetto al LES di riferimento. Gli autori notano che la variabilità relativamente bassa del campo turbolento in questa regione porta a contributi trascurabili dei termini stocastici nella loro configurazione diagnostica attuale. Il modello fatica anche vicino alla superficie, probabilmente a causa dell'incapacità della rappresentazione stocastica di catturare le fluttuazioni specifiche a piccola scala risolte dal LES.

Significato e Affermazioni
Gli autori affermano che questo lavoro fornisce un quadro di modellazione strutturato che apre nuove prospettive per la dinamica oceanica in un contesto stocastico.

• Giustificazione Fisica: A differenza degli approcci fenomenologici, i termini aggiuntivi nelle equazioni stocastiche sono derivati rigorosamente dalle leggi di conservazione e dal calcolo stocastico, fornendo una chiara giustificazione fisica per la loro forma e il loro ruolo (ad esempio, analoghi al lavoro baropicinale nella LES comprimibile).
• Coerenza Energetica: Il quadro offre una via per migliorare la coerenza energetica dei modelli verticali a scala subgriglia utilizzati nei Modelli di Circolazione Generale Oceanica (OGCM), tenendo esplicitamente conto delle incertezze fisiche e delle approssimazioni, in particolare riguardo all'accoppiamento tra dinamica e termodinamica.
• Convezione Penetrativa: La capacità dei termini stocastici di riprodurre qualitativamente gli stress di Reynolds mancanti alla base dello strato misto suggerisce un potenziale per rappresentare gli effetti della convezione penetrativa senza il costo computazionale di solver non idrostatici completi.
• Direzioni Future: Il documento suggerisce modestamente che, sebbene l'approccio diagnostico attuale sia promettente, il lavoro futuro dovrebbe concentrarsi sulla modellazione prognostica delle caratteristiche del rumore (ad esempio, tramite modelli basati sull'energia cinetica turbolenta, TKE) e sull'inclusione di forzanti di pressione superficiale casuali per catturare meglio la dinamica all'interno dello strato misto e vicino alla superficie.