pith. sign in

arxiv: 1907.03611 · v1 · pith:ORINLANVnew · submitted 2019-07-08 · 🧮 math.AP

The Cauchy problem for the standard one pressure system of two fluid flows with energy equations

M. Colombeau This is my paper

Pith reviewed 2026-05-25 01:03 UTC · model grok-4.3

classification 🧮 math.AP
keywords two-fluid flowsCauchy problemRadon measuresweak solutionsapproximate solutionsODE reductionenergy equationsshock tube
0
0 comments X

The pith

Approximate solutions to the two-fluid flow equations with energy reduce to ODEs and converge weakly to Radon measure solutions.

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

The paper establishes existence of solutions to the Cauchy problem for the standard one-pressure two-fluid model with energy equations by constructing a specific family of approximate solutions. These approximations are built so the original PDE system reduces to a closed system of ordinary differential equations while preserving bounds that allow passage to a limit via weak compactness. The resulting limit objects are Radon measures that satisfy the equations in the natural weak sense. A reader would care because the construction simultaneously supplies a rigorous existence proof and a practical numerical method, tested on the Toumi shock tube where it reproduces results from independent schemes.

Core claim

We construct with full rigorous mathematical proof a family of approximate solutions to the Cauchy problem for the standard system of two fluid flows with energy equations and we pass to the limit by weak compactness to obtain Radon measures that satisfy the equations in a natural weak sense. Our method provides a convergent numerical method for the numerical calculation of these Radon measures by reducing the system of partial differential equations in the case of these approximate solutions to a system of ordinary differential equations.

What carries the argument

The family of approximate solutions that reduce the PDE system to a closed system of ODEs while satisfying the a priori bounds needed for weak compactness.

If this is right

  • The limit Radon measures satisfy the equations in the natural weak sense.
  • The construction supplies a convergent numerical method by reducing the PDEs to ODEs.
  • Numerical results on the Toumi shock tube problem agree with those from other methods.
  • An independent standard scheme with splittings and vanishing viscosity reproduces exactly the same numerical solution.

Where Pith is reading between the lines

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

  • The ODE reduction may permit explicit analysis of long-time behavior or stability properties that remain inaccessible in the original PDE setting.
  • The appearance of Radon measures indicates that concentrations or discontinuities are intrinsic to the model and cannot be avoided by classical solution concepts.
  • The same approximation strategy could be tested on related hyperbolic systems with energy equations where standard weak solutions are also unavailable.

Load-bearing premise

The family of approximate solutions can be constructed so that they satisfy the necessary a priori bounds and reduce the original PDE system to a closed system of ODEs while preserving the structure needed for weak compactness to apply.

What would settle it

A direct check showing that the limit Radon measures fail to satisfy the weak form of the equations, or that the numerical outputs on the Toumi shock tube disagree with the solutions produced by independent schemes using splittings and vanishing viscosity.

Figures

Figures reproduced from arXiv: 1907.03611 by M. Colombeau.

