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

This study utilizes direct numerical simulations and a population balance equation framework to reveal that in turbulent channel flows, cohesive aggregates undergo a size-dependent circulation where larger aggregates form in the channel center and migrate to the wall to break, while smaller fragments are transported away to grow, thereby balancing local aggregation and breakup dynamics.

The Gist

Imagine a busy, chaotic highway (the turbulent channel flow) filled with thousands of tiny, sticky marbles (the cohesive particles). In this highway, the wind is so strong that it constantly throws the marbles around. But because the marbles are sticky, when they bump into each other, they sometimes stick together to form bigger clumps. However, the wind is also strong enough to smash these clumps apart again.

This paper is a detailed investigation into how these sticky marbles behave when they are trapped between two walls (like a pipe or a riverbed), rather than floating freely in open space.

Here is the story of what the researchers found, broken down into simple concepts:

1. The Big Picture: A Never-Ending Cycle

Think of the sticky marbles as a population of "families."

• Aggregation: When two small families bump into each other, they merge into a bigger family.
• Breakup: When a big family gets hit too hard by the wind, it shatters back into smaller families.

The researchers ran computer simulations to see what happens when you make the marbles stickier and stickier. They found that the system doesn't just get bigger or smaller; it settles into a steady state where the total number of "people" (mass) stays the same, but the size of the families keeps changing.

2. The "Global" Balance vs. The "Local" Chaos

If you were to stand on a bridge looking down at the entire highway, you would see a perfect balance: for every big family that breaks apart, a small family merges to replace it. The total number of marbles is constant.

However, the researchers discovered that this balance is a lie if you look at just one specific spot on the road.

• Near the walls: It's a demolition zone. The wind is turbulent and the walls create friction. Big clumps crash into the wall and shatter into tiny pieces. Here, you have a surplus of small marbles and a shortage of big ones.
• In the center of the channel: It's a construction zone. The wind is smoother here. Small marbles drift in, bump into each other, and stick together to form massive clumps. Here, you have a surplus of big clumps and a shortage of small ones.

So, locally, the system is out of balance. The walls are destroying families, and the center is building them.

3. The Secret Ingredient: The "Elevator" (Transport)

If the walls are destroying families and the center is building them, why doesn't the whole highway run out of big clumps or fill up with tiny ones?

The answer is movement. The researchers found that the clumps act like passengers on a giant, invisible elevator system that moves them up and down the channel.

• The Upward Journey: Tiny, lonely marbles (created by the smashing near the walls) get swept away from the wall and drift toward the center of the channel.
• The Growth Phase: Once they reach the calm center, they have a party. They bump into each other and grow into huge, heavy clumps.
• The Downward Journey: These heavy clumps are now too big to stay in the center. They get pushed back toward the walls.
• The Crash: When they hit the turbulent zone near the wall, they smash apart, creating a fresh batch of tiny marbles, and the cycle starts all over again.

4. The "Statistical Circulation"

The most exciting discovery is that this isn't just a random mess; it's a perfectly organized loop.

Imagine a conveyor belt in a factory:

1. Raw Materials (Small particles) are dumped near the wall.
2. They get shipped to the factory floor (the center) where they are assembled into products (large aggregates).
3. The finished products are shipped back to the loading dock (the wall).
4. At the dock, the products are recycled (broken down) back into raw materials.

This loop happens continuously. The "size" of the particle determines where it goes. Small ones go to the center to grow; big ones go to the wall to break.

Why Does This Matter?

This isn't just about sticky marbles. This process happens everywhere in nature and engineering:

• Rivers and Oceans: Mud and silt (which are sticky) clump together, sink, and get stirred up again. Understanding this helps us predict how sediment moves and how rivers change shape.
• Pollution: Oil spills or plastic waste often clump together. Knowing how they break and reform helps in cleanup efforts.
• Industry: In factories making paint, ink, or medicine, particles often clump. If you don't understand this cycle, your product might be the wrong size or texture.

The Takeaway

The paper teaches us that in a turbulent environment with sticky particles, you cannot understand the whole system just by looking at the average. You have to look at the journey of the particles.

The system works like a giant, self-sustaining cycle: Break near the wall $\rightarrow$ Drift to the center $\rightarrow$ Grow in the center $\rightarrow$ Drift back to the wall $\rightarrow$ Break again. It is a beautiful, chaotic dance of destruction and creation that keeps the system in perfect equilibrium.

Technical Summary

1. Problem Statement

Cohesive particles suspended in turbulent flows undergo complex processes of aggregation (flocculation), breakup, and restructuring. While significant research exists on these phenomena in homogeneous isotropic turbulence (HIT), there is a critical gap in understanding them in wall-bounded flows (such as channels or pipes). In wall-bounded flows, turbulence intensity and shear rates are highly inhomogeneous, particularly near the walls.

Previous studies have often relied on simplified parameterizations or focused only on initial transients. A key unresolved question is whether the rates of aggregation and breakup balance locally in the wall-normal direction, or if spatial imbalances exist that drive net transport of aggregates. This study aims to resolve these dynamics using high-fidelity simulations to understand the statistical circulation of aggregate mass in both physical space and size space.

