Inertial focusing of neutrally buoyant spherical particle in shallow microchannels
Pith reviewed 2026-05-07 15:09 UTC · model grok-4.3
The pith
An explicit formula predicts the inertial lift force on neutrally buoyant spherical particles up to 35 percent of channel height at low Reynolds numbers.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using immersed boundary method simulations of 2D planar Poiseuille flow, the study derives an explicit formula for the lift force that accurately predicts the force for particle radius to height ratios from 0.03 to 0.35 and particle Reynolds numbers up to 1. The formula remains valid even as near-wall lift forces decrease at higher Reynolds numbers within this range. Predictions of particle trajectories using this model match published experimental observations.
What carries the argument
The explicit formula for the lift force obtained from immersed boundary method calculations in planar Poiseuille flow.
If this is right
- The formula enables direct estimation of particle migration paths in microfluidic devices without repeated simulations.
- Increased wall slip length reduces near-wall lift force and shifts particle equilibrium positions closer to the walls.
- The model applies to particles up to 35 percent of channel height, expanding the usable size range beyond earlier limits.
- Trajectory predictions from the derived model align with experimental data on particle migration.
Where Pith is reading between the lines
- If the formula holds in three-dimensional rectangular channels, it could simplify design of inertial focusing devices for separating particles by size.
- The observed drop in near-wall lift at higher Reynolds numbers within the stated range indicates that inertial corrections become necessary even at nominally low Reynolds numbers.
- The same simulation approach could be tested on non-spherical particles or different flow profiles to check whether a similar explicit formula emerges.
Load-bearing premise
The two-dimensional immersed boundary method simulations accurately represent the three-dimensional lift forces experienced by particles in actual microchannels.
What would settle it
Experimental measurements of equilibrium positions or lift forces for particles with radius-to-height ratio 0.35 at particle Reynolds number 1 that differ substantially from the formula's predictions.
Figures
read the original abstract
This study investigates the lift force acting on a finite-size, neutrally buoyant spherical particle suspended in a liquid while flowing through a shallow channel at low Reynolds numbers. Using an immersed boundary method, we calculate the lift force for particle radius-to-channel height ratios spanning \(0.03 \leq a/H \leq 0.35\) in 2D planar Poiseuille flows. We propose an explicit formula that accurately predicts the lift force for particles as large as \(a/H = 0.35\) and remains valid for particle Reynolds number \(Re_p \leq 1\), despite a reduction in near-wall lift force at higher \(Re_p\). The influence of slip boundary conditions is also explored, demonstrating that increased slip length reduces near-wall lift force and shifts the particle equilibrium position closer to the wall. Predictions of the particle trajectory from the derived model are in good agreement to the published experimental data. These findings offer a practical framework for estimating the migration of large particles in microfluidic devices.(This article has been accepted for publication in Physics of Fluids. After publication, it will be available via the AIP Publishing website.)
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses 2D immersed-boundary simulations of planar Poiseuille flow to compute the lift force on neutrally buoyant spheres for 0.03 ≤ a/H ≤ 0.35. It proposes an explicit formula asserted to remain accurate up to a/H = 0.35 and Re_p ≤ 1 (despite near-wall lift reduction at higher Re_p), examines the effect of slip boundary conditions on equilibrium position, and reports that integrated trajectories agree with published experimental data.
Significance. If the 2D reduction is shown to be quantitatively reliable, the explicit formula would provide a practical, ready-to-use expression for inertial migration of large particles in shallow microfluidic channels. The reported experimental agreement is a strength, but the absence of 3D benchmarks leaves the central predictive claim dependent on an unverified dimensional approximation.
major comments (2)
- [Simulation method and results sections] The central claim rests on 2D IBM results in planar Poiseuille flow being representative of 3D spherical particles. For a/H = 0.35 the particle nearly spans the channel height; 3D flow around the equator will modify both pressure and viscous contributions to lift. No comparison to 3D simulations, to the small-a/H asymptotic lift expressions of Ho & Leal or Schonberg & Hinch, or to any 3D benchmark is described, making the validity range a/H ≤ 0.35 load-bearing and unverified.
- [Derivation of the explicit lift-force formula] The explicit formula is stated to be 'derived from simulations.' The manuscript must clarify whether its coefficients were obtained by least-squares fit to the same IBM data later used for validation; if so, the formula is a correlation rather than an independent prediction, and the reported accuracy must be re-evaluated on held-out cases or independent 3D data.
minor comments (2)
- [Abstract and validation section] The abstract states 'good agreement' with experiments; quantitative metrics (RMS error in equilibrium position or lift coefficient) should be reported in the main text or a table.
- [Results] Notation for the proposed formula (e.g., definition of all non-dimensional groups and the precise functional form) should be collected in one place for reproducibility.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive review. We address the two major comments point by point below. Revisions have been made to clarify the nature of the explicit formula and to discuss the limitations of the 2D approximation, while acknowledging that direct 3D validation lies outside the present study.
read point-by-point responses
-
Referee: The central claim rests on 2D IBM results in planar Poiseuille flow being representative of 3D spherical particles. For a/H = 0.35 the particle nearly spans the channel height; 3D flow around the equator will modify both pressure and viscous contributions to lift. No comparison to 3D simulations, to the small-a/H asymptotic lift expressions of Ho & Leal or Schonberg & Hinch, or to any 3D benchmark is described, making the validity range a/H ≤ 0.35 load-bearing and unverified.
Authors: We agree that the 2D planar simulations cannot fully capture three-dimensional flow features around the particle equator, particularly when a/H approaches 0.35. The 2D model is motivated by the shallow-channel limit (height ≪ width) in which the mid-plane flow is approximately planar; however, we recognize that this approximation becomes less accurate for large particles. In the revised manuscript we have added an explicit limitations paragraph that (i) cites the Ho & Leal and Schonberg & Hinch asymptotics for small a/H, (ii) notes the absence of 3D benchmarks, and (iii) qualifies the stated validity range accordingly. The experimental trajectory agreement remains the primary external check, but we do not claim it substitutes for 3D lift-force validation. revision: partial
-
Referee: The explicit formula is stated to be 'derived from simulations.' The manuscript must clarify whether its coefficients were obtained by least-squares fit to the same IBM data later used for validation; if so, the formula is a correlation rather than an independent prediction, and the reported accuracy must be re-evaluated on held-out cases or independent 3D data.
Authors: The coefficients were obtained by a least-squares fit to the full set of 2D IBM lift-force data. We have revised the text to state explicitly that the formula is an empirical correlation fitted to the simulated cases. No held-out subset was reserved; instead, the independent validation we provide is the integration of the fitted force into particle trajectories that are compared against published experimental data. We have added a sentence noting that the accuracy figures therefore reflect in-sample performance of the correlation within the 2D parameter space. revision: yes
- Direct three-dimensional simulation benchmarks or comparisons to 3D asymptotic lift expressions for a/H up to 0.35
Circularity Check
No significant circularity; derivation relies on external experimental validation
full rationale
The paper computes lift forces via 2D immersed-boundary simulations across the stated a/H range, constructs an explicit formula from those results, and validates resulting trajectory predictions against independent published experimental data. No load-bearing self-citations, self-definitional reductions, or fitted inputs called predictions on the same dataset appear in the abstract or described chain. The central claim is therefore self-contained against an external benchmark rather than tautological.
Axiom & Free-Parameter Ledger
free parameters (1)
- coefficients in the proposed explicit lift-force formula
axioms (2)
- domain assumption The flow is two-dimensional planar Poiseuille flow.
- domain assumption The particle is neutrally buoyant, rigid, and spherical.
Reference graph
Works this paper leans on
-
[1]
experimentally finds the particle is force-free at zeq = 0.25H (zeq is the equilibrium location from the particle center to the tube wall) for a/H = 0.15 and shifts towards the tube cente as a/H increases to 0 .39; Staben et al. [40] observed in numerical simulations that particle streamwi se slip velocity increases with the particle to channel size ra tio...
-
[2]
experimentally finds the particle force-free location at zeq/H ∼ 0.3 for a/H = 0.05 and shifts towards the center as a/H increases to 0.45 in a square duct. Additionally, through numerical sim- ulations, they find the lift-force is proportional to a3 far from the wall while it is proportional to a6 in the near-wall region for 0.22 ≤ a/H ≤ 0.38. To clarify t...
work page 1962
-
[3]
combine the coefficients obtained by V asseur and Cox [ 15] and Cherukat and McLaughlin
-
[4]
to generalized expressions for the lift force of small part icles (a/H ≤ 0.05) and finite Re ≤
-
[5]
[37] propose a modification of expressions for the coefficients reported by Asmolov et al
For larger particles ( a/H ≤ 0.15), Nizkaya et al. [37] propose a modification of expressions for the coefficients reported by Asmolov et al. [35] by fitting to their own simulation data. The key advantage of the models proposed by Asmolov et al. [35] and Nizkaya et al. [37] is their simplicity and accuracy in predicting the lift force, makin g them well-sui...
-
[6]
whereas about 30 percents smaller than simulation result o f Asmolov et al. [35]. We note that the lattice Boltzmann method (LBM) is employed in both A smolov et al. [35] and Nizkaya et al. [37], in which there are about 80 lattices placed in the wall norm al direction corresponding to 1 ∼ 2 lattices in the nearest gap between the sphere and wall. Whe n c...
-
[7]
[37] (red dashed curve), and Hood et al
(blue dash-dotted curve), Nizkaya et al. [37] (red dashed curve), and Hood et al. [43] (magenta circles), along with LBM results from Asmolov et al. [35] (blue plus symbols) and Nizkaya et al. [37] (red squares). Insets provide magnified views highlightin g local differences between numerical and model results, particularly where cL is small. tions (zp) is...
-
[8]
extended the parameter space to higher Reynolds numbers ( Re = O(1) and Re = O(103), respectively), and both studies retained the same a4 scaling as in Ho and Leal [ 14]. Fig. 5a presents cL representing the a4 scaling of the lift force FL, based on numerical sim- ulations for particle sizes ranging from a/H = 0.03 to 0 .35 and the model prediction for sm...
-
[9]
[39] corresponds to a tube, and V asseur and Cox [15] is for a/H ≪ 1
correspond to a square channel on symmetry planes, Karnis et al. [39] corresponds to a tube, and V asseur and Cox [15] is for a/H ≪ 1. Model predictions of Asmolov et al. [35] and Nizkaya et al. [37] are for 2D Poiseuille flow. B. Effect of Rep on the validity of the lift force model In the mixed scaling law model proposed by Hood et al. [43], the flow velo...
-
[10]
J. M. Martel and M. Toner, Annual review of biomedical eng ineering 16, 371 (2014)
work page 2014
- [11]
- [12]
-
[13]
D. Di Carlo, D. Irimia, R. G. Tompkins, and M. Toner, Proce edings of the National Academy of Sciences 104, 18892 (2007)
work page 2007
- [14]
- [15]
-
[16]
S. S. Shevkoplyas, T. Y oshida, L. L. Munn, and M. W. Bitens ky, Analytical chemistry 77, 933 (2005)
work page 2005
-
[17]
P . Han, S. Y osinski, Z. A. Kobos, R. Chaudhury, J. S. Lee, T . M. Fahmy, and M. A. Reed, ACS nano 14, 8646 (2020)
work page 2020
-
[18]
K. Jang, Y . Tanaka, J. Wakabayashi, R. Ishii, K. Sato, K. M awatari, M. Nilsson, and T. Kitamori, Biomicrofluidics 6, 044117 (2012)
work page 2012
-
[19]
T. Roebroek, W. V an Roy, S. Roth, L. Chacon Orellana, Z. L uo, Y . El Jerrari, C. Arnett, K. Claes, S. Ha, K. Jans, R. Labie, Z. Lin, M. Obst, D. Origuella, V . Pham , F. V an Bellinghen, A. V . Giriku- mar Krishna, E. V aezzadeh, W. V anhove, J. V an Duppen, G. Wang, S. Willems, P . Peumans, M. Jaya- pala, and T. Stakenborg, Analytical Chemistry 0, null...
- [20]
-
[21]
R. C. Jeffrey and J. Pearson, Journal of Fluid Mechanics 22, 721 (1965)
work page 1965
- [22]
- [23]
- [24]
-
[25]
J. A. Schonberg and E. Hinch, Journal of Fluid Mechanics 203, 517 (1989)
work page 1989
-
[26]
J. Feng, H. H. Hu, and D. D. Joseph, Journal of Fluid Mecha nics 261, 95 (1994)
work page 1994
-
[27]
E. S. Asmolov, Journal of fluid mechanics 381, 63 (1999)
work page 1999
-
[28]
M. Han, C. Kim, M. Kim, and S. Lee, Journal of rheology 43, 1157 (1999)
work page 1999
- [29]
-
[30]
Z. Y u, N. Phan-Thien, and R. I. Tanner, Journal of Fluid M echanics 518, 61 (2004)
work page 2004
-
[31]
B. H. Y ang, J. Wang, D. D. Joseph, H. H. Hu, T.-W. Pan, and R. Glowinski, Journal of Fluid mechanics 540, 109 (2005). 26
work page 2005
- [32]
-
[33]
A. A. S. Bhagat, S. S. Kuntaegowdanahalli, and I. Papaut sky, Physics of Fluids 20 (2008)
work page 2008
-
[34]
Y . W. Kim and J. Y . Y oo, Journal of Micromechanics and Microengineering 18, 065015 (2008)
work page 2008
-
[35]
X. Shao, Z. Y u, and B. Sun, Physics of Fluids 20 (2008)
work page 2008
- [36]
- [37]
- [38]
- [39]
-
[40]
N. Nakagawa, T. Y abu, R. Otomo, A. Kase, M. Makino, T. Ita no, and M. Sugihara-Seki, Journal of fluid mechanics 779, 776 (2015)
work page 2015
- [41]
-
[42]
H. T. Kazerooni, W. Fornari, J. Hussong, and L. Brandt, P hysical review fluids 2, 084301 (2017)
work page 2017
- [43]
-
[44]
E. S. Asmolov, A. L. Dubov, T. V . Nizkaya, J. Harting, and O. I. Vinogradova, Journal of fluid mechanics 840, 613 (2018)
work page 2018
-
[45]
S. Nakayama, H. Y amashita, T. Y abu, T. Itano, and M. Sugi hara-Seki, Journal of Fluid Mechanics 871, 952 (2019)
work page 2019
-
[46]
T. V . Nizkaya, A. S. Gekova, J. Harting, E. S. Asmolov, an d O. I. Vinogradova, Physics of Fluids 32 (2020)
work page 2020
-
[47]
J. Su, X. Zheng, and G. Hu, Physics of Fluids 35 (2023)
work page 2023
- [48]
-
[49]
M. E. Staben, A. Z. Zinchenko, and R. H. Davis, physics of fluids 15, 1711 (2003)
work page 2003
-
[50]
M. E. Staben and R. H. Davis, International journal of mu ltiphase flow 31, 529 (2005)
work page 2005
-
[51]
D. Di Carlo, J. F. Edd, K. J. Humphry, H. A. Stone, and M. To ner, Physical review letters 102, 094503 (2009)
work page 2009
-
[52]
K. Hood, S. Lee, and M. Roper, Journal of Fluid Mechanics 765, 452 (2015)
work page 2015
- [53]
-
[54]
P . Cherukat and J. B. McLaughlin, Journal of Fluid Mecha nics 263, 1 (1994)
work page 1994
-
[55]
E. P . V an Der Poel, R. Ostilla-Mónico, J. Donners, and R. V erzicco, Computers & Fluids 116, 10 27 (2015)
work page 2015
-
[56]
M. D. de Tullio and G. Pascazio, Journal of Computationa l Physics 325, 201 (2016)
work page 2016
-
[57]
V . Spandan, D. Lohse, M. D. de Tullio, and R. V erzicco, Jo urnal of computational physics 375, 228 (2018)
work page 2018
-
[58]
Uhlmann, Journal of computational physics 209, 448 (2005)
M. Uhlmann, Journal of computational physics 209, 448 (2005)
work page 2005
-
[59]
Breugem, Journal of Computational Physics 231, 4469 (2012)
W.-P . Breugem, Journal of Computational Physics 231, 4469 (2012)
work page 2012
- [60]
-
[61]
G. Wang, M. Abbas, and É. Climent, Physical Review Fluid s 2, 084302 (2017)
work page 2017
- [62]
-
[63]
N. Vrancken, J. Li, S. Sergeant, G. V ereecke, G. Doumen, F. Holsteyns, C. Chen, H. Terryn, S. De Gendt, and X. Xu, Scientific reports 8, 11637 (2018)
work page 2018
-
[64]
T. V . Nizkaya, E. S. Asmolov, J. Harting, and O. I. Vinogr adova, Physical review fluids 5, 014201 (2020)
work page 2020
-
[65]
O. I. Vinogradova and A. V . Belyaev, Journal of Physics: Condensed Matter 23, 184104 (2011)
work page 2011
-
[66]
O. I. Vinogradova, Langmuir 11, 2213 (1995)
work page 1995
- [67]
- [68]
- [69]
-
[70]
P . Reschiglian, D. Melucci, G. Torsi, and A. Zattoni, Ch romatographia 51, 87 (2000)
work page 2000
-
[71]
P . Cherukat and J. B. Mclaughlin, Journal of Fluid Mecha nics 285, 407–407 (1995)
work page 1995
-
[72]
Brenner, Chemical engineering science 16, 242 (1961)
H. Brenner, Chemical engineering science 16, 242 (1961)
work page 1961
-
[73]
J. R. Smart and D. T. Leighton Jr, Physics of Fluids A: Flu id Dynamics 1, 52 (1989). 28
work page 1989
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.