Existence and regularity of solutions to parabolic-elliptic nonlinear systems

Titel: Existenz und Regularität von Lösungen zu parabolisch-elliptischen nichtlinearen Systemen

1. Problemstellung

Der Artikel untersucht ein gekoppeltes System aus partiellen Differentialgleichungen (PDEs) auf einem beschränkten Gebiet $\Omega \subset \mathbb{R}^N$ ($N \ge 3$) über den Zeitintervall $(0, T)$. Das System besteht aus einer parabolischen Gleichung für eine Variable $u$ und einer elliptischen Gleichung für eine Variable $\psi$:

$$\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{auf } \partial\Omega \times (0, T), \\ u(x, 0) = 0 & \text{in } \Omega. \end{cases}$$

Wichtige Merkmale des Modells:

• Kopplung: Die Gleichungen sind durch den Term $u M(x) \nabla \psi$ (in der ersten Gleichung) und $|u|^\theta$ (in der zweiten Gleichung) gekoppelt.
• Koeffizienten: Die Matrizen $A(x,t)$ und $M(x)$ sind messbar, beschränkt ($L^\infty$) und gleichmäßig elliptisch. Sie können unstetig sein.
• Quellen: $f(x,t)$ ist eine nichtnegative Funktion in $L^1(\Omega_T)$.
• Exponent: Der Parameter $\theta$ liegt im Intervall $(0, 2/N)$. Da $N > 2$, gilt $\theta < 1$, was eine sublineare Abhängigkeit in der elliptischen Gleichung bedeutet.
• Kontext: Das System ist eine modifizierte Version des Keller-Segel-Modells (Chemotaxis), jedoch mit diskontinuierlichen Koeffizienten, einem Quellterm $f$ und einer sublinearen Kopplung. Im Gegensatz zum klassischen Keller-Segel-Modell ist die zweite Gleichung elliptisch statt parabolisch.

2. Methodik

Der Beweis der Existenz und Regularität erfolgt durch eine Approximationsmethode (Regularisierung) in Kombination mit a-priori-Abschätzungen und Kompaktheitsargumenten.

1. Approximationsproblem:
   Es wird eine Folge von Lösungen $(u_n, \psi_n)$ für ein regularisiertes System betrachtet, bei dem die nichtlinearen Terme $u$ und $f$ durch Trunkierungsfunktionen $T_n(\cdot)$ (z.B. $T_n(s) = \max(-n, \min(s, n))$) beschränkt werden.

   • Die Existenz für das approximierende System wird mittels eines Fixpunktsatzes von Schauder nachgewiesen. Dabei wird ein Operator definiert, der eine Funktion $w$ auf die Lösung $v$ der ersten Gleichung abbildet, wobei die zweite Gleichung als elliptisches Problem mit dem Quellterm $|T_n(w)|^\theta$ gelöst wird.

2. A-priori-Abschätzungen:
   Je nach der Summierbarkeit (Integrabilität) der Daten $f \in L^m(\Omega_T)$ werden unterschiedliche Abschätzungen für die Folge $(u_n, \psi_n)$ hergeleitet:

   • $L^1$-Abschätzung: Durch Testen mit $T_h(u_n)$ wird gezeigt, dass $u_n$ in $L^\infty(0, T; L^1(\Omega))$ beschränkt ist.
   • $L^\infty$-Abschätzung: Für hinreichend glatte Daten ($m$ groß) wird die Beschränktheit von $u_n$ in $L^\infty(\Omega_T)$ bewiesen (unter Verwendung von Stampacchias Lemma und Gagliardo-Nirenberg-Ungleichungen).
   • Energieabschätzungen: Für $f \in L^m$ mit $m \ge \frac{2N+4}{N+4}$ werden Abschätzungen in $L^2(0, T; H^1_0(\Omega))$ und $L^{m^{**}}(\Omega_T)$ hergeleitet.
   • Entropie-Abschätzungen: Für singuläre Daten ($1 \le m < \frac{2N+4}{N+4}$) werden Abschätzungen für Trunkierungen $T_k(u_n)$ in $L^2(0, T; H^1_0(\Omega))$ und für die Gradienten in $L^q$-Räumen mit $q < \frac{N+2}{N+1}$ durchgeführt.

