Global solutions of the 3D inhomogeneous incompressible viscoelastic system without structure assumptions

Titolo: Soluzioni globali del sistema viscoelastico incomprimibile non omogeneo 3D senza ipotesi strutturali

Autori: Chengfei Ai, Mengxing Bei e Yong Wang.

---

1. Il Problema Studiato

Il lavoro si concentra sul problema di Cauchy per il sistema viscoelastico incomprimibile non omogeneo in tre dimensioni spaziali ($\mathbb{R}^3$). Il sistema è descritto dalle seguenti equazioni (1.1):

• Equazione di continuità: $\tilde{\rho}_t + u \cdot \nabla \tilde{\rho} = 0$ (conservazione della massa con densità variabile $\tilde{\rho}$).
• Equazione del momento: $\tilde{\rho}u_t + \tilde{\rho}u \cdot \nabla u + \nabla p = \mu\Delta u + c^2 \text{div}(\tilde{\rho}FF^\top)$ (conservazione della quantità di moto, dove $p$ è la pressione, $\mu$ la viscosità e $c$ la velocità dell'onda elastica).
• Equazione di evoluzione del gradiente di deformazione: $F_t + u \cdot \nabla F = \nabla u F$ (dove $F$ è il gradiente di deformazione).
• Vincolo di incomprimibilità: $\text{div } u = 0$.

Le condizioni al contorno all'infinito richiedono che la densità tenda a 1, la velocità a 0 e il tensore di deformazione $F$ alla matrice identità $I$ quando $|x| \to \infty$.

La sfida principale:
La maggior parte dei risultati precedenti sulla globalità delle soluzioni per sistemi viscoelastici richiede l'assunzione di specifiche identità strutturali (spesso chiamate strutture "div-curl"), come $\text{div}(\tilde{\rho}F^\top) = 0$ o condizioni di "curl-free" su $F$. Questo paper mira a dimostrare l'esistenza globale di soluzioni forti senza alcuna di queste ipotesi strutturali aggiuntive, un problema notevolmente più difficile, specialmente nel caso non omogeneo dove la densità $\tilde{\rho}$ non è costante.

2. Metodologia e Strategie Chiave

Gli autori sviluppano una nuova strategia che combina tecniche di trasformazione del sistema e analisi spettrale avanzata.

A. Trasformazione del Sistema e Tensore Effettivo

Poiché la densità non è costante, l'operatore di proiezione di Helmholtz non può eliminare direttamente il termine di pressione $\nabla p$. Per superare questo ostacolo:

1. Viene introdotto un tensore effettivo $G$ definito come $G := \tilde{\rho}FF^\top - I$.
2. Viene definita una nuova variabile di flusso effettivo $\mathcal{G} := \Lambda^{-1}\text{div } G$ (dove $\Lambda = \sqrt{-\Delta}$ è l'operatore frazionario).
3. Viene introdotta una trasformazione per la pressione: $\nabla \tilde{p} := \frac{1}{\rho+1}\nabla p$ (dove $\rho = \tilde{\rho}-1$).
   Queste trasformazioni permettono di riscrivere il sistema originale (1.1) in un sistema dissipativo (2.4)-(2.5) per le variabili $(u, G)$, isolando meglio le componenti dissipative e rendendo il sistema più gestibile per l'analisi energetica.

B. Analisi Spettrale e Decadimento Temporale

Il cuore della dimostrazione risiede nell'uso dell'analisi spettrale per ottenere tassi di decadimento temporale migliorati per il campo di velocità $u$.

• Il sistema linearizzato viene analizzato nel dominio di Fourier. La matrice del sistema lineare ha autovalori che determinano il comportamento asintotico.
• Gli autori dividono lo spazio delle frequenze in due sottodomini ($S_1$ per basse frequenze e $S_2$ per alte frequenze) per stimare i nuclei del semigruppo associato.
• A differenza di lavori precedenti che richiedevano condizioni iniziali nello spazio $H^{-1}$, questo lavoro sostituisce tale condizione con condizioni $L^1$ sui dati iniziali ($u_0$ e $G_0$).
• L'analisi spettrale permette di dimostrare che la velocità $u$ decade nel tempo con un tasso accelerato, essenziale per controllare i termini non lineari.

C. Stime Energetiche Pesate nel Tempo

Viene utilizzata una metodologia energetica con pesi temporali.

• Si definisce un'energia ausiliaria $E_a(t)$ che include un peso $(1+t)^\sigma$.
• Si combinano le stime energetiche standard (metodo energetico) con le stime di decadimento ottenute dall'analisi spettrale.
• L'obiettivo è chiudere le stime a priori assumendo che le norme $H^2$ delle soluzioni rimangano piccole, dimostrando poi che questa assunzione è valida globalmente nel tempo.

3. Risultati Principali

Il teorema principale (Teorema 1.1) stabilisce quanto segue:

• Ipotesi sui Dati Iniziali: I dati iniziali $(\tilde{\rho}_0, u_0, F_0)$ devono soddisfare:

  1. Una condizione di picco in norma $H^2$: $\|u_0\|_{H^2} + \|\tilde{\rho}_0 - 1\|_{H^2} + \|G_0\|_{H^2} \le \delta$ (con $\delta$ sufficientemente piccolo).
  2. Una condizione di integrabilità $L^1$: $\|u_0\|_{L^1} + \|G_0\|_{L^1} = M_0 < +\infty$.
     Nota: La costante $M_0$ non deve essere piccola, a differenza di quanto richiesto in lavori precedenti basati su condizioni $H^{-1}$.

• Esistenza Globale: Esiste una soluzione unica globale $(\tilde{\rho}, u, F)$ definita per tutto $t \ge 0$.

• Stime Uniformi e Decadimento: La soluzione soddisfa stime uniformi che garantiscono la limitatezza delle norme $H^2$ per la densità e il tensore di deformazione, e un decadimento temporale per la velocità:
  $$\sup_{0 \le t \le \infty} \left[ \|( \tilde{\rho}-1, F-I)(t) \|_{H^2}^2 + (1+t)^\sigma \|u(t)\|_{H^2}^2 \right] + \int_0^\infty (1+t)^\sigma \|\nabla u(t)\|_{H^2}^2 dt \lesssim \text{dati iniziali}$$
  dove $\sigma \in (1, \sigma_0)$ con $1 < \sigma_0 < 3/2$. Questo intervallo per $\sigma$ è ottimale.

4. Contributi Chiave e Innovazioni

1. Rimozione delle Ipotesi Strutturali: Il risultato più significativo è la dimostrazione della ben-postezza globale senza assumere le identità "div-curl" ($\text{div}(\tilde{\rho}F^\top)=0$ o simili), che erano considerate necessarie per gestire la complessità del tensore di deformazione $F$.
2. Sostituzione delle Condizioni $H^{-1}$ con $L^1$: Gli autori migliorano le condizioni sui dati iniziali sostituendo la richiesta di appartenenza a spazi di Sobolev negativi ($H^{-1}$) con condizioni di integrabilità $L^1$. Questo è un risultato tecnico importante perché le condizioni $L^1$ sono spesso più fisiche e meno restrittive per la regolarità iniziale.
3. Nuove Tecniche di Trasformazione: L'introduzione del tensore $G$ e la trasformazione della pressione permettono di gestire il termine di pressione in un sistema a coefficienti variabili (densità non costante), un problema che non si presenta nel caso omogeneo.
4. Ottimizzazione del Tasso di Decadimento: L'uso combinato di analisi spettrale e pesi temporali permette di ottenere un tasso di decadimento per la velocità che è sufficiente a controllare le non linearità del sistema non omogeneo.

5. Significato e Impatto

Questo lavoro rappresenta un passo avanti significativo nella teoria matematica dei fluidi viscoelastici.

• Generalità: Rimuove restrizioni artificiali (ipotesi strutturali) che limitavano la validità fisica dei modelli precedenti, rendendo il risultato applicabile a una classe più ampia di fluidi reali.
• Robustezza del Metodo: Dimostra che l'analisi spettrale, combinata con tecniche energetiche pesate, è uno strumento potente per trattare sistemi accoppiati non lineari complessi senza dipendere da strutture speciali.
• Estensione al Caso Non Omogeneo: Estende i risultati noti per i fluidi omogenei e comprimibili al caso più difficile e fisicamente rilevante dei fluidi non omogenei incomprimibili, affrontando la difficoltà aggiuntiva rappresentata dalla densità variabile e dal termine di pressione.

In sintesi, il paper risolve un problema aperto nella teoria delle equazioni alle derivate parziali per i fluidi complessi, fornendo un quadro teorico solido per l'esistenza globale di soluzioni forti in assenza di simmetrie o vincoli strutturali specifici.