2. Methodology

The authors employed Particle-Resolved Direct Numerical Simulations (PR-DNS) to investigate a fully developed turbulent channel flow laden with cohesive particles.

• Simulation Setup:

  • Domain: A channel of size $(6h, 2h, 3h)$ with a Reynolds number $Re = 2800$ (based on half-height $h$).
  • Particles: $N_p = 3841$ cohesive particles with a volume fraction $\phi = 0.015$.
  • Cohesion Model: Based on Vowinckel et al. (2019), characterized by the cohesive number $Co = \max|F_{coh}| / (\rho_f U^2 d^2)$.
  • Cases: Five simulations were performed with increasing cohesive strength: $Co = \{0, 0.24, 0.48, 1.0, 1.5\}$.
  • Initial Conditions: Particles were randomly distributed in a fully developed single-phase flow.

• Analytical Framework:

  • A Population Balance Equation (PBE) framework was applied to the simulation data.
  • The PBE tracks the evolution of aggregate number density $f(n, x, t)$, where $n$ is the number of constituent particles in an aggregate.
  • The equation includes terms for advection (transport), aggregation (source/sink), and breakup (source/sink).
  • The study analyzed the PBE in two ways:
    1. Bulk-averaged: Integrated over the entire domain to check global mass conservation.
    2. Plane-averaged: Averaged over streamwise and spanwise directions to resolve variations in the wall-normal direction ($y$).

3. Key Contributions

1. First PR-DNS of Resolved Cohesive Particles in Wall-Bounded Flow: The study provides the first fully resolved investigation of aggregation and breakup in a turbulent channel, moving beyond point-particle approximations.
2. Local Imbalance Discovery: The authors demonstrate that while aggregation and breakup balance globally (over the whole channel), they do not balance locally in the wall-normal direction.
3. Size-Dependent Transport Mechanism: The study quantifies the role of wall-normal advection in closing the PBE locally, revealing that aggregates of different sizes are transported differently by the flow.
4. Statistical Circulation Cycle: The work identifies a closed statistical cycle of aggregate mass involving production near walls, growth in the channel center, and breakup near walls.

4. Key Results

A. Global vs. Local Balance

• Global Balance: When integrated over the full channel height, the net source term ($Q = \text{Aggregation} + \text{Breakup}$) is zero for all aggregate sizes. This confirms that aggregation and breakup are globally balanced.
• Local Imbalance: When analyzed locally in the wall-normal direction ($y$), significant imbalances exist ($Q(n, y) \neq 0$).
  • Near the wall ($y \approx 0$): Breakup dominates for larger aggregates ($n > 3$), creating a net depletion of large aggregates and a net production of monomers ($n=1$).
  • Channel center: Aggregation dominates for intermediate and large sizes, leading to net growth.

B. Size-Dependent Dynamics

• Monomers ($n=1$): Produced near the wall via breakup and consumed in the center via aggregation.
• Large Aggregates ($n > 3$): Formed in the channel center via aggregation and destroyed near the wall via high shear and turbulence.
• Intermediate Sizes: Exhibit complex behavior depending on $Co$. At high cohesion, they behave more like monomers (unable to break in the center); at lower cohesion, they behave like larger aggregates.

C. The Statistical Circulation Cycle

By analyzing the $(n, y)$ space (size vs. wall-normal position), the authors identified a four-region cyclic process:

1. Region 1 (Near Wall, Small $n$): Aggregates are broken down, producing monomers/small clusters.
2. Region 2 (Moving Inward, Small $n$): These small aggregates are advected toward the channel center.
3. Region 3 (Channel Center, Growing $n$): In the lower turbulence intensity of the center, aggregates grow via collisions (aggregation) into larger sizes.
4. Region 4 (Moving Outward, Large $n$): Large aggregates are advected back toward the wall, where high shear stresses cause them to break, returning mass to Region 1.

This cycle creates a log-normal distribution of aggregate sizes for $Co \geq 0.48$, a feature previously observed but not mechanistically explained in this specific flow context.

D. Role of Advection

The study proves that the wall-normal advective flux ($\partial_y (v_n f)$) is the critical term required to close the PBE locally. Without this transport term, the local production/depletion rates would violate mass conservation. The direction and magnitude of this flux are size-dependent, driving the circulation.

5. Significance

• Modeling Implications: The findings challenge the validity of local equilibrium assumptions often used in low-order models (like standard PBEs without transport terms) for wall-bounded flows. Models must account for size-dependent transport to accurately predict aggregate size distributions.
• Environmental and Engineering Applications: This work provides a fundamental understanding of sediment transport, wastewater treatment, and industrial mixing processes where cohesive particles (like clay or flocs) are present in turbulent channels.
• Mechanistic Insight: It resolves the "where" and "how" of aggregate lifecycle in shear flows, showing that the channel acts as a reactor where particles continuously cycle between growth and destruction zones rather than reaching a static equilibrium state.

In conclusion, the paper establishes that the dynamics of cohesive particles in turbulent channels are governed by a spatially separated, size-dependent circulation loop, driven by the interplay of local turbulence shear, aggregation, and advection.