3. Grenzübergang ($n \to \infty$):

   • Kompaktheit: Mittels des Aubin-Lions-Lemmas (oder Korollar 4 aus [14]) wird die starke Konvergenz der Folge $(u_n)$ in geeigneten $L^p$-Räumen gezeigt.
   • Konvergenz der elliptischen Komponente: Da die rechte Seite der elliptischen Gleichung $|u_n|^\theta$ stark konvergiert, folgt die starke Konvergenz von $\psi_n$ und $\nabla \psi_n$ in $L^2$.
   • Entropielösungen: Für den Fall sehr singulärer Daten ($m < \frac{2N+4}{N+6}$), wo der Term $u \nabla \psi$ nicht mehr in $L^1$ liegt, wird der Grenzübergang im Sinne von Entropielösungen (nach dem Konzept von Boccardo, Gallouët, Orsina) durchgeführt. Dies erfordert das Testen mit Funktionen der Form $T_k(u - \phi)$.

3. Wichtige Beiträge und Ergebnisse

Die Hauptergebnisse des Papers hängen von der Integrabilitätsklasse $m$ der Quelldaten $f \in L^m(\Omega_T)$ ab:

• Fall 1: Hohe Regularität ($m > \frac{N+2}{2}$)

  • Es existiert eine stark reguläre Lösung.
  • $u \in C^0([0, T]; L^2(\Omega)) \cap L^2(0, T; H^1_0(\Omega)) \cap L^\infty(\Omega_T)$.
  • Die Lösung erfüllt die Gleichungen im klassischen schwachen Sinne.

• Fall 2: Mittlere Regularität ($\frac{2N+4}{N+4} \le m < \frac{N+2}{2}$)

  • Die Lösung $u$ ist nicht notwendigerweise beschränkt, aber sie besitzt eine verbesserte Summierbarkeit.
  • $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))$.
  • Hier gilt $m^{**} = \frac{(N+2)m}{N+2-2m}$. Dies entspricht den Regularitätssätzen von Stampacchia für parabolische Gleichungen.

• Fall 3: Geringe Regularität / Singuläre Daten ($1 \le m < \frac{2N+4}{N+4}$)

  • In diesem Bereich kann keine Lösung mit endlicher Energie ($H^1$) garantiert werden.
  • Es wird die Existenz einer Entropielösung bewiesen.
  • Für $1 < m < \frac{2N+4}{N+4}$ gilt: $u \in L^\infty(0, T; L^1(\Omega)) \cap L^{m^{**}}(\Omega_T) \cap L^{m^*}(0, T; W^{1, m^*}_0(\Omega))$.
  • Für $m=1$: $u \in L^\infty(0, T; L^1(\Omega)) \cap L^s(\Omega_T) \cap L^q(0, T; W^{1, q}_0(\Omega))$ für alle $s < 1^{**}$ und $q < 1^*$.
  • Die Entropie-Bedingung kompensiert das Fehlen von $L^1$-Integrabilität des drift-Terms $u \nabla \psi$.

Spezifische Regularitätsverbesserung:
Ein zentrales Ergebnis ist, dass trotz des nicht-regulären Drift-Terms $u M(x) \nabla \psi$ (der im Allgemeinen nur in $L^1$ liegen würde) die Kopplung mit der elliptischen Gleichung zu einer Verbesserung der Regularität führt. Die Lösung $u$ gehört zu besseren Räumen als bei einer reinen parabolischen Gleichung mit gleichem Quellterm. Dies wird als "Interplay" (Wechselspiel) zwischen den beiden Gleichungen bezeichnet.

4. Bedeutung und Fazit

• Theoretischer Fortschritt: Das Paper erweitert die Theorie der parabolisch-elliptischen Systeme auf den Fall von diskontinuierlichen Koeffizienten und sublinearer Kopplung ($\theta < 1$).
• Behandlung singulärer Daten: Es liefert einen vollständigen Existenzbeweis für Daten in $L^1$, was für Anwendungen in der Physik und Biologie (z.B. Chemotaxis mit massiven Quellen) relevant ist, wo Daten oft nur als Maße oder $L^1$-Funktionen vorliegen.
• Methodische Klarheit: Die Arbeit demonstriert effektiv, wie man bei Systemen, bei denen der Drift-Term die Standard-Energieabschätzungen zerstört, auf Entropielösungen zurückgreift, um Existenz und Regularität zu sichern.
• Vergleich mit Literatur: Die Ergebnisse korrelieren mit den bekannten $L^p$-Regularitätssätzen für parabolische Gleichungen (wie in [1] und [5] zitiert), zeigen aber, dass die spezifische Struktur des Keller-Segel-Typs (elliptische Rückkopplung) die Regularität der Lösung $u$ über das hinaushebt, was für isolierte parabolische Gleichungen mit gleichen Daten zu erwarten wäre.

Zusammenfassend stellt das Paper einen robusten analytischen Rahmen für die Analyse eines wichtigen nichtlinearen PDE-Systems mit realistischen, aber mathematisch anspruchsvollen Bedingungen (diskontinuierliche Koeffizienten, $L^1$-Daten) bereit.