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

Titel: Optimaal bilineair beheer beperkt tot het driedimensionale chemo-repulsiemodel met potentiële productie

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

---

1. Probleemstelling

Het artikel onderzoekt een gekoppeld stelsel van parabolische partiële differentiaalvergelijkingen (PDE's) dat het gedrag van cellen beschrijft die zich verplaatsen als reactie op een chemisch signaal (chemotaxis), specifiek in een chemo-repulsie scenario (afstoting). Het model is gedefinieerd in een begrensde domein $\Omega \subset \mathbb{R}^3$ over een tijdsinterval $(0, T)$.

Het systeem omvat twee variabelen:

• $u(t, x)$: De dichtheid van de cellen.
• $v(t, x)$: De concentratie van het chemische signaal.

De vergelijkingen zijn:

1. Voor de cel-dichtheid ($u$):
   $$\partial_t u - \Delta u = \nabla \cdot (u \nabla v) + r u - \mu u^p$$
   Hierin vertegenwoordigt $\nabla \cdot (u \nabla v)$ het repulsie-mechanisme (cellen bewegen weg van hoge concentraties). De termen $ru - \mu u^p$ beschrijven een logistische groei (proliferatie en competitie), waarbij $p > 1$.
2. Voor de chemische concentratie ($v$):
   $$\partial_t v - \Delta v + v = u^p + f v \mathbb{1}_{\Omega_c}$$
   De term $u^p$ staat voor de niet-lineaire productie van het chemische signaal door de cellen. De term $f v \mathbb{1}_{\Omega_c}$ is een bilineaire controle, waarbij $f(t, x)$ een controlefunctie is die werkt op een subdomein $\Omega_c \subseteq \Omega$. Deze controle fungeert als een proliferatie- of degradatiecoëfficiënt voor het chemische signaal.

Doel: Het bewijzen van het bestaan van globale zwakke oplossingen, het vaststellen van regulariteitscriteria voor sterke oplossingen, en het oplossen van een optimaal beheerprobleem waarbij men een kostenfunctionaal minimaliseert door de controle $f$ te variëren.

---

2. Methodologie

De auteurs hanteren een geavanceerde analysebenadering die bestaat uit de volgende stappen:

• Energie-ongelijkheden: Ze leiden formele energie-ongelijkheden af om de groei van de oplossingen te beheersen. Een cruciaal punt is de analyse van de interactie tussen de diffusie-term en de niet-lineaire productieterm $u^p$.
• Regulering en Grenswaarden: Om het bestaan van zwakke oplossingen te bewijzen, introduceren ze een familie van geregulariseerde problemen (met een parameter $\epsilon > 0$). Ze gebruiken de Leray-Schauder-vastpuntstelling om het bestaan van sterke oplossingen voor deze geregulariseerde systemen aan te tonen. Vervolgens laten ze zien dat deze oplossingen convergeren naar een zwakke oplossing van het originele probleem als $\epsilon \to 0$.
• Regulariteitscriteria (Bootstrapping): Ze ontwikkelen een niet-standaard "bootstrapping"-argument gebaseerd op lineaire diffusie-convectie-reactievergelijkingen met niet-reguliere coëfficiënten. Dit wordt gebruikt om aan te tonen dat onder bepaalde voorwaarden zwakke oplossingen in feite sterke oplossingen zijn (met hogere regulariteit).
• Optimaal Beheer: Voor het beheerprobleem gebruiken ze de Implicit Function Theorem in Banachruimten om de differentieerbaarheid van de afbeelding "controle-naar-toestand" (control-to-state mapping) te bewijzen. In plaats van een generieke Lagrange-multiplicator-benadering, berekenen ze de afgeleide expliciet via het geadjungeerde systeem (adjoint problem), opgelost via een dualiteitsargument.

---

3. Belangrijkste Bijdragen en Resultaten

A. Bestaan van Zwakke Oplossingen

De auteurs bewijzen het bestaan van globale-in-tijd zwakke oplossingen voor het 3D-systeem, afhankelijk van de exponent $p$:

• Geen logistieke term ($r=\mu=0$): Voor $1 < p \leq 5/3$ bestaat er een zwakke oplossing. De diffusie is hier voldoende om de energie te beheersen.
• Met logistieke term ($r, \mu > 0$): Voor $p > 5/3$ is de logistische term $ru - \mu u^p$ noodzakelijk om de energie-ongelijkheid te sluiten en het bestaan van een zwakke oplossing te garanderen. De logistische competitie fungeert als een rem op de explosieve groei veroorzaakt door de sterke niet-lineariteit.

B. Regulariteitscriteria en Sterke Oplossingen

Ze stellen een Serrin-type regulariteitscriterium op. Als de controle $f$ behoort tot $L^{5/2}(0, T; L^{5/2+}(\Omega_c))$ en de oplossing voldoet aan $u \in L^{5p/2}(Q) \cap L^{10/3}(Q)$, dan is de oplossing uniek en heeft het de volgende hoge regulariteit:

• $u, \nabla v \in L^5(Q)$ voor $p \leq 2$.
• $u, \nabla v \in L^{5(p-1)}(Q)$ voor $p > 2$.
  Deze regulariteit is essentieel voor de goedgesteldeheid van het lineaire systeem en het optimaal beheerprobleem.

C. Optimaal Bilineair Beheer

Het artikel behandelt een optimaal beheerprobleem met een kostenfunctionaal $J$ die de afwijking van de gewenste toestanden ($u_d, v_d$) en de kosten van de controle $f$ minimaliseert.

• Bestaan van een optimum: Er wordt bewezen dat er minstens één globaal optimaal oplossing bestaat binnen de toelaatbare verzameling.
• Differentieerbaarheid: De auteurs bewijzen dat de afbeelding van controle naar toestand Fréchet-differentieerbaar is.
• Noodzakelijke optimaliteitsvoorwaarden: Door het gebruik van het geadjungeerde systeem (dat de gradiënt van de gereduceerde kostenfunctionaal identificeert), leiden ze eerste-orde noodzakelijke optimaliteitsvoorwaarden af voor lokale optimale oplossingen. Dit resulteert in een variatiële ongelijkheid die de optimale controle $f^*$ karakteriseert.

---

4. Significatie en Impact

• Vulling van een Kennislacune: Tot nu toe waren de meeste resultaten voor optimaal beheer van chemotaxis-modellen beperkt tot 2D-domeinen of lineaire productie ($p=1$). Dit artikel vult een belangrijke lacune door de analyse uit te breiden naar 3D-domeinen met niet-lineaire productie ($p > 1$).
• Rol van de Logistieke Term: Het werk illustreert wiskundig precies wanneer de logistische term noodzakelijk is (namelijk voor $p > 5/3$) om de regulariteit en het bestaan van oplossingen te garanderen in drie dimensies.
• Methodologische Vooruitgang: De benadering om de optimaliteitsvoorwaarden af te leiden via de Implicit Function Theorem en een expliciet geadjungeerd systeem (in plaats van alleen abstracte Lagrange-multiplicatoren) biedt een directere route naar de gradiënt van de kostenfunctionaal. Dit is cruciaal voor het ontwikkelen van numerieke algoritmen (zoals gradiëntafdaalmethoden) voor het oplossen van dergelijke complexe biologische beheerproblemen.
• Toepassingsgebied: De resultaten zijn relevant voor het modelleren van biologische processen zoals wondheling, kankerinvasie en predator-prooi interacties, waarbij het mogelijk is om chemische signalen lokaal te manipuleren om het celgedrag te sturen.

Samenvattend biedt dit artikel een rigoureuze wiskundige basis voor het begrijpen en optimaliseren van complexe chemo-repulsie systemen in drie dimensies, met een sterke focus op de interactie tussen niet-lineariteit, regulariteit en controle-theorie.