Global existence of classical solutions for the multi-dimensional compressible Navier-Stokes-Poisson equations on solid balls for arbitrary spherically symmetric large initial data

Titolo: Esistenza globale di soluzioni classiche per le equazioni di Navier-Stokes-Poisson comprimibili multi-dimensionali su sfere solide per dati iniziali grandi sfericamente simmetrici.

1. Il Problema

Il lavoro affronta una questione aperta di lunga data nella fisica matematica e nell'analisi delle equazioni alle derivate parziali: l'esistenza di soluzioni classiche globali per le equazioni di Navier-Stokes-Poisson (NSP) comprimibili in tre dimensioni (e due dimensioni) per dati iniziali arbitrariamente grandi.

Le equazioni NSP governano il moto di stelle gassose o plasmi in un campo gravitazionale o elettrico auto-consistente. Il sistema è dato da:
$$\begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ \rho(u_t + u \cdot \nabla u) + \nabla P(\rho) = -\rho\nabla\Phi + \text{div}(\mu(\rho)Du) + \nabla(\lambda(\rho)\text{div}u), \\ \Delta\Phi = \kappa(\rho - \bar{\rho}), \end{cases}$$
dove $\rho$ è la densità, $u$ la velocità, $\Phi$ il potenziale (gravitazionale se $\kappa > 0$, elettrico se $\kappa < 0$), e i coefficienti di viscosità $\mu(\rho)$ e $\lambda(\rho)$ dipendono dalla densità.

Sfide principali:

• Dati iniziali grandi: La maggior parte dei risultati precedenti richiede perturbazioni piccole di uno stato stazionario.
• Singolarità al centro: In coordinate sferiche, i termini di Poisson e le derivate radiali introducono singolarità forti all'origine ($r=0$), specialmente in 3D dove il termine di Poisson scala come $x/r^2$.
• Viscosità degenerata: Il lavoro considera viscosità dipendenti dalla densità (modello BD), dove la viscosità può annullarsi nel vuoto (anche se qui si studiano soluzioni "lontane dal vuoto", la struttura della viscosità è cruciale).

2. Metodologia

Gli autori, Fan, Huang e Lei, sviluppano una serie di stime a priori sofisticate per estendere le soluzioni locali a soluzioni globali. La strategia si basa su:

• Simmetria Sferica: Il problema viene ridotto a una dimensione spaziale (radiale) su un dominio fisso $[0, R]$, trasformando le equazioni in coordinate di Lagrange per gestire meglio la densità.
• Relazione BD (Bresch-Desjardins): Sfruttano la relazione specifica tra i coefficienti di viscosità: $\mu(\rho) = \rho^\alpha$ e $\lambda(\rho) = (\alpha-1)\rho^\alpha$. Questa relazione garantisce l'esistenza di una disuguaglianza di entropia BD, che fornisce stime aggiuntive cruciali per il controllo della densità e della velocità.
• Velocità Effettiva: Introducono la variabile "velocità effettiva" $w = u + r^{N-1}(\rho^\alpha)_x$. Questa combinazione cancella parzialmente i termini di pressione e viscosità, permettendo di ottenere stime di regolarità migliori.
• Stime Pesate e Controllo della Singolarità:
  • Il cuore della dimostrazione risiede nel controllo della crescita della densità e del potenziale gravitazionale vicino alla singolarità centrale ($r=0$).
  • Gli autori derivano stime pesate della forma $\|\rho^{\alpha-1/2} r^\xi\|_{L^\infty}$, che permettono di controllare integrali singolari come $\int \frac{x}{r} dx$ o $\int (\frac{x}{r^2})^{2k} dx$.
  • Stabiliscono una relazione funzionale tra la coordinata di Lagrange $x$ e la coordinata Euleriana $r$ (es. $x \leq C r^\beta$), che è essenziale per eliminare le singolarità nei termini di accoppiamento.
• Integrabilità Alta: Costruiscono stime di integrabilità di ordine elevato ($L^{2k}$) per la velocità $u$ e la velocità effettiva $w$, sfruttando la densità dell'insieme $A_{set}$ (insieme di parametri admissibili) per scegliere esponenti critici che chiudono le stime.

