pith. sign in

arxiv: 2511.03956 · v2 · submitted 2025-11-06 · ⚛️ physics.flu-dyn · math-ph· math.MP

Thin gap approximations for microfluidic device design

Lingyun Ding , Terry Wang , Marcus Roper This is my paper

Pith reviewed 2026-05-18 01:38 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn math-phmath.MP
keywords Hele-Shaw flowmicrofluidicsthin gap approximationmethod of weighted residualsreduced modelingNavier-Stokesinertial microfluidicsorthogonal polynomials
0
0 comments X

The pith

Higher-order Hele-Shaw models capture non-parabolic flows in microfluidic devices

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

This paper derives a family of two-dimensional models for flow in thin gaps by applying the method of weighted residuals to the three-dimensional Navier-Stokes equations and expanding the velocity in orthogonal polynomials across the gap height. The leading term recovers the classic Hele-Shaw potential flow, but the next terms add corrections for non-parabolic velocity profiles and flow components perpendicular to the gap plane. A reader should care because microfluidic devices almost always have thin gaps, so these reduced models let engineers simulate and design devices much faster than with full three-dimensional calculations while keeping the essential physics. Numerical tests on real device shapes confirm that the corrections improve accuracy over the basic Hele-Shaw model.

Core claim

Using the Method of Weighted Residuals, the authors reinterpret the Hele-Shaw approximation as the leading term of an orthogonal polynomial expansion across the gap. Extending this expansion produces a new reduced two-dimensional model that includes non-parabolic gap-wise velocity profiles and out-of-plane flow effects. Substantial numerical evidence demonstrates that these approximate equations successfully model flows in actual microfluidic and inertial-microfluidic device geometries.

What carries the argument

Orthogonal polynomial expansion of the velocity field across the gap, obtained via the Method of Weighted Residuals applied to the Navier-Stokes equations.

If this is right

  • Engineers can use the extended model to predict flows in complex microfluidic geometries without solving the full three-dimensional equations.
  • Out-of-plane velocity components become accessible in the reduced two-dimensional framework.
  • The same approach works for both low-Reynolds-number and inertial microfluidic flows.
  • Device design cycles shorten because two-dimensional meshes replace three-dimensional ones.

Where Pith is reading between the lines

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

  • These corrections might allow the model to handle moderate gap aspect ratios where the basic Hele-Shaw approximation breaks down.
  • Integration with topology optimization could help discover optimal channel layouts for specific flow tasks.
  • Analogous expansions may improve reduced models for other thin-layer fluid problems outside microfluidics.

Load-bearing premise

The gap height remains small enough compared with lateral dimensions that a polynomial expansion in the gap direction converges to the true three-dimensional solution.

What would settle it

Running a full three-dimensional simulation of flow through a thin microfluidic T-junction and comparing the gap-averaged velocity and out-of-plane components against the reduced model's predictions; large mismatches at moderate Reynolds numbers would falsify the claim.

Figures

Figures reproduced from arXiv: 2511.03956 by Lingyun Ding, Marcus Roper, Terry Wang.

Figure 1
Figure 1. Figure 1: The Hele–Shaw approximation models thin-gap flows via 2D potential theory. (a) Streamlines around a model airfoil: top—experiment Werle (1973); bottom—2D potential flow computed in COMSOL 5.4. (b) Thin-gap geometry as a testbed for Saffman–Taylor instability, where a low-viscosity fluid (blue-dyed water) displaces a more viscous one (glycerol). (c) Example microfluidic device with thin gap: a centrifuge-on… view at source ↗
Figure 2
Figure 2. Figure 2: Validating a 2D approximation for the pressure-driven flow in a rectangular channel. (a) Flow profile at z = 0: exact solution (black curve), approximation from Eq. (2.3) (red curve), and approximation from Eq. (2.5) (blue dashed curve), for L = 1. (b) Relative error, measured by the L ∞ norm ∥u−uapp∥∞ ∥u∥∞ , as a function of aspect ratio L. Red curve: Eq. (2.3), blue dashed curve: Eq. (2.5). integrating t… view at source ↗
Figure 3
Figure 3. Figure 3: Coaxial flow. (a) Simulation geometry based on Anna et al. (2003): left— inlets, right- outlet. Outer inlets have mass flux Q, the inner inlet has flux 1. Inlet and outlet lengths are extended by 20 to ensure fully developed flow. (b) Dividing streamlines for Q = 1.7: classical Hele–Shaw (red), derived a = 6/5 approximation (blue dotted), and 3D simulation (black dashed). in which our optimal 2D approximat… view at source ↗
Figure 4
Figure 4. Figure 4: 2D approximation for finite Re flows in a centrifuge-on-a-chip geometry. Comparison of the separation bubble width in 3D simulation (black), 2D approximation (blue, Eq. (3.2a)), and unweighted residual approximation (red, Eq. (3.2b)). Right panels: Streamline patterns colored by velocity magnitude (red = larger, blue = smaller), comparing the 2D approximation (top) with the full 3D simulation (bottom). and… view at source ↗
Figure 5
Figure 5. Figure 5: Fidelity of the numerically computed 3D flow field to the Hele-Shaw (parabolic flow profile) for a centrifuge-on-a-chip flow. (a) r (defined in Eq. (4.1)) at Re = 200 across the centrifuge chamber. (b,c) u− profiles against z sampled on two different transects (locations shown in (a)). To compare shapes, all profiles are normalized so that u = 1 at z = 0 (b) shows a small r-region, where the Hele-Shaw appr… view at source ↗
read the original abstract