Figure 1
Figure 1. Figure 1: Comparison of the asymptotic solution constructed in this pa￾per (black +) and the result from the transport-averaging-correction scheme of [14] with δ = 0 (red, continuous line). One observes a perfect coincidence. INSERT FIGURE 2 [PITH_FULL_IMAGE:figures/full_fig_p020_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: The liquid flow rate (kg/s) at time t=0.03 (top left panel), t=0.06 (top right panel), t=0.12 (bottom panel, in a twice longer tube). One observes that the top value increases very slowly with time, and that its width increases proportionally to time: one interval=2.5 centimeter in each panel. In the three figures the system of ODEs (6-8) is solved by the explicit Eu￾ler order 1 scheme. Space which is 100 … view at source ↗
Figure 1
Figure 1. Figure 1: Comparison of the asymptotic solution constructed in this paper (black +) and the result from the transport correction scheme (red, continuous line) [PITH_FULL_IMAGE:figures/full_fig_p025_1.png] view at source ↗
read the original abstract

We construct with full rigorous mathematical proof a family of approximate solutions to the Cauchy problem for the standard system of two fluid flows with energy equations and we pass to the limit by weak compactness to obtain Radon measures that satisfy the equations in a natural weak sense. Our method provides a convergent numerical method for the numerical calculation of these Radon measures by reducing the system of partial differential equations in the case of these approximate solutions to a system of ordinary differential equations. We observe numerically on the standard Toumi shock tube problem that the Radon measures from our method agree with the numerical solutions previously obtained by other authors with various different numerical methods. In a subsequent numerical paper, using a standard confident scheme with splittings and vanishing viscosity (independent on the above construction), we observe exactly the numerical solution given by our mathematical proof.

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

0 major / 2 minor

Summary. The manuscript constructs a family of approximate solutions to the Cauchy problem for the standard one-pressure two-fluid flow system with energy equations. These approximations are obtained via a regularization that reduces the PDE system to a closed system of ODEs while preserving uniform a priori bounds. Weak compactness is then applied to pass to the limit, yielding Radon measures that satisfy the original equations in the distributional sense. Numerical results on the Toumi shock tube problem are shown to agree with solutions from other methods, and the approach is claimed to provide a convergent numerical scheme.

Significance. If the construction and compactness arguments hold, the result supplies a rigorous existence theory for measure-valued weak solutions to this hyperbolic system of conservation laws, which is relevant for two-phase flows where shocks and discontinuities arise. The explicit reduction to ODEs also yields a numerical method whose convergence is tied to the mathematical proof. The independent numerical confirmation via a different scheme (vanishing viscosity with splittings) adds credibility. Strengths include the self-contained construction and the direct derivation of bounds from the ODE structure rather than ad-hoc assumptions.

minor comments (2)
  1. The abstract refers to 'a subsequent numerical paper' using a 'standard confident scheme'; this cross-reference should be clarified with a citation or title once available, even if the current work is self-contained.
  2. Notation for the two-fluid variables (e.g., densities, velocities, energies) and the precise form of the one-pressure closure should be stated explicitly in the introduction for readers unfamiliar with the 'standard system'.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the positive recommendation to accept.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivation consists of an explicit construction of approximate solutions that reduce the PDE system to a closed ODE system while preserving uniform a priori bounds, followed by a direct application of weak compactness to obtain Radon-measure solutions satisfying the equations distributionally. The numerical agreement on the Toumi problem is presented as an independent observation against prior external methods, and the subsequent numerical paper is explicitly stated to be independent of the construction. No step reduces by definition to its own input, renames a fitted quantity as a prediction, or relies on a self-citation chain for the central existence result.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No explicit free parameters, axioms, or invented entities are described in the abstract; the construction relies on standard functional-analytic tools for hyperbolic systems whose precise invocation cannot be audited from the given text.

pith-pipeline@v0.9.0 · 5661 in / 1130 out tokens · 18777 ms · 2026-05-25T01:03:22.862404+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages · 1 internal anchor

  1. [1]

    Abreu, M

    E. Abreu, M. Colombeau, E. Panov. Weak asymptotic methods for scalar equations and systems. J. Math. Anal. Appl. 44,2, 2016, pp. 1203-1232

    work page 2016

  2. [2]

    C. S. Avelar, P. R. Ribeiro, K. Seperhrnoori. Deep water gas kick simulation. J. Pet. Sci. Eng. 67, 2009, pp. 13-22

    work page 2009

  3. [3]

    C. S. Avelar. Well contol modelling: a finite difference approach (por- tuguese). Master Thesis, State University of Campinas, 26/08/2008

    work page 2008

  4. [4]

    Bendiksen, D

    K.H. Bendiksen, D. Maines, R. Moe, S. Nuland. The dynamic two- fluid model OLGA. Theory and Application. SPE Prod. Eng. 6, 1991, pp. 171-180

    work page 1991

  5. [5]

    D. Bestion. The physical closure laws in the CATHARE code. Nucl. Eng. Des. 124, 1990, pp. 229-245

    work page 1990

  6. [6]

    E. M. Bezeira. Study of well control consideruing phase behabior of gas-liquid mixture (portuguese). Master Thesis, State University of Campinas, 07/04/2006

    work page 2006

  7. [7]

    H. Brezis. Functional Analysis, Sobolev spaces and Partial Differential Equations. Springer, New York, Dordrecht, Heidelberg, London. 2010. 21

    work page 2010

  8. [8]

    K. I. Choi. A dynamic lagrange-euler numerical model for biphasic flows in oil wells. Master Thesis, CEPETRO, 01/03/1996

    work page 1996

  9. [9]

    Colombeau

    M. Colombeau. A method of projection of delta waves in a Godunov scheme and application to pressureless fluid dynamics. SIAM J. Nu- mer. Anal. 48, 5, 2010, pp. 1900-1919

    work page 2010

  10. [10]

    Colombeau

    M. Colombeau. Irregular shock wave solutions as continuations of the analytic solutions. Appl. Anal. 94,9,2015, pp.1800-1820

    work page 2015

  11. [11]

    Colombeau

    M. Colombeau. Weak asymptotic methods for 3-D selfgravitating pres- sureless fluids. Application to the creation and evolution of solar sys- tems from the fully nonlinear Euler-Poisson equations. J. of Mathe- matical Physics 56, 061506, 20 pages, 2015

    work page 2015

  12. [12]

    Colombeau

    M. Colombeau. Approximate solutions to the initial value problem for some compressible flows. Zeitschrift fur Angewandte Mathematik und Physik. 66, 2015,5, pp. 2575-2599

    work page 2015

  13. [13]

    Colombeau

    M. Colombeau. Asymptotic study of the initial value problem to a standard one pressure model of multifluid flows in nondivergence form. J. Diff. Eq. 260, 2016,1, pp. 197-217

    work page 2016

  14. [14]

    A numerical approximation for the standard one pressure system of two fluid flows with energy equations

    M. Colombeau. A numerical approximation of the standard one pres- sure system of two fluid flows with energy equations under its four versions. arXiv.org. 1808.08467

    work page internal anchor Pith review Pith/arXiv arXiv

  15. [15]

    S. Evje, T. Flatten. Hybrid Flux-splitting Schemes for a common two fluid model. J. Comput. Physics 192, 2003, pp. 175-210

    work page 2003

  16. [16]

    K.T. Joseph. Boundary layers in approximate solutions. Trans. Amer. Math. Soc. 314, 1989, pp. 709-726

    work page 1989

  17. [17]

    K.T. Joseph. P.L. Sadchev. Exact solutions for some nonconservative hyperbolic systems. Intern. J. Nonlinear Mechanics 38,2003, 9, pp. 1377-1386

    work page 2003

  18. [18]

    K.T. Joseph. Generalized solutions to a Cauchy problem for a noncon- servative hyperbolic system. J. Math. Anal. Appl. 2007, 1997, pp.361- 387

    work page 2007

  19. [19]

    Joseph, Ph LeFloch

    K.T. Joseph, Ph LeFloch. Singular limits for the Riemann problem CR Math Acad. Sci. Paris 344, 2007,1, pp.59-64. Analytical Approaches to mutidimensional balance laws. pp. 143-171. Nova Sci. Publi. New York 2006. 22

    work page 2007

  20. [20]

    Joseph, M.R

    K.T. Joseph, M.R. Sahoo. Vanishing viscosity appoach to a system of conservation laws admitting δ”-waves. Comm. Pure Appl. Analysis 12, 2013, 5, pp.2091-2118

    work page 2013

  21. [21]

    K.T. Joseph. Exact solution of a nonconservative system in elastody- namics. Electron. J. Diff. Eq. 2015, art. 259, 7p

    work page 2015

  22. [22]

    B.L. Keyfitz. Properties of Nonhyperbolic Models for Incompressible Two Phase Flows. Report 77204-3476 of University of Houston, Fields Institute and University of Ohio

    work page

  23. [23]

    B.L. Keyfitz. Hold that light! Modelling of traffic flow by Differential Equations. In ”‘Six Themes on variation”’, Robert Hardt editor, Stu- dent Mathematical Library, vol. 26, AMS, Providence, Rhode Island, USA, 2004, pp. 127-153

    work page 2004

  24. [24]

    B.L. Keyfitz. Admissibility conditions for shocks in conservation laws that change type. Siam J. Math. Anal. 22, 5, 1991, pp. 1284-1292

    work page 1991

  25. [25]

    B.L. Keyfitz. Change of type in simple models for two phase flows . In Viscous Profiles and Numerical Methods for shock waves”’ ed. M. Shearer, 1991, pp 84-104, SIAM, Philadelphy

    work page 1991

  26. [26]

    B.L. Keyfitz. Change of type in three phase flows: a simple analog. J. Diff. Eq. 80, 1989, pp. 280-305

    work page 1989

  27. [27]

    B.L. Keyfitz. Conservation laws that change type and porous medium flow:a review. In ”‘Modeling and Analysis of Diffusive and Advective Processes in Geosciences”’pp. 122-145. eds W.E. Fitzgibbon and M.F. Wheeler, SIAM, Philadelphia, 1992

    work page 1992

  28. [28]

    B.L. Keyfitz. Multiphase saturation equations, change of type and inaccessible regions. In Proceedings of the Oberwolfach conference on porous media. pp. 103-116. eds J. Douglas, C.J. Van Duijn and U. Hornung. Birhkhauser, 1993

    work page 1993

  29. [29]

    B.L. Keyfitz. A geometric theory of conservation laws which change type. Zeitschrift fur Angewandte Mathematik und Mechanik, 75, pp. 571-581, 1995

    work page 1995

  30. [30]

    Keyfitz, M

    B.L. Keyfitz, M. Lopes-Filho. A geometric study of shocks in equations that change type. J. of Dyn and Diff. Eqs, 6,3, 1994, pp. 351-393. 23

    work page 1994

  31. [31]

    Keyfitz, R

    B.L. Keyfitz, R. Sanders, M. Sever. Lack of hyperbolicity in the two- fluid model for two-phase incompressible flows. Discrete and continu- ous dynamical systems, serie B,3,4,2003,pp. 541-563

    work page 2003

  32. [32]

    Kunzinger, G

    M. Kunzinger, G. Rein, R. Steinbauer, G. Teschl. Global weak solution of the relativistic Vlassov-Klein Gordon system. Comm. Math. Pkys. 238, 2003, 1-2, pp. 367-378

    work page 2003

  33. [33]

    Larsen, E

    M. Larsen, E. Hustvedt, P. Hedne, T. Straume. Petra: a novel com- puter code for simulation of slug flow. In SPI Annual Technical Con- ference and Exhibition. SPI 38841, 1997, pp. 1-12

    work page 1997

  34. [34]

    P.D. Lax. Computational Fluid Dynamics. J. Sci. Comput. 31, 2007, pp. 185-193

    work page 2007

  35. [35]

    P.D. Lax. Mathematics and Physics. Bull. A.M.S. 45, 1, 2008, pp. 135-152

    work page 2008

  36. [36]

    C. H. O. Loiola. Compositional well control simulation (portuguese). Master Thesis, CEPETRO, State university of Campinas,22/09/2015

    work page 2015

  37. [37]

    Martinez-Ferrer, T

    P.J. Martinez-Ferrer, T. Flatten, S.T. Munkejord. On the effect of temperature and velocity relaxation in two-phase flow models M2AN 46, 2012, pp. 411-442

    work page 2012

  38. [38]

    A. Morin,T. Flatten, S.V. Munkejord. A Roe scheme for a compress- ible six-equation two fluid model. International Journal for Numerical Methods in fluids, 2013,72, pp. 478-504

    work page 2013

  39. [39]

    S. T. Munkejord, S. Evje, T. Flatten. A Musta scheme for a non- conservative two-fluid model. SIAM J. Sci. Comput. 31, 4, 2009, pp. 2587-2622

    work page 2009

  40. [40]

    Riemann Solvers and Numerical Methods for Fluid Dynamics

    E.Toro. Riemann Solvers and Numerical Methods for Fluid Dynamics. Springer Verlag, 1999

    work page 1999

  41. [41]

    Toumi, A

    I. Toumi, A. Kumbaro. An approximate linearized Riemann solver for a two fluid model. J. Comput. Physics 66, 1986, pp. 62-82

    work page 1986

  42. [42]

    JSI Report IJS-DP-8841, Josef Stefan Insti- tute, Ljubljana, Slovenia, 2004

    WAHA3 Code Manual. JSI Report IJS-DP-8841, Josef Stefan Insti- tute, Ljubljana, Slovenia, 2004

    work page 2004

  43. [43]

    H. B. Stewart, B. Wendroff. Two-phase flows: Models and Methods. J. Comput. Phys. 56, 1984, pp. 363-409. 24 Figure 1. Comparison of the asymptotic solution constructed in this paper (black +) and the result from the transport correction scheme (red, continuous line). 0 500 1000 1500 2000 2500 3000 3500 4000 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 gas volume fract...

    work page 1984