pith. sign in

arxiv: 2604.21064 · v1 · submitted 2026-04-22 · ⚛️ physics.flu-dyn

Particle-resolved simulations of settling particles: A methodology for long time-integration intervals

M. Moriche , M. Garc\'ia-Villalba , M. Uhlmann This is my paper

Pith reviewed 2026-05-09 22:39 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords settling particlesparticle-resolved simulationsmoving reference framedilute suspensionscollective effectsgravity-driven flowfluid-particle interaction
0
0 comments X

The pith

A moving reference frame removes vertical periodicity to enable long-time simulations of settling particles.

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

This paper presents a method for simulating dilute suspensions of particles settling under gravity that avoids the strong vertical correlations created by periodic boundaries. The approach drops vertical periodicity and shifts to a moving reference frame, so that particles settle through a fixed computational domain while the inlet velocity is adjusted accordingly. This setup permits efficient long integrations for single particles and for many-particle cases, recovering previously observed collective effects over an interval of roughly 600 particle diameters divided by the settling velocity. The only changes needed are updates to the inlet speed or fluid viscosity and the time step, which can be done in existing solvers without altering their core code. The configuration opens direct access to questions such as how still fluid meets the leading particle layer and what turbulent structures remain after a swarm has passed.

Core claim

The central claim is that a moving reference frame with non-periodic vertical boundaries produces physically equivalent settling dynamics to the laboratory frame, allowing converged simulations over intervals of approximately 600 D/U_g in many-particle cases without artificial corrections or periodicity-induced correlations.

What carries the argument

The moving reference frame, which translates the computational domain so that particles fall through it at a controlled rate while the inlet free-stream velocity is tuned to maintain physical equivalence.

Load-bearing premise

That adjusting the inlet velocity or viscosity and the time step in the moving frame produces settling dynamics identical to the laboratory frame without new numerical artifacts or changes to collective particle behavior.

What would settle it

A side-by-side comparison of particle velocity statistics, cluster size distributions, or fluid turbulence spectra between the moving-frame run and an equivalent fixed-frame simulation with a long vertical domain that would directly test whether the reported 600 D/U_g interval matches without discrepancies.

Figures

Figures reproduced from arXiv: 2604.21064 by M. Garc\'ia-Villalba, M. Moriche, M. Uhlmann.

