Systems of Nonlocal Conservation Laws with Memory and Their Zero Retention Limit

1. Énoncé du problème

L'article traite de l'analyse mathématique et de l'approximation numérique de systèmes de lois de conservation hyperboliques non locales qui intègrent des effets de mémoire (non-localité à la fois dans l'espace et dans le temps).

Le système gouvernant pour $N$ équations couplées est donné par :
$$\partial_t U^k + \partial_x \left( f^k(U^k) \nu^k((\Theta \ast U)_k) \right) = 0, \quad (t, x) \in (0, T) \times \mathbb{R}$$
où :

• $U = (U^1, \dots, U^N)$ est le vecteur des inconnues.
• Le flux dépend d'un opérateur non local $(\Theta \ast U)_k$, qui est une convolution impliquant à la fois un noyau spatial $\mu$ et un noyau de mémoire temporelle $\Gamma$.
• Plus précisément, $(\Theta \ast U)_k(t, x) = \sum_{j} \int_0^t \int_{\mathbb{R}} U^j(\tau, \xi) \mu_{j,k}(x-\xi) \Gamma_{j,k}(t-\tau) \, d\xi \, d\tau$.

Défis clés :

1. Absence de propriété de semi-groupe : La convolution temporelle détruit la propriété de semi-groupe standard trouvée dans les lois de conservation locales, rendant les techniques classiques de bien-posé inapplicables.
2. Flux discontinu : Le terme de convolution temporelle manque souvent de régularité spatiale, rendant le flux discontinu par rapport à la variable spatiale $x$.
3. Limites asymptotiques : Il est nécessaire de comprendre le comportement des solutions lorsque le noyau de mémoire rétrécit (la limite « de mémoire à sans mémoire ») et de savoir si les schémas numériques restent compatibles avec cette limite (compatibilité asymptotique).

2. Méthodologie

Les auteurs emploient une combinaison de techniques analytiques et de méthodes numériques de volumes finis :

