pith. sign in

arxiv: 2603.20153 · v2 · submitted 2026-03-20 · 🧮 math.AP

Global solutions to cross-diffusion systems with independent advections in one dimension

Pith reviewed 2026-05-15 07:58 UTC · model grok-4.3

classification 🧮 math.AP
keywords cross-diffusionadvectionglobal existenceone dimensionvanishing viscositycompensated compactnessYoung measure
0
0 comments X

The pith

Global solutions exist for one-dimensional cross-diffusion systems with independent advections for all pressure exponents and arbitrary data.

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

The paper constructs global-in-time solutions to a system of two cross-diffusing species in one space dimension, where each species moves according to Darcy's law under a power-law pressure that depends on total density plus independent advection fields. Earlier methods required extra structural conditions on the advections or the initial data, especially when the pressure exponent exceeds one. The new construction proceeds by passing to the limit in a vanishing-viscosity regularization. The decisive step is showing that oscillations in one species density and in the pressure gradient remain correlated inside the Young measure, which lets the product term pass to the limit using only three entropy-entropy flux pairs.

Core claim

Solutions are recovered as the vanishing-viscosity limit of regularized equations. The central technical point is that any oscillations present in the approximating sequence satisfy a correlation between the first species density and the pressure gradient; this correlation identifies the limit of their product inside the compensated-compactness framework. The argument closes with only three entropy pairs and applies uniformly to every pressure exponent greater than zero and to arbitrary initial data.

What carries the argument

Vanishing-viscosity approximation together with a compensated-compactness argument that exploits the correlation of oscillations between species density and pressure gradient inside the Young measure, implemented with three entropy-entropy flux pairs.

If this is right

  • Global solutions exist for any initial data in one space dimension.
  • The result holds for every pressure exponent alpha greater than zero.
  • No structural assumptions on the advection fields are required.
  • The construction works uniformly in both the sublinear and superlinear pressure regimes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The oscillation-correlation property may be reusable in other one-dimensional systems whose symmetry is broken by transport terms.
  • Fewer entropy pairs suffice when the advection structure forces a specific dependence between density and pressure gradient oscillations.
  • Direct numerical checks of the product limit in oscillating regimes could test whether the correlation persists beyond the analytic setting.

Load-bearing premise

Oscillations of the species density and the pressure gradient remain correlated in the Young measure so that their product limit can be identified without knowing the full measure.

What would settle it

A concrete sequence of vanishing-viscosity solutions in which the weak limit of u times the pressure gradient differs from the product of the individual weak limits would show that the correlation step fails and block passage to the limit.

read the original abstract

We consider cross-diffusion systems describing evolution of two species $u$ and $v$ moving according to Darcy's law with the pressure law $p(s) = \frac{1}{\alpha-1} s^{\alpha-1}$ where $s=u+v$. One of the most challenging questions in the field is the construction of solutions to the problem in the presence of additional advection fields, without imposing any artificial structure on the fields or the initial conditions. Although advection arises naturally in these models, it breaks the symmetry of the system and prevents application of techniques developed in recent years. Here, we provide a new approach to construct solutions in one space dimension that works in a unified manner for all pressure exponents $\alpha \in (0,\infty)$ and for arbitrary initial data. In~particular, in the regime $\alpha > 1$, this yields the first existence result of its kind, obtained without any structural assumptions. We construct the solutions as a limit of a vanishing viscosity approximation $(u_{\varepsilon}, v_{\varepsilon})$. The main challenge is to identify the limit of $u_{\varepsilon} \, \partial_x p(s_{\varepsilon})$, where $s_{\varepsilon} = u_{\varepsilon} + v_{\varepsilon}$. The key new insight is that possible oscillations of $u_\varepsilon$ and $\partial_x p(s_\varepsilon)$ are correlated, simplifying the Young measure analysis in the compensated compactness argument and allowing identification of the limit. Somewhat surprisingly, in contrast to the theory of $2\times2$ hyperbolic systems, the argument relies on only three entropy-entropy flux pairs. This is particularly useful for $\alpha>2$, where it is unclear whether additional entropies are available.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proves global existence of weak solutions to a one-dimensional cross-diffusion system for two species with independent advection fields and pressure p(s) = 1/(α-1)s^{α-1} for all α > 0 and arbitrary nonnegative initial data. The proof proceeds by constructing solutions as limits of a vanishing-viscosity regularization and applying compensated compactness to pass to the limit in the nonlinear advection term u ∂_x p(s), exploiting a correlation between oscillations of u and ∂_x p(s) that allows the use of only three entropy-entropy flux pairs.