Figure 1
Figure 1. Figure 1: a) Sketch of the particles settling in the unbounded domain and b) representation [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Sketch of the possible scenarios for values of [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of the solution of the illustrative example in the a) moving and b) laboratory [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of a) vertical velocity expressed as a particle Reynolds number (labora [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: a) Lateral and b) top view of the computational domain and the initial position of [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Time history of a) Πs m as they are computed (moving frame) and b) vertical velocity in the laboratory frame and width of the distribution. In a) the thin gray lines indicate change of moving frame, the black dashed line the last correction, the red dashed lines lines the limits of the computational domain and the green lines the minimum and maximum thresholds used in the correction algorithm. The vertical… view at source ↗
Figure 8
Figure 8. Figure 8: Visualization of a slab of the problem of a) the Q-criterion and b) the vertical velocity. [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Contour plot of particle concentration as a function of time and the vertical coordinate. [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: a) Vertical velocity of the particle ensemble (black line) and standard deviation of [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗
read the original abstract

We present a methodology for simulating dilute suspensions of particles settling under gravity, with the main purpose of overcoming limitations of triply periodic configurations, mainly the strong vertical correlation that hinders the study of cluster dynamics. The current approach removes vertical periodicity and employs a moving reference frame, enabling efficient simulations of both single- and many-particle cases. We illustrate the method with two examples of increasing complexity: a single particle in the steady vertical regime, and a many-particle case at a parametric point where collective effects were previously observed and recovered here. A converged, free-of-corrections time interval of approximately $600 D/U_g$ is simulated in the many-particle case, representing the first simulation of this kind to date. New physical insights can be explored thanks to this new configuration, for example the effect of still fluid on the first layer of particles encountered by the fluid, or the turbulent character of the flow after a swarm of particles has passed by. Finally, the method only requires parameter tuning, allowing implementation within existing solvers without changes to their core formulation: for a standard configuration with an imposed free stream velocity at the inlet, only the input velocity (or the viscosity of the fluid) and the time step need to be updated.

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 / 1 minor

Summary. The manuscript presents a methodology for particle-resolved simulations of dilute settling particle suspensions that removes vertical periodicity via a moving reference frame. This enables long time-integration intervals. The approach is illustrated on a single-particle steady settling case and a many-particle case at parameters where collective effects were previously observed and are claimed to be recovered here. A converged simulation interval of approximately 600 D/U_g is reported for the many-particle case, stated as the first of its kind. The method requires only updating the inlet free-stream velocity (or viscosity) and time step in existing solvers.

Significance. If the moving-frame inlet adjustment is shown to preserve the correct non-dimensional groups and collective dynamics without introducing artifacts, the methodology would overcome a major limitation of triply periodic configurations for studying long-term cluster formation and post-swarm turbulence in sedimentation. This could open new avenues for exploring dilute suspension behavior over extended times.

major comments (3)
  1. [Abstract] Abstract: the claim that collective effects observed in prior work are recovered in the many-particle case is not supported by any quantitative comparison (e.g., particle velocity distributions, fluid fluctuation spectra, or cluster statistics) to reference data or an equivalent lab-frame run. This is load-bearing for validating that the inlet adjustment preserves the physics.
  2. [Abstract] Abstract: the assertion of a converged, free-of-corrections interval of ~600 D/U_g lacks any supporting convergence study, error metrics, or time-series analysis to demonstrate that the simulation has reached a statistically stationary state without residual transients or numerical drift.
  3. [Abstract] Abstract: the equivalence between the moving-frame inlet velocity/viscosity adjustment and the lab-frame settling problem is assumed but not verified; no demonstration is given that the uniform inlet condition and any retuning preserve wake interactions and collective effects at long times.
minor comments (1)
  1. [Abstract] The abstract would benefit from explicitly stating the Galileo number, density ratio, and other non-dimensional parameters for both the single- and many-particle examples to allow direct comparison with prior literature.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful review and constructive comments on our manuscript. We address each major comment point by point below. We agree that additional quantitative support will strengthen the validation of the methodology and propose specific revisions to the manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that collective effects observed in prior work are recovered in the many-particle case is not supported by any quantitative comparison (e.g., particle velocity distributions, fluid fluctuation spectra, or cluster statistics) to reference data or an equivalent lab-frame run. This is load-bearing for validating that the inlet adjustment preserves the physics.

    Authors: We acknowledge that the manuscript currently presents the recovery of collective effects primarily through qualitative consistency with prior observations at the chosen parametric point. To strengthen this, the revised manuscript will include quantitative comparisons, such as particle velocity distributions and basic cluster statistics extracted from the simulation data, benchmarked against the referenced prior work. This will explicitly demonstrate preservation of the key non-dimensional groups and collective dynamics under the moving-frame inlet adjustment. revision: yes

  2. Referee: [Abstract] Abstract: the assertion of a converged, free-of-corrections interval of ~600 D/U_g lacks any supporting convergence study, error metrics, or time-series analysis to demonstrate that the simulation has reached a statistically stationary state without residual transients or numerical drift.

    Authors: The current text asserts convergence over this interval based on the observed stabilization of bulk quantities, but we agree that explicit supporting analysis is required. In the revision, we will add time-series plots of average particle settling velocity and fluid fluctuation intensity, along with error metrics (e.g., running averages and variance stabilization) to confirm that the simulation reaches a statistically stationary regime after an initial transient, with no residual drift over the reported interval. revision: yes

  3. Referee: [Abstract] Abstract: the equivalence between the moving-frame inlet velocity/viscosity adjustment and the lab-frame settling problem is assumed but not verified; no demonstration is given that the uniform inlet condition and any retuning preserve wake interactions and collective effects at long times.

    Authors: Equivalence for the single-particle case is verified by exact matching of the steady settling velocity to the analytical value. For the many-particle case, the formulation preserves relative velocities, Reynolds number, and Froude number by construction through the inlet adjustment, ensuring wake interactions remain unaltered. We will expand the revised manuscript with a short dedicated subsection showing that short-time flow structures and particle-pair statistics match expectations from equivalent lab-frame configurations, thereby confirming that long-time collective effects are not affected by the uniform inlet condition. revision: yes

Circularity Check

0 steps flagged

No circularity: methodology relies on explicit parameter adjustment without self-referential derivation

full rationale

The paper describes a moving-frame simulation technique that adjusts inlet velocity (or viscosity) and time step to enable long-time integration of settling particles. No derivation chain, first-principles result, or prediction is presented that reduces by construction to its own inputs or fitted quantities. The central claim of physical equivalence is an explicit modeling assumption of the boundary-condition change rather than a derived outcome; no equations or uniqueness theorems are invoked that loop back to self-citation or ansatz smuggling. The work is self-contained as a methodological description implementable in existing solvers, with no load-bearing self-references or renaming of known results.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The method rests on the assumption that a moving reference frame with adjusted inlet velocity preserves the physics of gravity-driven settling. It uses standard incompressible Navier-Stokes solvers with no new entities postulated.

free parameters (2)
  • input velocity (or viscosity)
    Adjusted to impose the correct relative flow in the moving frame for the target settling regime.
  • time step
    Updated to maintain numerical stability under the new boundary conditions.
axioms (1)
  • domain assumption The flow field in the moving reference frame is dynamically equivalent to the inertial-frame settling problem when inlet velocity is properly chosen.
    Invoked to justify that the configuration reproduces physical settling without periodicity artifacts.

pith-pipeline@v0.9.0 · 5528 in / 1323 out tokens · 32011 ms · 2026-05-09T22:39:26.770976+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

117 extracted references · 117 canonical work pages

  1. [1]

    Richardson and W.N

    J.F. Richardson and W.N. Zaki. The sedimentation of a suspension of uniform spheres under conditions of viscous flow. Chem. Eng. Sci. 1954. doi:https://doi.org/10.1016/0009-2509(54)85015-9

    work page doi:10.1016/0009-2509(54)85015-9 1954

  2. [2]

    Maxey,M. R. and Riley,J. J. , title =. Phys. Fluids , volume =. 1983 , doi =

    work page 1983

  3. [3]

    Interaction between particle clusters and particle--induced turbulence , author=. 2002. doi:10.1016/S0142-727X(02)00159-5

    work page doi:10.1016/s0142-727x(02)00159-5 2002

  4. [4]

    Influence of particle rotation on the interaction between particle clusters and particle--induced turbulence , author=. Int. J. Heat Fluid Flow , volume=. 2004 , publisher=

    work page 2004

  5. [5]

    Jenny, M. and Du. Instabilities and transition of a sphere falling or ascending freely in a. J. Fluid Mech. , volume=

    work page

  6. [6]

    , journal=

    Uhlmann, M. , journal=. An immersed boundary method with direct forcing for the simulation of particulate flows. 2005. doi:10.1016/j.jcp.2005.03.017

    work page doi:10.1016/j.jcp.2005.03.017 2005

  7. [7]

    2011 , publisher=

    Multiphase flows with droplets and particles , author=. 2011 , publisher=

    work page 2011

  8. [8]

    and Doychev, T

    Uhlmann, M. and Doychev, T. , journal=. Sedimentation of a dilute suspension of rigid spheres at intermediate. doi:10.1017/jfm.2014.330 , year=

    work page doi:10.1017/jfm.2014.330 2014

  9. [9]

    The motion of a single heavy sphere in ambient fluid: a benchmark for interface--resolved particulate flow simulations with significant relative velocities , author=. 2014. doi:https://doi.org/10.1016/j.ijmultiphaseflow.2013.10.010

    work page doi:10.1016/j.ijmultiphaseflow.2013.10.010 2014

  10. [10]

    , year =

    Doychev, T. , year =. The dynamics of finite--size settling particles , doi =

    work page

  11. [11]

    Numerical study of the sedimentation of spheroidal particles , author=. Int. J. Multiph. Flow , volume=

    work page

  12. [12]

    de Tullio, G

    M.D. A moving--least--squares immersed boundary method for simulating the fluid–-structure interaction of elastic bodies with arbitrary thickness , journal =. 2016 , issn =. doi:https://doi.org/10.1016/j.jcp.2016.08.020 , url =

    work page doi:10.1016/j.jcp.2016.08.020 2016

  13. [13]

    and Picano, F

    Fornari, W. and Picano, F. and Brandt, L. , year=. Sedimentation of finite-size spheres in quiescent and turbulent environments , volume=. doi:10.1017/jfm.2015.698 , journal=

    work page doi:10.1017/jfm.2015.698 2015

  14. [14]

    Columnar structure formation of a dilute suspension of settling spherical particles in a quiescent fluid , author =. Phys. Rev. Fluids , volume =. 2016 , doi =

    work page 2016

  15. [15]

    and Ardekani, M

    Fornari, W. and Ardekani, M. N. and Brandt, L. , journal=. Clustering and increased settling speed of oblate particles at finite

    work page

  16. [16]

    Uhlmann and J

    M. Uhlmann and J. Derksen and A. Wachs and L.-P. Wang and M. Moriche , title =. Modelling approaches and computational methods for particle-laden turbulent flows , publisher =

    work page

  17. [17]

    2023 , publisher=

    Modeling approaches and computational methods for particle-laden turbulent flows , editor=. 2023 , publisher=

    work page 2023

  18. [18]

    A single oblate spheroid settling in unbounded ambient fluid: a benchmark for simulations in steady and unsteady wake regimes , author=. Int. J. Multiph. Flow , volume=. 2021 , doi=

    work page 2021

  19. [19]

    Sedimentation of inertial monodisperse suspensions of cubes and spheres , author =. Phys. Rev. Fluids , volume =. 2021 , doi =

    work page 2021

  20. [20]

    and Xu, X

    Lu, J. and Xu, X. and Zhong, S. and Ni, R. and Tryggvason, G. , journal=. The dynamics of suspensions of prolate spheroidal particles --

    work page

  21. [21]

    On the clustering of low-aspect-ratio oblate spheroids settling in ambient fluid , author=. J. Fluid Mech. , volume=

    work page

  22. [22]

    Prolate spheroids settling in a quiescent fluid: clustering, microstructures and collisions , author=. J. Fluid Mech. , volume=

    work page

  23. [23]

    Numerical methods for multiphase flows , author=. Int. J. Multiph. Flow , volume=

    work page

  24. [24]

    Zaidi, A. A. and Tsuji, T. and Tanaka, T. , journal=. Direct numerical simulation of finite sized particles settling for high. 2014 , doi =

    work page 2014

  25. [25]

    Shape effects on the local dynamics of suspensions of spheroidal particles , author=. Phys. Fluids , volume=

    work page

  26. [26]

    Low efficiency of large volcanic eruptions in transporting very fine ash into the atmosphere , author=. Sci. Rep. , volume=

    work page

  27. [27]

    Wake-induced oscillatory paths of bodies freely rising or falling in fluids , author=. Annu. Rev. Fluid Mech. , volume=

    work page

  28. [28]

    Powder Tech

    On the motion of non-spherical particles at high Reynolds number , author=. Powder Tech. , volume=

    work page

  29. [29]

    Orientations and aspect ratios of falling snow , author=. Geophys. Res. Lett. , volume=

    work page

  30. [30]

    Shape effects on turbulent modulation by large nearly neutrally buoyant particles , author=. J. Fluid Mech. , volume=

    work page

  31. [31]

    Acta Mech

    Turbulence modulation by large ellipsoidal particles: concentration effects , author=. Acta Mech. , volume=

    work page

  32. [32]

    Shape-dependence of particle rotation in isotropic turbulence , author=. Phys. Fluids , volume=

    work page

  33. [33]

    Measurements of the coupling between the tumbling of rods and the velocity gradient tensor in turbulence , author=. J. Fluid Mech. , volume=

    work page

  34. [34]

    and Girolami, L

    Wachs, A. and Girolami, L. and Vinay, G. and Ferrer, G. , journal=

    work page

  35. [35]

    A fast procedure for computing the distance between complex objects in three-dimensional space , author=. IEEE J. Robotics Automation , volume=

    work page

  36. [36]

    2007 , publisher=

    Atmospheric and space flight dynamics , author=. 2007 , publisher=

    work page 2007

  37. [37]

    Marcus, G. G. and Parsa, S. and Kramel, S. and Ni, R. and Voth, G.A. , journal=. Measurements of the solid-body rotation of anisotropic particles in

    work page

  38. [38]

    Ellipsoidal particles transport and deposition in turbulent channel flows , author=. Int. J. Multiphase Flow , volume=

    work page

  39. [39]

    Fluctuations and instability in sedimentation , author=. Annu. Rev. Fluid. Mech. , volume=

    work page

  40. [40]

    A physical introduction to suspension dynamics , author=

    work page

  41. [41]

    The instability of a dispersion of sedimenting spheroids , author=. J. Fluid Mech. , volume=

    work page

  42. [42]

    Collective effects in settling of spheroids under steady-state sedimentation , author=. Phys. Rev. Lett. , volume=

    work page

  43. [43]

    and Koch, D

    Shin, M. and Koch, D. L. and Subramanian, G. , journal=. Structure and dynamics of dilute suspensions of finite-

    work page

  44. [44]

    Orientation statistics and settling velocity of ellipsoids in decaying turbulence , author=. Atmos. Res. , volume=

    work page

  45. [45]

    Qi, G. Q. and Nathan, G. J. and Kelso, R. M. , journal=

    work page

  46. [46]

    Powder Tech

    The influence of aspect ratio on distributions of settling velocities and orientations of long fibres , author=. Powder Tech. , volume=

    work page

  47. [47]

    and Wu, T.-H

    Liu, Y. and Wu, T.-H. and Guo, R.-S. and Lee, Y.-H. and Qi, D. , journal=. Dynamics of sedimentation of flexible fibers in moderate

    work page

  48. [48]

    Orientation and rotation of inertial disk particles in wall turbulence , author=. J. Fluid Mech. , volume=

    work page

  49. [49]

    Computation of particle settling speed and orientation distribution in suspensions of prolate spheroids , author=. J. Eng. Math. , volume=

    work page

  50. [50]

    Simulation of finite-size fibers in turbulent channel flows , author=. Phys. Rev. E , volume=

    work page

  51. [51]

    Drag reduction in turbulent channel flow laden with finite-size oblate spheroids , author=

    work page

  52. [52]

    Finite-length effects on dynamical behaviour of rod-like particles in wall-bounded turbulent flow , author=. Int. J. Multiphase Flow , volume=

    work page

  53. [53]

    Rotation of nonspherical particles in turbulent channel flow , author=. Phys. Rev. Lett. , volume=

    work page

  54. [54]

    Fall velocity of multi-shaped clasts , author=. J. Volcanology Geothermal Res. , volume=

    work page

  55. [55]

    Acta Mech

    Anisotropic particles in turbulence: status and outlook , author=. Acta Mech. , volume=

    work page

  56. [56]

    Dynamics of prolate ellipsoidal particles in a turbulent channel flow , author=. Phys. Fluids , volume=

    work page

  57. [57]

    On the orientation of ellipsoidal particles in a turbulent shear flow , author=. Int. J. Multiphase Flow , volume=

    work page

  58. [58]

    Orientation, distribution, and deposition of elongated, inertial fibers in turbulent channel flow , author=. Phys. Fluids , volume=

    work page

  59. [59]

    Acta Mech

    Rotation statistics of fibers in wall shear turbulence , author=. Acta Mech. , volume=

    work page

  60. [60]

    Rotation rate of rods in turbulent fluid flow , author=. Phys. Rev. Lett. , volume=

    work page

  61. [61]

    Nonlinear disturbance growth during sedimentation in dilute fibre suspensions , author=. J. Fluid Mech. , volume=

    work page

  62. [62]

    The aerodynamic behaviour of cylindrical and spheroidal particles when settling under gravity , author=. J. Aerosol Sci. , volume=

    work page

  63. [63]

    Sedimentary Geology , volume=

    Grains in motion: A review , author=. Sedimentary Geology , volume=. 2005 , publisher=

    work page 2005

  64. [64]

    Dedicated vertical wind tunnel for the study of sedimentation of non-spherical particles , author=. Rev. Sci. Instrum. , volume=

    work page

  65. [65]

    Powder Tech

    On the drag of freely falling non-spherical particles , author=. Powder Tech. , volume=

    work page

  66. [66]

    Orientation dynamics of small, triaxial--ellipsoidal particles in isotropic turbulence , author=. J. Fluid Mech. , volume=

    work page

  67. [67]

    Tumbling of small axisymmetric particles in random and turbulent flows , author=. Phys. Rev. Lett. , volume=

    work page

  68. [68]

    Orientation statistics of small particles in turbulence , author=. New J. Phys. , volume=

    work page

  69. [69]

    Effect of solid body degrees of freedom on the path instabilities of freely falling or rising flat cylinders , author=. J. Fluids Struct. , volume=

    work page

  70. [70]

    Numerical simulation of the dynamics of freely falling discs , author=. Phys. Fluids , volume=

    work page

  71. [71]

    Effect of inclination on the transition scenario in the wake of fixed disks and flat cylinders , author=. J. Fluid Mech. , volume=

    work page

  72. [72]

    and Lee, C

    Zhong, H. and Lee, C. and Su, Z. and Chen, S. and Zhou, M. and Wu, J. , journal=

    work page

  73. [73]

    and Su, Z

    Lee, C. and Su, Z. and Zhong, H. and Chen, S. and Zhou, M. and Wu, J. , journal=

    work page

  74. [74]

    Influence of aspect ratio on tumbling plates , author=. J. Fluid Mech. , volume=

    work page

  75. [75]

    Simulation Methods for Particulate Flows and Concentrated Suspensions , author=. Annu. Rev. Fluid Mech. , volume=

    work page

  76. [76]

    2016 , school=

    Path instabilities of spheres, spheroids and bubbles , author=. 2016 , school=

    work page 2016

  77. [77]

    Chrust, M. , year=. Etude num

    work page

  78. [78]

    Free falling and rising of spherical and angular particles , author=. Phys. Fluids , volume=

    work page

  79. [79]

    Oscillatory motion and wake instability of freely rising axisymmetric bodies , author=. J. Fluid Mech. , volume=

    work page

  80. [80]

    Dynamics of axisymmetric bodies rising along a zigzag path , author=. J. Fluid Mech. , volume=

    work page

Showing first 80 references.