pith. sign in

arxiv: 2604.20512 · v1 · submitted 2026-04-22 · ⚛️ physics.flu-dyn

Nonisothermal global-pressure exactness in fractured multiphase flow with evolving fracture aperture

Christian Tantardini , Fernando Alonso-Marroquin This is my paper

Pith reviewed 2026-05-09 23:31 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords global pressure formulationmultiphase Darcy flownonisothermal flowfractured mediacapillary pressureexact equivalenceevolving apertureheat transport
0
0 comments X

The pith

Exact equivalence between global-pressure and phase-pressure multiphase flow holds only when a mobility-weighted capillary one-form closes on the saturation-temperature space.

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

The paper establishes the precise conditions for a global pressure to drive the total flux identically to the set of phase pressures in nonisothermal multiphase flow. It focuses on temperature-dependent mobilities and capillary pressures inside fractured media whose aperture can change. The governing structure is the closure of a weighted capillary one-form over the joint saturation-temperature variables, which recovers the usual saturation-only compatibility rules plus a new cross-term that appears solely when temperature varies. When the condition holds, the global-pressure model is mathematically identical to the phase-pressure version; when it fails, the total flux cannot be expressed as the gradient of any scalar pressure. This distinction matters for efficient simulation of heat-carrying flows in reservoirs and geothermal systems, where temperature changes and fracture dynamics can move the system across exact and approximate regimes.

Core claim

Equivalence between global-pressure and phase-pressure formulations is controlled by closure of the mobility-weighted capillary one-form on the augmented saturation-temperature state space. This closure supplies both the classical compatibility conditions inside the saturation sector and an additional mixed saturation-temperature condition required only in the nonisothermal case. The criterion is inserted into a reduced matrix-fracture model that incorporates heat transport, matrix-fracture thermal exchange, and evolving fracture aperture, with benchmarks confirming three regimes: globally exact, exact on fixed-temperature slices but not across the full space, and fully nonexact.

What carries the argument

The mobility-weighted capillary one-form on the saturation-temperature state space, whose closure decides whether the total flux is exactly the gradient of a single scalar pressure.

If this is right

  • Thermal forcing alone can drive a fractured system from an exact regime into a nonexact one even if saturation compatibility is satisfied.
  • Aperture evolution changes the trajectory through saturation-temperature space and can therefore switch the model between exact and nonexact behavior.
  • When full closure fails, a temperature-slice least-squares projection still supplies a conservative scalar pressure together with explicit defect measures.
  • The same one-form structure unifies the isothermal compatibility conditions with their nonisothermal extension inside a single reduced model.

Where Pith is reading between the lines

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

  • The closure test could be applied to other state variables such as composition or salinity to obtain analogous exactness criteria for compositional flows.
  • In practice, measuring local temperature gradients alongside saturation profiles would indicate whether a global-pressure simulation remains accurate without additional projection.
  • Dynamic aperture laws that depend on effective stress could be rewritten to track the path through state space and flag intervals where exactness is lost.

Load-bearing premise

The temperature-dependent mobility and capillary-pressure functions supplied by the user must make the mobility-weighted capillary one-form closed on the saturation-temperature domain.

What would settle it

A direct numerical comparison, for constitutive data violating the mixed saturation-temperature closure, between the total flux computed from phase pressures and the flux implied by the gradient of any candidate global pressure.

Figures

Figures reproduced from arXiv: 2604.20512 by Christian Tantardini, Fernando Alonso-Marroquin.

