Optimal Bilinear control restricted to the three-dimensional chemo-repulsion model with potential production

Titel: Optimale bilineare Steuerung eingeschränkt auf das dreidimensionale Chemo-Repulsions-Modell mit potenzieller Produktion

Autoren: Francisco Guillén-González, Exequiel Mallea-Zepeda, María A. Rodríguez-Bellido, und Élder J. Villamizar-Roa

1. Problemstellung

Das Papier untersucht ein gekoppeltes System parabolischer partieller Differentialgleichungen (PDEs), das ein Chemo-Repulsions-Modell in einem dreidimensionalen beschränkten Gebiet $\Omega \subset \mathbb{R}^3$ beschreibt. Das System modelliert die Bewegung von Zellen (Dichte $u$), die sich von einem chemischen Signal (Konzentration $v$) abstößt, wobei die Zellen das Signal selbst produzieren.

Das System wird durch folgende Gleichungen in $Q = (0, T) \times \Omega$ definiert:

1. Zelldichte ($u$):
   $$\partial_t u - \Delta u = \nabla \cdot (u \nabla v) + r u - \mu u^p$$
   Der Term $\nabla \cdot (u \nabla v)$ beschreibt die Repulsion (im Gegensatz zur Anziehung bei Chemotaxis). Der Term $r u - \mu u^p$ ist ein logistischer Wachstumsterm (für $r, \mu \ge 0$), der Konkurrenz und Sättigung modelliert. Für den Fall ohne logistisches Wachstum ($r=\mu=0$) wird nur der Diffusionsterm betrachtet.

2. Chemische Konzentration ($v$):
   $$\partial_t v - \Delta v + v = u^p + f v \mathbf{1}_{\Omega_c}$$
   Hier beschreibt $u^p$ die nichtlineare Produktion des Signals durch die Zellen. Der Term $f v \mathbf{1}_{\Omega_c}$ stellt eine bilineare Steuerung dar, wobei $f(t,x)$ als Kontrollfunktion auf einem Teilgebiet $\Omega_c \subseteq \Omega$ wirkt und als Proliferations- oder Degradationskoeffizient für das chemische Signal fungiert.

Das System ist mit Anfangsbedingungen und homogenen Neumann-Randbedingungen (kein Fluss) versehen.

Ziel: Das Hauptziel ist die Untersuchung der Existenz und Regularität von Lösungen für dieses gesteuerte System in 3D und die Formulierung sowie Analyse eines zugehörigen optimalen Steuerungsproblems.

2. Methodik und mathematischer Rahmen

Die Autoren verwenden eine Kombination aus analytischen Techniken der nichtlinearen PDE-Theorie und der Optimalsteuerungstheorie in Banachräumen:

• Existenz schwacher Lösungen:

  • Es werden formale Energieabschätzungen durchgeführt, um die Rolle der Diffusion im Vergleich zur nichtlinearen Produktion $u^p$ zu analysieren.
  • Für den Fall $1 < p \le 5/3$ wird gezeigt, dass die Diffusion ausreicht, um die Energie zu kontrollieren (ohne logistischen Term).
  • Für $p > 5/3$ ist der logistische Term ($-\mu u^p$) notwendig, um die Energieabschätzungen zu schließen und das "Blow-up" zu verhindern.
  • Die Existenz wird durch eine Regularisierungsmethode (Einführung eines Parameters $\varepsilon$) und den Leray-Schauder-Fixpunktsatz bewiesen, gefolgt von Kompaktheitsargumenten (Aubin-Lions-Lemma).

• Regularitätskriterium (Starke Lösungen):

  • Um ein Optimalsteuerungsproblem zu lösen, sind stärkere Lösungen erforderlich. Die Autoren leiten ein Serrin-artiges Regularitätskriterium her.
  • Unter der Annahme, dass die Steuerung $f \in L^{5/2}(0, T; L^{5/2+}(\Omega_c))$ und die Lösung $u$ bestimmte Integrabilitätsbedingungen erfüllt ($u \in L^{5p/2} \cap L^{10/3}$), wird gezeigt, dass schwache Lösungen zu starken Lösungen werden ($u, \nabla v \in L^5(Q)$).
  • Ein nicht-standardisiertes "Bootstrapping"-Argument wird verwendet, um die Regularität der Konvektions- und Reaktionsterme zu verbessern.

