pith. sign in

arxiv: 2509.08034 · v1 · submitted 2025-09-09 · ⚛️ physics.flu-dyn

Towards enhanced mixing of a high viscous miscible blob in porous media

Mijanur Rahaman , Jiten C. Kalita , Satyajit Pramanik This is my paper

Pith reviewed 2026-05-18 17:32 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords miscible displacementviscous fingeringporous mediamixing enhancementPeclet numbermobility ratioblob deformationnumerical simulation
0
0 comments X

The pith

A viscous circular blob in porous media mixes most at intermediate flow speeds and viscosity contrasts due to its initial shape.

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

The paper uses high-order numerical simulations to study how a highly viscous miscible circular blob deforms and spreads when displaced by a less viscous fluid in a homogeneous porous medium with no-flux boundaries. It identifies three main patterns—comet shapes, lump shapes, and viscous fingers—that emerge as the Peclet number and log-mobility ratio change. The key finding is that deformation, spreading, and mixing do not increase steadily with higher speeds or larger viscosity differences; instead they peak at moderate values because the blob begins with curvature that shapes how the interface evolves over time. This matters for applications where controlling how fluids mix inside rock affects efficiency or safety.

Core claim

Simulations for Pe up to 3000 and R up to 7 reveal comet-shape, lump-shape, and viscous fingering patterns. The deformation, spreading, and mixing of the blob vary non-ideally with both Pe and R because of the blob's initial curvature, so enhanced mixing occurs at intermediate values of Pe and R and points to an optimal mixing condition.

What carries the argument

The initial curvature of the circular blob, which produces the non-monotonic dependence of mixing efficiency on Peclet number and log-mobility ratio.

If this is right

  • Mixing efficiency can be increased by operating at moderate rather than extreme flow rates or viscosity contrasts.
  • The blob transitions through comet, lump, and fingering shapes as Pe and R rise.
  • Results apply directly to tuning displacement processes in oil recovery, CO2 sequestration, and contaminant cleanup.

Where Pith is reading between the lines

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

  • Non-circular starting shapes could move or remove the optimal mixing window entirely.
  • Heterogeneous real-world porous media might shift the intermediate values that give best mixing.
  • Targeted lab experiments with curved initial blobs at controlled Pe and R could confirm the curvature mechanism.

Load-bearing premise

The non-monotonic mixing behavior is caused by the blob's initial curvature rather than by details of the numerical scheme, domain size, or boundary conditions.

What would settle it

Running the same displacement with an initially square blob instead of a circular one and finding that mixing no longer peaks at intermediate Pe and R would undermine the curvature explanation.

Figures

Figures reproduced from arXiv: 2509.08034 by Jiten C. Kalita, Mijanur Rahaman, Satyajit Pramanik.

