pith. machine review for the scientific record.
sign in

arxiv: 2512.20846 · v3 · submitted 2025-12-23 · ⚛️ physics.flu-dyn · cond-mat.soft

Droplet Breakup Against an Isolated Obstacle

David J. Meer , Shivnag Sista , Mark D. Shattuck , Corey S. O'Hern , Eric R. Weeks This is my paper

Pith reviewed 2026-05-16 19:40 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn cond-mat.soft
keywords droplet breakupmicrofluidic chambercapillary numberobstacle collisiondiscrete element methodnondimensional number
0
0 comments X

The pith

A nondimensional breakup number Bk based on capillary number predicts whether droplets break or slip past an isolated obstacle.

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

The authors introduce a breakup number Bk that scales with the capillary number to forecast the outcome of droplet collisions with a circular obstacle in a thin microfluidic chamber. Experiments and simulations demonstrate that droplets rarely break for Bk much below 1 but break consistently for Bk much above 1, with the chance of breakup rising sharply near Bk equals 1. They further show that Bk grows with the fourth power of three of a symmetry measure S for the collision geometry. This framework explains why higher flow speeds, larger droplets, lower surface tension, and more head-on impacts all promote breakup, and why thicker chambers increase the breakup rate. If accurate, it offers a parameter-free way to design flows that either preserve or fragment droplets at obstacles.

Core claim

The paper establishes that droplet breakup against an obstacle is governed by a nondimensional breakup number Bk proportional to the capillary number Ca. As Bk increases, the probability of breakup shifts from zero for Bk much less than one to one for Bk much greater than one, with a rapid change near Bk equals one. Additionally, Bk scales with the collision symmetry S raised to the four-thirds power, which indicates that the critical symmetry for breakup is set by a length scale comparable to the obstacle radius R.

What carries the argument

The breakup number Bk, a nondimensional quantity proportional to the capillary number that also incorporates collision symmetry to determine breakup probability.

If this is right

  • Droplets never break when Bk is much less than 1 regardless of other parameters.
  • Droplets always break when Bk is much greater than 1.
  • The breakup probability changes rapidly around Bk equals 1.
  • Bk scales as S to the power 4/3, linking breakup directly to collision symmetry.
  • The minimum symmetry needed for breakup is set by a distance h about equal to the obstacle radius R.

Where Pith is reading between the lines

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

  • Changing the chamber height alters the effective Bk, suggesting height as a control knob for breakup in microfluidic devices.
  • This scaling may extend to other deformable particle systems like cells or bubbles interacting with obstacles.
  • Testing the Bk transition in different fluid viscosities or obstacle shapes could confirm the generality of the model.

Load-bearing premise

The quasi-two-dimensional chamber and deformable particle model in the simulations capture all relevant breakup and slipping physics, including chamber height effects, without extra parameters or three-dimensional corrections.

What would settle it

An experiment in a three-dimensional chamber where breakup probability does not show a sharp transition near Bk equals 1 would disprove the claim that the model suffices without 3D corrections.

Figures

Figures reproduced from arXiv: 2512.20846 by Corey S. O'Hern, David J. Meer, Eric R. Weeks, Mark D. Shattuck, Shivnag Sista.

Figure 2
Figure 2. Figure 2: ) Minimizing rij subject to this constraint on sij prevents unphysical breakup events. Specifically, in the limit dmin-sep → 0, daughter droplets form with an unrealistically small number of vertices. Conversely, as dmin-sep → D0, it becomes increasingly difficult to identify a vertex pair whose shortest connecting arc length exceeds the threshold. Consequently, dmin-sep must fall within a small range. If … view at source ↗
read the original abstract

We describe combined experiments and simulations of droplet breakup during flow-driven interactions with a circular obstacle in a quasi-two-dimensional microfluidic chamber. Due to a lack of in-plane confinement, the droplets can also slip past the obstacle without breaking. Droplets are more likely to break when they have a higher flow velocity, larger size (relative to the obstacle radius R), smaller surface tension, and for head-on collisions with the obstacle. We also observe that droplet-obstacle collisions are more likely to result in breakup when the height of the sample chamber is increased. We define a nondimensional breakup number Bk ~ Ca, where Ca is the Capillary number, that accounts for changes in the likelihood of droplet break up with variations in these parameters. As Bk increases, we find in both experiments and discrete element method (DEM) simulations of the deformable particle model that the behavior changes from droplets never breaking (Bk << 1) to always breaking for Bk >> 1, with a rapid change in the probability of droplet breakup near Bk = 1. We also find that Bk ~ S^(4/3), where S characterizes the symmetry of the collision, which implies that the minimum symmetry required for breakup is controlled by a characteristic distance h ~ R.

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

2 major / 2 minor

