Global existence and optimal decay for a three-dimensional penalized Navier--Stokes system with biharmonic damping

Sintesi Tecnica: Esistenza Globale e Decadimento Ottimale per un Sistema di Navier–Stokes Penalizzato Tridimensionale con Smorzamento Biharmonico

Enunciato del Problema
Questo articolo investiga un sistema parabolico tridimensionale che funge da approssimazione iperviscosa e penalizzata delle equazioni di Navier–Stokes incomprimibili. Il modello affronta le sfide matematiche delle equazioni di Navier–Stokes classiche (1.1), in particolare il problema aperto della regolarità globale e dell'unicità per dati iniziali arbitrari in tre dimensioni. Il sistema proposto (1.3) combina tre meccanismi:

1. Diffusione viscosa classica ($-\nu\Delta u$).
2. Smorzamento biharmonico (iperviscoso) ($\beta\Delta^2 u$), che fornisce una dissipazione di ordine superiore.
3. Penalizzazione della divergenza ($-\frac{1}{\varepsilon}\nabla \text{div} u$), che rilassa il vincolo di incompressibilità stretta ($\text{div} u = 0$) imponendo una quasi-incompressibilità al limite $\varepsilon \to 0$.

Fondamentalmente, il termine di convezione non lineare è modificato utilizzando una correzione antisimmetrica di tipo Temam ($N(u) = (u \cdot \nabla)u + \frac{1}{2}u \text{div} u$). Questa correzione compensa la debole comprimibilità introdotta dalla penalizzazione, ripristinando la proprietà di cancellazione $L^2$ del termine di trasporto, che viene tipicamente persa nelle formulazioni comprimibili.

Metodologia
L'analisi si basa su una combinazione di metodi energetici, teoria dei semigruppi e tecniche di scissione di Fourier (Fourier splitting), con un focus centrale sulla derivazione di stime che siano uniformi rispetto al parametro di penalizzazione $\varepsilon > 0$.

• Soluzioni Deboli (Teorema 1.3): L'esistenza globale di soluzioni deboli per dati iniziali $L^2$ arbitrari è stabilita tramite un'approssimazione di tipo Galerkin di Friedrichs. Gli autori proiettano l'equazione su sfere di Fourier utilizzando operatori di cut-off $J_n$ anziché il proiettore di Leray. Questa distinzione è critica; l'uso del proiettore di Leray annullerebbe il termine di penalizzazione, alterando il problema limite. Stime a priori uniformi sono derivate dall'identità energetica, e la compattezza è ottenuta tramite il lemma di Aubin–Lions–Simon.
• Soluzioni Forti (Teorema 1.4): Per dati iniziali piccoli in $H^2(\mathbb{R}^3)$, l'esistenza globale di soluzioni forti è provata utilizzando un argomento di punto fisso sulla formulazione mite associata al semigruppo lineare $S_\varepsilon(t) = e^{tL_\varepsilon}$. L'operatore lineare $L_\varepsilon = \nu\Delta - \beta\Delta^2 + \frac{1}{\varepsilon}\nabla\text{div}$ genera un semigruppo con proprietà di regolarizzazione (smoothing). La prova si basa su una stima energetica $H^2$ uniforme, dove i termini non lineari sono assorbiti dalla dissipazione viscosa e biarmonica.
• Stime di Decadimento (Teorema 1.6): Ottimi tassi di decadimento nel tempo lungo sono derivati utilizzando il metodo di scissione di Fourier (Schonbek). L'analisi separa le basse e le alte frequenze. Il comportamento a bassa frequenza è controllato tramite la formula di Duhamel e la specifica decomposizione della non linearità, mentre le alte frequenze sono gestite attraverso disuguaglianze energetiche differenziate e interpolazione di Gagliardo–Nirenberg.

Contributi Chiave e Risultati

1. Esistenza Globale Debole: Il saggio prova l'esistenza globale di soluzioni deboli per qualsiasi dato iniziale $u_0 \in L^2(\mathbb{R}^3)$. Le soluzioni soddisfano specifiche classi di regolarità, tra cui $L^\infty(0, T; L^2)$, $L^2(0, T; H^2)$ e $L^4(0, T; L^4)$.
2. Soluzioni Forti Globali con Piccolezza Uniforme: Per dati iniziali sufficientemente piccoli in $H^2(\mathbb{R}^3)$, il sistema ammette una soluzione forte globale unica. Un contributo significativo è che la soglia di piccolezza $\delta_0$ e le associate stime di energia sono indipendenti da $\varepsilon$. Questa uniformità garantisce che i risultati si applichino all'intera famiglia di sistemi penalizzati, non solo a un singolo caso fissato.
3. Tassi di Decadimento Ottimali: Per dati iniziali in $L^1(\mathbb{R}^3) \cap H^2(\mathbb{R}^3)$, gli autori derivano stime di decadimento ottimali:
   $$\|u(t)\|_{L^2} \leq C(1+t)^{-3/4}$$
   $$\|\nabla^k u(t)\|_{L^2} \leq C(1+t)^{-3/4 - k/2}, \quad k=1, 2$$
   Questi tassi corrispondono al decadimento asintotico dell'equazione del calore classica in tre dimensioni. L'analisi conferma che, sebbene il termine biarmonico domini la dissipazione ad alta frequenza, il decadimento a lungo termine è governato dal comportamento del flusso di calore delle basse frequenze.
4. Uniformità in $\varepsilon$: Una caratteristica definente dell'analisi è che tutte le costanti e le stime rimangono uniformi al tendere di $\varepsilon \to 0$. Il termine di penalizzazione $\frac{1}{\varepsilon}\|\text{div} u\|^2$ è trattato come un termine di dissipazione non negativo che viene o mantenuto o scartato, mai limitato superiormente da una quantità che coinvolge $1/\varepsilon$.

Significatività e Rivendicazioni
Il saggio sostiene di fornire una "fondazione analitica stabile e rigorosa" per l'approssimazione penalizzata dei flussi incomprimibili. Incorporando lo smorzamento biarmonico e la correzione di Temam, il modello raggiunge risultati di ben l'esistenza (well-posedness) migliorati rispetto alle equazioni di Navier–Stokes classiche, specificamente riguardo alla disponibilità di stime energetiche $L^2$ incondizionate e di una maggiore regolarità (smoothing).

Gli autori osservano che le loro stime uniformi suggeriscono che, per $\varepsilon \to 0$, le soluzioni forti convergano a un sistema limite che risolve un sistema di Navier–Stokes iperviscoso incomprimibile. Tuttavia, essi dichiarano esplicitamente che l'identificazione del sistema limite e la prova dei tassi di convergenza in $\varepsilon$ sono direzioni lasciate a lavori futuri. Inoltre, sebbene la regolarità globale per dati arbitrari sia nota per il sistema iperviscoso incomprimibile (grazie a Lions), gli autori non rivendicano un tale risultato per il sistema penalizzato con dati grandi, notando che la struttura comprimibile antisimmetrica modifica le stime energetiche di ordine superiore in modo tale da impedire una diretta estensione del risultato per dati grandi.

Il lavoro si distingue dagli studi precedenti sulle equazioni di Navier–Stokes smorzate (ad esempio, con termini di assorbimento) concentrandosi su un meccanismo di smorzamento di ordine superiore e un metodo di penalizzazione che elimina il termine di pressione, semplificando così la formulazione debole rimuovendo il vincolo di divergenza sulle funzioni di test.