Significance. If the key oscillation-correlation argument holds, the result is significant because it provides the first existence theorem for α > 1 without any structural assumptions on the advection fields or initial data, unifying the treatment across all pressure exponents. The reduction to three entropy pairs is a technical simplification that may have broader applicability in 1D cross-diffusion models.

major comments (2)
  1. [Compensated compactness / limit identification of u_ε ∂_x p(s_ε)] In the compensated-compactness argument for identifying lim u_ε ∂_x p(s_ε) (the main technical step after the vanishing-viscosity construction), the claim that oscillations of u_ε and ∂_x p(s_ε) are correlated because of the specific pressure law and 1D advection structure must be verified explicitly. The sketch indicates that only three entropy-entropy flux pairs suffice, but the explicit computation showing why the Young measure permits passage to the limit (without additional structural hypotheses) is not visible in sufficient detail; this identification is load-bearing for the global-solution claim, especially when α > 1.
  2. [Vanishing-viscosity approximation and a-priori estimates] The uniform a-priori estimates and compactness for the approximating system (u_ε, v_ε) are stated to hold for arbitrary initial data, but the passage to the limit in the cross terms requires confirming that the oscillation correlation remains valid uniformly in ε for all advection fields; if the correlation relies on hidden regularity or sign conditions on the advection velocities, the result would not cover the stated generality.
minor comments (2)
  1. [Entropy pairs] Clarify the precise definition of the entropy-entropy flux pairs used in the three-pair argument and state explicitly which pairs are employed for the different ranges of α.
  2. [Introduction] In the introduction, add a short comparison table or explicit citations distinguishing the new result from prior works that required structural assumptions on the advection for α > 1.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the two major points below and will revise the manuscript to provide the requested explicit computations and clarifications while preserving the stated generality.

read point-by-point responses
  1. Referee: [Compensated compactness / limit identification of u_ε ∂_x p(s_ε)] In the compensated-compactness argument for identifying lim u_ε ∂_x p(s_ε) (the main technical step after the vanishing-viscosity construction), the claim that oscillations of u_ε and ∂_x p(s_ε) are correlated because of the specific pressure law and 1D advection structure must be verified explicitly. The sketch indicates that only three entropy-entropy flux pairs suffice, but the explicit computation showing why the Young measure permits passage to the limit (without additional structural hypotheses) is not visible in sufficient detail; this identification is load-bearing for the global-solution claim, especially when α > 1.

    Authors: We agree that an explicit verification strengthens the argument. The correlation arises because the 1D continuity equations for u_ε and s_ε together with the specific power-law pressure imply that the Young measure ν satisfies ∫ (λ - ū) (p'(s) - p'(s̄)) dν = 0 for the relevant test functions, which is verified by direct computation on the three entropy pairs: the basic L^1 entropy, the α-dependent power entropy (valid for α > 1), and the third pair obtained by multiplying the advection equation by a suitable function of s. This computation shows that the measure is supported on the graph u = ū(s), allowing passage to the limit in u ∂_x p(s) without further hypotheses. We will insert the full calculation in the revised Section 4, confirming uniformity in ε from the ε-independent a priori bounds. revision: yes

  2. Referee: [Vanishing-viscosity approximation and a-priori estimates] The uniform a-priori estimates and compactness for the approximating system (u_ε, v_ε) are stated to hold for arbitrary initial data, but the passage to the limit in the cross terms requires confirming that the oscillation correlation remains valid uniformly in ε for all advection fields; if the correlation relies on hidden regularity or sign conditions on the advection velocities, the result would not cover the stated generality.

    Authors: The a priori estimates (L^∞ bounds on s_ε, BV bounds on the fluxes, and entropy dissipation) are obtained from the vanishing-viscosity system and depend only on the L^1 norm of the initial data and the L^∞ norm of the advection fields; no sign conditions or extra regularity on the velocities are used. The oscillation correlation is a purely algebraic consequence of the 1D structure and the pressure law once these bounds are available, and therefore holds uniformly in ε for arbitrary bounded advection fields. We will add a clarifying paragraph after the statement of the a priori estimates to make this independence explicit. revision: yes