3. Contributi Chiave

1. Risposta Affermativa al Problema Aperto: Forniscono la prima prova dell'esistenza globale di soluzioni classiche per dati iniziali grandi (non piccoli) in un dominio multidimensionale reale (non solo 1D) con caratteristiche di alta dimensionalità.
2. Gestione della Singolarità Centrale: Risolvono il problema tecnico della singolarità all'origine, che è più forte in 3D ($x/r^2$) rispetto al 2D ($x/r$). Dimostrano che, grazie alla struttura BD e alla simmetria sferica, è possibile controllare questi termini senza richiedere restrizioni eccessive sui dati iniziali.
3. Distinzione tra Casi 2D e 3D:
   • 2D: Le stime sono ottenute per $\alpha \in (1/2, 1]$ e $\gamma > 2-\alpha$.
   • 3D: Richiedono condizioni più stringenti a causa della singolarità più forte: $\alpha \in (5/6, 1]$ e $\gamma > 4/3$ (per stelle gassose, $\kappa=1$). Per i plasmi ($\kappa=-1$), la condizione su $\gamma$ è meno restrittiva ($\gamma > 1$).
4. Caso $\alpha = 1$: Estendono i risultati al caso di viscosità lineare ($\mu=\rho$), che è fisicamente rilevante, dimostrando l'esistenza globale per $\gamma > 1$ in 2D e $4/3 < \gamma < 3$ in 3D.

4. Risultati Principali

Il teorema principale stabilisce l'esistenza e l'unicità di una soluzione classica sfericamente simmetrica $(\rho, u, \Phi)$ definita per ogni tempo $t \in [0, \infty)$, con le seguenti proprietà:

• Limiti Superiori e Inferiori della Densità: Esistono costanti $C(T)$ tali che $(C(T))^{-1} \leq \rho(x,t) \leq C(T)$. Questo garantisce che la soluzione non sviluppi vuoti (vacuum) né singolarità di densità (blow-up) in tempo finito.
• Regolarità: La soluzione appartiene agli spazi di Sobolev appropriati ($H^3$ per la densità, $H^4$ per la velocità, ecc.), confermando che si tratta di soluzioni classiche.
• Condizioni sui Parametri:
  • Per $N=2$: $\alpha \in (1/2, 1]$, $\gamma > 2-\alpha$ (se $\alpha < 1$) o $\gamma > 1$ (se $\alpha=1$).
  • Per $N=3$: $\alpha \in (5/6, 1]$, con vincoli specifici su $\gamma$ che dipendono da $\alpha$ e dal segno di $\kappa$ (gravità vs plasma).

5. Significatività

Questo lavoro rappresenta un progresso fondamentale nella teoria delle equazioni di fluidi comprimibili:

• Superamento dei limiti precedenti: Molti risultati precedenti richiedevano dati piccoli o assumevano la presenza di un "nucleo duro" (hard core) per evitare la singolarità centrale. Questo lavoro rimuove tali ipotesi, trattando il caso di una sfera solida che include l'origine.
• Applicabilità Astrofisica e Fisica dei Plasmi: I risultati sono direttamente applicabili alla modellazione di stelle gassose (dove la gravità è dominante) e plasmi, fornendo una base matematica rigorosa per l'evoluzione a lungo termine di questi sistemi con condizioni iniziali realistiche e grandi.
• Metodologia Innovativa: La tecnica di combinare le stime di entropia BD con stime pesate specifiche per la singolarità centrale e l'uso della velocità effettiva offre un nuovo paradigma per affrontare problemi con viscosità degenerata e termini di campo a lungo raggio in domini multidimensionali.

In sintesi, il paper risolve un problema aperto di decenni, dimostrando che le equazioni di Navier-Stokes-Poisson ammettono soluzioni globali classiche per dati grandi in domini sferici, superando le barriere poste dalle singolarità centrali attraverso un'analisi tecnica estremamente raffinata.