Singular limits for non-isentropic compressible rotating fluids

1. Il Problema Studiato

L'articolo indaga il comportamento asintotico (limiti singolari) delle equazioni di Navier-Stokes-Korteweg per fluidi comprimibili, non isoentropici e rotanti, in un dominio tridimensionale illimitato confinato tra due piani paralleli ($\Omega = \mathbb{R}^2 \times (0, 1)$).

Il sistema fisico è governato da tre parametri adimensionali che tendono a valori limite estremi contemporaneamente:

• Numero di Mach ($Ma$) basso: $Ma \propto \varepsilon \to 0$. Indica che la velocità del fluido è molto inferiore alla velocità del suono, portando verso l'incomprimibilità.
• Numero di Rossby ($Ro$) basso: $Ro \propto \varepsilon \to 0$. Indica una rapida rotazione del sistema, dove le forze di Coriolis dominano, portando verso un comportamento bidimensionale.
• Numero di Reynolds ($Re$) alto: $Re \propto 1/\varepsilon \to \infty$. Indica che gli effetti viscosi diventano trascurabili rispetto alle forze d'inerzia, portando verso l'inviscidità.

Il modello include anche l'effetto capillare (tensione superficiale) tramite il tensore di Korteweg, controllato da un coefficiente $\kappa$. L'analisi distingue due casi fondamentali basati su un parametro $\alpha$:

• Caso $\alpha = 1$: L'effetto capillare è presente ma trattato in modo tale da non generare termini di ordine superiore critici nella scala di $\varepsilon$ (o meglio, il termine di capillarità scala in modo da permettere l'uso di stime di dispersione).
• Caso $\alpha = 0$: L'effetto capillare è dominante nella scala di $\varepsilon$, richiedendo un'analisi diversa basata su stime di errore.

2. Metodologia

Gli autori adottano un approccio analitico rigoroso basato sulla teoria delle equazioni alle derivate parziali non lineari, articolato in due fasi principali a seconda del caso di $\alpha$:

A. Stime Uniformi e Analisi delle Onde Acustiche ($\alpha = 1$)

1. Stime di Energia: Vengono derivati limiti uniformi per le soluzioni forti $(q_\varepsilon, \mathbf{u}_\varepsilon, \theta_\varepsilon)$ (dove $q$ e $\theta$ sono le variazioni di densità e temperatura) in spazi di Sobolev $H^s$. Queste stime sono indipendenti da $\varepsilon$ su un intervallo di tempo finito $[0, T]$.
2. Analisi delle Onde Acustiche: Per gestire i termini che divergono come $1/\varepsilon$ (legati alla compressibilità e alla rotazione), il sistema viene riscritto in termini di onde acustiche.
3. Teorema di RAGE: Viene applicato il teorema di RAGE (Ruelle-Agmon-Rage) per dimostrare la dispersione delle onde acustiche. Questo teorema, applicato all'operatore associato alle onde acustiche, garantisce che le componenti oscillanti ad alta frequenza (onde sonore) si disperdano e convergano a zero in senso debole, lasciando solo la componente geostrofica (lenta).
4. Compattezza: Si utilizzano teoremi di immersione di Sobolev e il teorema di Ascoli-Arzela per estrarre sottosuccessioni convergenti.

B. Stime di Errore ($\alpha = 0$)

1. Soluzione di Riferimento: Si assume l'esistenza di una soluzione regolare $(\mathbf{w}, \pi)$ per il sistema limite (equazioni di Navier-Stokes incompressibili 2D).
2. Costruzione di Soluzioni Approssimate: Si costruisce una soluzione approssimata basata sul sistema limite.
3. Analisi dell'Errore: Si definiscono le variabili di errore come la differenza tra la soluzione del sistema completo e quella approssimata. Attraverso un'analisi energetica dettagliata delle equazioni differenziali per l'errore, si dimostra che l'errore è dell'ordine di $O(\varepsilon)$ in norme di Sobolev appropriate. Questo metodo evita la necessità di stime di dispersione complesse, affidandosi invece alla stabilità della soluzione limite.

3. Risultati Chiave

Caso $\alpha = 1$ (Teorema 2)