Summary. The manuscript reports combined experiments and quasi-2D discrete-element-method (DEM) simulations of droplet breakup against a circular obstacle in a microfluidic chamber. The authors introduce a dimensionless breakup number Bk proportional to the capillary number Ca that incorporates relative droplet size and collision symmetry S; they show that breakup probability transitions from near zero (Bk ≪ 1) to near unity (Bk ≫ 1) around Bk ≈ 1 in both data sets, and report an empirical scaling Bk ~ S^{4/3}.

Significance. If the Bk collapse holds across the reported parameter ranges, including chamber height, the result supplies a compact predictive metric for breakup versus slip in confined droplet-obstacle flows, with direct relevance to microfluidic design and multiphase transport. The reported agreement between experiment and DEM on the Bk dependence is a positive feature, provided the quasi-2D model adequately captures the height-dependent physics.

major comments (2)
  1. [Abstract / Bk definition] Abstract and definition of Bk: Bk is defined proportionally to Ca (and stated to account for changes with velocity, size, tension, and symmetry) yet the text separately reports that breakup probability increases with chamber height h. Because Ca itself contains no h dependence and the DEM is quasi-2D, the manuscript must show either that Bk contains an explicit h factor or that data for multiple h values still collapse onto the same Bk = 1 transition; otherwise the claimed universality is not demonstrated.
  2. [DEM simulations section] DEM model and comparison with experiment: The deformable-particle DEM is quasi-two-dimensional and therefore omits gap-flow lubrication and meniscus curvature that scale with h. If the simulated transition remains fixed at Bk ≈ 1 while the experimental 50 % breakup threshold shifts with h, the quantitative match asserted in the abstract cannot hold; the paper should either introduce an effective-height correction or restrict the comparison to fixed h.
minor comments (2)
  1. [Methods / Bk definition] The proportionality constant relating Bk to Ca (including any dependence on relative size) should be stated explicitly, together with the precise definition of the symmetry parameter S.
  2. [Results / Figures] Error bars, number of realizations, and fitting procedure for the probability-versus-Bk curves and the S^{4/3} scaling should be reported so that the sharpness of the Bk = 1 transition can be assessed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which correctly highlight the need to address the chamber-height dependence in both the Bk definition and the DEM-experiment comparison. We respond point by point below and will revise the manuscript to incorporate an explicit h dependence, demonstrate data collapse, and clarify model limitations.

read point-by-point responses

  1. Referee: [Abstract / Bk definition] Abstract and definition of Bk: Bk is defined proportionally to Ca (and stated to account for changes with velocity, size, tension, and symmetry) yet the text separately reports that breakup probability increases with chamber height h. Because Ca itself contains no h dependence and the DEM is quasi-2D, the manuscript must show either that Bk contains an explicit h factor or that data for multiple h values still collapse onto the same Bk = 1 transition; otherwise the claimed universality is not demonstrated.

    Authors: We agree that the present Bk definition omits explicit h dependence despite the reported increase in breakup probability with chamber height. In the revised manuscript we will augment Bk with a multiplicative factor (h/R)^α, where the exponent α is chosen to be consistent with the observed S^{4/3} scaling (which already implies a characteristic length comparable to h). We will then show that experimental datasets obtained at several distinct h values collapse onto a single transition near Bk ≈ 1. The abstract will be updated to state that Bk now incorporates the chamber-height effect. revision: yes

  2. Referee: [DEM simulations section] DEM model and comparison with experiment: The deformable-particle DEM is quasi-two-dimensional and therefore omits gap-flow lubrication and meniscus curvature that scale with h. If the simulated transition remains fixed at Bk ≈ 1 while the experimental 50 % breakup threshold shifts with h, the quantitative match asserted in the abstract cannot hold; the paper should either introduce an effective-height correction or restrict the comparison to fixed h.

    Authors: The referee is correct that the quasi-2D DEM omits lubrication and meniscus effects that vary with h. We will therefore introduce an effective-height correction into the Bk used for the DEM runs, calibrated directly from the experimental h dependence. With this correction the simulated breakup transition remains at Bk ≈ 1 and matches the experimental threshold across the studied h range. A new paragraph in the DEM section will explicitly discuss the quasi-2D approximation and the rationale for the correction, while the abstract claim of quantitative agreement will be qualified accordingly. revision: yes

Circularity Check

1 steps flagged

Bk defined to collapse parameters; transition near 1 is normalized rather than independently predicted

specific steps
  1. fitted input called prediction [Abstract]
    "We define a nondimensional breakup number Bk ~ Ca, where Ca is the Capillary number, that accounts for changes in the likelihood of droplet break up with variations in these parameters. As Bk increases, we find in both experiments and discrete element method (DEM) simulations of the deformable particle model that the behavior changes from droplets never breaking (Bk << 1) to always breaking for Bk >> 1, with a rapid change in the probability of droplet breakup near Bk = 1."

    Bk is explicitly constructed to account for the listed parameters (velocity via Ca, relative size, tension, symmetry S). Normalizing the group so the observed transition occurs near unity makes the 'rapid change near Bk=1' a direct consequence of the fitting constants chosen for collapse, rather than a parameter-free prediction.

full rationale

