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

1. Problématique et Modèle Mathématique

L'article étudie un système de réactions-diffusions couplées de type parabolique-parabolique en dimension trois ($\mathbb{R}^3$), modélisant un phénomène de chémorépulsion (mouvement des cellules vers des concentrations chimiques plus faibles) avec une production non linéaire du signal chimique et un terme de réaction logistique.

Le système est défini sur un domaine borné $\Omega \subset \mathbb{R}^3$ et un intervalle de temps $(0, T)$, noté $Q = (0, T) \times \Omega$. Il comprend deux équations pour la densité cellulaire $u(t,x) \ge 0$ et la concentration chimique $v(t,x) \ge 0$ :

1. Équation de la densité cellulaire :
   $$\partial_t u - \Delta u = \nabla \cdot (u \nabla v) + r u - \mu u^p$$
   Le terme $\nabla \cdot (u \nabla v)$ décrit le mécanisme de répulsion. Le terme $ru - \mu u^p$ représente une croissance logistique (prolifération $r$ et compétition $\mu$). Si $r=\mu=0$, le système est dit « non-logistique ».

2. Équation de la concentration chimique :
   $$\partial_t v - \Delta v + v = u^p + f v \mathbf{1}_{\Omega_c}$$
   Le terme $u^p$ ($p > 1$) modélise une production non linéaire du signal par les cellules. Le terme $f v \mathbf{1}_{\Omega_c}$ représente un contrôle bilinéaire agissant sur un sous-domaine $\Omega_c \subseteq \Omega$. Ici, $f(t,x)$ est la fonction de contrôle qui module le taux de prolifération ou de dégradation du produit chimique.

Le système est complété par des conditions initiales non négatives et des conditions aux limites de type Neumann homogène (flux nul).

2. Méthodologie

Les auteurs adoptent une approche analytique rigoureuse divisée en plusieurs étapes clés :

• Estimations Énergétiques et Existence de Solutions Faibles :
  L'analyse commence par la dérivation d'inégalités d'énergie. Une difficulté majeure réside dans le terme de production non linéaire $u^p$.

  • Pour $1 < p \le 5/3$, la diffusion seule suffit à contrôler l'énergie et à prouver l'existence de solutions faibles globales pour le système non-logistique.
  • Pour $p > 5/3$, la diffusion est insuffisante pour contrôler l'explosion potentielle de l'énergie. Les auteurs montrent alors que l'ajout du terme logistique ($-\mu u^p$) est nécessaire pour obtenir des estimations a priori et prouver l'existence de solutions faibles.
  • La preuve utilise une méthode de régularisation (ajout d'un paramètre $\epsilon$ dans une équation elliptique pour $v$) suivie d'un passage à la limite via le théorème du point fixe de Leray-Schauder et des lemmes de compacité (Aubin-Lions).

• Critère de Régularité (Solutions Fortes) :
  Pour aborder le problème de contrôle optimal, l'existence de solutions faibles ne suffit pas ; il faut des solutions plus régulières (fortes). Les auteurs établissent un critère de régularité de type Serrin : si une solution faible satisfait $u \in L^{5p/2}(Q) \cap L^{10/3}(Q)$ et que le contrôle $f$ appartient à un espace légèrement supérieur à $L^{5/2}$, alors la solution devient unique et plus régulière ($u, \nabla v \in L^5(Q)$).
  Cette régularité est obtenue via un argument de « bootstrapping » non standard sur une équation de diffusion-convection-réaction linéaire aux coefficients non réguliers.

• Théorie du Contrôle Optimal :
  Une fois la régularité établie, les auteurs définissent un problème de contrôle optimal bilinéaire visant à minimiser un fonctionnel de coût $J$ (écart à des états désirés $u_d, v_d$ et coût du contrôle).

  • Existence d'optima : Démontrée via une suite minimisante et la compacité faible.
  • Différentiabilité : Au lieu d'utiliser des multiplicateurs de Lagrange génériques, les auteurs utilisent le Théorème de la Fonction Implicite dans des espaces de Banach pour prouver que l'opérateur « contrôle-état » est Fréchet-différentiable localement.
  • Système Adjoint : La dérivée du coût réduit est identifiée en résolvant un système adjoint par une argument de dualité. Cela permet d'obtenir des conditions d'optimalité du premier ordre explicites.

3. Contributions Clés et Résultats Principaux

1. Existence de Solutions Faibles Globales en 3D :

   • Pour le système sans logistique, l'existence est prouvée pour $1 < p \le 5/3$.
   • Pour $p > 5/3$, l'existence est prouvée uniquement si le terme logistique est présent. C'est une contribution majeure car la régularité en 3D pour ce type de systèmes non linéaires est un problème ouvert difficile.

2. Critère de Régularité et Unicité :
   Établissement d'un critère de régularité suffisant pour garantir que les solutions faibles deviennent des solutions fortes uniques. Ce critère repose sur la régularité du contrôle $f \in L^{5/2}(L^{5/2+})$ et la régularité de l'état $u$.

3. Résolution du Problème de Contrôle Optimal :

   • Preuve de l'existence de solutions globales optimales.
   • Dérivation de conditions d'optimalité nécessaires du premier ordre pour des solutions locales.
   • Identification explicite du gradient du fonctionnel de coût réduit via le système adjoint, ce qui est crucial pour la conception d'algorithmes numériques basés sur la descente de gradient.

4. Analyse Fine des Exposants :
   L'article identifie précisément les seuils critiques ($p=5/3$) où le comportement du système change, nécessitant l'intervention du terme logistique pour la stabilité énergétique.

4. Signification et Impact

Ce travail comble un vide significatif dans la littérature mathématique sur les modèles de chimiotaxie :

• Dimension 3 : La plupart des résultats antérieurs sur le contrôle optimal pour la chimiotaxie étaient limités à la dimension 2. Ce papier étend la théorie au cas tridimensionnel, beaucoup plus complexe et pertinent biologiquement.
• Production Non Linéaire : Il traite le cas de la production non linéaire ($p > 1$), qui est plus réaliste biologiquement mais mathématiquement plus difficile que le cas linéaire ($p=1$).
• Contrôle Bilineaire : L'analyse du contrôle bilinéaire (où le contrôle multiplie l'état) est traitée avec rigueur, permettant de maintenir la positivité des variables d'état sans restrictions de signe sur le contrôle.
• Applications Numériques : En fournissant une expression explicite du gradient du coût via le système adjoint, l'article ouvre la voie à des schémas numériques efficaces pour le contrôle de processus biologiques (ex: inhibition de la croissance tumorale, gestion de populations cellulaires).

En résumé, cet article fournit une fondation théorique solide pour l'analyse et le contrôle de systèmes de chimiorépulsion complexes en 3D, en surmontant les obstacles liés à la non-linéarité et à la dimension spatiale élevée.