Figure 1
Figure 1. Figure 1: FIG. 1. Benchmark A: line diagnostics on the augmented state space ∆ [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Benchmark A: representative state-space defect maps on fixed-temperature slices of ∆ [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Benchmark B: fixed-aperture fractured thermal-front benchmark. (a) Time histories of the maximum local defects in [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Benchmark C: fractured thermal-front benchmark with evolving aperture. (a) Time histories of the maximum local [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Benchmark C: aperture-feedback metrics on sepa [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
read the original abstract

Global-pressure formulations recast multiphase Darcy flow in terms of a single pressure driving the total flux. Their exact equivalence to phase-pressure formulations, however, holds only when the constitutive data satisfy the compatibility conditions required for a total-differential structure and its generalized nonisothermal extension. In this work, we derive the corresponding exactness criterion for temperature-dependent mobilities and capillary pressures. We show that equivalence is governed by the closure of a mobility-weighted capillary one-form on the augmented state space of saturation and temperature. This yields both the classical compatibility conditions within the saturation sector and a distinct mixed saturation--temperature condition that arises only in the nonisothermal setting. We then incorporate this structure into a reduced matrix--fracture model with heat transport, matrix--fracture thermal exchange, and evolving fracture aperture. Numerical benchmarks recover the three regimes predicted by the theory: globally exact, exact on each fixed-temperature slice but not on the full saturation--temperature space, and fully nonexact. In fractured systems, thermal forcing alone can drive transitions between these regimes, while aperture evolution changes the path through state space. When exactness fails, a least-squares projection performed independently on each fixed-temperature slice provides a conservative scalar-pressure surrogate together with quantitative defect diagnostics. The resulting framework unifies nonisothermal exactness theory, fractured-flow dynamics, and conservative reduced closure within a single global-pressure formulation.

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 / 3 minor

Summary. The paper derives an exactness criterion for global-pressure formulations of nonisothermal multiphase Darcy flow, showing that equivalence to phase-pressure formulations holds if and only if a mobility-weighted capillary one-form closes on the augmented (saturation, temperature) state space. This produces the classical saturation-sector compatibility conditions together with one new mixed saturation-temperature condition. The criterion is incorporated into a reduced matrix-fracture model that includes heat transport, matrix-fracture thermal exchange, and evolving fracture aperture (treated as path-dependent forcing that does not enlarge the state space). Numerical benchmarks are stated to recover the three predicted regimes (globally exact, exact on fixed-temperature slices, and fully nonexact), with a least-squares projection on each temperature slice offered as a conservative surrogate when exactness fails.

Significance. If the derivation is correct, the work supplies a parameter-free mathematical criterion that unifies nonisothermal exactness theory with fractured multiphase flow and supplies a practical defect diagnostic plus surrogate closure. The differential-form approach, the explicit identification of the new mixed condition, and the clean separation of aperture evolution as path-dependent forcing are genuine strengths. The numerical recovery of regime transitions driven by thermal forcing alone is a useful illustration of the theory's predictive power.

minor comments (3)
  1. The abstract and introduction refer to 'three predicted regimes' without enumerating them; a short explicit list (globally exact, slice-exact, nonexact) would improve immediate readability.
  2. The numerical benchmarks section states that the three regimes are recovered but provides no quantitative error norms, L2 residuals, or direct side-by-side comparison tables against the phase-pressure reference solution; adding such metrics would strengthen the validation without altering the central claim.
  3. Notation for the mobility-weighted capillary one-form is introduced in the derivation but its explicit coordinate expression on the (S,T) plane is not repeated in the fractured-model section; a single boxed equation restating the closure condition would aid cross-referencing.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The referee's summary accurately captures the derivation of the mobility-weighted capillary one-form closure condition on the saturation-temperature state space, the resulting mixed compatibility condition, the incorporation into the reduced matrix-fracture model with evolving aperture, and the numerical recovery of the three exactness regimes. No major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper constructs a mobility-weighted capillary one-form directly from the input constitutive maps (temperature-dependent mobilities and capillary pressures) and requires its closure on the saturation-temperature state space. The resulting compatibility conditions, including the new mixed saturation-temperature condition, are mathematical consequences of dω = 0 rather than inputs or fitted quantities. The reduced fractured model is assembled to preserve this closed structure, with aperture evolution entering only as a path-dependent forcing term that does not enlarge the state space or alter the closure criterion. No load-bearing self-citation, ansatz smuggling, or renaming of known results appears in the derivation; the exactness criterion is therefore an independent output of the differential-form requirement applied to the given constitutive data.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard domain assumptions for multiphase Darcy flow and the mathematical requirement that a one-form be closed for exactness; no free parameters or new postulated entities are introduced in the abstract.

axioms (2)
  • domain assumption Multiphase flow obeys Darcy's law with temperature-dependent mobilities and capillary pressures.
    This is the standard physical model invoked for the system.
  • standard math Exact equivalence of global-pressure and phase-pressure formulations holds precisely when the mobility-weighted capillary one-form is closed.
    This follows from the theory of differential forms and total differentials on the state space.

pith-pipeline@v0.9.0 · 5550 in / 1406 out tokens · 26952 ms · 2026-05-09T23:31:46.943191+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

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

  1. [1]

    S. E. Buckley and M. C. Leverett, Transactions of the AIME146, 107 (1942)

    work page 1942

  2. [2]

    M. C. Leverett, Transactions of the AIME142, 152 (1941)

    work page 1941

  3. [3]

    Chavent and J

    G. Chavent and J. Jaffr´ e,Mathematical Models and Fi- nite Elements for Reservoir Simulation: Single Phase, Multiphase and Multicomponent Flows through Porous Media(North-Holland, Amsterdam, 1986)

    work page 1986

  4. [4]

    Amaziane and M

    B. Amaziane and M. Jurak, Comptes Rendus M´ ecanique 336, 600 (2008)

    work page 2008

  5. [5]

    Amaziane, M

    B. Amaziane, M. Jurak, and A. ˇZ. Keko, Journal of Dif- ferential Equations250, 1685 (2011)

    work page 2011

  6. [6]

    Holden, SIAM Journal on Applied Mathematics50, 667 (1990)

    L. Holden, SIAM Journal on Applied Mathematics50, 667 (1990)

    work page 1990

  7. [7]

    Chavent, Applicable Analysis88, 1527 (2009)

    G. Chavent, Applicable Analysis88, 1527 (2009)

    work page 2009

  8. [8]

    di Chiara Roupert, G

    R. di Chiara Roupert, G. Chavent, and G. Sch¨ afer, Jour- nal of Computational Physics229, 4762 (2010)

    work page 2010

  9. [9]

    Sch¨ afer, R

    G. Sch¨ afer, R. di Chiara Roupert, A. H. Alizadeh, and M. Piri, Advances in Water Resources145, 103731 (2020)

    work page 2020

  10. [10]

    Chen and R

    Z. Chen and R. E. Ewing, Journal of Computational Physics132, 362 (1997)

    work page 1997

  11. [11]

    Global Buckley--Leverett for Multicomponent Flow in Fractured Media: Isothermal Equation-of-State Coupling and Dynamic Capillarity

    C. Tantardini and F. Alonso-Marroquin, Global buckley– leverett for multicomponent flow in fractured media: Isothermal equation-of-state coupling and dynamic cap- illarity (2025), preprint, arXiv:2511.06233 [physics.flu- dyn]

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  12. [12]

    S. M. Hassanizadeh and W. G. Gray, Water Resources Research29, 3389 (1993)

    work page 1993

  13. [13]

    H. Y. She and B. E. Sleep, Water Resources Research34, 2587 (1998)

    work page 1998

  14. [14]

    Jurak, A

    M. Jurak, A. Koldoba, A. Konyukhov, and L. Pankratov, Comptes Rendus M´ ecanique347, 920 (2019)

    work page 2019

  15. [15]

    D. C. Standnes and P. Fotland, Water Resources Re- search57, e2021WR029887 (2021)

    work page 2021

  16. [16]

    C. T. Miller, K. Bruning, C. L. Talbot, J. E. McClure, and W. G. Gray, Water Resources Research55, 6825 (2019)

    work page 2019

  17. [17]

    Berre, F

    I. Berre, F. Doster, and E. Keilegavlen, Transport in Porous Media130, 215 (2019)

    work page 2019

  18. [18]

    Pruess and T

    K. Pruess and T. N. Narasimhan, SPE Journal25, 14 (1985)

    work page 1985

  19. [19]

    Ahmed, J

    E. Ahmed, J. Jaffr´ e, and J. E. Roberts, Mathematics and Computers in Simulation137, 49 (2017)

    work page 2017

  20. [20]

    Y. Wang, S. de Hoop, D. Voskov, D. Bruhn, and G. Bertotti, Advances in Water Resources154, 103985 17 (2021)

    work page 2021

  21. [21]

    D. Y. Ding, Transport in Porous Media143, 169 (2022)

    work page 2022

  22. [22]

    P. A. Witherspoon, J. S. Y. Wang, K. Iwai, and J. E. Gale, Water Resources Research16, 1016 (1980)

    work page 1980

  23. [23]

    Ghassemi, A

    A. Ghassemi, A. Nygren, and A. Cheng, Geothermics37, 525 (2008)

    work page 2008

  24. [24]

    S. N. Pandey, A. Chaudhuri, and S. Kelkar, Geothermics 65, 17 (2017)

    work page 2017

  25. [25]

    Bisdom, G

    K. Bisdom, G. Bertotti, and H. M. Nick, Journal of Geo- physical Research: Solid Earth121, 4045 (2016)

    work page 2016

  26. [26]

    Bear,Dynamics of Fluids in Porous Media(American Elsevier, New York, 1972)

    J. Bear,Dynamics of Fluids in Porous Media(American Elsevier, New York, 1972)

    work page 1972

  27. [27]

    Weishaupt, V

    K. Weishaupt, V. Joekar-Niasar, and R. Helmig, Water Resources Research57, e2020WR028772 (2021)

    work page 2021

  28. [28]

    Hamidi, M

    S. Hamidi, M. Massoudi, and T. Antoun, International Journal of Heat and Mass Transfer135, 641 (2019)

    work page 2019

  29. [29]

    Heinzeet al., Earth-Science Reviews251, 104786 (2024)

    T. Heinzeet al., Earth-Science Reviews251, 104786 (2024)

    work page 2024

  30. [30]

    Heinze, S

    T. Heinze, S. Hamidi, and B. Galvan, Geothermics65, 10 (2017)

    work page 2017

  31. [31]

    S. N. Pandey, A. Chaudhuri, S. Kelkar, S. Karra, and A. Ghassemi, Geothermics65, 17 (2017)

    work page 2017