Existence result for a 2 x 2 system of conservation laws with discontinuous flux and applications

1. Problème Étudié

L'article s'intéresse à la résolution du problème de Cauchy pour un système de lois de conservation hyperbolique $2 \times 2$ en dimension un, caractérisé par un flux discontinu par rapport à la variable spatiale $x$. Le système est de la forme :

$$\begin{cases} \partial_t \rho + \partial_x \left( c(x) V\left(\frac{q}{\rho} - p(\rho)\right) \rho \right) = 0 \\ \partial_t q + \partial_x \left( c(x) V\left(\frac{q}{\rho} - p(\rho)\right) q \right) = 0 \end{cases}$$

où $U = (\rho, q)$ représente les variables conservatives (densité et quantité de mouvement).

Caractéristiques principales du modèle :

• Flux Discontinu : Le coefficient $c(x)$ est une fonction en escalier (constante par morceaux), discontinue en un nombre fini de points fixes $\xi_i$. Cela modélise des changements abrupts dans les propriétés de la route (capacité, loi de vitesse).
• Type Temple : Le système appartient à la classe des systèmes de Temple, où les courbes de choc et de détente pour le champ non-linéaire coïncident.
• Absence d'hypothèse de monotonie : Contrairement à de nombreux travaux antérieurs, l'article n'impose aucune hypothèse de monotonie sur l'application $x \mapsto f(x, U)$.
• Vitesse caractéristique : La deuxième vitesse caractéristique est autorisée à s'annuler, ce qui rend le cadre non-caractéristique plus général et plus difficile à traiter.
• Loi de pression : Le choix d'une pression logarithmique $p(\rho) = V_{ref} \ln(\rho/R_{ref})$ est crucial pour éviter l'apparition de vide ($\rho > 0$) tout en permettant une analyse rigoureuse.

2. Méthodologie

La preuve d'existence repose sur la méthode de tracking des fronts d'ondes (wave-front tracking), une technique classique pour les systèmes hyperboliques, mais adaptée ici aux spécificités des flux discontinus.

Étapes clés de la méthode :

1. Approximation des données initiales : Les données initiales $U_0$ sont approchées par des fonctions constantes par morceaux.
2. Constructeur de Solveur de Riemann Adapté :
   • En dehors des points de discontinuité de $c(x)$, le solveur de Riemann standard est utilisé.
   • Aux points de discontinuité $\xi_i$, un solveur de Riemann adapté est construit. Ce solveur détermine les états intermédiaires $\hat{U}$ et $\check{U}$ en satisfaisant les conditions de Rankine-Hugoniot stationnaires.
   • Innovation majeure : Les auteurs introduisent un nouveau système de coordonnées basé sur les invariants de Riemann, noté $\tilde{W} = (\tilde{w}_1, w_2)$, où $\tilde{w}_1$ est une version "adaptée" de l'invariant $w_1$ tenant compte du coefficient $c$. Cette transformation est essentielle car le domaine physique standard n'est pas invariant sous l'action du solveur de Riemann aux interfaces.
3. Estimations Uniformes :
   • L'existence globale de la solution approchée $U_n$ est garantie par une borne uniforme sur le nombre de fronts d'ondes. Bien que le nombre d'ondes puisse augmenter lors des interactions, cela ne se produit qu'aux discontinuités de $c$, et l'augmentation est contrôlée (au plus +1).
   • Contrôle de la Variation Totale (TV) : Une difficulté majeure est que la variation totale de la solution $TV(U_n(t))$ ou même celle des invariants de Riemann standards $TV(W(U_n(t)))$ peut augmenter après une interaction. Pour contourner cela, les auteurs prouvent que la variation totale exprimée dans les coordonnées adaptées $\tilde{W}$, notée $TV(\tilde{W}(U_n(t)))$, reste bornée uniformément dans le temps.
4. Convergence et Limites :
   • Grâce au théorème de Helly et aux estimations de variation totale sur $\tilde{W}$, on extrait une sous-suite convergente vers une fonction limite $U$.
   • La limite est prouvée être une solution faible du problème.
5. Condition d'Entropie :
   • Pour sélectionner la solution physique unique, une condition d'entropie de type Kruzhkov est introduite.
   • Cette condition inclut un terme compensatif $N_\kappa$ aux points de discontinuité de $c$, qui rend compte de la dissipation d'entropie supplémentaire générée par l'interface. La forme de ce terme est déduite de la construction du solveur de Riemann adapté.

3. Résultats Principaux

Le résultat central est l'énoncé du Théorème 2.8 :

• Existence Globale : Pour des données initiales de grande amplitude (dans un espace $L^\infty$ avec variation totale bornée dans les coordonnées adaptées), il existe une solution globale d'entropie au sens de la définition 2.6.
• Stabilité : La solution satisfait des estimations de stabilité $L^1$ et une borne sur la variation totale dans les coordonnées adaptées :
  $$TV(c S(w_1(U(t)))) + TV(w_2(U(t))) \leq TV(c S(w_1(U_0))) + TV(w_2(U_0))$$
  où $S$ est une fonction strictement croissante liée au flux.
• Unicité (Partielle) : Bien que l'unicité ne soit pas prouvée dans cet article (c'est une question ouverte pour ce type de système avec flux discontinu), la condition d'entropie proposée est conçue pour sélectionner la solution physiquement pertinente.

4. Contributions Techniques et Innovations

1. Coordonnées Invariants Adaptées ($\tilde{W}$) : L'introduction de la coordonnée $\tilde{w}_1(U, c) = c S(w_1(U))$ est la contribution technique la plus significative. Elle permet de restaurer la propriété d'invariance du domaine et de contrôler la variation totale là où les coordonnées classiques échouent.
2. Solveur de Riemann aux Interfaces : La construction explicite d'un solveur de Riemann pour les discontinuités de flux qui respecte les contraintes physiques (domaine invariant relaxé) et les conditions d'entropie.
3. Condition d'Entropie avec Terme Compensatif : La formulation d'une inégalité d'entropie qui intègre explicitement un terme de bord aux points de discontinuité de $c$, justifié par la structure du solveur de Riemann.
4. Généralité : Le résultat s'applique sans hypothèse de monotonie sur $c(x)$ et permet à la vitesse caractéristique de s'annuler, couvrant des cas plus généraux que la littérature précédente sur les systèmes de Temple.

5. Signification et Applications

Modélisation du Trafic Routier :
L'application principale discutée est la modélisation du trafic routier sur des routes hétérogènes.

• Le système (2.1) généralise le modèle Aw-Rascle-Zhang (ARZ) et le modèle GARZ (Generalized ARZ).
• Le coefficient $c(x)$ représente des changements collectifs abrupts de la capacité de la route ou de la loi de vitesse (ex: passage d'une autoroute à une route nationale, zones de travaux, péages).
• Le modèle permet de capturer des phénomènes complexes comme les ondes de choc stationnaires et les changements de régime de flux dus à ces hétérogénéités spatiales.

Impact Théorique :
Ce travail comble un vide dans la théorie des systèmes de conservation lois avec flux discontinu, en particulier pour les systèmes $2 \times 2$ non-symétriques de type Temple. Il offre un cadre robuste pour l'analyse mathématique et le développement de schémas numériques pour des problèmes d'ingénierie réels où les paramètres spatiaux varient de manière discontinue.

En résumé, cet article fournit une preuve rigoureuse d'existence pour une classe importante de systèmes hyperboliques avec flux discontinu, en surmontant les difficultés liées à la perte de contrôle de la variation totale grâce à une transformation intelligente des coordonnées et à une condition d'entropie adaptée.