Figure 1
Figure 1. Figure 1: FIG. 1: Schematic of the displacement of a circular blob in a two-dimensional homogeneous [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Rescaled longitudinal- and transverse-averaged concentration, [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Spatio-temporal distribution of a circular blob at different times, [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Phase plot in the [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Spatial distribution of a circular blob of radius [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Spatial distribution of a circular blob of radius [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Spatial distribution of a circular blob of radius [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Spatial distributions of a circular blob of radius [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Effects of log-mobility ratio on the evolution of variances rescaled with their initial value [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Effects of log-mobility ratio on the evolution of (a) area covered by the sample rescaled [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Variation of (a) [PITH_FULL_IMAGE:figures/full_fig_p023_11.png] view at source ↗
read the original abstract

In this study, we investigate the rectilinear displacement and deformation of a highly viscous, miscible circular blob influenced by a less viscous fluid within a homogeneous porous medium featuring physically realistic no-flux boundaries. We utilize a fourth-order accurate compact finite difference scheme for the spatial discretization of the nonlinear partial differential equations that govern this phenomenon. The resulting semi-discrete equations are then integrated using the second-order Crank-Nicolson (CN) method. We conduct numerical simulations for a P\'eclet number ($Pe \leq 3000$) and a log-mobility ratio $0 \leq R \leq 7$, which reveal three distinct pattern formations: comet-shape, lump-shape, and viscous fingering instability. Our results demonstrate that the deformation, spreading, and mixing of the blob vary non-ideally with both $Pe$ and $R$, a behavior attributed to the blob's initial curvature. Consequently, enhanced mixing can be achieved at intermediate values of $Pe$ and $R$, suggesting the existence of an optimal mixing condition. These findings have significant implications for fields such as oil recovery, CO$_2$ sequestration, pollution remediation, and chromatography separation.

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

3 major / 2 minor

Summary. The manuscript numerically studies the rectilinear displacement of a high-viscosity miscible circular blob by a less viscous fluid in a homogeneous porous medium with no-flux boundaries. Employing a fourth-order compact finite-difference spatial discretization and second-order Crank-Nicolson time integration, the authors simulate the governing PDEs over Pe ≤ 3000 and log-mobility ratio 0 ≤ R ≤ 7. They identify three pattern regimes (comet-shape, lump-shape, viscous fingering) and report that blob deformation, spreading, and mixing exhibit non-monotonic dependence on both Pe and R, which they attribute to the initial curvature; this leads to the claim of enhanced mixing at intermediate parameter values.

Significance. If the reported non-monotonic trends prove robust, the work would usefully demonstrate how initial geometry can produce an optimal mixing window in miscible porous-media flows, with direct relevance to oil recovery, CO2 sequestration, and remediation. The direct numerical approach (standard scheme, no fitted parameters, systematic variation of Pe and R) is a methodological strength that avoids circularity.

major comments (3)
  1. [Numerical Methods] Numerical Methods section: no grid-convergence study, resolution tests, or domain-size sensitivity analysis is reported. At Pe = 3000 sharp fronts form; without these checks the non-monotonic dependence on Pe and R (and the claimed optimal mixing regime) cannot be confidently distinguished from numerical artifacts or boundary-induced recirculation.
  2. [Results] Results section: the central attribution of non-ideal mixing behavior to the blob’s initial curvature is not supported by a control simulation using a flat initial interface. Without this comparison it remains unclear whether the observed pattern transitions and mixing enhancement arise from curvature or from other factors such as the no-flux boundaries or the chosen mobility-ratio range.
  3. [Abstract and Results] Abstract and Results: the existence of an “optimal mixing condition” at intermediate Pe and R is asserted without quantitative mixing metrics (e.g., mixing time, scalar dissipation rate, or efficiency with uncertainty estimates). The non-monotonic claim therefore lacks a clear, reproducible definition that would allow readers to verify the optimum.
minor comments (2)
  1. [Figures and Results] Figure captions and text should explicitly define the quantitative measures used to quantify “mixing” and “spreading” so that the non-ideal dependence can be assessed independently.
  2. [Numerical Methods] The manuscript would benefit from a brief statement of the computational domain aspect ratio and any tests confirming that no-flux boundaries do not artificially influence the long-time blob evolution.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below, indicating where we agree and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Numerical Methods] Numerical Methods section: no grid-convergence study, resolution tests, or domain-size sensitivity analysis is reported. At Pe = 3000 sharp fronts form; without these checks the non-monotonic dependence on Pe and R (and the claimed optimal mixing regime) cannot be confidently distinguished from numerical artifacts or boundary-induced recirculation.

    Authors: We agree that explicit numerical verification is essential, especially at high Péclet numbers where sharp fronts develop. Although our fourth-order compact finite-difference scheme with Crank-Nicolson time integration is designed for accuracy, we did not report a dedicated grid-convergence or domain-size study in the original manuscript. In the revised version we will add resolution tests for representative high-Pe cases (including Pe = 3000) together with domain-size sensitivity checks to confirm that the reported non-monotonic trends and pattern regimes are free of numerical artifacts or boundary effects. revision: yes

  2. Referee: [Results] Results section: the central attribution of non-ideal mixing behavior to the blob’s initial curvature is not supported by a control simulation using a flat initial interface. Without this comparison it remains unclear whether the observed pattern transitions and mixing enhancement arise from curvature or from other factors such as the no-flux boundaries or the chosen mobility-ratio range.

    Authors: The comet-shape and lump-shape regimes arise directly from the interaction of the curved interface with the rectilinear flow; these morphologies are not observed in conventional flat-interface displacements. Nevertheless, we recognize that an explicit flat-interface control simulation would isolate the curvature effect more rigorously and rule out influences from the no-flux boundaries or the R range. We will therefore perform and include such a control comparison in the revised manuscript. revision: yes

  3. Referee: [Abstract and Results] Abstract and Results: the existence of an “optimal mixing condition” at intermediate Pe and R is asserted without quantitative mixing metrics (e.g., mixing time, scalar dissipation rate, or efficiency with uncertainty estimates). The non-monotonic claim therefore lacks a clear, reproducible definition that would allow readers to verify the optimum.

    Authors: Mixing is quantified in the results via the decay of concentration variance and the growth of iso-concentration contour length, both of which display non-monotonic dependence on Pe and R; the optimal regime is identified where the variance decay rate is largest. To improve clarity and reproducibility we will revise the abstract and results sections to state these metrics explicitly, add scalar-dissipation-rate diagnostics, and include uncertainty estimates derived from the simulations. revision: partial

Circularity Check

0 steps flagged

Numerical simulations of governing PDEs with independent input parameters and no self-referential reductions

full rationale

The paper solves the standard nonlinear PDEs for miscible displacement in porous media via a fourth-order compact finite-difference spatial discretization and second-order Crank-Nicolson time stepping. Pe and R are supplied as independent control parameters that are systematically varied; the resulting pattern transitions and non-monotonic mixing trends are direct outputs of these integrations rather than quantities fitted or redefined from the same data. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior work by the same authors appear in the derivation chain. The attribution of non-ideal behavior to initial curvature is an interpretive remark on the simulation results, not a definitional equivalence or statistical forcing. The study is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard continuum description of miscible displacement in porous media together with the numerical discretization; no new physical constants or entities are introduced.

axioms (2)
  • domain assumption The flow obeys Darcy's law and the concentration field satisfies a convection-diffusion equation with constant dispersion coefficients.
    Invoked when setting up the governing PDEs for the rectilinear displacement problem.
  • domain assumption The porous medium is homogeneous and isotropic with no-flux boundaries on the lateral walls.
    Stated in the problem description and used to close the computational domain.

pith-pipeline@v0.9.0 · 5746 in / 1438 out tokens · 31097 ms · 2026-05-18T17:32:34.032728+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.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We utilize a fourth-order accurate compact finite difference scheme for the spatial discretization of the nonlinear partial differential equations... Crank-Nicolson (CN) method... Pe ≤ 3000 and 0 ≤ R ≤ 7... three distinct pattern formations: comet-shape, lump-shape, and viscous fingering instability.

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.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. HOC simulations of miscible viscous fingering of a finite slice: A new insight

    physics.flu-dyn 2026-04 unverdicted novelty 4.0

    Simulations show permeable transverse boundaries increase solute mass in miscible viscous fingering, driving stronger instabilities and larger mixing lengths than other boundaries.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages · cited by 1 Pith paper

  1. [1]

    While the pseudo-spectral method is known for its accuracy, it necessitates periodic boundary conditions

    implemented an alternating-direction implicit (ADI) technique combined with a Hartley-based pseudo-spectral method for similar computations. While the pseudo-spectral method is known for its accuracy, it necessitates periodic boundary conditions. A higher-order compact (HOC) finite difference method offers a way to overcome the restriction of boundary con...

    work page

  2. [2]

    In the absence of viscosity contrast between the displacing and displaced fluids, the blob is carried along with the displacing fluid. Thus, in the moving frame of reference, the dynamics of the circular blob can be described by solving a 2D diffusion equation subjected to no-flux boundary conditions and an appropriate initial condition. To better underst...

    work page 2000

  3. [3]

    G. M. Homsy, Viscous fingering in porous media, Annual Review of Fluid Mechanics19, 271 (1987)

    work page 1987

  4. [4]

    H. E. Huppert and J. A. Neufeld, The fluid mechanics of carbon dioxide sequestration, Annual Review of Fluid Mechanics46, 255 (2014)

    work page 2014

  5. [5]

    De Wit, Y

    A. De Wit, Y . Bertho, and M. Martin, Viscous fingering of miscible slices, Physics of Fluids17(2005). 24

    work page 2005

  6. [6]

    P. G. Saffman and G. I. Taylor, The penetration of a fluid into a porous medium or hele-shaw cell containing a more viscous liquid, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences245, 312 (1958)

    work page 1958

  7. [7]

    A. Riaz, M. Hesse, H. A. Tchelepi, and F. M. Orr, Onset of convection in a gravitationally unstable diffusive boundary layer in porous media, Journal of Fluid Mechanics548, 87 (2006)

    work page 2006

  8. [8]

    R. Maes, G. Rousseaux, B. Scheid, M. Mishra, P. Colinet, and A. De Wit, Experimental study of dispersion and miscible viscous fingering of initially circular samples in hele-shaw cells, Physics of Fluids22(2010)

    work page 2010

  9. [9]

    Mishra, M

    M. Mishra, M. Martin, and A. De Wit, Differences in miscible viscous fingering of finite width slices with positive or negative log-mobility ratio, Physical Review E78, 066306 (2008)

    work page 2008

  10. [10]

    Pramanik and M

    S. Pramanik and M. Mishra, Fingering instability and mixing of a blob in porous media, Physical Review E94, 043106 (2016)

    work page 2016

  11. [11]

    Pramanik, A

    S. Pramanik, A. De Wit, and M. Mishra, Viscous fingering and deformation of a miscible circular blob in a rectilinear displacement in porous media, Journal of Fluid Mechanics782, R2 (2015)

    work page 2015

  12. [12]

    C. T. Tan and G. M. Homsy, Simulation of nonlinear viscous fingering in miscible displacement, Physics of Fluids31, 1330 (1988)

    work page 1988

  13. [13]

    Li and B

    J. Li and B. Rivière, High order discontinuous galerkin method for simulating miscible flooding in porous media, Computational Geosciences19, 1251 (2015)

    work page 2015

  14. [14]

    Li and B

    J. Li and B. Rivière, Numerical modeling of miscible viscous fingering instabilities by high-order methods, Transport in Porous Media113, 607 (2016)

    work page 2016

  15. [15]

    J. J. Hidalgo, J. Fe, L. C. Felgueroso, and R. Juanes, Scaling of convective mixing in porous media, Physical Review Letters109, 264503 (2012)

    work page 2012

  16. [16]

    M. N. Islam and J. Azaiez, Fully implicit finite difference pseudo-spectral method for simulating high mobility-ratio miscible displacements, International Journal for Numerical Methods in Fluids47, 161 (2005)

    work page 2005

  17. [17]

    M. M. Gupta, R. M. Manohar, and J. H. Stephenson, A single cell high order scheme for the convection-diffusion equation with variable coefficients, International Journal for Numerical Methods in Fluids4, 641 (1984)

    work page 1984

  18. [18]

    S. C. R. Dennis and J. D. Hudson, Compact h4 finite-difference approximations to operators of Navier- Stokes type, Journal of Computational Physics85, 390 (1989). 25

    work page 1989

  19. [19]

    E. C. Gartland, Jr., Discrete weighted mean approximation of a model convection-diffusion equation, SIAM Journal on Scientific and Statistical Computing3, 460 (1982)

    work page 1982

  20. [20]

    B. J. Noye and H. H. Tan, A third order semi-implicit method for solving the one-dimensional convection–diffusion equation, International Journal for Numerical Methods in Engineering26, 1615 (1988)

    work page 1988

  21. [21]

    B. J. Noye and H. H. Tan, Finite difference methods for solving the two-dimensional advection– diffusion equation, International Journal for Numerical Methods in Fluids9, 75 (1989)

    work page 1989

  22. [22]

    M. Li, T. Tang, and B. Fornberg, A compact fourth-order finite difference scheme for the steady in- compressible navier-stokes equations, International Journal for Numerical Methods in Fluids20, 1137 (1995)

    work page 1995

  23. [23]

    N. S. Wilkes and C. P. Thompson, An evaluation of higher-order upwind differencing for elliptic flow problems, Numerical Methods in Laminar and Turbulent Flow , 248 (1983)

    work page 1983

  24. [24]

    Yanwen, F

    M. Yanwen, F. Dexun, T. Kobayashi, and N. Taniguchi, Numerical solution of the incompressible navier–stokes equations with an upwind compact difference scheme, International Journal for Numer- ical Methods in Fluids30, 509 (1999)

    work page 1999

  25. [25]

    Sesterhenn, A characteristic-type formulation of the navier–stokes equations for high order upwind schemes, Computers & fluids30, 37 (2000)

    J. Sesterhenn, A characteristic-type formulation of the navier–stokes equations for high order upwind schemes, Computers & fluids30, 37 (2000)

    work page 2000

  26. [26]

    Lax and B

    P. Lax and B. Wendroff,Systems of conservation laws, Tech. Rep. (1959)

    work page 1959

  27. [27]

    R. H. Pletcher, J. C. Tannehill, and D. Anderson,Computational fluid mechanics and heat transfer (CRC press, 2012)

    work page 2012

  28. [28]

    J. D. Anderson,Computational fluid dynamics: the basics with applications(McGraw-Hill New York, 2002)

    work page 2002

  29. [29]

    R. J. MacKinnon and G. F. Carey, Analysis of material interface discontinuities and superconvergent fluxes in finite difference theory, Journal of Computational Physics75, 151 (1988)

    work page 1988

  30. [30]

    W. F. Spotz and G. F. Carey,High-order compact finite difference methods with applications to viscous flows(Texas Institute for Computational and Applied Mathematics, University of Texas, Austin, 1994)

    work page 1994

  31. [31]

    W. F. Spotz and G. F. Carey, High-order compact scheme for the steady stream-function vorticity equations, International Journal for Numerical Methods in Engineering38, 3497 (1995)

    work page 1995

  32. [32]

    J. C. Kalita, D. C. Dalal, and A. K. Dass, A class of higher order compact schemes for the unsteady two- dimensional convection–diffusion equation with variable convection coefficients, International Journal 26 for Numerical Methods in Fluids38, 1111 (2002)

    work page 2002

  33. [33]

    J. C. Kalita,HOC Schemes for Incompressible Viscous Flows: Application and Development, Ph.D. thesis (2001)

    work page 2001

  34. [34]

    M. M. Gupta and J. C. Kalita, A new paradigm for solving Navier–Stokes equations: streamfunction– velocity formulation, Journal of Computational Physics207, 52 (2005)

    work page 2005

  35. [35]

    Nishida, A

    A. Nishida, A. Fujii, and Y . Oyanagi, Lis: Library of iterative solvers for linear systems (2022)., URL: http://www.ssisc.org/lis/

    work page 2022

  36. [36]

    Sharma, S

    V . Sharma, S. Pramanik, and M. Mishra, Dynamics of a highly viscous circular blob in homogeneous porous media, Fluids2, 32 (2017)

    work page 2017

  37. [37]

    Kretz, P

    V . Kretz, P. Berest, J. P. Hulin, and D. Salin, An experimental study of the effects of density and viscosity contrasts on macrodispersion in porous media, Water Resources Research39(2003)

    work page 2003

  38. [38]

    J. C. Bacri, D. Salin, and R. Woumeni, Three-dimensional miscible viscous fingering in porous media, Physical Review Letters67, 2005 (1991). 27

    work page 2005