Pith. sign in

REVIEW 3 major objections 4 minor 39 references

Continuum of coupled Wasserstein gradient flows

T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read A continuum of coupled Wasserstein gradient flows gives the entropic relaxation of the AHT scheme an unconditional convergence guarantee that the original algorithm lacks.

desk verdict A substantial and mostly rigorous paper that proves existence, stability, and long-time convergence for a continuum of coupled Wasserstein gradient flows; the main caveat is a skipped compactness argument in the stability lemma that deserves a careful referee. read the letter →

arxiv 2411.13969 v1 pith:AMWTAOPC submitted 2024-11-21 math.AP

classification math.AP MSC 35K5135A1535A35
keywords multiphaseporousmediaflowfiberedWassersteingradientminimizingmovementschemeAngenent–Haker–Tannenbaumentropicoptimaltransportweaksolutionslong-timeasymptoticsdiscrete-to-continuouslimit
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper studies a system of drift-diffusion equations for a potentially infinite number of incompressible phases occupying a porous domain, each phase labelled by $y \in Y$ and all sharing a fixed volume-occupancy profile $\mu(x)$. It interprets this system as a gradient flow in the space of couplings between $\mu$ and $\nu$, measured with a fibered Wasserstein distance that transports mass separately inside every $y$-fiber. The main results are existence of weak solutions by a minimizing-movement construction, stability as the phase distribution $\nu$ approaches a continuum, and long-time convergence: every weak solution converges to the unique minimizer of the energy, which is exactly the objective of entropy-regularized optimal transport. Since that energy is the entropic version of the Angenent–Haker–Tannenbaum (AHT) scheme's cost, the paper establishes that this relaxation is unconditionally convergent, in contrast to the classical AHT scheme, which can stall at suboptimal stationary maps.

What carries the argument

The carrying object is the fibered Wasserstein distance $W_F$, which for couplings with the same second marginal $\nu$ disintegrates each coupling into fiber measures on $X$ and integrates their squared $W_2$ distances over $Y$. The construction proceeds by a minimizing-movement (JKO-type) scheme with steps $E_\tau(\rho\mid\bar\rho)=\frac{1}{2\tau}W_F(\rho,\bar\rho)^2+E(\rho)$; each step is equivalent to a convex program whose optimality conditions yield the discrete pressure $\Pi^n$ and the transport relation $\rho^{n-1}_y = (\mathrm{id}+\tau\nabla_X(V-\Pi^n+\kappa\log r^n))_\sharp\,\rho^n_y$. The pressure $\Pi$ is determined by the elliptic equation $\nabla\cdot(\mu\nabla\Pi)=\kappa\Delta\mu+\nabla\cdot(\int_Y \nabla_X V\, r\,d\nu)$, acting as the Lagrange multiplier for the volume constraint. The technical bottleneck is that both density and pressure converge only weakly, so the product $\rho\nabla\Pi$ is handled through strong $L^2$ convergence of vertical averages of the density, which gives the joint convergence of the product needed in the weak formulation. Long-time convergence is closed by showing that vanishing dissipation $I(\rho,\Pi)=\int|\nabla_X(V-\Pi+\kappa\log r)|^2\,d\rho$ forces $\rho$ to be the unique diagonal-scaling minimizer of $E$.

What would settle it

Take $X=Y=[0,1]$, $\mu=\nu$ uniform, $V(x,y)=|x-y|^2/2$, and $\kappa>0$; run the minimizing-movement scheme at decreasing $\tau$ from the flipped initial plan $\rho_0=(\mathrm{id},1-\mathrm{id})_\sharp\mu$. The theorem predicts $E(\rho_t)-E(\rho^*) \to 0$ along the whole trajectory; observing that the trajectory approaches any stationary point other than the entropic plan $\rho^*$, or that the energy gap plateaus above the discretization error as $\tau \to 0$, would falsify Lemma 16.

Watch

Extended reading notes

Core claim