Over 125 years ago, Henry Selby Hele-Shaw realized that the depth-averaged flow in thin gap geometries can be closely approximated by two-dimensional (2D) potential flow, in a surprising marriage between the theories of viscous-dominated and inviscid flows. Hele-Shaw approximation allows visualization of potential flows over 2D airfoils and also undergirds important discoveries in the dynamics of interfacial instabilities and convection, yet it has found little use in modeling flows in microfluidic devices, although these devices often have thin gap geometries. Here, we derive a Hele-Shaw approximation for the flow in the kinds of thin gap geometries created within microfluidic devices. Using the Method of Weighted Residuals (MWR), we reinterpret the Hele-Shaw approximation as the leading term of an orthogonal polynomial expansion that can be systematically extended to higher-order corrections. The resulting leading-order equation coincides with the previously derived 2D approximations, but our derivation is shorter and more direct. By extending the expansion beyond leading order, we obtain a new reduced model that captures non-parabolic gap-wise velocity profiles and out-of-plane flow effects. We provide substantial numerical evidence showing that approximate equations can successfully model real microfluidic and inertial-microfluidic device geometries. By reducing three-dimensional (3D) flows to 2D models, our validated model will allow for accelerated device modeling and design.

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 derives a thin-gap reduced model for microfluidic flows by applying the Method of Weighted Residuals to the 3D Navier-Stokes equations, interpreting the classical Hele-Shaw approximation as the leading term of an orthogonal polynomial expansion across the gap height. Higher-order terms are introduced to capture non-parabolic gap-wise velocity profiles and out-of-plane flow components. The authors present numerical comparisons demonstrating that the resulting 2D models accurately reproduce flows in realistic microfluidic and inertial-microfluidic device geometries.

Significance. If the numerical evidence is robust, the work supplies a systematic, extensible framework for reducing 3D thin-gap flows to 2D models with controlled corrections. This could materially accelerate iterative design of microfluidic devices by replacing repeated full 3D simulations with cheaper 2D solves while retaining accuracy beyond the classical parabolic Hele-Shaw limit. The absence of free parameters in the derivation and the explicit connection to orthogonal expansions are positive features.

major comments (2)
  1. [§4] §4 (numerical results for inertial-microfluidic geometries): the central claim that the higher-order MWR model remains accurate for inertial-microfluidic cases (Re ~ 10–100) rests on the assertion that a fixed low-order polynomial basis continues to span the gap-wise profiles once convective inertia distorts them. The manuscript should report the truncation residual or L2 error of the weighted residual for at least one inertial case at a location where secondary flow is expected; without this metric it is unclear whether the basis choice remains adequate or whether additional terms are required.
  2. [§3.2] §3.2 (MWR formulation): the weighting functions and polynomial basis are chosen to satisfy no-slip at the walls under the thin-gap scaling. When inertial terms are retained, the reduced Reynolds number based on gap height is no longer ≪1; the paper should state the maximum reduced Re for which the reported truncation order is guaranteed to control the error, or demonstrate convergence with increasing polynomial degree for an inertial benchmark.
minor comments (2)
  1. Figure captions should explicitly label which curve corresponds to leading-order, first correction, and full 3D reference solutions.
  2. The abstract states that the leading-order equation 'coincides with previously derived 2D approximations'; a brief side-by-side comparison of the final PDE with one standard reference (e.g., the depth-averaged Stokes or Hele-Shaw form) would help readers verify the claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address each major comment below and indicate the revisions that will be incorporated into the next version of the manuscript.

