A new Timestep Criterion for the Simulation of Immiscible Two-Phase Flow with IMPES Solvers
Pith reviewed 2026-05-22 12:09 UTC · model grok-4.3
The pith
A CFL timestep criterion based on computed velocity derivatives allows stable IMPES simulations of two-phase flows without any parameter adjustments.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The new timestep criterion uses the Courant-Friedrichs-Lewy condition applied to numerically computed velocity derivatives to adapt the timestep size regardless of whether pressure drop or capillary forces dominate the flow. When used with the iterative IMPES solver, it achieves efficiency and robustness in various flow scenarios with compressible and incompressible fluids, reaches expected stationary states even with discontinuous porous media parameters, requires fewer iterations than the Coats criterion with similar accuracy on test problems, and yields better non-wetting phase mass conservation in air compression cases.
What carries the argument
The novel CFL-based timestep criterion that relies on numerically computed velocity derivatives to determine adaptive timestep sizes independent of flow regime.
If this is right
- It requires fewer time iterations than the Coats criterion while keeping comparable accuracy on the Buckley-Leverett problem.
- It reaches the expected stationary states in gravity-capillary equalization with known solutions.
- It achieves significantly better non-wetting phase mass conservation in air compression examples.
- It handles cases with discontinuous porosity, permeabilities, and capillary pressure functions without instability.
Where Pith is reading between the lines
- If the velocity derivatives can be computed reliably, this criterion could be applied to other numerical schemes for multiphase flows.
- Simulations of reservoir flows might become more automated since no manual timestep parameter tuning is needed.
- Improved mass conservation could lead to more trustworthy long-term predictions in compressible fluid scenarios.
Load-bearing premise
The numerically computed velocity derivatives remain sufficiently accurate and stable even when porosity, permeabilities, and capillary pressure functions are discontinuous.
What would settle it
Running the method on a test case with extreme discontinuities in capillary pressure that causes the simulation to become unstable or lose significant mass would show the criterion does not work as claimed.
read the original abstract
We present an iterative IMPES solver and a novel timestep criterion for the simulation of immiscible two-phase flow involving compressible fluid phases. The novel timestep criterion uses the Courant-Friedrichs-Lewy (CFL) condition and employs numerically computed velocity derivatives to adapt the timestep size, regardless of the dominant flow characteristics. The solver combined with this timestep criterion demonstrates both efficiency and robustness across a range of flow scenarios, including pressure drop dominated and capillary dominated flows with compressible and incompressible fluid phases, without the need to adjust any numerical parameters. Furthermore, it successfully reaches the expected stationary states in a case involving discontinuous porous media parameters such as porosity, permeabilities, and capillary pressure function. Comparison with the established Coats timestep criterion reveals that our approach requires fewer time iterations while maintaining comparable accuracy on the Buckley-Leverett problem and a gravity-capillary equalization example with a known stationary state. Additionally, in an example with air compression, the new timestep criterion leads to a significantly improved non-wetting phase mass conservation compared to the Coats criterion.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces an iterative IMPES solver paired with a novel timestep criterion for immiscible two-phase flow with compressible phases. The criterion adapts the timestep via the CFL condition using on-the-fly numerical derivatives of velocity, independent of whether pressure drop or capillary forces dominate. Tests on the Buckley-Leverett problem, gravity-capillary equalization to a known stationary state, and an air-compression case are reported to require fewer iterations than the Coats criterion while achieving comparable accuracy and markedly better non-wetting-phase mass conservation; the method is also shown to reach the expected stationary state on a discontinuous porous-medium example without any parameter tuning.
Significance. If substantiated, the parameter-free adaptive timestep would address a persistent practical limitation of IMPES schemes in heterogeneous and compressible two-phase flow, offering both computational efficiency and improved conservation properties. The explicit demonstration of robustness across pressure-dominated, capillary-dominated, and discontinuous-media regimes without case-specific tuning is a clear strength.
major comments (3)
- [timestep criterion section] § on the new timestep criterion: the CFL bound is constructed from numerically computed velocity derivatives, yet the manuscript provides no description of the finite-difference stencil, any limiting or smoothing, or upwinding applied to these derivatives. This detail is load-bearing for the central robustness claim, especially given the discontinuous-media test.
- [discontinuous media test] Discontinuous porous-media example: the paper states that the expected stationary state is reached despite jumps in porosity, permeability, and capillary-pressure function, but supplies no diagnostic on the sign or magnitude of the computed velocity derivatives near the interfaces. Without such evidence or explicit stabilization, the assertion of parameter-free reliability across all regimes rests on an unexamined assumption.
- [air compression example] Air-compression comparison: the claim of 'significantly improved' non-wetting-phase mass conservation relative to the Coats criterion is central to the efficiency/robustness narrative, yet the manuscript presents only a high-level statement without tabulated mass-error values, iteration counts, or convergence plots that would allow quantitative verification.
minor comments (2)
- [Abstract] The abstract is concise but would benefit from a single sentence summarizing the mathematical form of the new CFL expression.
- [results figures] Figure captions for the Buckley-Leverett and equalization cases should explicitly state the grid resolution and the exact definition of the reported error norms.
Simulated Author's Rebuttal
We thank the referee for the thorough and constructive review. The comments highlight important areas where additional detail and quantitative support will strengthen the manuscript. We address each major comment below and will incorporate the suggested clarifications and data in a revised version.
read point-by-point responses
-
Referee: [timestep criterion section] § on the new timestep criterion: the CFL bound is constructed from numerically computed velocity derivatives, yet the manuscript provides no description of the finite-difference stencil, any limiting or smoothing, or upwinding applied to these derivatives. This detail is load-bearing for the central robustness claim, especially given the discontinuous-media test.
Authors: We agree that the current manuscript lacks sufficient detail on the computation of the numerical velocity derivatives. In the revised version we will add an explicit subsection describing the finite-difference stencil (a compact central difference with upwinding based on the local flow direction), the mild smoothing applied to suppress oscillations at sharp fronts, and the simple limiting procedure used to prevent unphysical negative timesteps. These additions will directly support the robustness claim for the discontinuous-media case. revision: yes
-
Referee: [discontinuous media test] Discontinuous porous-media example: the paper states that the expected stationary state is reached despite jumps in porosity, permeability, and capillary-pressure function, but supplies no diagnostic on the sign or magnitude of the computed velocity derivatives near the interfaces. Without such evidence or explicit stabilization, the assertion of parameter-free reliability across all regimes rests on an unexamined assumption.
Authors: We acknowledge that diagnostics on the velocity derivatives near the material interfaces would provide valuable evidence. The revised manuscript will include additional figures showing the spatial distribution and magnitude of the computed velocity derivatives across the discontinuity, together with a short discussion of how the CFL adaptation itself supplies the necessary stabilization without parameter tuning. revision: yes
-
Referee: [air compression example] Air-compression comparison: the claim of 'significantly improved' non-wetting-phase mass conservation relative to the Coats criterion is central to the efficiency/robustness narrative, yet the manuscript presents only a high-level statement without tabulated mass-error values, iteration counts, or convergence plots that would allow quantitative verification.
Authors: We agree that the quantitative comparison is currently insufficient. In the revision we will add a dedicated table reporting the final non-wetting-phase mass error, average iterations per time step, and total number of time steps for both criteria, along with convergence plots of mass error versus simulation time. These data will allow direct verification of the improvement. revision: yes
Circularity Check
Timestep criterion derived from standard CFL plus numerical derivatives; no reduction to inputs by construction
full rationale
The paper introduces a timestep criterion that directly applies the standard CFL condition using on-the-fly numerically computed velocity derivatives. This construction is independent of the target simulation outcomes or fitted parameters; the reported improvements in iteration count, mass conservation, and robustness on discontinuous-media cases are shown via direct comparison to the Coats criterion rather than being forced by the criterion's definition. No self-citation load-bearing steps, ansatz smuggling, or self-definitional reductions appear in the derivation chain. The approach remains self-contained against external benchmarks such as the Buckley-Leverett problem and known stationary states.
Axiom & Free-Parameter Ledger
axioms (1)
- standard math The Courant-Friedrichs-Lewy condition supplies a sufficient stability limit for explicit or semi-implicit transport steps when velocity derivatives are available.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/ArithmeticFromLogic.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The novel timestep criterion uses the Courant-Friedrichs-Levy (CFL) condition and employs numerically computed velocity derivatives to adapt the timestep size... generalized characteristic wave velocity criterion (Eq. 42)
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
discontinuous porous media parameters such as porosity, permeabilities, and capillary pressure function
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.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.