Circularity Check

0 steps flagged

Direct limit argument with independent compensated compactness analysis

full rationale

The paper constructs global solutions as limits of a vanishing-viscosity approximation (u_ε, v_ε) and identifies lim u_ε ∂x p(s_ε) via compensated compactness in one dimension. The central step asserts that oscillations of u_ε and ∂x p(s_ε) are correlated due to the explicit pressure p(s) = 1/(α-1)s^{α-1} and the 1D advection structure, allowing the limit to be identified with only three entropy-entropy flux pairs. This correlation is presented as a new verifiable mathematical property of the system rather than a self-definition, fitted parameter, or reduction to prior self-citations. The argument is self-contained against external benchmarks (Young measure theory and compensated compactness) with no load-bearing self-referential definitions or renamings. No steps reduce by construction to the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The argument relies on standard existence of Young measures for bounded sequences in L^∞ and on the theory of compensated compactness for entropy pairs; no new free parameters or invented entities are introduced.

axioms (2)
  • standard math Existence of Young measures for bounded sequences in L^∞(Ω)
    Invoked to represent weak limits of oscillating sequences u_ε and ∂x p(s_ε)
  • standard math Compensated compactness lemma for entropy-entropy flux pairs
    Used to pass to the limit in the nonlinear term after identifying the correlation

pith-pipeline@v0.9.0 · 5606 in / 1276 out tokens · 41528 ms · 2026-05-15T07:58:10.760092+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