read point-by-point responses

  1. Referee: [§4] §4 (numerical results for inertial-microfluidic geometries): the central claim that the higher-order MWR model remains accurate for inertial-microfluidic cases (Re ~ 10–100) rests on the assertion that a fixed low-order polynomial basis continues to span the gap-wise profiles once convective inertia distorts them. The manuscript should report the truncation residual or L2 error of the weighted residual for at least one inertial case at a location where secondary flow is expected; without this metric it is unclear whether the basis choice remains adequate or whether additional terms are required.

    Authors: We agree that an explicit truncation residual or L2 error of the weighted residual for an inertial case would strengthen the validation. In the revised manuscript we will add this diagnostic for a representative inertial-microfluidic geometry at a station where secondary flow is observed, confirming that the residual remains small for the chosen polynomial degree. revision: yes

  2. Referee: [§3.2] §3.2 (MWR formulation): the weighting functions and polynomial basis are chosen to satisfy no-slip at the walls under the thin-gap scaling. When inertial terms are retained, the reduced Reynolds number based on gap height is no longer ≪1; the paper should state the maximum reduced Re for which the reported truncation order is guaranteed to control the error, or demonstrate convergence with increasing polynomial degree for an inertial benchmark.

    Authors: We accept the suggestion to clarify the range of applicability. The derivation is performed under the thin-gap scaling with reduced Reynolds number Re_δ = Re δ². In the revised §3.2 we will state the regime Re_δ ≪ 1 in which the present truncation order is expected to control the error, together with a short numerical demonstration of convergence under mesh refinement in polynomial degree for one inertial benchmark. revision: yes

Circularity Check

0 steps flagged

Derivation from 3D Navier-Stokes via standard MWR is self-contained with no circular reductions

full rationale

The paper begins from the three-dimensional viscous flow equations in thin-gap geometries and applies the Method of Weighted Residuals using an orthogonal polynomial expansion across the gap height. The leading-order term recovers the classical Hele-Shaw approximation as a direct consequence of the truncation, while higher-order terms systematically extend the same weighted-residual procedure to capture non-parabolic profiles and out-of-plane effects. No parameters are fitted to data and then relabeled as predictions, no self-citations form the load-bearing justification for the core steps, and the numerical evidence is presented as independent validation rather than part of the derivation chain. The thin-gap scaling and polynomial basis choice are standard assumptions justified by the geometry (height much smaller than lateral dimensions), not smuggled in via prior self-work or defined circularly. This is a normal, non-circular reduction of 3D equations to 2D models.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on the standard incompressible Navier-Stokes equations for viscous flow together with the geometric assumption that the gap height is small compared with lateral dimensions; no new free parameters, ad-hoc constants, or postulated entities are introduced in the abstract.

axioms (2)
  • domain assumption The gap height is much smaller than the lateral dimensions of the device, allowing depth-averaging and polynomial expansion across the gap.
    This thin-gap premise is invoked at the start of the derivation to justify reducing the 3D problem to a 2D model with corrections.
  • standard math The flow is governed by the incompressible Navier-Stokes equations at low to moderate Reynolds number.
    Standard governing equations for viscous microfluidic flow are assumed without re-derivation.

pith-pipeline@v0.9.0 · 5777 in / 1603 out tokens · 38304 ms · 2026-05-18T01:38:35.092583+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