• Convergenza: Le soluzioni del sistema comprimibile, rotante e capillare convergono (in senso debole e debole-* in spazi di Sobolev) verso un sistema limite.
• Sistema Limite: Il sistema limite è descritto dalle equazioni di Eulero bidimensionali incompressibili (o una loro variante con struttura simile).
• Struttura Geostrofica: La soluzione limite soddisfa l'equazione di equilibrio geostrofico: $\mathbf{e}_3 \times \mathbf{u} + \nabla \sigma = 0$, dove $\sigma = q + \theta$.
• Dimensionalità: Il flusso diventa bidimensionale ($\mathbf{u}$ e $\sigma$ sono indipendenti dalla coordinata verticale $x_3$) e incomprimibile ($\text{div } \mathbf{u} = 0$).
• Densità e Temperatura: La densità e la temperatura convergono a 1 (stato di riferimento), con fluttuazioni che si annullano nel limite.

Caso $\alpha = 0$ (Teorema 3)

• Convergenza Forte: Viene dimostrata una convergenza forte delle soluzioni.
• Stima dell'Errore: È stato provato che esiste una costante $K$ indipendente da $\varepsilon$ tale che:
  $$\sup_{t \in [0, T^*]} \left( \|q_\varepsilon - \varepsilon^2 \pi\|_{H^{s+1}} + \|\mathbf{u}_\varepsilon - \mathbf{w}\|_{H^s} + \|\theta_\varepsilon - \varepsilon \pi\|_{H^s} \right) \leq K \varepsilon$$
• Sistema Limite: Anche in questo caso, il sistema converge alle equazioni di Navier-Stokes incompressibili bidimensionali, ma la presenza del termine capillare dominante ($\alpha=0$) modifica la struttura dell'errore e la scala delle fluttuazioni di densità e temperatura rispetto al caso $\alpha=1$.

4. Contributi Principali

1. Trattamento Unificato con Capillarità: Il lavoro è uno dei primi a trattare sistematicamente il limite singolare triplo (basso Mach, basso Rossby, alto Reynolds) per fluidi non isoentropici includendo esplicitamente gli effetti capillari (equazioni di Korteweg).
2. Distinzione dei Regimi Capillari: L'articolo distingue chiaramente tra due regimi fisici ($\alpha=1$ e $\alpha=0$) relativi alla scala del coefficiente capillare, fornendo metodi di prova diversi (dispersione vs. stima di errore) per ciascuno.
3. Riduzione Dimensionale Rigorosa: Dimostra matematicamente come la combinazione di rapida rotazione e bassa comprimibilità porti alla degenerazione di un sistema 3D complesso in un sistema 2D incompressibile, anche in presenza di gradienti di temperatura e forze capillari.
4. Analisi di Dati Iniziali "Ill-Prepared": Il lavoro affronta il caso di dati iniziali che non soddisfano le condizioni di equilibrio geostrofico (dati "ill-prepared"), mostrando come le onde acustiche e le oscillazioni di Rossby si disperdano rapidamente, lasciando il sistema a evolvere secondo la dinamica lenta.

5. Significato Scientifico

Questo studio è significativo per la fluidodinamica matematica e geofisica per diverse ragioni:

• Modellistica Atmosferica e Oceanica: Fornisce una giustificazione teorica rigorosa per l'uso di modelli bidimensionali incompressibili (come le equazioni di Eulero o Navier-Stokes 2D) per descrivere fenomeni su larga scala nell'atmosfera e negli oceani, dove la rotazione terrestre e la bassa velocità del flusso sono dominanti.
• Fluidi Complessi: L'inclusione degli effetti capillari è cruciale per la comprensione di fluidi in prossimità di punti critici, transizioni di fase, o in microfluidica, dove le forze di superficie giocano un ruolo fondamentale.
• Validazione Numerica: I risultati forniscono un benchmark teorico per la validazione di schemi numerici che simulano flussi comprimibili ad alta velocità di rotazione, confermando che tali schemi devono recuperare correttamente la dinamica 2D incompressibile nel limite.
• Avanzamento Teorico: Estende la teoria dei limiti singolari (sviluppata precedentemente per fluidi isoentropici o senza capillarità) a contesti più realistici e fisicamente complessi, risolvendo le difficoltà analitiche poste dall'accoppiamento tra non-linearità, compressibilità, rotazione e tensione superficiale.

In sintesi, il paper conferma che, nonostante la complessità delle equazioni di Navier-Stokes-Korteweg non isoentropiche, nel limite di Mach basso, Rossby basso e Reynolds alto, la fisica del sistema si semplifica drasticamente, riducendosi a una dinamica bidimensionale incompressibile governata da leggi di conservazione della vorticità.