60 extracted references · 60 canonical work pages

  1. [1]

    N. D. Alikakos.L p bounds of solutions of reaction-diffusion equations.Comm. Partial Differential Equations, 4(8):827–868, 1979

  2. [2]

    Ambrosio, N

    L. Ambrosio, N. Fusco, and D. Pallara.Functions of bounded variation and free discontinuity problems. Oxford Mathematical Monographs. The Clarendon Press, Oxford University Press, New York, 2000

  3. [3]

    N.J.Armstrong, K.J.Painter, andJ.A.Sherratt.Acontinuumapproachtomodellingcell-celladhesion. J. Theoret. Biol., 243(1):98–113, 2006

  4. [4]

    D. G. Aronson and P. Bénilan. Régularité des solutions de l’équation des milieux poreux dansRN.C. R. Acad. Sci. Paris Sér. A-B, 288(2):A103–A105, 1979

  5. [5]

    J. M. Ball. A version of the fundamental theorem for Young measures. InPDEs and continuum models of phase transitions (Nice, 1988), volume 344 ofLecture Notes in Phys., pages 207–215. Springer, Berlin, 1989

  6. [6]

    Bansaye, A

    V. Bansaye, A. Moussa, and F. Muñoz Hernández. Stability of a cross-diffusion system and approxi- mation by repulsive random walks: a duality approach.J. Eur. Math. Soc. (JEMS), 27(9):3889–3928, 2025

  7. [7]

    Bertsch, R

    M. Bertsch, R. Dal Passo, and M. Mimura. A free boundary problem arising in a simplified tumour growth model of contact inhibition.Interfaces Free Bound., 12(2):235–250, 2010. 56 JAKUB SKRZECZKOWSKI

  8. [8]

    Bertsch, M

    M. Bertsch, M. E. Gurtin, and D. Hilhorst. On interacting populations that disperse to avoid crowding: the case of equal dispersal velocities.Nonlinear Anal., 11(4):493–499, 1987

  9. [9]

    Bertsch, M

    M. Bertsch, M. E. Gurtin, D. Hilhorst, and L. A. Peletier. On interacting populations that disperse to avoid crowding: preservation of segregation.J. Math. Biol., 23(1):1–13, 1985

  10. [10]

    Brenier, C

    Y. Brenier, C. De Lellis, and L. Székelyhidi, Jr. Weak-strong uniqueness for measure-valued solutions. Comm. Math. Phys., 305(2):351–361, 2011

  11. [11]

    Bubba, B

    F. Bubba, B. Perthame, C. Pouchol, and M. Schmidtchen. Hele-Shaw limit for a system of two reaction- (cross-)diffusion equations for living tissues.Arch. Ration. Mech. Anal., 236(2):735–766, 2020

  12. [12]

    J. A. Carrillo, S. Fagioli, F. Santambrogio, and M. Schmidtchen. Splitting schemes and segregation in reaction cross-diffusion systems.SIAM J. Math. Anal., 50(5):5695–5718, 2018

  13. [13]

    J.A.Carrillo, T.Lorenzi, andF.R.Macfarlane.Spatialsegregationacrosstravellingfrontsinindividual- based and continuum models for the growth of heterogeneous cell populations.Bull. Math. Biol., 87(6):Paper No. 77, 40, 2025

  14. [14]

    J. A. Carrillo, H. Murakawa, M. Sato, H. Togashi, and O. Trush. A population dynamics model of cell- cell adhesion incorporating population pressure and density saturation.J. Theoret. Biol., 474:14–24, 2019

  15. [15]

    J. A. Carrillo, Y. Salmaniw, and J. Skrzeczkowski. Well-posedness of aggregation-diffusion systems with irregular kernels.arXiv preprint arXiv:2406.09227, to appear in Ann. Inst. H. Poincaré C Anal. Non Linéaire, 2024

  16. [16]

    J. A. Carrillo, J. Skrzeczkowski, and J. Warnett. The Stein-log-Sobolev inequality and the exponen- tial rate of convergence for the continuous Stein variational gradient descent method.arXiv preprint arXiv:2412.10295, 2024

  17. [17]

    J. A. Carrillo and J. L. Vázquez. Fine asymptotics for fast diffusion equations.Comm. Partial Differ- ential Equations, 28(5-6):1023–1056, 2003

  18. [18]

    Chaudhuri, P

    N. Chaudhuri, P. Gwiazda, and E. Zatorska. Analysis of the generalized Aw-Rascle model.Comm. Partial Differential Equations, 48(3):440–477, 2023

  19. [19]

    G.-Q. Chen, M. Slemrod, and D. Wang. Vanishing viscosity method for transonic flow.Arch. Ration. Mech. Anal., 189(1):159–188, 2008

  20. [20]

    X. Chen, E. S. Daus, and A. Jüngel. Global existence analysis of cross-diffusion population systems for multiple species.Arch. Ration. Mech. Anal., 227(2):715–747, 2018

  21. [21]

    Conti, G

    S. Conti, G. Dolzmann, and S. Müller. The div-curl lemma for sequences whose divergence and curl are compact inW −1,1.C. R. Math. Acad. Sci. Paris, 349(3-4):175–178, 2011

  22. [22]

    Darcy.Les fontaines publiques de la ville de Dijon, volume 2

    H. Darcy.Les fontaines publiques de la ville de Dijon, volume 2. 1856

  23. [23]

    David, T

    N. David, T. Dębiec, M. Mandal, and M. Schmidtchen. A degenerate cross-diffusion system as the inviscid limit of a nonlocal tissue growth model.SIAM J. Math. Anal., 56(2):2090–2114, 2024. 57

  24. [24]

    David and M

    N. David and M. Schmidtchen. On the incompressible limit for a tumour growth model incorporating convective effects.Comm. Pure Appl. Math., 77(5):2613–2650, 2024

  25. [25]

    R. J. DiPerna. Convergence of approximate solutions to conservation laws.Arch. Rational Mech. Anal., 82(1):27–70, 1983

  26. [26]

    Dębiec, B

    T. Dębiec, B. Perthame, M. Schmidtchen, and N. Vauchelet. Incompressible limit for a two-species model with coupling through Brinkman’s law in any dimension.J. Math. Pures Appl. (9), 145:204–239, 2021

  27. [27]

    Elbar and F

    C. Elbar and F. Santambrogio. A cross-diffusion system with independent drifts and fast diffusion.arXiv preprint arXiv:2510.07937, 2025

  28. [28]

    Elbar and J

    C. Elbar and J. Skrzeczkowski. On the inviscid limit connecting Brinkman’s and Darcy’s models of tissue growth with nonlinear pressure.J. Math. Fluid Mech., 27(2):Paper No. 28, 14, 2025

  29. [29]

    L. C. Evans.Weak convergence methods for nonlinear partial differential equations, volume 74 ofCBMS Regional Conference Series in Mathematics. Conference Board of the Mathematical Sciences, Washing- ton, DC; by the American Mathematical Society, Providence, RI, 1990

  30. [30]

    Falcó, R

    C. Falcó, R. E. Baker, and J. A. Carrillo. A nonlocal-to-local approach to aggregation-diffusion equa- tions.SIAM Rev., 67(2):353–372, 2025

  31. [31]

    Feireisl, P

    E. Feireisl, P. Gwiazda, A. Świerczewska-Gwiazda, and E. Wiedemann. Dissipative measure-valued solutions to the compressible Navier-Stokes system.Calc. Var. Partial Differential Equations, 55(6):Art. 141, 20, 2016

  32. [32]

    B. Gess. Optimal regularity for the porous medium equation.J. Eur. Math. Soc. (JEMS), 23(2):425–465, 2021

  33. [33]

    H. P. Greenspan. On the growth and stability of cell cultures and solid tumors.J. Theoret. Biol., 56(1):229–242, 1976

  34. [34]

    M. E. Gurtin and A. C. Pipkin. A note on interacting populations that disperse to avoid crowding. Quart. Appl. Math., 42(1):87–94, 1984

  35. [35]

    Gwiazda, B

    P. Gwiazda, B. Perthame, and A. Świerczewska Gwiazda. A two-species hyperbolic-parabolic model of tissue growth.Comm. Partial Differential Equations, 44(12):1605–1618, 2019

  36. [36]

    M. A. Herrero and M. Pierre. The Cauchy problem forut = ∆um when0< m <1.Trans. Amer. Math. Soc., 291(1):145–158, 1985

  37. [37]

    Kim and A

    I. Kim and A. R. Mészáros. On nonlinear cross-diffusion systems: an optimal transport approach.Calc. Var. Partial Differential Equations, 57(3):Paper No. 79, 40, 2018

  38. [38]

    Lions, B

    P.-L. Lions, B. Perthame, and P. E. Souganidis. Existence and stability of entropy solutions for the hyperbolic systems of isentropic gas dynamics in Eulerian and Lagrangian coordinates.Comm. Pure Appl. Math., 49(6):599–638, 1996. 58 JAKUB SKRZECZKOWSKI

  39. [39]

    Lions, B

    P.-L. Lions, B. Perthame, and E. Tadmor. Kinetic formulation of the isentropic gas dynamics and p-systems.Comm. Math. Phys., 163(2):415–431, 1994

  40. [40]

    Matthes and C

    D. Matthes and C. Parsch. Convergence to equilibrium for cross diffusion systems with nonlocal inter- action.Calc. Var. Partial Differential Equations, 65(3):Paper No. 74, 2026

  41. [41]

    A. R. Mészáros and G. Parker. On a cross-diffusion system with independent drifts and no self-diffusion: The existence of totally mixed solutions.arXiv preprint arXiv:2504.18484, 2025

  42. [42]

    A. R. Mészáros and G. Parker. Existence theory for a cross-diffusion system with independent drifts: Mixing dynamics.arXiv preprint arXiv:2603.18770, 2026

  43. [43]

    Mogilner and L

    A. Mogilner and L. Edelstein-Keshet. A non-local model for a swarm.J. Math. Biol., 38(6):534–570, 1999

  44. [44]

    Murakawa and H

    H. Murakawa and H. Togashi. Continuous models for cell–cell adhesion.Journal of Theoretical Biology, 374:1–12, 2015

  45. [45]

    F. Murat. Compacité par compensation.Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 5(3):489–507, 1978

  46. [46]

    Pedregal.Parametrized measures and variational principles, volume 30 ofProgress in Nonlinear Differential Equations and their Applications

    P. Pedregal.Parametrized measures and variational principles, volume 30 ofProgress in Nonlinear Differential Equations and their Applications. Birkhäuser Verlag, Basel, 1997

  47. [47]

    Perthame, F

    B. Perthame, F. Quirós, and J. L. Vázquez. The Hele-Shaw asymptotics for mechanical models of tumor growth.Arch. Ration. Mech. Anal., 212(1):93–127, 2014

  48. [48]

    Perthame and J

    B. Perthame and J. Skrzeczkowski. Fast reaction limit with nonmonotone reaction function.Comm. Pure Appl. Math., 76(7):1495–1527, 2023

  49. [49]

    P. I. Plotnikov. Passage to the limit with respect to viscosity in an equation with a variable direction of parabolicity.Differentsialnye Uravneniya, 30(4):665–674, 734, 1994

  50. [50]

    B. C. Price and X. Xu. Global existence theorem for a model governing the motion of two cell popula- tions.Kinet. Relat. Models, 13(6):1175–1191, 2020

  51. [51]

    Rindler.Calculus of variations

    F. Rindler.Calculus of variations. Universitext. Springer, Cham, 2018

  52. [52]

    Roubíček.Nonlinear partial differential equations with applications, volume 153 ofInternational Series of Numerical Mathematics

    T. Roubíček.Nonlinear partial differential equations with applications, volume 153 ofInternational Series of Numerical Mathematics. Birkhäuser/Springer Basel AG, Basel, second edition, 2013

  53. [53]

    Shigesada, K

    N. Shigesada, K. Kawasaki, and E. Teramoto. Spatial segregation of interacting species.J. Theoret. Biol., 79(1):83–99, 1979

  54. [54]

    Skrzeczkowski

    J. Skrzeczkowski. Fast reaction limit and forward-backward diffusion: a Radon-Nikodym approach.C. R. Math. Acad. Sci. Paris, 360:189–203, 2022

  55. [55]

    Tadmor and T

    E. Tadmor and T. Tao. Velocity averaging, kinetic formulations, and regularizing effects in quasi-linear PDEs.Comm. Pure Appl. Math., 60(10):1488–1521, 2007

  56. [56]

    L. Tartar. Compensated compactness and applications to partial differential equations. InNonlinear analysis and mechanics: Heriot-Watt Symposium, Vol. IV, volume 39 ofRes. Notes in Math., pages 136–212. Pitman, Boston, Mass.-London, 1979. 59

  57. [57]

    C. M. Topaz and A. L. Bertozzi. Swarming patterns in a two-dimensional kinematic model for biological groups.SIAM J. Appl. Math., 65(1):152–174, 2004

  58. [58]

    J. L. Vázquez. Behaviour of the velocity of one-dimensional flows in porous media.Trans. Amer. Math. Soc., 286(2):787–802, 1984

  59. [59]

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

  60. [60]

    L. C. Young. Generalized surfaces in the calculus of variations.Ann. of Math. (2), 43:84–103, 1942. Jakub Skrzeczkowski:St John’s College, University of Oxford, St Giles, Oxford, OX1 3JP, United Kingdom & Mathematical Institute, University of Oxford, Woodstock Road, Oxford, OX2 6GG, United Kingdom Email address:jakub.skrzeczkowski@maths.ox.ac.uk