Résumé technique : Distances de $p$-Wasserstein sur les réseaux et convergence de la 3D vers la 1D

1. Énoncé du problème

Le papier traite de la modélisation mathématique des réseaux de transport de gaz à l'aide de la théorie du transport optimal. Plus précisément, il étudie deux problèmes interconnectés :

1. Transport dynamique sur des graphes métriques : Étendre la formulation dynamique de la distance de Wasserstein (Benamou-Brenier) aux graphes métriques représentant des réseaux de gaz. Cela implique de gérer la conservation de la masse sur les arêtes et de définir des conditions de couplage appropriées aux sommets (jonctions de conduites). Les auteurs considèrent deux scénarios physiques distincts :
   • Conditions de Kirchhoff : La masse est conservée globalement aux nœuds sans stockage (jonctions de conduites classiques).
   • Stockage de masse aux nœuds : La masse peut être stockée aux nœuds, modélisée par des équations de continuité supplémentaires et des lois de Kirchhoff généralisées.
2. Réduction de dimension (3D vers 1D) : Les auteurs cherchent à justifier rigoureusement la pratique courante de l'ingénierie consistant à modéliser les conduites de gaz comme des entités unidimensionnelles (graphes métriques) plutôt que comme des cylindres tridimensionnels. Ils étudient le comportement asymptotique de la distance de Wasserstein statique sur un domaine de réseau 3D lorsque le diamètre de la conduite ($\varepsilon$) converge vers zéro, prouvant que le coût de transport 3D converge vers le coût de transport du graphe métrique 1D.

2. Méthodologie

Formulation dynamique sur les graphes métriques

Les auteurs généralisent l'équation de continuité aux graphes métriques. Au lieu d'un champ de vecteurs vitesse, ils utilisent la quantité de mouvement (flux de masse) pour assurer la linéarité dans l'équation de continuité et la convexité dans le problème d'optimisation.

• Cas du stockage de masse : Ils définissent une équation de continuité où la dérivée temporelle de la densité du nœud ($\gamma_v$) est égale au flux net entrant dans le nœud depuis les arêtes adjacentes. Cela conduit à une distance de Wasserstein dynamique définie via une fonctionnelle d'action impliquant une fonction perspective $h(a,b) = |a|^p/b^{p-1}$, pénalisant les densités négatives et les flux non nuls à densité nulle.
• Cas de Kirchhoff : Ils utilisent une approche par espace quotient où les arêtes sont collées aux nœuds. L'équation de continuité est formulée à l'aide de fonctions tests continues sur l'ensemble du graphe, imposant implicitement la loi de Kirchhoff (la somme des flux est égale à zéro).
• Écoulements de gradient : Le papier examine comment des modèles d'écoulement de gaz spécifiques apparaissent comme des écoulements de gradient dans ces espaces de Wasserstein. Notamment, le modèle ISO3 pour l'écoulement de gaz est identifié comme un écoulement de gradient $W_{K,3}$ d'une fonctionnelle d'énergie spécifique, tandis que les équations de diffusion-dérive correspondent à des écoulements de gradient $W_{K,2}$.

Analyse de convergence 3D vers 1D

Pour prouver la convergence de la distance de Wasserstein d'un domaine 3D ($N_\varepsilon$) vers un graphe métrique 1D ($N_0$), les auteurs emploient les étapes suivantes :

• Définition de la fonction de coût : Ils définissent une fonction de coût $c_\varepsilon(x,y)$ sur le domaine 3D comme l'infimum de l'intégrale d'action sur des courbes contraintes de rester dans le réseau 3D. Ce coût incorpore une fonction indicatrice $\iota_{N_\varepsilon}$ pour pénaliser les chemins sortant du réseau.
• Plans de transport optimaux : Ils analysent les plans de transport optimaux ($\pi_\varepsilon$). Un défi majeur identifié est la non-unicité potentielle des applications optimales en raison des géodésiques bifurquant aux jonctions de réseaux (où plusieurs chemins de longueur égale existent).
• Monotonie cyclique : La preuve repose largement sur la caractérisation des plans de transport optimaux comme des ensembles $c$-cycliquement monotones. En établissant que les fonctions de coût $c_\varepsilon$ convergent vers le coût 1D $c_0$ et en utilisant la stabilité de la $c$-monotonie cyclique, ils démontrent la convergence des coûts et des plans de transport optimaux.

