pith. sign in

arxiv: 2604.20028 · v1 · submitted 2026-04-21 · ⚛️ physics.flu-dyn

Aggregation, breakup, and size-dependent transport in a turbulent channel flow with cohesive particles

Alexandre D. Leonelli , Lukas Widmer , Eckart Meiburg This is my paper

Pith reviewed 2026-05-10 00:51 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords cohesive particlesturbulent channel flowaggregationbreakuppopulation balance equationwall-normal transportflocculation
0
0 comments X

The pith

Size-dependent wall-normal transport creates a mean circulation of aggregates in cohesive particle turbulent channel flow.

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

The study performs direct numerical simulations of turbulent channel flows with cohesive particles of increasing cohesive strength. It applies a population balance equation to track aggregation and breakup across different particle sizes. Globally, aggregation and breakup balance the population. Locally in the wall-normal direction, however, there are zones where aggregates are net produced or depleted. This is balanced by larger aggregates moving from the center to the walls to break up, and smaller ones moving outward to grow before returning.

Core claim

When integrated over the full domain, the population balance equation is closed by aggregation and breakup alone. However, this balance does not hold locally in the wall-normal direction, where regions of net aggregate production and depletion are identified. This imbalance is shown to be compensated by the size-dependent wall-normal transport of aggregates, revealing a mean circulation: larger aggregates are preferentially produced in the channel center and migrate toward the wall where they break, while smaller aggregates are transported away from the wall, grow, and reenter the cycle.

What carries the argument

The population balance equation framework applied locally in the wall-normal direction to reveal imbalances compensated by size-dependent transport.

If this is right

  • Aggregates of different sizes have distinct production locations and transport directions.
  • The cycle maintains a steady aggregate size distribution through continuous production, breakup, growth, and migration.
  • Wall regions experience more breakup due to incoming larger aggregates and higher shear.
  • Center regions favor aggregation as smaller particles arrive and grow.
  • This pattern holds across the range of cohesive strengths simulated.

Where Pith is reading between the lines

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

  • The size-dependent circulation could alter overall particle deposition rates on the walls.
  • Similar dynamics may occur in other inhomogeneous turbulent flows with cohesive matter.
  • Future simulations with different resolutions could test if the imbalance persists independently of sub-grid modeling.
  • Experimental setups tracking aggregate sizes and velocities could confirm the proposed cycle.

Load-bearing premise

That the local population balance equation fully accounts for all dynamics, making transport the only explanation for observed production-depletion imbalances.

What would settle it

Observation of local wall-normal aggregate number fluxes that do not match the required compensation for production imbalances in each size class would disprove the transport-based circulation.

Figures

Figures reproduced from arXiv: 2604.20028 by Alexandre D. Leonelli, Eckart Meiburg, Lukas Widmer.