• Optimale Steuerung:

  • Das Steuerungsproblem minimiert ein Kostenfunktional $J$, das die Abweichung von gewünschten Zuständen ($u_d, v_d$) und die Kosten der Steuerung $f$ misst.
  • Die Differentiierbarkeit der Abbildung "Steuerung-zu-Zustand" wird mittels des Satzes über die implizite Funktion in Banachräumen bewiesen. Dies erfordert die Wohlgestelltheit des linearisierten Systems.
  • Die Adjungierten (Dual-)Probleme werden mittels eines Dualitätsarguments gelöst, um den Gradienten des reduzierten Kostenfunktionals explizit zu identifizieren.

3. Hauptbeiträge und Ergebnisse

Die Arbeit liefert mehrere signifikante theoretische Durchbrüche:

1. Existenz globaler schwacher Lösungen in 3D:

   • Es wird bewiesen, dass für den nicht-logistischen Fall ($r=\mu=0$) globale schwache Lösungen existieren, wenn $1 < p \le 5/3$.
   • Für den Fall $p > 5/3$ wird gezeigt, dass der Zusatz des logistischen Terms ($ru - \mu u^p$) notwendig und hinreichend ist, um die Existenz globaler schwacher Lösungen zu garantieren. Dies füllt eine Lücke in der Literatur, da 3D-Chemotaxis-Modelle oft nur für lineare Produktion ($p=1$) oder in 2D gelöst wurden.

2. Reguläritätskriterium und Eindeutigkeit:

   • Es wird ein Kriterium etabliert, das schwache Lösungen in starke Lösungen überführt. Dies ist entscheidend, da die Optimalsteuerungstheorie auf differenzierbaren Abbildungen basiert, die nur für starke Lösungen definiert sind.
   • Die Eindeutigkeit der starken Lösungen wird unter diesen Regularitätsannahmen bewiesen.

3. Existenz globaler optimaler Lösungen:

   • Unter Verwendung eines Minimierungsfolgen-Arguments wird die Existenz mindestens einer globalen optimalen Lösung für das gesteuerte Problem nachgewiesen.

4. Notwendige Optimalitätsbedingungen erster Ordnung:

   • Im Gegensatz zu früheren Arbeiten, die oft auf generischen Lagrange-Multiplikatoren basierten, berechnen die Autoren die Ableitung der Zustandsvariablen nach der Steuerung explizit über den Satz über die implizite Funktion.
   • Durch die Lösung des adjungierten Systems wird der Gradient des reduzierten Kostenfunktionals identifiziert.
   • Daraus werden notwendige Optimalitätsbedingungen erster Ordnung für lokale optimale Lösungen hergeleitet (Variationsungleichung).

4. Signifikanz und Bedeutung

• Mathematische Tiefe: Die Arbeit adressiert die komplexe Wechselwirkung zwischen nichtlinearer Produktion ($u^p$), Repulsion und bilinearer Steuerung in drei Dimensionen, was mathematisch anspruchsvoller ist als die üblichen 2D- oder linearen Fälle.
• Brücke zur Numerik: Die explizite Identifikation des Gradienten des Kostenfunktionals ist von großer praktischer Bedeutung. Sie ermöglicht die Entwicklung von numerischen Algorithmen auf Gradientenabstiegsbasis (z. B. für die Steuerung von Zellpopulationen oder chemischen Prozessen), was in der Chemotaxis-Forschung bisher wenig untersucht wurde.
• Erweiterung des Zustandsraums: Die Arbeit erweitert den bekannten Rahmen der optimalen Steuerung für Chemotaxis-Modelle von linearen Produktionstermen ($p=1$) auf den nichtlinearen Fall ($p>1$) in 3D, was biologisch realistischer ist (z. B. bei Sättigungseffekten).
• Methodischer Fortschritt: Die Kombination aus Energieabschätzungen, Regularitätskriterien und der Anwendung des Satzes über die implizite Funktion in spezifischen Banachräumen ($X_u, X_v$) bietet einen robusten Rahmen für die Analyse ähnlicher gekoppelter PDE-Systeme mit Steuerung.

Zusammenfassend liefert dieses Papier eine vollständige theoretische Grundlage für die Analyse und Steuerung von 3D-Chemo-Repulsionsmodellen mit nichtlinearer Produktion und bilinearer Steuerung, von der Existenz schwacher Lösungen bis hin zu notwendigen Bedingungen für optimale Steuerungen.