3. Contributions clés et résultats

Résultats théoriques sur les graphes métriques

• Existence et unicité (Well-posedness) : Le papier établit l'existence de courbes absolument continues dans l'espace $p$-Wasserstein sur les graphes métriques, les reliant à des solutions faibles de l'équation de continuité avec une énergie $p$ finie (Théorème 1).
• Identification de l'écoulement de gradient : Il confirme que le modèle d'écoulement de gaz ISO3 est un écoulement de gradient par rapport à la distance $W_{K,3}$, fournissant une interprétation variationnelle du modèle physique.

Résultats de convergence (3D vers 1D)

• Convergence des coûts : Le résultat principal (Proposition 1) stipule que pour des séquences de mesures de probabilité absolument continues convergeant faiblement vers des mesures sans atome sur le graphe 1D, le coût de transport optimal $OT(\mu_\varepsilon, \nu_\varepsilon, c_\varepsilon)$ converge vers $OT(\mu_0, \nu_0, c_0)$ quand le diamètre de la conduite $\varepsilon \to 0$.
• Convergence des plans : Le Théorème 3 prouve que si le plan de transport optimal pour la limite 1D est unique, la séquence des plans optimaux 3D converge faiblement vers cette limite unique. Si le plan limite n'est pas unique, chaque sous-suite convergente des plans 3D converge vers un plan 1D optimal.
• Stabilité avec la topologie : Le Théorème 4 fournit des bornes sur la façon dont le coût de transport change lorsque la topologie du réseau est modifiée (par exemple, en supprimant ou en ajoutant des arêtes), montrant que la différence de coût est bornée par la différence maximale des fonctions de coût pondérées par la masse transportée.

Illustrations numériques

Le papier inclut des exemples numériques utilisant une version discrétisée dans l'espace du problème. Ceux-ci illustrent que les géodésiques dans le réseau peuvent bifurquer et s'effondrer aux jonctions, confirmant visuellement la difficulté théorique d'établir des applications de transport uniques en présence de structures de branchement.

4. Signification et affirmations

Le papier affirme fournir un fondement mathématique rigoureux pour l'utilisation des graphes métriques dans la modélisation des réseaux de gaz.

• Justification des modèles 1D : En prouvant la convergence de la distance de Wasserstein de la 3D vers la 1D, les auteurs justifient l'argument de "moyennage" utilisé en ingénierie pour traiter les conduites comme des entités 1D, spécifiquement dans le contexte du transport optimal et des écoulements de gradient.
• Cadre unifié : Le travail unifie la formulation dynamique des distances de Wasserstein avec la physique des réseaux de gaz, liant explicitement le modèle ISO3 aux écoulements de gradient dans l'espace $W_{K,3}$.
• Gestion de la non-unicité : Les auteurs reconnaissent et abordent explicitement le problème de la non-unicité des plans de transport causée par les géodésiques de branchement, un phénomène souvent négligé dans les contextes euclidiens standards. Ils suggèrent que le levage espace-temps (Section 3.5) peut être nécessaire pour dériver des applications de transport uniques dans ces exemples "bien élevés" mais présentant des branchements.

Les auteurs restent modestes concernant la convergence de la formulation dynamique ; ils déclarent que la convergence prouvée de la distance statique sert de "point de départ pour des investigations ultérieures" concernant la convergence de la formulation dynamique et des écoulements de gradient de Wasserstein.