• Cadre des solutions entropiques : Ils définissent un concept de solution entropique adapté aux systèmes non locaux pour assurer l'unicité, car les solutions faibles ne sont pas uniques en raison de la non-linéarité.
• Approximation par volumes finis : Un schéma de volumes finis convergent d'ordre un est proposé. Le schéma utilise une formule de marche pour mettre à jour les moyennes de cellules, intégrant le terme de convolution non local via des approximations discrètes des noyaux.
• Stabilité et bornes BV : Les auteurs dérivent des estimations uniformes de Variation Bornée (BV) pour les solutions numériques. Crucialement, ils prouvent que ces bornes sont indépendantes de la norme $L^\infty$ du noyau temporel $\Gamma$, dépendant uniquement de sa norme $L^1$. Ceci est essentiel pour l'analyse asymptotique.
• Estimations d'entropie relative : Pour analyser la convergence lorsque le noyau de mémoire rétrécit, ils utilisent une estimation d'entropie relative (doublement des variables de type Kruzhkov) adaptée aux lois de conservation non locales.
• Analyse asymptotique : Ils étudient la limite lorsque le support du noyau temporel $\Gamma_\delta$ rétrécit vers un delta de Dirac (c'est-à-dire $\delta \to 0$), supprimant efficacement l'effet de mémoire.

3. Contributions clés

A. Théorie de bien-posé

• Existence et unicité : L'article établit l'existence et l'unicité des solutions entropiques pour le système couplé avec des flux non linéaires généraux et des noyaux réguliers généraux.
• Dépendance continue : Ils prouvent la dépendance continue par rapport aux données initiales, montrant que l'opérateur solution est lipschitzien dans la norme $L^1$.
• Généralisation : Contrairement aux travaux antérieurs restreints à des noyaux spécifiques ou des flux linéaires, ce cadre gère des systèmes couplés $N$ généraux avec des flux non linéaires $f^k$ et des noyaux réguliers généraux.

B. Analyse asymptotique (Limite de mémoire à sans mémoire)

• Convergence vers la limite non locale spatiale : Les auteurs prouvent que lorsque le rayon du noyau temporel $\delta \to 0$, la solution entropique du système non local espace-temps converge vers la solution entropique du système purement non local spatial (où la mémoire est supprimée).
• Taux de convergence : Ils dérivent une estimation d'erreur explicite d'ordre $O(\sqrt{\delta})$ pour la convergence des solutions exactes.

C. Compatibilité asymptotique des schémas numériques

• Compatibilité du schéma : L'article prouve que le schéma de volumes finis proposé est asymptotiquement compatible. Cela signifie que lorsque à la fois la taille de maillage $\Delta x \to 0$ et le rayon de mémoire $\delta \to 0$, la solution numérique converge vers la solution entropique de la limite sans mémoire (purement spatiale).
• Estimations d'erreur : L'erreur totale pour le schéma numérique dans la limite sans mémoire est estimée à $O(\sqrt{\delta} + \sqrt{\Delta x})$.
• Nouveauté : Il s'agit du premier travail fournissant une telle convergence et analyse asymptotique pour les schémas de volumes finis appliqués aux lois de conservation non locales avec mémoire.

D. Assouplissement des contraintes géométriques

• L'analyse est menée sans imposer de restrictions géométriques (telles que la convexité) sur les noyaux de convolution spatiale ou temporelle, ce qui constitue une amélioration significative par rapport à la littérature existante sur les lois non locales uniquement dans l'espace.

4. Résultats clés

1. Théorème 1 (Unicité) : L'unicité des solutions entropiques est établie avec une estimation de stabilité dépendante des données initiales et des constantes de Lipschitz des flux et des fonctions non locales.
2. Théorème 3 (Convergence) : Les approximations par volumes finis convergent vers la solution entropique unique du problème de Cauchy.
3. Théorème 6 (Limite de la solution exacte) : Lorsque $\delta \to 0$, $\|U_\delta(T) - U(T)\|_{L^1} = O(\sqrt{\delta})$, où $U_\delta$ est la solution avec mémoire et $U$ est la solution sans mémoire.
4. Théorème 7 (Limite numérique) : La solution numérique $U^\Delta_\delta$ converge vers la solution sans mémoire $U$ avec un taux $O(\sqrt{\delta} + \sqrt{\Delta x})$, à condition que le noyau temporel ait un premier moment fini.
5. Expériences numériques :
   • Des simulations sur un modèle de type système de Keyfitz-Kranzer confirment les taux de convergence théoriques.
   • Les taux de convergence observés dans les expériences étaient légèrement meilleurs ($\alpha > 0.5$) que la borne théorique $O(\sqrt{\delta})$, validant la robustesse du schéma.
   • Des visualisations démontrent la transition des dynamiques « avec mémoire » vers les dynamiques « sans mémoire » lorsque $\delta$ diminue.

5. Importance et perspectives futures

Importance :

• Combler une lacune théorique : Ce travail comble une lacune critique dans la théorie des lois de conservation hyperboliques en traitant rigoureusement des systèmes avec à la fois une non-localité spatiale et temporelle.
• Modélisation pratique : Les résultats sont directement applicables à la modélisation de phénomènes avec mémoire atténuée, tels que l'écoulement du trafic avec les temps de réaction des conducteurs, la viscoélasticité et le transport de gaz dans les milieux poreux.
• Fiabilité numérique : En prouvant la compatibilité asymptotique, les auteurs garantissent que les simulations numériques de modèles basés sur la mémoire approximeront correctement les modèles plus simples, sans mémoire, dans la limite, empêchant ainsi les artefacts numériques lorsque les effets de mémoire sont négligeables.

Perspectives futures et problèmes ouverts :

• Non-localité uniquement dans le temps : Le cas où la non-localité spatiale est supprimée (mémoire purement temporelle) reste ouvert pour les flux non linéaires, car le système se comporte comme une loi locale avec flux discontinu, perdant les bornes BV spatiales.
• Limite entièrement locale : Le passage complet du système non local espace-temps vers un système entièrement local (suppression de $\mu$ et $\Gamma$) pour $N > 1$ avec des flux non linéaires reste un problème ouvert et difficile en raison du mal-posé des systèmes couplés locaux correspondants.
• Extensions : Les travaux futurs traiteront des termes sources, de noyaux moins réguliers et de l'analyse d'erreur pour la transition non locale vers locale.