pith. sign in

arxiv: 2603.28981 · v2 · submitted 2026-03-30 · 🧮 math.NA · cs.NA· physics.flu-dyn

A bounded-interval multiwavelet formulation with conservative finite-volume transport for one-dimensional Buckley--Leverett waterflooding

Christian Tantardini This is my paper

Pith reviewed 2026-05-14 01:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.flu-dyn MSC 65M0865M1276S05
keywords Buckley-Leverett equationmultiwaveletfinite-volume methodconservative transportwaterfloodinghyperbolic conservation lawssaturation profilesBerea benchmark
0
0 comments X

The pith

A hybrid finite-volume scheme with bounded-interval multiwavelet reconstruction solves the one-dimensional Buckley-Leverett equation while preserving shocks and conservation.

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

The paper develops a hybrid method for the Buckley-Leverett equation that performs the nonlinear saturation update with a conservative finite-volume scheme using monotone fluxes. This step ensures correct entropy-admissible shock transport, while the state is represented and reconstructed in a bounded-interval multiwavelet basis to supply a hierarchical multiresolution description. The approach maintains conservation and monotonicity at discontinuities. Validation on a Berea benchmark shows close agreement with reference solutions in saturation histories, spatial profiles, front locations, and error measures. The multiwavelet part tracks the underlying finite-volume state with essentially exact fidelity.

Core claim

The central claim is that a bounded-interval multiwavelet formulation can be paired with a conservative finite-volume transport step for the deterministic one-dimensional Buckley-Leverett equation. The finite-volume scheme handles the entropy-admissible shocks via monotone numerical fluxes, while the multiwavelet basis provides the evolving state representation and reconstruction. This hybrid strategy preserves the correct shock-compatible transport mechanism and simultaneously yields a hierarchical multiresolution description of the saturation field. Numerical tests on the Berea benchmark confirm excellent agreement with reference profiles across probe histories, spatial distributions, and

What carries the argument

Bounded-interval multiwavelet reconstruction inserted around a conservative finite-volume transport step with monotone fluxes. It carries the argument by enabling hierarchical representation of the saturation field without violating conservation or shock conditions.

If this is right

  • The formulation supplies a reliable first step toward native multiwavelet transport solvers for porous-media flow.
  • The method maintains correct shock speeds and global conservation in nonlinear hyperbolic conservation laws.
  • It enables simultaneous multiresolution analysis of saturation fronts during waterflooding.
  • The hybrid structure supports extension to higher-dimensional or heterogeneous reservoir models.

Where Pith is reading between the lines

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

  • The multiwavelet coefficients could enable adaptive refinement focused on the moving front, reducing cost in large simulations.
  • Similar hybrids may apply to other hyperbolic problems with shocks, such as Euler equations or traffic flow models.
  • The approach could later incorporate spatially varying permeability to handle realistic heterogeneous media.

Load-bearing premise

Inserting the bounded-interval multiwavelet reconstruction around the finite-volume transport step does not introduce non-monotone artifacts or violate the entropy condition at shocks.

What would settle it

If the computed saturation profiles, shock front locations, or probe histories deviate significantly from the reference Buckley-Leverett solution on the Berea benchmark or display oscillations, the central claim would be falsified.

Figures