39 extracted references · 39 canonical work pages

  1. [1]

    flow focusing

    Anna, Shelley L , Bontoux, Nathalie & Stone, Howard A 2003 Formation of dispersions using “flow focusing” in microchannels . Appl. Phys. Lett. 82 , 364--366

    work page 2003

  2. [2]

    Journal of Fluid Mechanics 372 , 25--44

    Balsa, Thomas F 1998 Secondary flow in a hele-shaw cell . Journal of Fluid Mechanics 372 , 25--44

    work page 1998

  3. [3]

    Journal of Applied Analysis 6 (2), 259--282

    Belinsky, Rachel 2000 Integrals of legendre polynomials and solution of some partial differential equations . Journal of Applied Analysis 6 (2), 259--282

    work page 2000

  4. [4]

    Proceedings of the Royal Society A 481 (2311), 20240613

    Booth, DJ , Griffiths, IM & Howell, PD 2025 The motion of a bubble in a non-uniform hele-shaw flow . Proceedings of the Royal Society A 481 (2311), 20240613

    work page 2025

  5. [5]

    Journal de Mathématiques Pures et Appliquées 13 (2), 377--424

    Boussinesq, Joseph 1868 Mémoire sur l'influence des frottements dans les mouvements réguliers des fluids . Journal de Mathématiques Pures et Appliquées 13 (2), 377--424

    work page

  6. [6]

    Flow, Turbulence and Combustion 1 (1), 27--34

    Brinkman, Hendrik C 1949 A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles . Flow, Turbulence and Combustion 1 (1), 27--34

    work page 1949

  7. [7]

    Transport in porous media 44 , 325--335

    Chen, Zhangxin , Lyons, Stephen L & Qin, Guan 2001 Derivation of the forchheimer law via homogenization . Transport in porous media 44 , 325--335

    work page 2001

  8. [8]

    Journal of fluid mechanics 552 , 83--97

    Chevalier, Christophe , Amar, Martine Ben , Bonn, Daniel & Lindner, Anke 2006 Inertial effects on saffman--taylor viscous fingering . Journal of fluid mechanics 552 , 83--97

    work page 2006

  9. [9]

    Proceedings of the Royal Society of London

    Cornish, R Jo 1928 Flow in a pipe of rectangular cross-section . Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 120 (786), 691--700

    work page 1928

  10. [10]

    Journal of Statistical Physics 131 , 941--967

    Desvillettes, Laurent , Golse, Fran c ois & Ricci, Valeria 2008 The mean-field limit for solid particles in a navier-stokes flow . Journal of Statistical Physics 131 , 941--967

    work page 2008

  11. [11]

    Lab on a Chip 9 (21), 3038--3046

    Di Carlo, Dino 2009 Inertial microfluidics . Lab on a Chip 9 (21), 3038--3046

    work page 2009

  12. [12]

    Journal of Fluid Mechanics 970 , A27

    Ding, Lingyun 2023 Shear dispersion of multispecies electrolyte solutions in the channel domain . Journal of Fluid Mechanics 970 , A27

    work page 2023

  13. [13]

    Finlayson, Bruce A 2013 The method of weighted residuals and variational principles\/ . SIAM

    work page 2013

  14. [14]

    Physics of Fluids 9 (11), 3267--3274

    Gondret, Philippe & Rabaud, Marc 1997 Shear instability of two-fluid parallel flow in a hele--shaw cell . Physics of Fluids 9 (11), 3267--3274

    work page 1997

  15. [15]

    Tamkang journal of mathematics 41 (3), 217--243

    Guta, Lemi & Sundar, S 2010 Navier-stokes-brinkman system for interaction of viscous waves with a submerged porous structure . Tamkang journal of mathematics 41 (3), 217--243

    work page 2010

  16. [16]

    Annual review of fluid mechanics 19 (1), 271--311

    Homsy, George M 1987 Viscous fingering in porous media . Annual review of fluid mechanics 19 (1), 271--311

    work page 1987

  17. [17]

    Journal of Fluid Mechanics 64 (3), 449--476

    Howells, ID 1974 Drag due to the motion of a newtonian fluid through a sparse random array of small fixed rigid objects . Journal of Fluid Mechanics 64 (3), 449--476

    work page 1974

  18. [18]

    Biomicrofluidics 5 (2)

    Hur, Soojung Claire , Mach, Albert J & Di Carlo, Dino 2011 High-throughput size-based rare cell enrichment using microscale vortices . Biomicrofluidics 5 (2)

    work page 2011

  19. [19]

    Physical Review Fluids 5 (6), 064101

    Inamdar, Tanmay C , Wang, Xiaojia & Christov, Ivan C 2020 Unsteady fluid-structure interactions in a soft-walled microchannel: A one-dimensional lubrication model for finite reynolds number . Physical Review Fluids 5 (6), 064101

    work page 2020

  20. [20]

    Physics of fluids 16 (8), 3051--3062

    Lauga, Eric , Stroock, Abraham D & Stone, Howard A 2004 Three-dimensional flows in slowly varying planar geometries . Physics of fluids 16 (8), 3051--3062

    work page 2004

  21. [21]

    Physics of Fluids 33 (1)

    Li, Di , Song, Le , Zhang, Cheng , Yu, Liandong & Xuan, Xiangchun 2021 A depth-averaged model for newtonian fluid flows in shallow microchannels . Physics of Fluids 33 (1)

    work page 2021

  22. [22]

    Journal of Fluid Mechanics 608 , 43--70

    Lin, Hao , Storey, Brian D & Santiago, Juan G 2008 A depth-averaged electrokinetic flow model for shallow microchannels . Journal of Fluid Mechanics 608 , 43--70

    work page 2008

  23. [23]

    Physics of Fluids 14 (2), 902--905

    Martin, J , Rakotomalala, N & Salin, D 2002 Gravitational instability of miscible fluids in a hele-shaw cell . Physics of Fluids 14 (2), 902--905

    work page 2002

  24. [24]

    Microfluidics and Nanofluidics 21 , 1--11

    Pai \`e , P , Che, J & Di Carlo, D 2017 Effect of reservoir geometry on vortex trapping of cancer cells . Microfluidics and Nanofluidics 21 , 1--11

    work page 2017

  25. [25]

    Journal of the Royal Society Interface 7 (46), 851--862

    Pepper, Rachel E , Roper, Marcus , Ryu, Sangjin , Matsudaira, Paul & Stone, Howard A 2010 Nearby boundaries create eddies near microscopic filter feeders . Journal of the Royal Society Interface 7 (46), 851--862

    work page 2010

  26. [26]

    Physics of Fluids 14 (3), 922--929

    Plourabou \'e , Franck & Hinch, E John 2002 Kelvin--helmholtz instability in a hele-shaw cell . Physics of Fluids 14 (3), 922--929

    work page 2002

  27. [27]

    Mechanics research communications 36 (1), 10--21

    Qian, Shizhi & Bau, Haim H 2009 Magneto-hydrodynamics based microfluidics . Mechanics research communications 36 (1), 10--21

    work page 2009

  28. [28]

    Comptes Rendus de l'Acad \'e mie des Sciences-Series IIB-Mechanics 329 (5), 337--342

    Ruyer-Quil, Christian 2001 Inertial corrections to the darcy law in a hele--shaw cell . Comptes Rendus de l'Acad \'e mie des Sciences-Series IIB-Mechanics 329 (5), 337--342

    work page 2001

  29. [29]

    Proceedings of the Royal Society of London

    Saffman, Philip Geoffrey & Taylor, Geoffrey Ingram 1958 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 Sciences 245 (1242), 312--329

    work page 1958

  30. [30]

    Journal of applied physics 101 (7)

    Salmon, Jean-Baptiste & Ajdari, Armand 2007 Transverse transport of solutes between co-flowing pressure-driven streams for microfluidic studies of diffusion/reaction processes . Journal of applied physics 101 (7)

    work page 2007

  31. [31]

    Stone, Howard A , Stroock, Abraham D & Ajdari, Armand 2004 Engineering flows in small devices: microfluidics toward a lab-on-a-chip . Annu. Rev. Fluid Mech. 36 , 381--411

    work page 2004

  32. [32]

    Science 295 (5555), 647--651

    Stroock, Abraham D , Dertinger, Stephan KW , Ajdari, Armand , Mezic, Igor , Stone, Howard A & Whitesides, George M 2002 Chaotic mixer for microchannels . Science 295 (5555), 647--651

    work page 2002

  33. [33]

    Journal of Fluid Mechanics 226 , 125--148

    Tsay, Ruey-Yug & Weinbaum, Sheldon 1991 Viscous flow in a channel with periodic cross-bridging fibres: exact solutions and brinkman approximation . Journal of Fluid Mechanics 226 , 125--148

    work page 1991

  34. [34]

    Annual Review of Fluid Mechanics 5 (Volume 5, 1973), 361--386

    Werle, H 1973 Hydrodynamic flow visualization . Annual Review of Fluid Mechanics 5 (Volume 5, 1973), 361--386

    work page 1973

  35. [35]

    Transport in porous media 1 , 3--25

    Whitaker, Stephen 1986 Flow in porous media i: A theoretical derivation of darcy's law . Transport in porous media 1 , 3--25

    work page 1986

  36. [36]

    Fluid Dynamics Research 47 (1), 015506

    Yuan, Qingwang & Azaiez, Jalel 2014 Inertial effects of miscible viscous fingering in a hele-shaw cell . Fluid Dynamics Research 47 (1), 015506

    work page 2014

  37. [37]

    Physics of Fluids 15 (12), 3829--3836

    Zeng, Jie , Yortsos, Yannis C & Salin, Dominique 2003 On the brinkman correction in unidirectional hele-shaw flows . Physics of Fluids 15 (12), 3829--3836

    work page 2003

  38. [38]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...

    work page

  39. [39]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

    work page