Global solutions to cross-diffusion systems with independent advections in one dimension
Pith reviewed 2026-05-15 07:58 UTC · model grok-4.3
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.
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
- 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.
Referee Report
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)
- [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.
- [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)
- [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 α.
- [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
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
-
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
-
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
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
axioms (2)
- standard math Existence of Young measures for bounded sequences in L^∞(Ω)
- standard math Compensated compactness lemma for entropy-entropy flux pairs
Reference graph
Works this paper leans on
-
[1]
N. D. Alikakos.L p bounds of solutions of reaction-diffusion equations.Comm. Partial Differential Equations, 4(8):827–868, 1979
work page 1979
-
[2]
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
work page 2000
-
[3]
N.J.Armstrong, K.J.Painter, andJ.A.Sherratt.Acontinuumapproachtomodellingcell-celladhesion. J. Theoret. Biol., 243(1):98–113, 2006
work page 2006
-
[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
work page 1979
-
[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
work page 1988
-
[6]
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
work page 2025
-
[7]
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
work page 2010
-
[8]
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
work page 1987
-
[9]
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
work page 1985
-
[10]
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
work page 2011
- [11]
-
[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
work page 2018
-
[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
work page 2025
-
[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
work page 2019
- [15]
- [16]
-
[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
work page 2003
-
[18]
N. Chaudhuri, P. Gwiazda, and E. Zatorska. Analysis of the generalized Aw-Rascle model.Comm. Partial Differential Equations, 48(3):440–477, 2023
work page 2023
-
[19]
G.-Q. Chen, M. Slemrod, and D. Wang. Vanishing viscosity method for transonic flow.Arch. Ration. Mech. Anal., 189(1):159–188, 2008
work page 2008
-
[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
work page 2018
- [21]
-
[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]
-
[24]
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
work page 2024
-
[25]
R. J. DiPerna. Convergence of approximate solutions to conservation laws.Arch. Rational Mech. Anal., 82(1):27–70, 1983
work page 1983
- [26]
-
[27]
C. Elbar and F. Santambrogio. A cross-diffusion system with independent drifts and fast diffusion.arXiv preprint arXiv:2510.07937, 2025
-
[28]
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
work page 2025
-
[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
work page 1990
- [30]
-
[31]
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
work page 2016
-
[32]
B. Gess. Optimal regularity for the porous medium equation.J. Eur. Math. Soc. (JEMS), 23(2):425–465, 2021
work page 2021
-
[33]
H. P. Greenspan. On the growth and stability of cell cultures and solid tumors.J. Theoret. Biol., 56(1):229–242, 1976
work page 1976
-
[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
work page 1984
-
[35]
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
work page 2019
-
[36]
M. A. Herrero and M. Pierre. The Cauchy problem forut = ∆um when0< m <1.Trans. Amer. Math. Soc., 291(1):145–158, 1985
work page 1985
- [37]
- [38]
- [39]
-
[40]
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
work page 2026
- [41]
- [42]
-
[43]
A. Mogilner and L. Edelstein-Keshet. A non-local model for a swarm.J. Math. Biol., 38(6):534–570, 1999
work page 1999
-
[44]
H. Murakawa and H. Togashi. Continuous models for cell–cell adhesion.Journal of Theoretical Biology, 374:1–12, 2015
work page 2015
-
[45]
F. Murat. Compacité par compensation.Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 5(3):489–507, 1978
work page 1978
-
[46]
P. Pedregal.Parametrized measures and variational principles, volume 30 ofProgress in Nonlinear Differential Equations and their Applications. Birkhäuser Verlag, Basel, 1997
work page 1997
-
[47]
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
work page 2014
-
[48]
B. Perthame and J. Skrzeczkowski. Fast reaction limit with nonmonotone reaction function.Comm. Pure Appl. Math., 76(7):1495–1527, 2023
work page 2023
-
[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
work page 1994
-
[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
work page 2020
-
[51]
Rindler.Calculus of variations
F. Rindler.Calculus of variations. Universitext. Springer, Cham, 2018
work page 2018
-
[52]
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
work page 2013
-
[53]
N. Shigesada, K. Kawasaki, and E. Teramoto. Spatial segregation of interacting species.J. Theoret. Biol., 79(1):83–99, 1979
work page 1979
-
[54]
J. Skrzeczkowski. Fast reaction limit and forward-backward diffusion: a Radon-Nikodym approach.C. R. Math. Acad. Sci. Paris, 360:189–203, 2022
work page 2022
-
[55]
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
work page 2007
-
[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
work page 1979
-
[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
work page 2004
-
[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
work page 1984
-
[59]
J. L. Vázquez.The porous medium equation. Oxford Mathematical Monographs. The Clarendon Press, Oxford University Press, Oxford, 2007. Mathematical theory
work page 2007
-
[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
work page 1942
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.