A new comparison principle for discrete Volterra equations with an application to convex sweeping processes with infinite delays

1. Problématique et Contexte

L'article s'intéresse à l'analyse mathématique et numérique d'un nouveau type de processus de balayage (sweeping process) avec retards infinis, motivé par la modélisation de systèmes biologiques multi-cellulaires possédant une « mémoire adhésive ».

• Contexte physique : Les cellules adhèrent à leurs positions passées via des liaisons transitoires, créant une force de mémoire modélisée par un noyau de convolution. Contrairement aux modèles classiques où la dynamique est locale (dérivée temporelle), ici l'état présent dépend de l'histoire entière du système.
• Le défi géométrique : Le processus est contraint à rester dans un ensemble convexe mobile $C(t)$. La particularité majeure de ce modèle est que la projection de la position moyenne (pondérée par le noyau de mémoire) sur l'ensemble contraint $C(t)$ n'est pas infinitésimale. La distance de projection est de l'ordre de $O(1)$ (et non $O(h)$ comme dans le processus de Moreau classique), ce qui brise les estimations de régularité usuelles et complique l'obtention de la compacité nécessaire à l'existence de solutions.
• Le problème mathématique : Il s'agit d'établir l'existence de solutions pour des équations de Volterra discrètes avec des noyaux non exponentiels (décroissance algébrique possible) et d'utiliser ces résultats pour prouver la convergence vers un processus de balayage continu avec mémoire.

2. Méthodologie

Les auteurs développent une approche originale combinant analyse fonctionnelle et analyse numérique, reposant sur trois piliers méthodologiques :

A. Un nouveau principe de comparaison pour les équations de Volterra

Les méthodes classiques (transformée de Laplace, décomposition spectrale de la résolvante) échouent ici car elles ne fournissent pas de bornes uniformes par rapport au pas de discrétisation $h$, surtout pour des noyaux non exponentiels.

• Innovation : Les auteurs introduisent un argument sans résolvante basé sur un principe de comparaison.
• Technique clé : Ils construisent une sur-solution (super-solution) composée de deux parties :
  1. Un terme dominant (ordre principal) gérant la source.
  2. Un correcteur de couche initiale (initial layer corrector) qui compense spécifiquement les erreurs introduites par la « queue » (tail) du noyau de convolution.
• Ce correcteur permet de dominer le terme de droite de l'équation de Volterra sans avoir besoin de connaître la décomposition spectrale explicite de la résolvante, garantissant ainsi des bornes $L^\infty$ uniformes en $h$.

B. Estimations d'énergie et compacité

Pour le processus de balayage retardé, les auteurs reformulent le problème comme un flot gradient retardé (minimisation d'une énergie fonctionnelle).

• Ils établissent une estimation de décroissance d'énergie qui inclut un terme de dissipation positif.
• Cette estimation permet de déduire des bornes uniformes en $H^1$ sur les solutions discrètes, condition sine qua non pour appliquer le théorème d'Ascoli-Arzelà et obtenir la compacité des suites de solutions approchées.

C. Analyse de convergence vers le continu

En utilisant les bornes $H^1$ obtenues via le principe de comparaison, ils prouvent que les solutions discrètes convergent vers une solution continue satisfaisant une inclusion intégrale de type Moreau :
$$X(t) - \int_0^\infty \varrho(s) X(t-s) ds \in -N_{C(t)}(X(t))$$
où $N_{C(t)}$ est le cône normal à l'ensemble convexe.

3. Résultats Principaux

1. Théorème de comparaison uniforme (Théorèmes 1 et 2) :

   • Établissement d'un principe de comparaison pour les équations de Volterra continues et discrètes avec des noyaux non négatifs, monotones et à décroissance polynomiale (ou plus générale).
   • Preuve de l'existence de bornes $L^\infty(0, T)$ uniformes par rapport au pas de temps $h$, indépendamment de la décomposition spectrale.

2. Estimations d'énergie pour les processus retardés (Proposition 3) :

   • Démonstration que le processus de balayage discret satisfait une loi de dissipation d'énergie.
   • Preuve que la distance entre le point projeté et le point non contraint reste bornée (contrairement au cas classique où elle tend vers 0), mais que cela suffit pour contrôler la dynamique globale.

3. Convergence pour des ensembles circulaires mobiles (Théorème 6) :

   • Preuve rigoureuse de la convergence des solutions discrètes vers une solution continue pour le cas spécifique d'ensembles convexes mobiles (disques).
   • Utilisation cruciale du principe de comparaison pour obtenir les bornes $H^1$ nécessaires à la compacité.

4. Résultats numériques :

   • Simulation du processus avec des contraintes circulaires et en forme de stade (capsule) suivant des courbes de Lissajous.
   • Observation d'une convergence d'ordre 2 en $L^2$ lorsque les conditions initiales sont compatibles, et d'ordre 1 sinon.
   • Validation numérique que le produit « projection-déplacement » (terme clé dans les estimations d'énergie) reste borné lorsque $h \to 0$.

4. Contributions Clés

• Nouveau principe de comparaison : La contribution la plus significative est la construction d'un correcteur de couche initiale permettant de contourner l'absence de décomposition spectrale uniforme pour les noyaux non exponentiels. Cela ouvre la voie à l'analyse de systèmes avec mémoire à longue portée (décroissance algébrique).
• Modélisation de la projection non infinitésimale : L'article formalise mathématiquement le fait que, dans les systèmes à mémoire, la projection sur la contrainte peut être macroscopique ($O(1)$), ce qui nécessite de nouvelles techniques d'estimation par rapport au processus de Moreau classique.
• Existence pour les processus de balayage avec mémoire infinie : Première preuve d'existence pour ce type de processus discret avec des noyaux généraux, comblant un vide dans la littérature où la compacité était difficile à établir sans dérivée temporelle explicite.

5. Signification et Perspectives

• Impact théorique : Ce travail étend considérablement la théorie des processus de balayage et des équations de Volterra en fournissant des outils robustes pour les systèmes non locaux avec contraintes géométriques. Il démontre que la compacité peut être récupérée même en l'absence de régularité différentielle classique, grâce à la structure dissipative et aux nouvelles bornes de comparaison.
• Applications biologiques : Le modèle est directement applicable à la modélisation de la migration cellulaire dans des environnements structurés, où l'adhésion aux positions passées joue un rôle crucial.
• Perspectives :
  • Extension des résultats à des ensembles convexes généraux (au-delà des disques mobiles).
  • Traitement de contraintes prox-régulières (non convexes).
  • Analyse de la limite où le paramètre de rigidité des liaisons tend vers zéro, reliant le processus retardé au processus de Moreau classique.
  • Développement d'algorithmes numériques plus efficaces pour le cas multi-cellulaire.

En résumé, cet article propose une avancée majeure en combinant une analyse fine des équations intégro-différentielles discrètes avec une application géométrique complexe, offrant ainsi un cadre rigoureux pour l'étude des systèmes dynamiques à mémoire longue sous contraintes.