The paper constructs Bk ~ Ca to incorporate velocity, size, tension, and symmetry effects so that breakup probability data collapses. The reported rapid change near Bk=1 follows from this normalization choice, matching pattern 2. However, the specific S^(4/3) scaling and DEM-experiment agreement retain independent empirical content from the raw data, preventing full reduction to tautology. No self-citation load-bearing or self-definitional steps appear in the provided text.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the definition of Bk as proportional to Ca and the empirical observation of the S^(4/3) scaling. No new physical entities are introduced. The quasi-2D approximation and the validity of the DEM model for breakup are taken as given.

free parameters (1)
  • proportionality constant in Bk ~ Ca
    The exact prefactor relating Bk to Ca is not numerically specified and is treated as order-one.
axioms (2)
  • domain assumption Quasi-two-dimensional flow in the microfluidic chamber adequately represents the breakup dynamics
    The abstract notes lack of in-plane confinement and the effect of chamber height, but treats the 2D model as sufficient.
  • domain assumption The discrete element method deformable particle model reproduces the experimental breakup behavior
    Agreement between experiments and simulations is stated without detailed validation metrics.

pith-pipeline@v0.9.0 · 5539 in / 1664 out tokens · 69859 ms · 2026-05-16T19:40:57.461866+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

6 extracted references · 6 canonical work pages

  1. [1]

    Bibette, L

    1 J. Bibette, L. Calderon and P. Poulin,Rep. Prog. Phys., 1999, 62, 969–1033. 2 W. C. K. Poon, A. T. Brown, S. O. L. Direito, D. J. M. Hodgson, L. L. Nagard, A. Lips, C. E. MacPhee, D. Marenduzzo, J. R. Royer, A. F. Silva, J. H. J. Thijssen and S. Titmuss,Soft Matter, 2020,16, 8310–8324. 3 P. Tabeling,Introduction to Microfluidics, Oxford University Press...

    work page 1999

  2. [2]

    Selimovic, F

    4 S. Selimovic, F. Gobeaux and S. Fraden,Lab Chip, 2010,10, 1696–1699. 5 J. J. M. Janssen, A. Boon and W. G. M. Agterof,Fluid Mech. Transp. Phenom., 1994,40, 1929–1939. 6 D. R. Link, S. L. Anna, D. A. Weitz and H. A. Stone,Physical Review Letters, 2004,92, 054503. 7 T. Fu, Y. Ma, D. Funfschilling and H. Z. Li,Chem. Eng. Sci., 2011,66, 4184–4195. 8 W. L. C...

    work page 2010

  3. [3]

    9 D. H. Wankawala and P. K. Mondal,S ¯adhan¯a Acad. Proc. Eng. Sci., 2023,48, 274–274. 10 D. H. Wankawala, A. Q. Shen and P. K. Mondal,Proc. R. Soc. A, 2025,481, 20240594. 11 A. Tiribocchi, M. Durve, M. Lauricella, A. Montessori, D. Marenduzzo and S. Succi,Phys. Fluids, 2023,35, 063106. 12 S. S. Dhruvkumar H. Wankawala and P. K. Mondal, Europhysics Letter...

    work page 2023

  4. [4]

    Chen and B

    20 S. Chen and B. Guo,Water Resour. Res., 2023,59, e2023WR034664. 21 J.-W. Liu, K.-H. Wei, S.-W. Xu, J. Cui, J. Ma, X.-L. Xiao, B.-D. Xi and X.-S. He,Sci. Total Environ., 2008,756, 144142. 22 S. S. Datta, H. Chiang, T. S. Ramakrishnan and D. A. Weitz, Phys. Rev. Lett., 2013,111, 064501. 23 S. S. Datta and D. A. Weitz,Europhys. Lett., 2013,101, 14002. 24 S...

    work page 2023

  5. [5]

    Shukla, N

    37 I. Shukla, N. Kofman, G. Balestra, L. Zhu and F. Gallaire, Journal of Fluid Mechanics, 2019,874, 1021–1040. 38 F. P. Bretherton,Journal of Fluid Mechanics, 1961,10, 166–188. 39 G. TAYLOR and P. G. SAFFMAN,The Quarterly Journal of Mechanics and Applied Mathematics, 1959,12, 265–279. 40 K. W. Desmond, P. J. Young, D. Chen and E. R. Weeks,Soft Matter, 201...

    work page 2019

  6. [6]

    42 X. Shi, M. P. Brenner and S. R. Nagel,Science, 1994,265, 219–222. 43 N. Bremond, A. R. Thiam and J. Bibette,Phys. Rev. Lett., 2008, 100, 024501. 44 A. Lai, N. Bremond and H. A. Stone,J. Fluid Mech., 2009,632, 97–107. 45 A. G. Yiotis, L. Talong and D. Salin,Phys. Rev. E, 2013,87, 033001. 46 E. Amstad, S. S. Datta and D. A. Weitz,Lab Chip, 2014,14, 705–7...

    work page 1994