The central claim is that the coupled system (1)–(3) is a well-posed gradient flow whose longtime limit solves entropic optimal transport. More precisely, for fixed marginals $\mu$ and $\nu$, positive regularization $\kappa$, and a finite-energy initial coupling, the minimizing-movement sequence built with the fibered Wasserstein distance $W_F$ yields a weak solution $(\rho_t, \Pi_t)$ of the coupled PDEs; a key step is showing that vertical averages of the density converge strongly enough to pass the product of pressure gradient and density to the limit, jointly over infinitely many phases. Any such weak solution converges as $t \to \infty$ to the unique minimizer of $E(\rho)=\int V\,d\rho+\kappa H(\rho)$, namely the entropic optimal transport plan between $\mu$ and $\nu$. The paper also shows that this convergence property is stable under approximation of the second marginal by finitely many phases, recovering multiphase porous-media models in the discrete-to-continuous limit. In the unregularized limit $\kappa=0$, by contrast, suboptimal stationary configurations persist, so relaxation alone does not fix the AHT scheme; some regularization is necessary.

Load-bearing premise

The porosity profile $\mu$ must be bounded away from zero and infinity and have finite Fisher information $\int_X |\nabla\sqrt{\mu}|^2\,dx < \infty$; this regularity is what gives the horizontal $H^1$ control on the density and $L^2$ control on the pressure gradient, so the proof does not cover discontinuous porosities such as the piecewise-constant bottleneck tested numerically.

Editorial extensions

If this is right

  • Every weak solution of (1)–(3) converges as $t \to \infty$ to the unique entropic optimal transport plan between $\mu$ and $\nu$, so the entropic relaxation of the AHT scheme is unconditionally convergent.
  • The continuum model is the discrete-to-continuous limit of finite multiphase porous-media flows: trajectories for $\nu$ made of finitely many Dirac masses converge to those for a general probability measure $\nu$.
  • The minimizing-movement scheme gives a constructive existence proof and a numerical algorithm; the experiments indicate linear convergence of the relative energy gap, opening a route to novel entropic optimal transport solvers on geometric domains.
  • For $\kappa=0$, suboptimal stationary configurations exist, so merely relaxing AHT to the space of plans without regularization does not guarantee convergence.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The numerically observed linear decay of $E(\rho_t)-E(\rho^*)$ suggests a quantitative dissipation inequality; proving one would turn the asymptotic guarantee into a complexity bound for entropy-regularized optimal transport computed by minimizing movements.
  • The bottleneck experiment with a discontinuous porosity $\mu$ indicates that the Fisher-information hypothesis may be relaxable to piecewise-smooth densities; a proof would extend the results to layered or fractured porous media.
  • Because the invariant measure of the flow is the entropic optimal transport plan, the dynamics can be read as an annealing of a Schrödinger bridge; this may connect the convergence proof to quantitative rates for iterative proportional fitting methods such as Sinkhorn's algorithm.
  • The same fibered-gradient-flow mechanism may apply to multi-marginal or dynamic transport problems where AHT-style projection algorithms stall, as long as a fiber metric with strong vertical-average compactness can be defined.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies the stratified drift-diffusion system (1)-(3), in which a continuum of phases y ∈ Y evolves under a joint pointwise volume constraint determined by a fixed first marginal μ. The dynamics are interpreted as a gradient flow in the fibered Wasserstein metric for the entropic optimal transport energy (4). The main results are: existence of global weak solutions via the JKO/minimizing-movement scheme (Section 4, Lemma 10); stability of weak solutions under joint variation of the marginals, including a discrete-to-continuous limit (Section 5, Lemma 13); and long-time convergence of every weak solution to the unique minimizer of the energy, i.e. the entropic optimal transport plan between μ and ν (Section 6, Lemma 16). The paper also shows that in the unregularized limit κ = 0, stationary suboptimal configurations persist, so regularization is essential for unconditional convergence, and it presents numerical experiments illustrating convergence, stability, and a bottleneck test with discontinuous μ.

Significance. If the main claims hold, this is a valuable contribution: it gives the first unconditional convergence guarantee for an entropic relaxation of the Angenent-Haker-Tannenbaum scheme and provides a rigorous variational framework for multiphase porous media flows with a continuum of phases. The paper is largely self-contained, the chain of estimates is detailed, and the energy functional is not fitted to the conclusions: the minimizer structure is cited from an independent source and the entropy parameter κ, time step τ, and discretization parameters are natural model parameters. The numerical section is qualitative and uses small academic examples, but it does support the theoretical statements and even probes regimes outside the stated assumptions. The main weakness is that Lemma 13, which is load-bearing for the stability and asymptotic-convergence results, contains an asserted but not proved compactness step for the varying-marginal case; this needs to be repaired before the central claims can be regarded as fully established.

