pith. sign in

arxiv: 2604.26189 · v1 · submitted 2026-04-29 · ❄️ cond-mat.soft · physics.flu-dyn

Viscous Settling of Bravais Unit-Cells

Sebastian B\"urger , Harshit Joshi , S Ganga Prasath , Rahul Chajwa , Rama Govindarajan This is my paper

Pith reviewed 2026-05-07 12:58 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.flu-dyn
keywords viscous settlingBravais latticeporosityStokes flowsedimentationpower-law scalingwall effectsFaxen correction
0
0 comments X

The pith

Bravais lattice units settle with speed proportional to porosity raised to the 0.3 power in unbounded fluid.

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

The paper examines how Bravais lattice unit-cells with controllable porosity fall through viscous fluid under Stokes flow. In container experiments the settling speed scales with solid fraction to the power 0.43 independent of lattice shape. Analytic and numerical work shows that distant walls produce this exponent; applying Faxen's correction removes the wall contribution and yields an intrinsic exponent of 0.30 for unbounded domains. The resulting power-law offers a route to estimating sedimentation rates of porous aggregates.

Core claim

Porous Bravais lattice unit-cells exhibit a power-law dependence of Stokes settling speed on solid fraction, with exponent 0.43 in confined experiments that corrects to 0.30 when extrapolated to unbounded domains using boundary corrections; the exponent is independent of lattice shape.

What carries the argument

The power-law relationship between settling speed U and solid fraction φ, obtained from Stokes flow around lattice units and extrapolated via Faxen's wall correction.

If this is right

  • Settling speed for these regular porous shapes can be predicted from porosity alone.
  • The scaling is the same across different Bravais lattices, indicating shape independence within this class.
  • Wall effects must be subtracted from confined experiments to recover the intrinsic scaling.
  • The relation supplies a building block for estimating bulk sedimentation fluxes of complex objects.

Where Pith is reading between the lines

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

  • If porosity dominates over detailed geometry, similar exponents may describe irregular natural aggregates.
  • Repeating the experiments with varied container aspect ratios would test the robustness of the Faxen extrapolation.
  • The 0.30 exponent could be used to close simple models of particle flux in clouds or ocean layers.

Load-bearing premise

Faxen's boundary correction accurately represents the influence of distant container walls on the settling speed.

What would settle it

Direct numerical simulation or measurement of the same lattice units settling in a domain large enough that wall effects are negligible, to test whether the exponent converges to 0.30.

Figures

Figures reproduced from arXiv: 2604.26189 by Harshit Joshi, Rahul Chajwa, Rama Govindarajan, Sebastian B\"urger, S Ganga Prasath.

