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

1. Il Problema Studiato

Il lavoro si concentra sull'analisi matematica e sul controllo ottimo di un sistema di equazioni alle derivate parziali (PDE) paraboliche-paraboliche che modella la chemorepulsione in un dominio tridimensionale limitato $\Omega \subset \mathbb{R}^3$. Il sistema descrive l'interazione tra la densità cellulare $u(t,x)$ e la concentrazione chimica $v(t,x)$.

Il modello considerato include:

• Termine di chemorepulsione: $\nabla \cdot (u \nabla v)$, che descrive il movimento delle cellule lontano da concentrazioni chimiche elevate.
• Produzione non lineare del segnale chimico: Il termine $u^p$ nell'equazione per $v$, con $p > 1$.
• Reazione logistica: Nel caso "logistico", un termine di crescita e competizione $ru - \mu u^p$ nell'equazione per $u$.
• Controllo bilineare: Una funzione di controllo $f(t,x)$ che agisce su un sottodominio $\Omega_c \subseteq \Omega$ nell'equazione per $v$, modellata come $f v \mathbb{1}_{\Omega_c}$. Questo permette di rappresentare l'iniezione o l'estrazione di sostanze chimiche mantenendo la non-negatività delle variabili di stato.

Il sistema è soggetto a condizioni al contorno di tipo Neumann (flusso nullo) e condizioni iniziali non negative.

2. Metodologia

Gli autori adottano un approccio rigoroso basato sull'analisi funzionale e sulla teoria delle equazioni paraboliche non lineari:

• Stime Energetiche e Esistenza di Soluzioni Deboli:

  • Vengono derivate disuguaglianze energetiche fondamentali. Si dimostra che per $1 < p \le 5/3$, la diffusione è sufficiente a controllare i termini non lineari e garantire l'esistenza di soluzioni deboli globali senza termine logistico.
  • Per $p > 5/3$, la produzione non lineare diventa troppo forte; è necessario introdurre il termine logistico ($ru - \mu u^p$) per ottenere stime a priori che garantiscano l'esistenza globale.
  • L'esistenza è provata tramite un metodo di regolarizzazione (introducendo un parametro $\epsilon$) e utilizzando il teorema del punto fisso di Leray-Schauder, seguito da un passaggio al limite.

• Criterio di Regolarità (Soluzioni Forti):

  • Viene stabilito un criterio di regolarità di tipo Serrin. Si dimostra che se la soluzione debole soddisfa condizioni di integrabilità specifiche (in particolare $u \in L^{5p/2}(Q) \cap L^{10/3}(Q)$), allora la soluzione diventa "forte" (più regolare).
  • Questo risultato è ottenuto attraverso un argomento di bootstrapping non standard su un'equazione lineare di diffusione-convezione-reazione con coefficienti non regolari.
  • La regolarità richiesta per il controllo è $f \in L^{5/2}(0, T; L^{5/2+}(\Omega_c))$.

• Teoria del Controllo Ottimo:

  • Viene formulato un problema di controllo ottimo bilineare con un funzionale di costo che penalizza la deviazione dallo stato desiderato e il costo del controllo.
  • Per caratterizzare le condizioni di ottimalità, gli autori evitano l'uso generico dei moltiplicatori di Lagrange. Invece, utilizzano il Teorema della Funzione Implicita negli spazi di Banach per dimostrare che l'operatore "controllo-stato" è Fréchet-differenziabile.
  • Viene poi risolto un sistema aggiunto (duale) per identificare esplicitamente il gradiente del funzionale di costo ridotto.

3. Contributi Chiave e Risultati Principali

1. Esistenza Globale di Soluzioni Deboli in 3D:

   • Per il sistema senza termine logistico, è stata provata l'esistenza di soluzioni deboli globali per $1 < p \le 5/3$.
   • Per il caso $p > 5/3$, l'esistenza è garantita solo se è presente il termine logistico di competizione, che stabilizza l'energia del sistema.
   • Questo colma un vuoto nella letteratura, che fino ad ora aveva trattato principalmente casi lineari ($p=1$) o bidimensionali.

2. Criterio di Regolarità e Unicità:

   • È stato stabilito un criterio sufficiente affinché una soluzione debole diventi una soluzione forte e unica. La regolarità ottenuta è $u, \nabla v \in L^5(Q)$ per $p \le 2$ e $u, \nabla v \in L^{5(p-1)}(Q)$ per $p > 2$.
   • Questa regolarità è critica per la ben-postezza del sistema linearizzato e per l'applicazione del teorema della funzione implicita.

3. Esistenza di Soluzioni Ottimali Globali:

   • È stata dimostrata l'esistenza di almeno una soluzione globale ottimale per il problema di controllo, utilizzando un argomento di successione minimizzante e la semicontinuità inferiore debole del funzionale di costo.

4. Condizioni di Ottimalità del Primo Ordine:

   • Grazie alla differenziabilità dell'operatore controllo-stato e alla risoluzione del sistema aggiunto, gli autori hanno derivato le condizioni necessarie di ottimalità del primo ordine per soluzioni locali.
   • Il gradiente del funzionale di costo ridotto è identificato esplicitamente come:
     $$J'_0(f^*)[f] = \int_0^T \int_{\Omega_c} \left( \gamma_f \frac{\text{sgn}(f^*)|f^*|^{3/2}}{\|f^*\|_{L^{5/2+}}} + v^* \eta \right) f \, dx dt$$
     dove $\eta$ è la seconda componente della soluzione del sistema aggiunto.

4. Significato e Implicazioni

• Avanzamento Teorico: Questo lavoro estende significativamente la teoria dei modelli di chemotassi/chemorepulsione al caso tridimensionale con produzione non lineare ($p>1$), un regime matematicamente molto più complesso e meno compreso rispetto al caso lineare o 2D.
• Gestione della Non Linearità: Dimostra come il termine logistico possa essere essenziale per la stabilità energetica in regimi di produzione chimica forte ($p > 5/3$).
• Applicabilità Numerica: La derivazione esplicita del gradiente del funzionale di costo (tramite il sistema aggiunto) è un risultato cruciale per lo sviluppo di schemi numerici basati su metodi di discesa del gradiente. Questo apre la strada a simulazioni numeriche efficienti per il controllo di processi biologici complessi.
• Metodologia Robusta: L'uso combinato di stime energetiche, criteri di regolarità non standard e il teorema della funzione implicita offre un quadro metodologico solido per affrontare problemi di controllo ottimo in sistemi PDE non lineari accoppiati in dimensioni superiori.

In sintesi, il paper fornisce una trattazione completa, dall'esistenza delle soluzioni alla caratterizzazione delle condizioni di ottimalità, per un modello di chemorepulsione 3D complesso, ponendo le basi teoriche per future applicazioni di controllo in biologia matematica.