major comments (3)
  1. [Section 5, Lemma 13] The proof of the strong convergence of the vertical averages in (95) and (96) is not supplied: after stating that the argument is 'similar' to Lemma 10, the text says that the variation of the marginals μn and νn 'can be factored into the proof without major issues', but no details are given. This is a load-bearing point: Lemma 16 uses Lemma 13 to pass to the limit in the cluster-point argument for long-time convergence, and Lemma 11 and Lemma 14 rely on the product convergence (96) for lower semicontinuity of the dissipation. In particular, the averages ωn are taken against νn, which is only known to converge weak*, so the Aubin-Lions compactness used for fixed ν in Lemma 10 does not transfer automatically. Please provide a complete proof of (95) and (96), or state and prove the needed variant of Aubin-Lions for varying reference measures.
  2. [Lemma 12, Eq. (83)] The statement (83) reads W2(ρt1,ρt2) ≤ W2_F(ρt1,ρt2) ≤ E(ρ0)(t2-t1), but the proof at the end of Lemma 12 derives W_F^2 ≤ (t2-t1) E(ρ0), so the correct inequality is W_F ≤ sqrt(E(ρ0)(t2-t1)). The printed bound is false in general and also has the wrong scaling in both E(ρ0) and time. Please correct the statement and check all places that use Lemma 12 for equicontinuity; the corrected square-root bound is sufficient for that purpose.
  3. [Section 1.5 and Section 8, Eq. (116)] The standing assumption (15) that μ is bounded away from 0 and ∞ and has finite Fisher information excludes piecewise-constant porosities, which are a standard case in porous media modeling. The numerical bottleneck experiment in Section 8 uses exactly such a discontinuous μ and says that the observed convergence 'may indicate that the Fisher information bound hypothesis could be further relaxed', but no proof is offered. This is a limitation of scope rather than an internal inconsistency, but it should be stated explicitly: the numerical experiment is outside the hypotheses of the theorems, and if a relaxation is intended, it should be proved or explicitly labeled as a conjecture.
minor comments (4)
  1. [Lemma 10, Eq. (72)] The sentence after (72) says 'weak* convergence of ωτ dx dt to ωτ dx dt'; the second factor should be ω dx dt. Also verify the displayed equation (72) itself, which appears to have a typo in the final equality.
  2. [Lemma 11] The proof refers to 'the regularity properties outlined in Lemma 12' and to Lemma 13 before those lemmas are stated. Please reorder the presentation or add explicit forward references so the logical dependence is clear.
  3. [Throughout] There are several typographical issues, e.g. 'Radon-Nykodym' for 'Radon-Nikodym', 'defing' for 'defining', and 'sufficiently' for 'sufficiently'. These do not affect the mathematics but should be corrected.
  4. [Section 8] The flipped-initialization experiment with N → ∞ is described as having initial entropy tending to infinity, so it is outside the assumptions of the existence/stability theorems; this should be stated clearly in the text, not only implicitly.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular derivation: the convergence theorem is built from self-contained PDE estimates; the one overlapping-author citation (Lemma 4) is an independent, standard entropic-OT characterization, and Lemma 13's asserted compactness step is a proof gap, not a circular input.

full rationale

The paper's derivation chain is predominantly self-contained. Existence (Lemma 10) is proven from the minimizing-movement scheme via JKO estimates, Aubin-Lions compactness, and the horizontal Fisher-information bound (49); no fitted quantity is renamed as a prediction. The long-time asymptotics (Lemma 16) follow a standard LaSalle-type argument: cluster points produce a limit solution (Lemma 13), summability of dissipation forces I=0 on the limit (Lemma 14), and Lemma 15 identifies I=0 with the unique minimizer. The only result imported from prior work is Lemma 4 on the structure of minimizers of E, cited to [6]; one author (Schmitzer) overlaps with the present paper, but this is an independently published, parameter-free characterization of entropic OT minimizers with stated assumptions that do not include the convergence theorem, so it is real evidence and not a load-bearing unverified self-citation. The stability lemma (Lemma 13) contains an asserted compactness step for vertical averages when marginals vary ('this can be factored into the proof without major issues') that is not written out; this is an omitted proof or potential gap, but it is not circular, because the convergence is not assumed as an input anywhere and the surrounding estimates are proved. Likewise the apparent dimension error in Lemma 12 (WF ≤ E0(t2−t1) instead of the square-root bound) is a typo-level correctness issue, not a circular step. The Fisher-information assumption (15) is a stated restriction, and the bottleneck experiment explicitly lies outside it; this limits generality but does not smuggle in the conclusion. No self-definitional, fitted-input, or renaming circularity is present.