Figure 1
Figure 1. Figure 1: FIG. 1: ( view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: (a) Schematic of the setup of the PIV-experiment with laser-sheet bisecting the object, revealing the flows view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Normalized terminal velocity obtained in the view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Analogous to Fig. 3, but with wall effects removed for each data point in the bounded domain by using view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Schematic of the setup used for the Stokesian view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Velocity field lines on the PIV plane (yz plane) near falling SC structure in the bounded container, obtained view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Excess terminal velocity of the SC structure for different lattice spacing view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Influence of connecting rods on the normalized view at source ↗
read the original abstract

We study experimentally and theoretically the Stokesian settling of a well-known class of porous shapes: Bravais lattice unit-cells, whose porosity we vary controllably by changing their lattice spacing. In our experiments, conducted in a square cuboidal container with its long-axis aligned along gravity, we find that the settling speed U and the solid fraction {\phi} of these lattice units obey a power-law relationship U $\propto$ {\phi}^{\gamma} , with an exponent {\gamma} = 0.43 independent of their shape. To understand the observed scaling exponent, we analytically and numerically investigate the settling of the simple cubic structure under different approximations. We find that the walls of the container, though far from the sinking object, have a defining effect. Our Stokesian boundary integral simulations show that the Faxen's boundary correction captures the wall-effects accurately and enables us to discount the wall-effect from the experimental data, yielding a power-law exponent {\gamma} = 0.30 for settling in an unbounded domain. The power-law relating sinking speed and porosity is a step towards predictively understanding the sedimentation fluxes of complex objects in the clouds and the oceans. However, the applicability of this universal scaling to irregular and biologically richer aggregates found in nature remains an open direction.

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 paper experimentally studies the Stokesian settling of Bravais lattice unit-cells (simple cubic, BCC, FCC) with controllable porosity in a square container, reporting that settling speed U scales as U ∝ φ^γ with γ = 0.43 independent of lattice type. Boundary-integral simulations for the simple-cubic case are used to validate Faxén's wall correction, which is then applied uniformly to extrapolate all experimental data to an unbounded domain, yielding γ = 0.30. The work positions this scaling as a step toward predicting sedimentation of complex porous objects.

Significance. If the shape-independent scaling and the uniform applicability of the Faxén correction hold, the result supplies a compact relation for sedimentation fluxes of lattice-like aggregates with relevance to marine snow and cloud microphysics. The controllable experimental porosity variation and the Stokesian simulation approach for at least one geometry constitute clear strengths.

major comments (2)
  1. [§4.2] §4.2 (boundary-integral simulations): The Faxén correction is validated numerically only for the simple-cubic lattice. No equivalent simulations are reported for BCC or FCC structures, yet data from all three lattices are combined in the single power-law fit that produces the headline unbounded-domain exponent γ = 0.30. Because far-field wall interactions depend on object extent and internal permeability (which differ by lattice symmetry), the correction factor may vary; if the variation is comparable to the shift from 0.43 to 0.30, the extrapolated universal scaling is not yet secured.
  2. [§3] §3 (experimental results): The power-law exponent γ = 0.43 is stated to be independent of shape, but the manuscript provides neither per-lattice fit uncertainties, R² values, nor a statistical test for the null hypothesis that the exponents differ across SC, BCC, and FCC data sets. This weakens the central claim that the scaling is shape-independent before any wall correction is applied.
minor comments (2)
  1. [Abstract] Abstract: No error bars, data-exclusion criteria, or quantitative fit diagnostics are mentioned for the reported exponents.
  2. [Notation] Notation: The solid fraction is written both as φ and as {phi}; adopt a single symbol throughout.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We address each of the major comments in detail below.

read point-by-point responses

  1. Referee: [§4.2] §4.2 (boundary-integral simulations): The Faxén correction is validated numerically only for the simple-cubic lattice. No equivalent simulations are reported for BCC or FCC structures, yet data from all three lattices are combined in the single power-law fit that produces the headline unbounded-domain exponent γ = 0.30. Because far-field wall interactions depend on object extent and internal permeability (which differ by lattice symmetry), the correction factor may vary; if the variation is comparable to the shift from 0.43 to 0.30, the extrapolated universal scaling is not yet secured.

    Authors: We appreciate this observation regarding the scope of our numerical validation. The boundary-integral simulations were indeed limited to the simple-cubic lattice to demonstrate the accuracy of the Faxén correction in our setup. We argue that the far-field wall correction primarily depends on the overall dimensions of the settling object and its distance to the container walls, rather than fine details of the internal lattice structure. Since the SC, BCC, and FCC unit-cells are designed with identical outer bounding boxes and comparable porosities in the experiments, the hydrodynamic 'effective radius' relevant for the Faxén formula should be similar across lattices. Internal permeability affects the drag but the wall correction is applied to the unbounded drag. Nevertheless, to fully address the concern, in the revised manuscript we will expand the discussion to include a theoretical justification for the uniform applicability and, if computational resources permit, add simulations for the BCC case as an example to verify the correction factor remains within a small percentage. revision: partial

  2. Referee: [§3] §3 (experimental results): The power-law exponent γ = 0.43 is stated to be independent of shape, but the manuscript provides neither per-lattice fit uncertainties, R² values, nor a statistical test for the null hypothesis that the exponents differ across SC, BCC, and FCC data sets. This weakens the central claim that the scaling is shape-independent before any wall correction is applied.

    Authors: We agree that including detailed statistical information would bolster the claim of shape-independent scaling. The current presentation combines the data for a single fit, but we have the individual datasets. In the revised manuscript, we will report the power-law fits performed separately for each lattice type, including the fitted exponents with standard errors or confidence intervals, the corresponding R² values, and a statistical test (e.g., a test for equality of slopes in log-log space or ANOVA on the residuals) to evaluate the null hypothesis that the exponents are the same across the three lattices. This will provide transparent evidence supporting or qualifying the independence from lattice symmetry. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in the derivation chain.

full rationale

The paper's core result is an empirical power-law U ∝ φ^γ (γ = 0.43) obtained by direct fitting to experimental settling-speed measurements across multiple Bravais lattices. The subsequent step applies the established Faxén wall correction (validated by independent boundary-integral simulations performed only for the simple-cubic case) to the same experimental data set, producing the extrapolated unbounded-domain exponent γ = 0.30. This correction is not derived from the experimental fit itself, nor does any equation in the manuscript reduce the reported exponents to a self-referential definition or to a parameter that was fitted and then relabeled as a prediction. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work appear in the provided text. The theoretical investigation of the simple-cubic structure is presented as a separate explanatory analysis rather than a closed loop that forces the experimental scaling. The derivation therefore remains self-contained against external hydrodynamic benchmarks and does not meet the criteria for any enumerated circularity pattern.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The scaling law is empirical and rests on the Stokes equations for low-Reynolds-number flow plus the applicability of Faxen's wall correction; no new entities are postulated.

free parameters (1)
  • power-law exponent gamma
    Fitted directly to experimental settling speeds versus solid fraction for both bounded and wall-corrected cases.
axioms (2)
  • domain assumption Stokes flow (zero Reynolds number) governs the settling
    Invoked for viscous fluid at the experimental scales.
  • domain assumption Faxen's correction formula accurately represents distant wall effects for the lattice objects
    Used to discount container influence and obtain the unbounded exponent.

pith-pipeline@v0.9.0 · 5537 in / 1425 out tokens · 101560 ms · 2026-05-07T12:58:44.339914+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 5 canonical work pages

  1. [1]

    Guazzelli and J

    E. Guazzelli and J. F. Morris,A Physical Introduction to Suspension Dynamics(Cambridge University Press, 2011)

    2011

  2. [2]

    Clift, J

    R. Clift, J. R. Grace, and M. E. Weber,Bubbles, Drops, and Particles(Courier Corporation, 2005)

    2005

  3. [3]

    S. Zhao, K. F. Kvale, L. Zhu, E. R. Zettler, M. Egger, T. J. Mincer, L. A. Amaral-Zettler, L. Lebreton, H. Nie- mann, R. Nakajima,et al., Nature641, 51 (2025)

    2025

  4. [4]

    A. L. Alldredge and C. Gotschalk, Lim- nology and Oceanography33, 339 (1988), https://aslopubs.onlinelibrary.wiley.com/doi/pdf/ 10.4319/lo.1988.33.3.0339

    work page doi:10.4319/lo.1988.33.3.0339 1988

  5. [5]

    M. H. DiBenedetto, Annual Review of Fluid Mechanics 58, 355–382 (2026)

    2026

  6. [6]

    P. W. Boyd, H. Claustre, M. Levy, D. A. Siegel, and T. Weber, Nature568, 327–335 (2019)

    2019

  7. [7]

    Nabika, M

    H. Nabika, M. Itatani, and I. Lagzi, Langmuir36, 481–497 (2019)

    2019

  8. [8]

    D. L. Mitchell, P. Rasch, D. Ivanova, G. McFarquhar, and T. Nousiainen, Geophysical Research Letters35 (2008), 10.1029/2008gl033552

    work page doi:10.1029/2008gl033552 2008

  9. [9]

    E. C. Laurenceau-Cornec, F. A. C. Le Moigne, M. Gal- linari, B. Moriceau, J. Toullec, M. H. Iversen, A. Engel, and C. L. De La Rocha, Limnology and Oceanography 65, 1264–1285 (2019)

    2019

  10. [10]

    Kim and S

    S. Kim and S. J. Karrila,Microhydrodynamics: Princi- ples and Selected Applications(Dover, 2005)

    2005

  11. [11]

    Happel and H

    J. Happel and H. Brenner,Low Reynolds Number Hydro- dynamics(Springer, New York, 1983)

    1983

  12. [12]

    Bear,Dynamics of fluids in porous media(Courier Corporation, 2013)

    J. Bear,Dynamics of fluids in porous media(Courier Corporation, 2013)

    2013

  13. [13]

    H. C. Brinkman, Applied Scientific Research1, 27 (1947)

    1947

  14. [14]

    Debye and A

    P. Debye and A. M. Bueche, The Journal of Chemical Physics16, 573–579 (1948)

    1948

  15. [15]

    V. M. Starov and V. G. Zhdanov, Colloids and Sur- faces A: Physicochemical and Engineering Aspects192, 363–375 (2001)

    2001

  16. [16]

    Stimson and G

    M. Stimson and G. B. Jeffery, Proceedings of the Royal Society of London. Series A111, 110 (1926)

    1926

  17. [17]

    Cichocki and K

    B. Cichocki and K. Hinsen, Physics of Fluids7, 285 (1995)

    1995

  18. [18]

    I. A. Lasso and P. D. Weidman, Physics of Fluids29, 3921 (1986)

    1986

  19. [19]

    Sangani and A

    A. Sangani and A. Acrivos, International Journal of Mul- tiphase Flow8, 343–360 (1982)

    1982

  20. [20]

    Alabrudzi´ nski, M

    S. Alabrudzi´ nski, M. L. Ekiel-Je˙ zewska, D. Chehata- G´ omez, and T. A. Kowalewski, Physics of Fluids21 (2009), 10.1063/1.3168615

    work page doi:10.1063/1.3168615 2009

  21. [21]

    Binder,Settling of Fractal Aggregates in Viscous Me- dia, Doctoral thesis, University of Erlangen-Nuremberg (2012)

    C. Binder,Settling of Fractal Aggregates in Viscous Me- dia, Doctoral thesis, University of Erlangen-Nuremberg (2012)

    2012

  22. [22]

    Kittel and J

    C. Kittel and J. F. Masi, Physics Today7, 18–19 (1954)

    1954

  23. [23]

    J. F. Brady and G. Bossis, Annual Review of Fluid Me- chanics20, 111 (1988)

    1988

  24. [24]

    Pozrikidis,Boundary Integral and Singularity Meth- ods for Linearized Viscous Flow(Cambridge University Press, 1992)

    C. Pozrikidis,Boundary Integral and Singularity Meth- ods for Linearized Viscous Flow(Cambridge University Press, 1992)

    1992

  25. [25]

    Hill and G

    R. Hill and G. Power, The Quarterly Journal of Mechan- ics and Applied Mathematics9, 313 (1956)

    1956

  26. [26]

    J. M. Crowley, Journal of Fluid Mechanics45, 151–159 (1971)

    1971

  27. [27]

    A. K. Townsend, Journal of Open Source Software9, 6011 (2024)

    2024

  28. [28]

    Pozrikidis,A Practical Guide to Boundary Element Methods(CRC Press, 2002)

    C. Pozrikidis,A Practical Guide to Boundary Element Methods(CRC Press, 2002)

    2002

  29. [29]

    Prosperetti and G

    A. Prosperetti and G. Tryggvason,Computational Meth- ods for Multiphase Flow(Cambridge University Press, 2009)

    2009

  30. [30]

    A. K. Townsend,The mechanics of suspensions, Ph.D. thesis, University College London (2017)

    2017

  31. [31]

    C. M. Flintrop, S. Ahmerkamp, N. Moradi, I. Kla- wonn, J. M¨ arz, C. H¨ orstmann, R. Kiko, A. Khalili, H.-P. Grossart, X. A. ´Alvarez Salgado, J. Ar´ ıstegui, and M. H. Iversen, (2023), 10.1101/2023.09.05.555549

    work page doi:10.1101/2023.09.05.555549 2023

  32. [32]

    Kiørboe, Scientia Marina65, 57–71 (2001)

    T. Kiørboe, Scientia Marina65, 57–71 (2001)

    2001

  33. [33]

    Crystal time-reversal symmetry breaking and spontaneous hall effect in collinear anti- ferromagnets,

    R. Chajwa, E. Flaum, K. D. Bidle, B. Van Mooy, and M. Prakash, Science386(2024), 10.1126/sci- ence.adl5767

    work page doi:10.1126/sci- 2024