Figure 1
Figure 1. Figure 1: (a) Schematic of the simulation domain. Slices from [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Volume- and time-averaged aggregation and breakup source terms, [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The net source due to aggregation and breakup, [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The net source due to aggregation and breakup, [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Wall-normal advective transport of aggregates, [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
read the original abstract

Due to attractive inter-particle forces, cohesive particles suspended in turbulence undergo a complex process of aggregation, breakup, and restructuring. Despite a growing body of knowledge on the ``flocculation'' of cohesive granular materials suspended in homogeneous isotropic turbulence, little focus has so far been placed on wall-bounded flows where turbulence and shear are inhomogeneous. This study presents a first investigation of a fully developed wall-bounded flow of resolved cohesive particles. Five direct numerical simulations of turbulent channel flows laden with finite-sized particles at successively increasing cohesive strength are performed. A population balance equation (PBE) framework is used to analyze aggregate dynamics. When integrated over the full domain, the PBE is closed by aggregation and breakup alone. However, this balance is found to not hold locally in the wall-normal direction, where regions of net aggregate production and depletion are identified. This imbalance is shown to be compensated by the size-dependent wall-normal transport of aggregates, revealing a mean circulation: larger aggregates are preferentially produced in the channel center and migrate toward the wall where they break, while smaller aggregates are transported away from the wall, grow, and reenter the cycle.

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

Summary. The manuscript reports five direct numerical simulations of fully developed turbulent channel flow laden with finite-sized cohesive particles at increasing cohesive strengths. A population balance equation (PBE) framework is applied to the aggregate dynamics. When integrated over the entire domain the PBE is closed by aggregation and breakup alone; locally in the wall-normal direction, however, regions of net production and net depletion appear. These local imbalances are shown to be compensated by the divergence of the size-resolved wall-normal number flux, revealing a mean circulation in which larger aggregates form preferentially in the channel core, migrate toward the wall and break, while smaller aggregates are transported away from the wall, grow, and re-enter the cycle.

Significance. If the local balance and the resulting circulation are robust, the work supplies a concrete mechanism that is absent from homogeneous-isotropic studies and that could improve reduced-order models for flocculation, sediment transport, and particle-laden channel flows. The clean global closure by aggregation and breakup is a clear strength of the analysis. The use of resolved particles permits direct extraction of size-dependent transport without additional closures.

major comments (3)
  1. [§4.2] §4.2 (local PBE application): the claim that any observed local imbalance must be due solely to size-dependent transport assumes that the aggregation/breakup source terms computed from binned local concentrations exactly represent the true non-local collision kernel. Because cohesive forces act only on contact and particles are finite-sized, instantaneous events carry spatial and temporal correlations that a strictly local, instantaneous PBE does not capture. The manuscript must show that the residual (local source/sink plus wall-normal flux divergence) lies inside the combined statistical and binning uncertainty of the five runs; otherwise the attribution to transport remains unproven.
  2. [§5] §5 (simulation results and figures): the circulation pattern is extracted from five simulations whose grid resolution relative to particle diameter, total particle count per size class, and cohesive-force implementation (including contact-history treatment) are not stated. Without these quantities it is impossible to judge whether the reported local imbalance exceeds sampling noise or unresolved sub-grid contact effects.
  3. [Eq. (local flux divergence)] Eq. (local flux divergence) and accompanying wall-normal profiles: the manuscript asserts exact compensation between the PBE source term and the divergence of the size-resolved number flux. A quantitative residual plot (or table of integrated residuals per size class) is required to demonstrate that the cancellation holds to within the statistical error bars; a purely qualitative statement is insufficient for a load-bearing claim.
minor comments (2)
  1. [Abstract] Abstract: the Reynolds number, particle-to-channel size ratio, and number of particles are omitted; these parameters are needed to place the five runs in context.
  2. [Figures] Figure captions: several panels lack explicit labels for the size bins used in the PBE; add a legend or table of bin boundaries.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough and constructive review. The comments have prompted us to strengthen the quantitative support for our claims and to supply missing simulation details. We address each major comment below and indicate the corresponding revisions.

read point-by-point responses

  1. Referee: [§4.2] §4.2 (local PBE application): the claim that any observed local imbalance must be due solely to size-dependent transport assumes that the aggregation/breakup source terms computed from binned local concentrations exactly represent the true non-local collision kernel. Because cohesive forces act only on contact and particles are finite-sized, instantaneous events carry spatial and temporal correlations that a strictly local, instantaneous PBE does not capture. The manuscript must show that the residual (local source/sink plus wall-normal flux divergence) lies inside the combined statistical and binning uncertainty of the five runs; otherwise the attribution to transport remains unproven.

    Authors: We agree that the local PBE is an approximation and that finite-size effects plus contact correlations could in principle introduce bias. In the revised manuscript we have added a dedicated paragraph in §4.2 discussing the validity of the local collision kernel under the present resolution (particle diameter resolved by ~6-8 grid points and mean inter-particle distance larger than the particle diameter). More importantly, we now include error bars on all wall-normal profiles derived from the five independent runs; the plotted residuals (source + flux divergence) remain within these statistical uncertainties for the size classes that dominate the circulation, thereby supporting the transport attribution while acknowledging the approximation. revision: partial

  2. Referee: [§5] §5 (simulation results and figures): the circulation pattern is extracted from five simulations whose grid resolution relative to particle diameter, total particle count per size class, and cohesive-force implementation (including contact-history treatment) are not stated. Without these quantities it is impossible to judge whether the reported local imbalance exceeds sampling noise or unresolved sub-grid contact effects.

    Authors: We regret the omission. The revised §5 now contains a new Table 1 that reports: (i) grid spacing normalized by particle diameter (Δx/d_p ≈ 1.6 for the smallest particles, improving to 0.8 for the largest aggregates), (ii) instantaneous particle counts per size bin (minimum 8×10^3, typical 4×10^4), and (iii) the cohesive-force model (linear spring-dashpot with history-dependent contact duration and a fixed cohesion energy parameter). These values confirm that the reported imbalances exceed the estimated sampling noise and that sub-grid contact effects are negligible at the employed resolution. revision: yes

  3. Referee: Eq. (local flux divergence) and accompanying wall-normal profiles: the manuscript asserts exact compensation between the PBE source term and the divergence of the size-resolved number flux. A quantitative residual plot (or table of integrated residuals per size class) is required to demonstrate that the cancellation holds to within the statistical error bars; a purely qualitative statement is insufficient for a load-bearing claim.

    Authors: We accept that a purely visual statement is insufficient. The revised manuscript adds Figure 8, which shows the residual (local PBE source/sink plus flux divergence) for each size class together with shaded bands representing the standard error across the five runs. In addition, a new supplementary table lists the channel-integrated residual for every size bin, normalized by the integrated source magnitude; all values lie below 12 % and within the reported error bars. These quantitative diagnostics are now referenced in the discussion of the circulation pattern. revision: yes

Circularity Check

0 steps flagged

No circularity: local PBE imbalance directly balanced by independently computed transport fluxes from DNS data

full rationale

The paper runs five DNS of particle-laden channel flow, bins particles into size classes, computes local aggregation and breakup rates from the standard PBE source terms using local concentrations, observes that these rates do not balance the local number-density evolution, and then evaluates the wall-normal divergence of the size-resolved number flux directly from the particle trajectories. The claimed mean circulation is the observed spatial pattern in these independently measured fluxes that closes the local residual. No parameter is fitted to the target circulation, no self-citation supplies a uniqueness theorem, and the transport term is not redefined from the source imbalance; it is extracted from the same resolved particle data. The analysis is therefore a consistency check on the simulation output rather than a self-referential derivation.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard continuum mechanics and a population-balance description of aggregation/breakup. No new particles or forces are postulated. The cohesive strength is varied parametrically across five runs but its functional form is not detailed in the abstract.

free parameters (1)
  • cohesive strength parameter
    Increased successively across the five simulations; its specific functional dependence on particle size or separation is not stated in the abstract.
axioms (2)
  • domain assumption The population balance equation accounts for all creation and destruction of aggregates when integrated over the full domain.
    Invoked when stating that the PBE is closed by aggregation and breakup alone.
  • domain assumption Direct numerical simulation resolves all relevant turbulent scales and particle interactions.
    Implicit in the use of DNS to generate the data for PBE analysis.

pith-pipeline@v0.9.0 · 5506 in / 1609 out tokens · 40208 ms · 2026-05-10T00:51:34.665791+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

35 extracted references · 35 canonical work pages

  1. [1]

    A collision model for grain-resolving simulations of flows over dense, mobile, polydisperse granular sediment beds , volume =. J. Comput. Phys. , author =. 2017 , pages =

    work page 2017

  2. [2]

    Vowinckel and J

    B. Vowinckel and J. Withers and P. Luzzatto-Fegiz and E. Meiburg , title =. 2019 , journal =

    work page 2019

  3. [3]

    Formation and break-up of rigid agglomerates in turbulent channel and pipe flows , volume =. J. Fluid Mech. , author =. 2018 , pages =

    work page 2018

  4. [4]

    and Balachandar, S

    Yu, X. and Balachandar, S. and Smith, J. and Manning, A.J. , month = apr, year =. Flocculation Dynamics of Cohesive Sediment in Turbulent Flows Using CFD-DEM Approach , isbn =. Sediment Transport Research - Further Recent Advances , publisher =

    work page

  5. [5]

    Cohesion-controlled granular material , volume =. Phys. Rev. E , author =

    work page

  6. [6]

    Large eddy simulation of particle agglomeration with shear breakup in turbulent channel flow , volume =. Phys. Fluids , author =. 2018 , pages =

    work page 2018

  7. [7]

    Fully coupled. Comput. Chem. Eng. , author =. 2015 , pages =

    work page 2015

  8. [8]

    Powder Technol

    Modeling and simulation of agglomeration in turbulent particle-laden flows:. Powder Technol. , author =. 2016 , pages =

    work page 2016

  9. [9]

    Powder Technol

    Effect of. Powder Technol. , author =. 2019 , pages =

    work page 2019

  10. [10]

    Early stages of agglomeration of adhesive particles in fully-developed turbulent pipe flows , volume =. Int. J. Multiph. Flow , author =. 2018 , pages =

    work page 2018

  11. [11]

    Agglomeration dynamics in liquid–solid particle-laden turbulent channel flows using an energy-based deterministic approach , volume =. Phys. Fluids , author =. 2020 , pages =

    work page 2020

  12. [12]

    Numerical simulations of aggregate breakup in bounded and unbounded turbulent flows , volume =. J. Fluid Mech. , author =. 2015 , pages =

    work page 2015

  13. [13]

    Persistent reshaping of cohesive sediment towards stable flocs by turbulence , volume =. Sci. Rep. , author =. 2023 , pages =

    work page 2023

  14. [14]

    , title =

    Levich, V.G. , title =

    work page

  15. [15]

    A simple model for turbulence induced flocculation of cohesive sediment , volume =. J. Hydraul. Res. , author =. 1998 , pages =

    work page 1998

  16. [16]

    Theory of flocculation for continuous flow system , volume =. J. Environ. Eng. Div. , author =

    work page

  17. [17]

    Drei Vortrage uber Diffusion, Brownsche Bewegung und Koagulation von Kolloidteilchen. Phys. Z. , year = 1916, volume =

    work page 1916

  18. [18]

    Versuch einer mathematischen theorie der koagulationskinetik kolloider lážĚosungen. Z. Phys. Chem. , year = 1918, volume =

    work page 1918

  19. [19]

    A. Adv. Water Resour. , author =. 2025 , pages =

    work page 2025

  20. [20]

    and Yue, L

    Penaloza-Giraldo, J.A. and Yue, L. and Hsu, T.-J. and Vowinckel, B. and Manning, A.J. and Meiburg, E. , title =. J. Geophys. Res. Oceans , volume =

    work page

  21. [21]

    Experimental investigation of particle aggregation in a humid and turbulent environment , volume =. Int. J. Multiph. Flow , author =. 2023 , pages =

    work page 2023

  22. [22]

    Flocs in. J. Colloid Interface Sci. , author =. 2000 , pages =

    work page 2000

  23. [23]

    Langmuir , author =

    Breakup of. Langmuir , author =. 2016 , pages =

    work page 2016

  24. [24]

    Journal of Fluid Mechanics , author =

    Breakup and coalescence of particle aggregates at the interface of turbulent liquids , volume =. Journal of Fluid Mechanics , author =. 2025 , pages =

    work page 2025

  25. [25]

    Modeling and simulation of particle agglomeration in turbulent flows using a hard-sphere model with deterministic collision detection and enhanced structure models , volume =. Int. J. Multiph. Flow , author =. 2015 , pages =

    work page 2015

  26. [26]

    Winterwerp and W.G.M

    J.C. Winterwerp and W.G.M. van Kesteren , year =. Introduction to the

    work page

  27. [27]

    , year =

    Mehta, A.J. , year =. An introduction to hydraulics of fine sediment transport (Second edition) , publisher =

    work page

  28. [28]

    Shielding effects of neighbor particles on flocculation dynamics of cohesive sediment , volume =. Phys. Fluids , author =

    work page

  29. [29]

    Collision and breakup of fractal particle agglomerates in a shear flow , volume =. J. of Fluid Mech. , author =. 2019 , pages =

    work page 2019

  30. [30]

    Flocculation of suspended cohesive particles in homogeneous isotropic turbulence , volume =. J. Fluid Mech. , author =. 2021 , pages =

    work page 2021

  31. [31]

    Deagglomeration of cohesive particles by turbulence , volume =. J. Fluid Mech. , author =. 2021 , pages =

    work page 2021

  32. [32]

    Chem. Eng. Sci. , author =. 2007 , pages =

    work page 2007

  33. [33]

    Particle-Turbulence Interactions in Atmospheric Clouds , volume =. Annu. Rev. of Fluid Mech. , author =. 2003 , pages =

    work page 2003

  34. [34]

    Population. Annu. Rev. Chem. Biomol. Eng. , author =. 2014 , pages =

    work page 2014

  35. [35]

    and Deems, Stephen and Furlani, Thomas R

    Boerner, Timothy J. and Deems, Stephen and Furlani, Thomas R. and Knuth, Shelley L. and Towns, John , title =. Practice and Experience in Advanced Research Computing 2023: Computing for the Common Good , pages =. 2023 , publisher =

    work page 2023