Assumptions & free parameters 5 free parameters · 9 assumptions · 0 invented entities

The analytical results depend only on standard mathematical objects; no new particles, fields, or forces are introduced. The numerical section has hand-chosen parameters, but none are fitted to data or used in the proofs. The central theorem is therefore not built on hidden fitted degrees of freedom.

free parameters (5)
  • κ (entropic regularization strength) = 0.01 in numerics; arbitrary κ>0 in theory
    Barely a free parameter in the analysis: fixed positive constant. Chosen by hand in Section 8.
  • τ (JKO time step) = 0.25 in numerics; tends to 0 in theory
    Discretization step of the minimizing movement scheme; not fitted, but affects numerical convergence artifacts.
  • M (discretization of X) = 64 or 256
    Grid size for numerical experiments, chosen in Section 8.
  • N (number of phases in numerics) = 4, 16, 64, 256
    Finite approximation of ν; central to the discrete-to-continuous stability illustrations.
  • δ, a, b in bottleneck porosity = δ=1/8, b=a/2, a normalizes ∫µ=1
    Hand-chosen to make a discontinuous µ example in Figures 7-8.
assumptions (9)
  • domain assumption X ⊂ R^d compact with connected interior and smooth boundary; Y compact
    Used throughout for compactness of probability measures and boundary conditions.
  • domain assumption µ absolutely continuous with density bounded away from 0 and ∞ and finite Fisher information (equation 15)
    Needed for horizontal H¹ control of r and L² control of ∇Π.
  • domain assumption V ∈ C²(X×Y) non-negative
    Gives smoothness for a priori estimates and dual potentials.
  • domain assumption ν arbitrary probability measure on Y
    The continuum of phases; central object of the paper.
  • domain assumption Initial coupling ρ0 has finite energy E(ρ0)
    Required for JKO energy estimates and positivity of iterates.
  • standard math Fibered Wasserstein metric WF is a metric with the representations in Lemma 1
    Cited to [32]; foundation of the minimizing movement scheme.
  • standard math Brenier's theorem: uniqueness and measurability of optimal transport maps for absolutely continuous marginals
    Used in Lemmas 1, 6, 8 to build optimal plans and dual potentials.
  • standard math Weighted Poincaré inequality and Aubin-Lions compactness criterion (Rossi-Savaré)
    Used in Lemma 10 for weak compactness of pressure and strong compactness of vertical averages.
  • domain assumption Twist condition: x ↦ ∇_X V(x,y) injective, for κ=0 stationarity example
    Assumed in Section 7; not needed for main κ>0 results.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Continuum of coupled Wasserstein gradient flows." pith.science (2026). https://pith.science/paper/AMWTAOPC

@misc{pith2026241113969,
  author       = {Pith},
  title        = {Pith review of: Continuum of coupled Wasserstein gradient flows},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AMWTAOPC}},
  note         = {Machine review of arXiv:2411.13969}
}
read the original abstract

