Existence and regularity of solutions to parabolic-elliptic nonlinear systems

1. Il Problema Studiato

Il paper si concentra sullo studio dell'esistenza e della regolarità delle soluzioni per un sistema di equazioni alle derivate parziali (PDE) di tipo parabolico-ellittico con coefficienti discontinui, definito in un dominio spaziotemporale $\Omega_T = \Omega \times (0, T)$, dove $\Omega \subset \mathbb{R}^N$ ($N > 2$) è limitato e aperto.

Il sistema è dato da:
$$\begin{cases} u_t - \text{div}(A(x, t)\nabla u) = -\text{div}(u M(x) \nabla \psi) + f(x, t) & \text{in } \Omega_T, \\ -\text{div}(M(x) \nabla \psi) = |u|^\theta & \text{in } \Omega_T, \\ u = \psi = 0 & \text{su } \partial\Omega \times (0, T), \\ u(x, 0) = 0 & \text{in } \Omega. \end{cases}$$

Caratteristiche principali del modello:

• Coefficienti: Le matrici $A(x,t)$ e $M(x)$ sono misurabili, limitate ($L^\infty$) e uniformemente ellittiche, ma non necessariamente continue.
• Dati: Il termine sorgente $f$ appartiene a $L^1(\Omega_T)$ (dati poco regolari).
• Non linearità: Il termine di accoppiamento $-\text{div}(u M(x) \nabla \psi)$ rappresenta un flusso convettivo. La seconda equazione è ellittica e dipende da $u$ in modo sublineare ($\theta \in (0, 2/N)$), a differenza del classico sistema di Keller-Segel che ha una dipendenza lineare.
• Contesto: Il sistema è una versione con coefficienti discontinui e sorgente esterna di una versione semplificata del modello di chemiotassi di Keller-Segel.

2. Metodologia

L'autore adotta un approccio basato su approssimazione e stime a priori, seguendo una strategia classica nell'analisi delle PDE non lineari con dati singolari:

1. Sistema Approssimante: Viene introdotto un sistema regolarizzato (3.1) dove i termini non lineari $u$ e $f$ sono troncati tramite l'operatore di troncamento $T_n(s) = \max(-n, \min(n, s))$. Questo garantisce l'esistenza di soluzioni regolari $(u_n, \psi_n)$ tramite il teorema del punto fisso di Schauder.
2. Stime A Priori: Vengono derivate stime energetiche e di sommabilità indipendenti da $n$ (o con dipendenza controllata) per le soluzioni approssimanti.
   • Si utilizzano funzioni di test specifiche (potenze di $u_n$, troncamenti $T_k(u_n)$, funzioni di Green $G_k$).
   • Si applicano disuguaglianze fondamentali come quella di Gagliardo-Nirenberg parabolica, Hölder, Young e il Lemma di Stampacchia per ottenere stime $L^\infty$ o $L^p$.
3. Passaggio al Limite:
   • Per dati più regolari ($f \in L^m$ con $m$ sufficientemente grande), si dimostra la convergenza forte in spazi di Lebesgue e la validità della soluzione nel senso debole classico.
   • Per dati singolari ($f \in L^1$), dove le stime energetiche standard falliscono, si ricorre alla teoria delle soluzioni di entropia.
   • Si utilizza un teorema di compattezza (Corollario 4 di [14], basato su Aubin-Lions) per estrarre sottosuccessioni convergenti.
   • Per trattare il termine non lineare $u_n \nabla \psi_n$, si sfrutta la convergenza quasi ovunque di $u_n$ e la regolarità di $\psi_n$ (ottenuta tramite teoremi di regolarità ellittica di Stampacchia).

3. Risultati Principali e Contributi Chiave

Il lavoro classifica l'esistenza e la regolarità della soluzione $u$ in base alla sommabilità del dato $f \in L^m(\Omega_T)$. Vengono definiti gli esponenti critici parabolici $m^*$ e $m^{**}$ (analoghi a quelli ellittici ma con $N$ sostituito da $N+2$).

A. Dati Regolari ($m > \frac{N+2}{2}$)

• Teorema 2.1: Se $f \in L^m$ con $m > \frac{N+2}{2}$, esiste una soluzione limitata ($u \in L^\infty(\Omega_T)$) con regolarità massima: $u \in C^0([0,T]; L^2) \cap L^2(0,T; H^1_0)$.
• Significato: L'interazione tra le due equazioni migliora la regolarità rispetto a equazioni paraboliche lineari con termini di trasporto singolari.

B. Dati Intermedi ($\frac{2N+4}{N+4} \le m < \frac{N+2}{2}$)

• Teorema 2.2: Se $f$ ha sommabilità inferiore, la soluzione non è più limitata, ma appartiene a spazi di Lebesgue e Sobolev specifici:
  $$u \in L^\infty(0, T; L^{\frac{Nm}{N+2-2m}}(\Omega)) \cap L^{m^{**}}(\Omega_T) \cap L^2(0, T; H^1_0(\Omega)).$$
• La condizione iniziale è soddisfatta in senso debole (tramite approssimazione).

C. Dati Singolari ($1 \le m < \frac{2N+4}{N+4}$)

In questo regime, il termine di trasporto $u M(x) \nabla \psi$ non appartiene necessariamente a $L^1$, rendendo impossibile la definizione di soluzione debole classica.

• Teorema 2.4 (Regime $m \in [\frac{2N+4}{N+6}, \frac{2N+4}{N+4})$): Esiste una soluzione debole con regolarità in $L^{m^*}(0,T; W^{1,m^*}_0)$.
• Teorema 2.6 (Regime $m \in [1, \frac{2N+4}{N+6})$): Per dati estremamente singolari ($f \in L^1$), si dimostra l'esistenza di una soluzione di entropia.
  • $u$ è una soluzione di entropia per la prima equazione.
  • $\psi$ è una soluzione debole per la seconda equazione.
  • Si ottengono stime di sommabilità: $u \in L^s(\Omega_T)$ e $\nabla u \in L^q(\Omega_T)$ per $s < 1^{**}$ e $q < 1^*$.

4. Significato e Innovazione

1. Generalizzazione del Modello Keller-Segel: Il lavoro estende i risultati noti per il sistema di Keller-Segel (spesso studiato con coefficienti lisci e dipendenza lineare) a un contesto con coefficienti discontinui e non linearità sublineare nel termine di sorgente ellittico.
2. Gestione dei Coefficienti Discontinui: Dimostra che la struttura ellittica della seconda equazione e la sublinearità ($\theta < 1$) permettono di ottenere stime di regolarità ($L^\infty$ o $L^p$) anche in presenza di coefficienti $A$ e $M$ solo misurabili e limitati, un risultato non banale per sistemi accoppiati.
3. Transizione tra Tipi di Soluzione: Il paper fornisce una mappa completa che collega la regolarità del dato $f$ al tipo di soluzione ottenibile (debole classica vs. entropia), chiarendo esattamente dove e perché le tecniche standard falliscono e dove è necessario ricorrere alla teoria dell'entropia.
4. Miglioramento della Regolarità: Viene evidenziato come l'accoppiamento parabolico-ellittico porti a un miglioramento della regolarità rispetto a equazioni paraboliche isolate con termini di trasporto dello stesso tipo (dove tipicamente si avrebbe solo $u \in L^1$ e $\log(1+|u|) \in H^1$).

In conclusione, il paper fornisce un quadro teorico rigoroso per l'analisi di sistemi di reazione-diffusione con chemiotassi in contesti matematici "sporchi" (coefficienti discontinui, dati $L^1$), stabilendo condizioni precise per l'esistenza e la regolarità delle soluzioni.