Figures reproduced from arXiv: 2603.28981 by Christian Tantardini.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the saturation history at the probe [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Spatial saturation profiles [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Temporal evolution of the most active dyadic de [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Additional transport diagnostics as functions of pore [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 2
Figure 2. Figure 2: Several features deserve attention. First, the [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
read the original abstract

We develop a hybrid conservative finite-volume / bounded-interval multiwavelet formulation for the deterministic one-dimensional Buckley--Leverett equation. Because Buckley--Leverett transport is a nonlinear hyperbolic conservation law with entropy-admissible shocks, the saturation update is performed by a conservative finite-volume scheme with monotone numerical fluxes, while the evolving state is represented and reconstructed in a bounded-interval multiwavelet basis. This strategy preserves the correct shock-compatible transport mechanism and simultaneously provides a hierarchical multiresolution description of the solution. Validation against reference Buckley--Leverett profiles for a Berea benchmark shows excellent agreement in probe saturation histories, spatial profiles, front-location diagnostics, and global error measures. The multiwavelet reconstruction also tracks the internal finite-volume state with essentially exact fidelity. The resulting formulation provides a reliable first step toward more native multiwavelet transport solvers for porous-media flow.

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 develops a hybrid method for the one-dimensional Buckley-Leverett equation in which the nonlinear hyperbolic transport step is performed exclusively with a conservative monotone finite-volume scheme (ensuring conservation and entropy-admissible shocks), while the saturation state is represented and reconstructed in a bounded-interval multiwavelet basis. Validation on a Berea benchmark reports close agreement with reference solutions in probe saturation histories, spatial profiles, front-location diagnostics, and global L1/L2 error measures, together with machine-precision recovery of the internal finite-volume cell averages by the multiwavelet reconstruction.

Significance. If the reported fidelity holds, the formulation supplies a direct, parameter-free route to embed hierarchical multiresolution representations inside standard conservative transport schemes for porous-media flow. The exact reconstruction property and the separation of the monotone update from the representation step are concrete strengths that could support subsequent adaptive or higher-dimensional extensions without re-deriving entropy fixes.

minor comments (2)
  1. [Abstract and Numerical results] The abstract and validation section describe agreement as 'excellent' and 'essentially exact' without quoting the concrete L1/L2 tolerances or front-location errors; adding a compact table of these quantities (with reference-solution details) would make the quantitative claim easier to assess.
  2. [Method formulation] The bounded-interval multiwavelet construction is introduced without an explicit reference to the precise wavelet family or the interval-extension technique employed; a short citation or one-sentence definition in §2 would improve accessibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of our manuscript and the recommendation for minor revision. The referee's summary accurately captures the key contributions of our hybrid conservative finite-volume and bounded-interval multiwavelet formulation for the Buckley-Leverett equation. Since no specific major comments were provided in the report, we have no additional points to address.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper constructs a hybrid method by applying a standard conservative monotone finite-volume scheme (with entropy-admissible fluxes) exclusively to the nonlinear Buckley-Leverett update, while restricting the bounded-interval multiwavelet basis to hierarchical representation and reconstruction only. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or renamed input; the reconstruction is reported to recover internal cell averages to machine precision, and validation metrics are compared against independent reference solutions. The central claim therefore rests on established conservation-law numerics and multiwavelet theory rather than any self-referential reduction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard mathematical properties of nonlinear hyperbolic conservation laws and the existence of monotone numerical fluxes; no new free parameters, ad-hoc constants, or invented physical entities are introduced in the abstract.

axioms (1)
  • domain assumption Buckley-Leverett transport is a nonlinear hyperbolic conservation law possessing entropy-admissible shocks.
    This premise directly motivates the choice of a conservative finite-volume scheme with monotone fluxes.

pith-pipeline@v0.9.0 · 5451 in / 1205 out tokens · 55275 ms · 2026-05-14T01:19:22.873826+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

24 extracted references · 24 canonical work pages

  1. [1]

    S. E. Buckley and M. C. Leverett,Mechanism of fluid displacement in sands, Transactions of the AIME146, 107–116 (1942)

    work page 1942

  2. [2]

    S. E. Buckley and M. C. Leverett,Mechanism of fluid displacements in sands, Transactions of the AIME146, 107–116 (1952), classical work on immiscible two-phase displacement theory using fractional flow analysis

    work page 1952

  3. [3]

    Belazreg, S

    L. Belazreg, S. M. Mahmood, and A. Aulia,Novel ap- proach for predicting water alternating gas injection re- covery factor, Journal of Petroleum Exploration and Pro- duction Technology9, 2893–2910 (2019), forecasting per- formance of immiscible WAG floods using analytical pre- diction tools

    work page 2019

  4. [4]

    E. F. Kaasschieter,Solving the buckley–leverett equation with gravity in a heterogeneous porous medium, Compu- tational Geosciences3, 23–48 (1999), hyperbolic limit and capillary regularisation discussion. 10

    work page 1999

  5. [5]

    R. J. LeVeque,Finite Volume Methods for Hyperbolic Problems, Cambridge Texts in Applied Mathematics (Cambridge University Press, Cambridge, 2002)

    work page 2002

  6. [6]

    K. R. Spayd,The buckley–leverett equation with dynamic capillary pressure, SIAM Journal on Applied Mathe- matics71, 1275–1300 (2011), shows how the Buck- ley–Leverett model becomes pseudo-parabolic when a dy- namic capillary term is included, and reduces to a hyper- bolic conservation law in the capillarity-free limit

    work page 2011

  7. [7]

    C. J. van Duijn, X. Cao, and I. S. Pop,Two-phase flow in porous media: Dynamic capillarity and heterogeneous media, Transport in Porous Media109, 333–357 (2015), analyzes two-phase flow with dynamic capillary effects in heterogeneous media, highlighting the transition between hyperbolic and regularised regimes

    work page 2015

  8. [8]

    X. Cao, I. S. Pop,et al.,Degenerate two-phase porous media flow model with dynamic capillarity, Journal of Differential Equations260, 2418–2456 (2016), provides mathematical analysis of a degenerate elliptic–parabolic (pseudo-parabolic) two-phase flow model including dy- namic capillary pressure, with existence/uniqueness re- sults

    work page 2016

  9. [9]

    C.-W. Shu,High order weno and dg methods for time- dependent convection-dominated problems, Journal of Computational Physics316, 598–658 (2016), survey of high-order finite volume/WENO and discontinuous Galerkin methods for hyperbolic conservation laws rel- evant to saturation transport

    work page 2016

  10. [10]

    Jiang and C.-W

    G.-S. Jiang and C.-W. Shu,Weighted essentially non- oscillatory (weno) methods, Journal of Computational Physics126, 202–228 (1996), foundational work on high- resolution WENO reconstruction for hyperbolic PDEs

    work page 1996

  11. [11]

    Harten,Multiresolution algorithms for the numeri- cal solution of hyperbolic conservation laws, Communica- tions on Pure and Applied Mathematics48, 1305–1342 (1995)

    A. Harten,Multiresolution algorithms for the numeri- cal solution of hyperbolic conservation laws, Communica- tions on Pure and Applied Mathematics48, 1305–1342 (1995)

    work page 1995

  12. [12]

    Cohen, S

    A. Cohen, S. M. Kaber, S. M¨ uller, and M. Postel,Fully adaptive multiresolution finite volume schemes for con- servation laws, Mathematics of Computation72, 183– 225 (2003)

    work page 2003

  13. [13]

    M¨ uller and Y

    S. M¨ uller and Y. Stiriba,Fully adaptive multiscale schemes for conservation laws employing locally varying time stepping, Journal of Scientific Computing30, 493– 531 (2007)

    work page 2007

  14. [14]

    M. J. Vuik and J. K. Ryan,Multiwavelet troubled- cell indicator for discontinuity detection of discontinuous galerkin schemes, Journal of Computational Physics270, 138–160 (2014)

    work page 2014

  15. [15]

    Gerhard, F

    N. Gerhard, F. Iacono, G. May, S. M¨ uller, and R. Sch¨ afer,A high-order discontinuous galerkin dis- cretization with multiwavelet-based grid adaptation for compressible flows, Journal of Scientific Computing62, 25–52 (2015)

    work page 2015

  16. [16]

    Huang and Y

    J. Huang and Y. Cheng,An adaptive multiresolution dis- continuous galerkin method with artificial viscosity for scalar hyperbolic conservation laws in multidimensions, SIAM Journal on Scientific Computing42, A2943–A2973 (2020)

    work page 2020

  17. [17]

    Pettersson and H

    P. Pettersson and H. A. Tchelepi,Stochastic galerkin framework with locally reduced bases for nonlinear two- phase transport in heterogeneous formations, Computer Methods in Applied Mechanics and Engineering310, 367–387 (2016)

    work page 2016

  18. [18]

    A. T. Corey,The interrelation between gas and oil relative permeabilities, Producers Monthly19, 38–41 (1954)

    work page 1954

  19. [19]

    R. H. Brooks and A. T. Corey,Hydraulic Properties of Porous Media, Hydrology Paper 3 (Colorado State Uni- versity, Fort Collins, Colorado, 1964)

    work page 1964

  20. [20]

    R. Bast, M. Bjorgve, R. Di Remigio, A. Durdek, L. Fre- diani, E. Fossgaard, G. Gerez, S. R. Jensen, J. Juselius, S. Lehtola, R. Monstad, and P. Wind, Mrcpp: Multires- olution computation program package (2023)

    work page 2023

  21. [21]

    Battistella, M

    E. Battistella, M. Bjorgve, R. Di Remigio, L. Frediani, G. Gerez, and S. R. Jensen, Vampyr: Very accurate mul- tiresolution python routines (2023)

    work page 2023

  22. [22]

    Vampyr repository (2024), accessed: 9 February 2024

    work page 2024

  23. [23]

    Bjørgve, C

    M. Bjørgve, C. Tantardini, S. R. Jensen, G. A. Gerez S., P. Wind, R. Di Remigio Eik˚ as, E. Dinvay, and L. Fredi- ani,Vampyr—a high-level python library for mathemati- cal operations in a multiwavelet representation, The Jour- nal of Chemical Physics160, 162502 (2024)

    work page 2024

  24. [24]

    Male,Pywaterflood: Well connectivity analysis through capacitance-resistance modeling, Journal of Open Source Software9, 6191 (2024)

    F. Male,Pywaterflood: Well connectivity analysis through capacitance-resistance modeling, Journal of Open Source Software9, 6191 (2024)

    work page 2024