We study a system of drift-diffusion PDEs for a potentially infinite number of incompressible phases, subject to a joint pointwise volume constraint. Our analysis is based on the interpretation as a collection of coupled Wasserstein gradient flows or, equivalently, as a gradient flow in the space of couplings under a `fibered' Wasserstein distance. We prove existence of weak solutions, long-time asymptotics, and stability with respect to the mass distribution of the phases, including the discrete to continuous limit. A key step is to establish convergence of the product of pressure gradient and density, jointly over the infinite number of phases. The underlying energy functional is the objective of entropy regularized optimal transport, which allows us to interpret the model as the relaxation of the classical Angenent-Haker-Tannenbaum (AHT) scheme to the entropic setting. However, in contrast to the AHT scheme's lack of convergence guarantees, the relaxed scheme is unconditionally convergent. We conclude with numerical illustrations of the main results.

Figures

Figures reproduced from arXiv: 2411.13969 by the authors.

Figure 1
Figure 1. shows an example of the first iterations in the minimizing scheme for a product initial￾ization and N = 4. Note how the accumulation of mass of species 1 at the left boundary and of species 4 at the right boundary makes up for the gap left by the other species. 0.0 0.5 1.0 1 2 3 4 4 species t = 0.0 0.0 0.5 1.0 t = 0.25 0.0 0.5 1.0 t = 0.5 0.0 0.5 1.0 t = 0.75 0.0 0.5 1.0 t = 1.0 [PITH_FULL_IMAGE:figures/full_fig_p0… view at source ↗
Figure 2
Figure 2. Comparison of discrete trajectories for different number of species N and product initialization. Note how the trajectories tend to a consistent limit as N increases. 0.0 0.5 1.0 1 2 3 4 4 species t = 0.0 0.0 0.5 1.0 t = 0.5 0.0 0.5 1.0 t = 1.0 0.0 0.5 1.0 t = 1.5 0.0 0.5 1.0 t = 2.0 [PITH_FULL_IMAGE:figures/full_fig_p029_2.png] view at source ↗
Figure 3
Figure 3. Discrete trajectory ¯ρτ for N = 4 and flipped initialization [PITH_FULL_IMAGE:figures/full_fig_p029_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Comparison of discrete trajectories for different number of species N and flipped initialization. Note how the trajectories again tend to a consistent limit, even when the initial conditions have unbounded entropy as N → ∞ in [PITH_FULL_IMAGE:figures/full_fig_p030_4.png]
Figure 5
Figure 5. Figure 5: Comparison of the final discrete iterate (corresponding to t = 10) and the minimizer of E. 0 2 4 6 8 10 t 10−2 10−1 100 101 ∆E M = 64, product init M = 64, flipped init M = 256, product init M = 256, flipped init [PITH_FULL_IMAGE:figures/full_fig_p031_5.png]
Figure 6
Figure 6. Figure 6: Evolution of the relative difference between E(ρt) and E(ρ∗). Finally, Figures 7 and 8 show a trajectory for a µ with discontinuous density, hence not satisfying our Fisher information bound assumption. Concretely, we use the density µ(x) := ( a for x ∈ [0, 1] \ [0.5 −…
Figure 7
Figure 7. Figure 7: Discrete trajectory ¯ρτ for µ featuring a bottleneck, N = 4 and prod￾uct initialization. We change the visualization with respect to previous figures to showcase better the effect of the bottleneck on the density of the different phases. Here the width the different ba…
Figure 8
Figure 8. Figure 8: Comparison of discrete trajectories for µ featuring a bottleneck, prod￾uct initialization and different number of species N. Note how the trajectories again tend to a consistent limit, even when µ has (in the continuum limit) un￾bounded Fisher information. Proof. This …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    Lectures in Mathematics ETH Z¨ urich

    Luigi Ambrosio, Nicola Gigli, and Giuseppe Savar´ e.Gradient flows in metric spaces and in the space of proba- bility measures. Lectures in Mathematics ETH Z¨ urich. Birkh¨ auser Verlag, Basel, second edition, 2008

  2. [2]

    Minimizing flows for the Monge-Kantorovich problem

    Sigurd Angenent, Steven Haker, and Allen Tannenbaum. Minimizing flows for the Monge-Kantorovich problem. SIAM J. Math. Anal. , 35(1):61–97, 2003

  3. [3]

    Aziz and A

    K. Aziz and A. Settari. Petroleum Reservoir Simulation . Elsevier Applied Science Publishers, Londres, 1979

  4. [4]

    Bear and Y

    J. Bear and Y. Bachmat. Introduction to modeling of transport phenomena in porous media . Kluwer Academic Publishers, Dordrecht, The Netherlands, 1990

  5. [5]

    Benamou, G

    J.-D. Benamou, G. Carlier, and M. Laborde. An augmented Lagrangian approach to Wasserstein gradient flows and applications. In Gradient flows: from theory to application , volume 54 of ESAIM Proc. Surveys, pages 1–17. EDP Sci., Les Ulis, 2016

  6. [6]

    Domain decomposition for entropy regularized optimal transport

    Mauro Bonafini and Bernhard Schmitzer. Domain decomposition for entropy regularized optimal transport. Numerische Mathematik , 149:819–870, 2021

  7. [7]

    Y. Brenier. Polar factorization and monotone rearrangement of vector-valued functions. Comm. Pure Appl. Math., 44(4):375–417, 1991

  8. [8]

    Optimal transport, convection, magnetic relaxation and generalized Boussinesq equations.Journal of Nonlinear Science , 19(5):547–570, 10 2009

    Yann Brenier. Optimal transport, convection, magnetic relaxation and generalized Boussinesq equations.Journal of Nonlinear Science , 19(5):547–570, 10 2009. CONTINUUM OF COUPLED W ASSERSTEIN GRADIENT FLOWS 35

Show all 39 references
  1. [9]

    Canc` es, T

    C. Canc` es, T. O. Gallou¨ et, M. Laborde, and L. Monsaingeon. Simulation of multiphase porous media flows with minimizing movement and finite volume schemes. European J. Appl. Math. , 30(6):1123–1152, 2019

  2. [10]

    Canc` es, T

    C. Canc` es, T. O. Gallou¨ et, and L. Monsaingeon. The gradient flow structure of immiscible incompressible two-phase flows in porous media. C. R. Acad. Sci. Paris S´ er. I Math., 353:985–989, 2015

  3. [11]

    Canc` es, T

    C. Canc` es, T. O. Gallou¨ et, and L. Monsaingeon. Incompressible immiscible multiphase flows in porous media: a variational approach. Anal. PDE, 10(8):1845–1876, 2017

  4. [12]

    Canc` es, D

    C. Canc` es, D. Matthes, and F. Nabet. A two-phase two-fluxes degenerate Cahn-Hilliard model as constrained Wasserstein gradient flow. Arch. Ration. Mech. Anal. , 233(2):837–866, 2019

  5. [13]

    Canc` es and A

    C. Canc` es and A. Zurek. A convergent finite volume scheme for dissipation driven models with volume filling constraint. Numer. Math. , 151:279–328, 2022

  6. [14]

    Canc` es and D

    C. Canc` es and D. Matthes. Construction of a two-phase flow with singular energy by gradient flow methods. J. Diff. Equ. , 344:439–470, 2023

  7. [15]

    Convergence of entropic schemes for optimal transport and gradient flows

    Guillaume Carlier, Vincent Duval, Gabriel Peyr´ e, and Bernhard Schmitzer. Convergence of entropic schemes for optimal transport and gradient flows. SIAM Journal on Mathematical Analysis , 49(2):1385–1418, 2017

  8. [16]

    A first-order primal-dual algorithm for convex problems with applications to imaging

    Antonin Chambolle and Thomas Pock. A first-order primal-dual algorithm for convex problems with applications to imaging. Journal of Mathematical Imaging and Vision , 40(1):120–145, May 2011

  9. [17]

    Y. Chen, T. T. Georgiou, and M. Pavon. On the relation between optimal transport and Schr¨ odinger bridges: a stochastic control viewpoint. Journal of Optimization Theory and Applications , 169:671–691, 2016

  10. [18]

    Z. Chen. Degenerate two-phase incompressible flow. I. Existence, uniqueness and regularity of a weak solution. J. Diff. Equ. , 171(2):203–232, 2001

  11. [19]

    Chen and R

    Z. Chen and R. E. Ewing. Fully discrete finite element analysis of multiphase flow in groundwater hydrology. SIAM J. Numer. Anal. , 34(6):2228–2253, 1997

  12. [20]

    Cordoba, D

    D. Cordoba, D. Faraco, and F. Gancedo. Lack of uniqueness for weak solutions of the incompressible porous media equation. Arch. Ration. Mech. An. , 200:725–746, 2011

  13. [21]

    M. Cuturi. Sinkhorn distances: Lightspeed computation of optimal transportation distances. In Advances in Neural Information Processing Systems 26 (NIPS 2013) , pages 2292–2300, 2013

  14. [22]

    D. E. Edmunds and B. Opic. Weighted Poincar´ e and Friedrichs inequalities.Journal of the London Mathematical Society, s2-47(1):79–96, 02 1993

  15. [23]

    Eymard, R

    R. Eymard, R. Herbin, and A. Michel. Mathematical study of a petroleum-engineering scheme. M2AN Math. Model. Numer. Anal. , 37(6):937–972, 2003

  16. [24]

    F¨ orster and L

    C. F¨ orster and L. Sz´ ekelyhidi. Piecewise constant subsolutions for the Muskat problem.Commun. Math. Phys. , 363:1051–1080, 2018

  17. [25]

    Wilfrid Gangbo and Robert J. McCann. The geometry of optimal transportation. Acta Math., 177(2):113–161, 1996

  18. [26]

    Gigli and F

    N. Gigli and F. Otto. Entropic Burgers’ equation via a minimizing movement scheme based on the Wasserstein metric. Calc. Var. Partial Differential Equations , 47(1-2):181–206, 2013

  19. [27]

    Multimarginal optimal transport with a tree- structured cost and the Schr¨ odinger bridge problem.SIAM Journal on Control and Optimization , 59(4):2428– 2453, 2021

    Isabel Haasler, Axel Ringh, Yongxin Chen, and Johan Karlsson. Multimarginal optimal transport with a tree- structured cost and the Schr¨ odinger bridge problem.SIAM Journal on Control and Optimization , 59(4):2428– 2453, 2021

  20. [28]

    Self-similarity in a thin film Muskat problem.SIAM Journal on Mathematical Analysis, 49(4):2790–2842, 2017

    Philippe Lauren¸ cot and Bogdan-Vasile Matioc. Self-similarity in a thin film Muskat problem.SIAM Journal on Mathematical Analysis, 49(4):2790–2842, 2017

  21. [29]

    Towards a mathematical theory of trajectory inference

    Hugo Lavenant, Stephen Zhang, Young-Heon Kim, and Geoffrey Schiebinger. Towards a mathematical theory of trajectory inference. Ann. Appl. Probab., 34(1A), 2024

  22. [30]

    Otto and W

    F. Otto and W. E. Thermodynamically driven incompressible fluid mixtures. J. Chem. Phys. , 107(23):10177– 10184, 1997

  23. [31]

    Evolution of microstructure in unstable porous media flow: a relaxational approach

    Felix Otto. Evolution of microstructure in unstable porous media flow: a relaxational approach. Communica- tions on Pure and Applied Mathematics: A Journal Issued by the Courant Institute of Mathematical Sciences , 52(7):873–915, 1999

  24. [32]

    Heterogeneous gradient flows in the topology of fibered optimal transport

    Jan Peszek and David Poyato. Heterogeneous gradient flows in the topology of fibered optimal transport. Cal- culus of Variations and Partial Differential Equations , 62(9):258, Oct 2023

  25. [33]

    Computational optimal transport.Foundations and Trends in Machine Learn- ing, 11(5–6):355–607, 2019

    Gabriel Peyr´ e and Marco Cuturi. Computational optimal transport.Foundations and Trends in Machine Learn- ing, 11(5–6):355–607, 2019

  26. [34]

    Tightness, integral equicontinuity and compactness for evolution problems in Banach spaces

    Riccarda Rossi and Giuseppe Savar´ e. Tightness, integral equicontinuity and compactness for evolution problems in Banach spaces. Ann. Sc. Norm. Super. Pisa Cl. Sci. (5) , 2(2):395–431, 2003. 36 CL ´EMENT CANC `ES, DANIEL MATTHES, ISMAEL MEDINA, AND BERNHARD SCHMITZER

  27. [35]

    Optimal Transport for Applied Mathematicians , volume 87 of Progress in Nonlinear Differential Equations and Their Applications

    Filippo Santambrogio. Optimal Transport for Applied Mathematicians , volume 87 of Progress in Nonlinear Differential Equations and Their Applications . Birkh¨ auser Boston, 2015

  28. [36]

    Sz´ ekelyhidi Jr

    L. Sz´ ekelyhidi Jr. Relaxation of the incompressible porous media equation. Ann. Sci. ´Ecole Norm. Sup. , 45(3):491–509, 2012

  29. [37]

    J. L. V´ azquez.The porous medium equation . Oxford Mathematical Monographs. The Clarendon Press Oxford University Press, Oxford, 2007. Mathematical theory

  30. [38]

    M. H. Vignal. Convergence of a finite volume scheme for an elliptic-hyperbolic system. RAIRO Mod´ el. Math. Anal. Num´ er., 30(7):841–872, 1996

  31. [39]

    C. Villani. Topics in Optimal Transportation, volume 58 of Graduate Studies in Mathematics . American Math- ematical Society, Providence, 2003. Univ. Lille, CNRS, Inria, UMR8524 – Laboratoire Paul Painlev ´e, F-59000 Lille Email address : clement.cances@inria.fr Department of ...

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.