Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation

Sintesi Tecnica: Formazione di Singolarità di Implosione nelle Equazioni di Navier–Stokes–Korteweg Compressibili in 3D

Enunciato del Problema
Questo articolo investiga la formazione di singolarità in tempo finito per il sistema Navier–Stokes–Korteweg (NSK) tridimensionale compressibile con coefficienti di viscosità e capillarietà dipendenti dalla densità. Il sistema è definito da:
$$\begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \nabla P = \text{div}(2\mu(\rho)D(u)) + \nabla(\lambda(\rho)\text{div} u) + \text{div} K, \end{cases}$$
dove i coefficienti soddisfano le specifiche relazioni di legge di potenza:
$$\mu(\rho) = \nu\rho^\alpha, \quad \lambda(\rho) = 2\nu(\alpha - 1)\rho^\alpha, \quad \kappa(\rho) = \varepsilon^2\alpha^2\rho^{2\alpha-3}.$$
Precedenti lavori di Gu–Huang–Meng–Zhou [19] e Huang–Lei–Zhou [21] hanno stabilito l'esistenza di soluzioni forti globali per dati iniziali arbitrariamente grandi (lontani dal vuoto) quando l'esponente $\alpha$ si trova in un intervallo specifico idoneo (ad esempio, $\alpha \in (\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}}, 1)$ per il 3D). La questione centrale affrontata qui è se le soluzioni esistano globalmente o esplodano in tempo finito quando $\alpha$ cade al di fuori di questo intervallo, specificamente per piccoli esponenti positivi ($\alpha < 1/2$). In questo regime, la viscosità volumetrica effettiva diventa negativa ($\lambda(\rho) + \frac{2}{3}\mu(\rho) < 0$), una caratteristica associata a mezzi acustici fuori equilibrio e materia attiva, che sfida le standard teorie di esistenza globale.

Metodologia
Gli autori utilizzano un approccio perturbativo basato sulla costruzione di profili di implosione auto-similari regolari derivati dalle equazioni di Euler compressibili. La metodologia procede attraverso diverse fasi chiave:

1. Riformulazione della Velocità Effettiva: Per gestire il tensore di Korteweg del terzo ordine, il sistema viene riformulato utilizzando una velocità effettiva $v = u + d\rho^{\alpha-2}\nabla\rho$ (dove $d = \nu - \sqrt{\nu^2-\varepsilon^2}$). Questa trasformazione riclassifica il sistema in una struttura parabolo-parabolica per $(\rho, v)$, rivelando un termine di diffusione degenere nascosto $d\Delta(\rho^\alpha)$ nell'equazione della massa.
2. Scaling Auto-simile: Il problema viene trasformato in coordinate auto-similari $(\tau, y)$ dove $\tau = -\log(T-t)/\Lambda$ e $y = x/(T-t)^{1/\Lambda}$. In queste coordinate, gli effetti viscosi e capillari appaiono come perturbazioni a decadimento esponenziale ($e^{-f_1\tau}$) a condizione che il parametro di scaling $\Lambda$ soddisfi una specifica condizione ($f_1 > 0$).
3. Analisi di Stabilità delle Perturbazioni: Il cuore della dimostrazione consiste nell'analizzare la stabilità non lineare del sistema attorno a un profilo di Euler auto-simile sfericamente simmetrico e regolare $(\bar{S}, \bar{V})$. Le equazioni di perturbazione vengono derivate e viene costruito un argomento di bootstrap per controllare la soluzione.
4. Stime di Energia e Decomposizione Spettrale:
   • Stime di Energia Pesata: Vengono stabilite stime di energia pesata di ordine superiore per controllare il decadimento delle perturbazioni. Gli autori utilizzano una funzione di peso specifica $\phi$ e una gerarchia di parametri per gestire l'interazione tra trasporto e dissipazione.
   • Argomento di Bootstrap: Viene progettato un rigoroso schema di bootstrap per chiudere le stime per le derivate di ordine inferiore (specificamente il limite inferiore della densità) e per le derivate di ordine superiore. Un'innovazione critica è l'uso sia della struttura di trasporto che di quella di dissipazione per stabilire un limite inferiore per la densità, superando la difficoltà per cui l'equazione della densità possiede una struttura di dissipazione oltre a quella di trasporto.
   • Controllo dei Modi Instabili: L'operatore linearizzato viene decomposto in sottospazi stabili e instabili. I componenti instabili, che corrispondono a uno spazio a dimensione finita, sono controllati selezionando coefficienti specifici dei dati iniziali. Questa selezione assicura che la soluzione rimanga su una varietà a codimensione finita che conduce all'esplosione.

Contributi Principali e Risultati
Il lavoro dimostra l'esistenza di dati iniziali regolari, con densità uniformemente separate dal vuoto, che evolvono in singolarità di implosione in tempo finito per il sistema NSK 3D quando $\alpha$ è sufficientemente piccolo ($\alpha < 1/2$).

• Teorema Principale (Teorema 1.1): Per parametri che soddisfano o:
  • (P1) $1 < \gamma < 1 + 2/\sqrt{3}$, $0 < \alpha < 1/2$, e $f_1 > 0$; oppure
  • (P2) $\gamma \in (1, 1 + 2/\sqrt{3}) \setminus \mathcal{J}$ (dove $\mathcal{J}$ è un insieme eccettuale numerabile) e $0 < \alpha < \alpha^*(\gamma)$,
    esistono dati iniziali $C^\infty$ $(\rho_0, u_0)$ tali che la corrispondente soluzione $(\rho, u)$ esiste su $[0, T)$ e sviluppa una singolarità di implosione al tempo $T$.
• Caratteristiche della Singolarità: Al tempo di esplosione $T$:
  • La densità diventa infinita all'origine: $\lim_{t \to T^-} \rho(t, 0) = +\infty$.
  • La velocità effettiva $v$ è illimitata in ogni vicinato dell'origine.
  • La soluzione converge, dopo riscalamento, al profilo di implosione auto-simile prescritto.
• Corollario 1.1: Il sistema originale sviluppa una singolarità in tempo finito, specificamente caratterizzata dall'esplosione della densità ($\lim_{t \to T^-} \rho(t, 0) = +\infty$). Tuttavia, l'esplosione della velocità fisica $u$ non può essere direttamente conclusa dall'esplosione della velocità effettiva, poiché il termine di gradiente $\nabla\rho$ potrebbe anch'esso diventare singolare vicino all'origine, e una possibile cancellazione nella relazione $v = u + d\rho^{\alpha-2}\nabla\rho$ non può essere esclusa.

Significatività
Il saggio sostiene di completare la teoria esistente di esistenza globale per il sistema NSK. Mentre i risultati precedenti hanno dimostrato la ben-posedness globale per un intervallo specifico di $\alpha$, questo lavoro identifica un meccanismo distinto di esplosione in tempo finito per il regime "small-$\alpha$". Questo regime è caratterizzato da una struttura di viscosità volumetrica effettiva negativa, un fenomeno rilevante per la fisica fuori equilibrio. I risultati dimostrano che anche con dati iniziali regolari lontani dal vuoto, l'interazione tra viscosità dipendente dalla densità, capillarietà e viscosità volumetrica negativa può portare a singolarità di implosione dove i campi di densità e di velocità effettiva esplodono essi stessi, piuttosto che solo le loro derivate (formazione di shock). La costruzione si basa sulla stabilità dei profili di Euler auto-similari, estendendo il meccanismo di implosione noto per le equazioni di Euler al contesto viscoso e capillare di NSK.