A suspension of heavy Kolmogorov-size spheres suppresses the inertial cascade in homogeneous and isotropic turbulence

Ludovico Foss\`a; Marco Edoardo Rosti

REVIEW 3 major objections 5 minor 72 references

A suspension of heavy Kolmogorov-size spheres suppresses the inertial cascade in homogeneous and isotropic turbulence

T0 review · 3 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper argues that heavy spheres as large as the Kolmogorov scale can suppress the inertial energy cascade: as their mass fraction rises, the carrier-flow spectrum shifts from the classical κ^{-5/3} law to κ^{-1}, while the nonlinear tra

desk verdict A substantial PR-DNS study showing a distinct cascade-suppression regime for heavy Kolmogorov-size spheres, but the missing single-phase control at matched Re_lambda is a real gap that a good referee should push on. read the letter →

arxiv 2602.05171 v2 pith:3DAXSUBL submitted 2026-02-05 physics.flu-dyn

Ludovico Fossà , Marco Edoardo Rosti This is my paper · ORCID

classification physics.flu-dyn MSC 76F0576F6576T20

keywords particle-ladenturbulenceKolmogorov-sizeparticleshomogeneousisotropicenergycascadesuppressionspectrumstructurefunctionsparticleclusteringpreferentialconcentration

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper tries to establish that a dilute suspension of spheres whose diameter equals the Kolmogorov length can, when the particles are dense enough, dismantle the usual turbulent energy cascade. Using particle-resolved simulations of homogeneous isotropic turbulence, the authors show that as the particle-to-fluid density ratio climbs from 100 to about 1500, the energy spectrum shifts from the classical κ^{-5/3} law to a κ^{-1} regime, the nonlinear energy flux nearly vanishes, and the scale-by-scale energy balance is instead carried by the fluid–solid interaction and viscous dissipation. The second-order structure function grows logarithmically with separation rather than as a power, which the authors read as velocity decorrelation produced by the particles acting as distributed momentum sinks. The work matters because Kolmogorov-size particles—too large for standard point-particle models and traditionally less studied—are common in natural and industrial flows, so showing they can short-circuit the cascade would change how such suspensions are modelled.

What carries the argument

The load-bearing object is the scale-by-scale kinetic-energy budget of the carrier flow, split into nonlinear transfer, viscous dissipation, large-scale forcing, and the fluid–particle interaction term produced by the immersed-boundary force. The paper’s arguments hinge on tracking how the nonlinear transfer weakens and the fluid–particle interaction strengthens as the mass fraction increases, and on the accompanying change of the energy spectrum from κ^{-5/3} to κ^{-1}. The particles themselves, at fixed diameter equal to the Kolmogorov length and volume fraction 10^{-3}, are the mechanical device: their inertia converts their relative slip into a direct drain of kinetic energy across a wid

What would settle it

Run the same triply periodic setup and forcing with no particles, tuned to the same Taylor-scale Reynolds number (about 107) and the same integral scale as the heaviest laden case. If that single-phase flow also develops a κ^{-1} spectrum, a logarithmic second-order structure function, and a nearly zero nonlinear energy flux, the paper’s central claim collapses. Conversely, if the unladen low-Reynolds flow keeps the κ^{-5/3} spectral shape and substantial nonlinear transfer, the particles are the cause.

Watch

Extended reading notes

Core claim

The central claim is that Kolmogorov-size spheres with high enough inertia act as a distributed energy sink that bypasses the inertial cascade. At mass fractions around 0.6, the kinetic energy spectrum of the carrier fluid follows E(k) ~ κ^{-1} over roughly a decade of wavenumbers (2 ≤ κ/κ_L < 70), instead of κ^{-5/3}; the spectral energy budget shows the nonlinear transfer term becoming negligible while the fluid-particle forcing term and viscous dissipation dominate; and the longitudinal structure function behaves as S_2 ~ log(r/η) for separations beyond the particle diameter, signalling a loss of velocity correlation. The paper also claims that increased particle inertia shifts small-scal

Load-bearing premise

The interpretation that particle inertia causes the cascade suppression assumes that the much lower Taylor-scale Reynolds number of the heaviest laden case (about 107, versus 149 unladen) does not by itself produce the same κ^{-1} spectrum and vanishing nonlinear flux; no single-phase simulation was run at the reduced Reynolds number.

Editorial extensions

If this is right

If heavy Kolmogorov-size spheres suppress the cascade, suspensions carrying such particles cannot be treated as a simple additive drag: the particles actively reshape the spectral energy distribution.
The collapse of the nonlinear flux implies that two-way coupled models must include the fluid–particle work at scales much larger than the particle diameter, not only at the particle scale.
The κ^{-1} spectrum and logarithmic second-order structure function give a testable signature: experiments with sub-Kolmogorov resolution should see this flat spectrum when mass loading is high.
Since clustering disappears at the highest inertia even while preferential sampling persists, cluster-formation criteria based on Stokes number alone need to be revised for dense Kolmogorov-size particles.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the κ^{-1} regime is a genuine inertial-range phenomenon rather than a low-Reynolds artifact, analogous suspensions at higher Reynolds number should show a progressively wider κ^{-1} plateau; this is testable and not claimed in the paper.
The same distributed-sink mechanism suggests a parallel with energy-bypassing flows: adding Kolmogorov-size heavy particles could remove energy before it reaches the dissipative scales, which might be exploitable for manipulating mixing or heat transfer.
The observed balance between axial and biaxial strain along the negative-R Vieillefosse tail may imply that heavy particles alter small-scale alignment statistics more broadly than the paper explicitly quantifies, for instance in enstrophy production.
Because the heaviest case approaches a random-Poisson spatial distribution while still preferentially sampling high-strain regions, the usual link between preferential concentration and clustering may break down for very inertial Kolmogorov-size particles.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The central result is genuine: for Kolmogorov-size spheres at mass fractions up to 0.6, the energy spectrum develops a κ^-1 range and the nonlinear flux nearly vanishes, with the balance taken over by the fluid-solid coupling and viscous dissipation. That is a distinct regime, not something I have seen for spheres. The simulations are large (2048^3, 74k particles) and the analysis is careful: the spectral budget, structure functions, and velocity-gradient invariants are all consistent with a cascade being bypassed. The authors also flag the Gibbs-artifact caveat on the κ^-4 sub-particle scaling, which I think is the right call; that part is not a clean physical result.

The main soft spot is the one the reader flagged: Re_lambda drops from 149 to 107 across the mass-fraction series, and there is no single-phase control at the lower Re. That makes the causal attribution to particles a bit less clean. I do not think the concern is fatal. A single-phase forced run at Re_lambda~107 would still show a -5/3 range in the spectrum and a nonlinear flux that is the dominant transfer; the near-zero Π in the M_p=0.6 case is a qualitative change that low-Re alone is unlikely to produce. But the control is cheap relative to the campaign and the authors should do it, or at least show a synthetic low-Re spectrum from a coarser run. The lack of error bars is a minor issue for a single simulation per case; runs at these costs rarely get ensemble error bars, and the trends are monotonic and internally consistent.

The clustering section is honest about the Stokes versus Reynolds ambiguity, and the collapse onto the random-Poisson limit at M_p=0.6 is nice. I would cite this paper for the cascade-suppression regime, and I would send it to peer review with a request for the matched-Re control. If that control comes back confirming the κ^-1 and the near-zero flux, this becomes a solid JFM paper.

Referee Report

3 major / 5 minor

Summary. This paper uses particle-resolved direct numerical simulations (immersed-boundary method, 2048^3 grid, 74,208 spheres, volume fraction 10^-3) of forced homogeneous isotropic turbulence at an unladen Taylor-scale Reynolds number of about 150 to study Kolmogorov-size (D = η) spherical particles with density ratios Ψ_p = 100–1500, corresponding to mass fractions M_p = 0.1–0.6. It reports that as particle inertia increases the carrier-flow energy spectrum departs from the canonical κ^-5/3 scaling and approaches κ^-1, the second-order structure function becomes logarithmic at super-particle separations, the spectral nonlinear flux is strongly suppressed while the fluid–solid interaction and viscous terms dominate the scale-by-scale budget, and particle clustering weakens for Ψ_p > 100. The authors interpret these observations as suppression of the inertial cascade by heavy Kolmogorov-size spheres.

Significance. If the causal interpretation is confirmed, this is a significant contribution to particle-laden turbulence: it identifies a regime distinct from larger-than-Kolmogorov particles (which provide a spectral shortcut) and from sub-Kolmogorov point particles, and it provides a high-quality benchmark PR-DNS dataset. The simulations are expensive and well executed, the diagnostics (spectral budget, structure functions, Voronoi clustering) are appropriate, and no free parameters are tuned to force the reported scalings. The paper also honestly discusses the possible Gibbs-origin of the high-wavenumber κ^-4 tail. However, because the Taylor-scale Reynolds number decreases from about 149 to 108 across the same cases, the central claim currently rests on comparisons that do not fully isolate particle inertia from the Reynolds-number reduction.

major comments (3)

[§3.1 (Table 1), §3.2–3.3] The main conclusion — that particle inertia suppresses the inertial cascade — is inferred by comparing all laden runs with the unladen case at Re_λ ≈ 149. Table 1 shows Re_λ decreases monotonically with M_p (149, 127, 118, 109, 108). No single-phase control at Re_λ ≈ 107 is reported. Since the inertial range is essentially absent and the bottleneck effect can flatten spectra at this Re_λ, the observed κ^-1 spectrum and the reduction of Π in Fig. 6 could be partly a low-Reynolds-number artifact. Please provide or cite a single-phase simulation at the same Re_λ (e.g. by reducing the forcing amplitude) and compare spectra, Π, and structure functions.
[Fig. 3(a) and §3.2] The κ^-1 scaling is claimed for 2 ≤ κ/κ_L < 70. With η/L ≈ 0.003, κ/κ_L = 70 corresponds to κη ≈ 1.3, well inside the dissipation range; the range also extends down to scales comparable to the ABC forcing wavenumber. At Re_λ ≈ 107 no well-defined inertial range exists, so a visual match to a −1 line over this range is not sufficient evidence. Please report a local logarithmic slope with uncertainty or a compensated plot, and test whether the same −1 slope appears in an unladen run at Re_λ ≈ 107.
[Fig. 6 and §3.3] The conclusion that the nonlinear flux is 'virtually zero' for M_p = 0.6 is central to the paper. The budget is time-averaged, but no averaging time (in eddy-turnover units), number of samples, or convergence check is reported. If the simulation time is too short, the near-zero Π could be statistical noise. Please provide the averaging interval and error estimates for ⟨Π(κ)⟩, particularly for the heaviest case, or show that the budget residual in (3.4) is small.

minor comments (5)

[§3.1, Table 1] The statement that ⟨L⟩ increases while ⟨K⟩ decreases is useful, but the definition of L via the spectral integral is sensitive to the low-wavenumber truncation; a brief note on this sensitivity would be helpful.
[Fig. 4(b) inset] The inset caption should clarify that the thin blue line is a guide of slope 1 in semi-logarithmic coordinates, i.e. S_2 ∼ log(r/η), rather than a plot of log(r/η) itself.
[Eq. (3.1), Table 2] The shell-averaged relative velocity in the definition of Re_p should specify whether the average is taken over the shell volume and over time; the current notation 〈·〉_sh is introduced only in words.
[§3.6, Table 3] The authors correctly note that the M_p = 0.6 Voronoi distribution collapses onto the RPP case. A sentence explaining why the lower and upper cross-over points are not reported for this case would improve readability (the table entry 'n.a.' is currently unexplained in the caption).
[Global] Since M_p is a deterministic function of Ψ_p at fixed Φ_p, the paper cannot separate density-ratio effects from mass-loading effects. This is a design choice, but the causal language would benefit from an explicit statement that the parameter sweep varies both quantities simultaneously.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central cascade-suppression claims are measured PR-DNS statistics, not fitted or definitionally derived quantities.

full rationale

The paper reports PR-DNS measurements, not a derived prediction. The central observations—κ^{-1} spectral scaling for M_p=0.6, the near-vanishing of ⟨Π⟩, S2∼log(r/η), strain-dominated Q/R statistics, and the weakening of clusters—are all computed from the velocity field and particle trajectories generated by (2.1)-(2.4), with no adjustable parameter fitted to produce them. The Stokes number uses the standard response time τ_p=Ψ_p D^2/(18ν_f), and the mass fraction is fixed by the imposed density ratio and volume fraction through (2.5). The spectral budget (3.2)-(3.4) is an exact Fourier-space rearrangement of the momentum equation, so the dominance of Π_fs and D_v in the heaviest case is a diagnostic, not a construction. The κ^{-1} claim is checked internally against the independently measured S2∼log(r/η) (inset of Fig. 4b) and compared with, not derived from, Olivieri et al. (2022a) and Chiarini et al. (2025); those self-citations are corroborative rather than load-bearing. The only genuine caveat—the absence of a single-phase control at the reduced Re_λ≈107, which the paper itself notes for clustering when it says 'separating the effect of each parameter is not straightforward'—is a confound affecting causal attribution, not a circular reduction: the laden spectra are not constructed from the unladen spectrum or from any function of Re_λ. No circular step is present.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new theoretical entities or fitted constitutive relations. Most assumptions are standard modeling choices for PR-DNS; the most consequential are the Reynolds-number confound and the unstated collision parameters.

free parameters (4)

κ^-1 scaling fit range = 2 ≤ κ/κ_L < 70
The wavenumber interval over which the heaviest-particle spectrum is identified as κ^-1 is chosen visually in Fig. 3(a); no automated fitting or uncertainty is reported. A different range would alter the exponent's credibility.
Shell outer diameter for Q-averaging = D_sh = 3D and 6D
Preferential sampling (Fig. 11) and joint PDFs (Fig. 12) depend on the volume-averaging shell thickness; the chosen values are based on Uhlmann & Chouippe (2017), but the results are sensitive to this choice.
Cluster cut-off distance = ⟨Vvoro⟩^{1/3}
In the clustering algorithm (§3.6), dense cells separated by more than the mean Voronoi spacing are excluded to prevent merging of distinct clusters; this threshold is standard but arbitrary.
Collision model parameters = not specified
The mass-spring-dashpot collision model parameters (stiffness, damping) are not reported; these parameters can affect particle trajectories and clustering statistics, especially at volume fraction 10^-3.

assumptions (5)

standard math The incompressible Navier-Stokes equations govern the carrier flow (Eq. 2.1).
The fundamental governing equations; universally accepted.
domain assumption The IBM of Hori et al. (2022) accurately captures fluid-solid coupling for D/η=1 and Re_p up to 4.3.
Section 2 relies on this method without in-paper validation for the present parameters; grid resolution is only ~6 grid points per particle diameter.
domain assumption The ABC forcing (Eq. 2.2) maintains a statistically steady, homogeneous and isotropic state in the laden cases.
The analysis assumes HIT; heavy particles may induce anisotropy or large-scale inhomogeneities that are not diagnosed.
domain assumption The Stokes drag-response time τ_p = Ψ D²/(18ν) characterizes Kolmogorov-size particles at Re_p > 1.
Table 2 uses this formula; it is strictly valid for point particles at Re_p << 1, which may not hold for the heaviest case.
domain assumption The collision model of Tsuji et al. (1993)/Costa et al. (2015) with unspecified parameters represents inter-particle contacts.
Collision parameters are not disclosed; clustering and Voro statistics could be sensitive to them.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A suspension of heavy Kolmogorov-size spheres suppresses the inertial cascade in homogeneous and isotropic turbulence." pith.science (2026). https://pith.science/paper/3DAXSUBL

@misc{pith2026260205171,
  author       = {Pith},
  title        = {Pith review of: A suspension of heavy Kolmogorov-size spheres suppresses the inertial cascade in homogeneous and isotropic turbulence},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3DAXSUBL}},
  note         = {Machine review of arXiv:2602.05171}
}

abstract

The effect of Kolmogorov-size spherical particles on homogeneous and isotropic turbulence is investigated using particle-resolved direct numerical simulations at an unladen Taylor-scale Reynolds number of $150$. Four monodisperse suspensions of particles with identical diameter and volume fraction $10^{-3}$ are considered, while the particle-to-fluid density ratio varies between $100$ and $1500$ and the mass fraction between $0.1$ and $0.6$. As particle inertia increases, the energy spectrum departs from the canonical Kolmogorov $\kappa^{-5/3}$ scaling and approaches a peculiar regime with $\kappa^{-1}$. In this limit, the nonlinear energy transfer is strongly suppressed and the kinetic energy balance is dominated by the fluid-solid interaction and the viscous dissipation. Consistently, the second-order structure function shows logarithmic scaling at separations larger than the particle diameter, indicating velocity decorrelation. Increasing particle inertia promotes axial strain and vortex compression in the vicinity of the particles and enhances the particle-fluid relative velocity. Particle clustering is maximum when the Stokes number based on the Kolmogorov time scale is $O(1)$ and weakens as the density ratio and the Stokes number increase, with the volume and the population of the clusters decreasing when inertia is enhanced. When clustering occurs, particles preferentially sample regions of high strain and low vorticity.

Figures

Figures reproduced from arXiv: 2602.05171 by the authors.

**Figure 1.** Effect of increasing particle-to-fluid density ratio and mass fraction on the turbulent dissipation rate, 𝜀. Each panel shows the dissipation rate with a logarithmic color scale in a particle-laden domain with different mass fractions: 𝑀𝑝 = 0.1 (𝑎), 0.2 (𝑏), 0.4 (𝑐), and 0.6 (𝑑). Insets in panels (a) and (d) provide a zoomed view of the dissipation in the vicinity of the particles (green spheres). Reynolds number is… view at source ↗

**Figure 2.** Effect of increasing particle-to-fluid density ratio and mass fraction on the enstrophy, E. Each panel shows the enstrophy with a logarithmic color scale in a particle-laden domain with different mass fractions: 𝑀𝑝 = 0.1 (a), 0.2 (b), 0.4 (c), and 0.6 (d). Insets in panels (a) and (d) provide a zoomed view of the enstrophy in the vicinity of the particles (green spheres). portion of the turbulent kinetic energy is t… view at source ↗

**Figure 3.** (𝑎) Kinetic energy spectra 𝐸ˆ(𝜅) and (𝑏) compensated energy spectra (𝜅/𝜅𝐿) 4 ⟨𝐸ˆ(𝜅)⟩, for different mass fraction 𝑀𝑝 of the suspension. All wavenumbers are scaled by 𝜅𝐿 = 2𝜋/𝐿. The black dashed line, the brown dash-dot line, and the black dotted line represent the (𝜅/𝜅𝐿) −5/3 , (𝜅/𝜅𝐿) −1 , and (𝜅/𝜅𝐿) −4 scaling, respectively. The black solid vertical line in (𝑎) denotes the wavenumber of the particle diameter 𝜅/𝜅𝐿 =… view at source ↗

Figures from the paper (9 more)

**Figure 4.** Figure 4: Structure functions of the fluid velocity ⟨𝑆𝑛⟩ with 𝑛 = 2, 4, 6. The thick solid curves represent the laden cases for 𝑀𝑝 = 0.2 (𝑎) and 𝑀𝑝 = 0.6 (𝑏), while the structure functions of the unladen flow field are also plotted for comparison with thin, dashed lines. The ver…

**Figure 5.** Figure 5: (𝑎) Skewness S and (𝑏) kurtosis K of the velocity field for increasing particle density 𝜌𝑝. The vertical line denotes the particle diameter normalised with the Kolmogorov length 𝐷/𝜂. The horizontal dashed line represent the values S = 0 in plot (𝑎) and K = 3 in plot (𝑏…

**Figure 6.** Figure 6: Scale-by-scale energy budget (3.4) for increasing particle inertia. The thick, solid curves represent the turbulent transport term (blue), particle forcing term (green) and dissipation term (magenta) of the cases 𝑀𝑝 = 0.1 (𝑎), 𝑀𝑝 = 0.2 (𝑏), 𝑀𝑝 = 0.4 (𝑐) and 𝑀𝑝 = 0.6 (𝑑…

**Figure 7.** Figure 7: Probability density function (PDF) of (𝑎) the second invariant 𝑄 (3.7a) and (𝑏) the third invariant 𝑅 (3.7b) of the velocity gradient tensor 𝜕𝑗𝑢𝑖 for increasing mass fraction 𝑀𝑝 (solid coloured curves) and in the unladen case (solid black curves). The values of 𝑄 and 𝑅…

**Figure 8.** Figure 8: Joint probability density function (JPDF) of the invariants 𝑄 and 𝑅 of the velocity gradient tensor 𝜕𝑗𝑢𝑖 for increasing solid mass fraction 𝑀𝑝. The invariants are normalized on the Kolmogorov time scale 𝜏𝜂. The five black curves in each panel represent isolines of cont…

**Figure 9.** Figure 9: (𝑎) PDF of the absolute value of the particle velocity |𝑈˜ 𝑖 | for increasing particle inertia (thick, coloured curves). The distributions overlap with the Gaussian curves (thin, dashed, black curves). The PDFs are rescaled with respect to the reference velocity (𝐹0𝐿) …

**Figure 10.** Figure 10: (𝑎) PDFs of the normalised Voronoi volumes Vvoro/⟨Vvoro⟩ for increasing particle-to-fluid density ratio 𝛹𝑝 and mass fraction 𝑀𝑝 (blue, orange and green) are compared to that obtained for a suspension of finite-size (f. s.) particles distributed via a random-Poisson pr…

**Figure 11.** Figure 11: Probability density function (PDF) of the second invariant of the velocity gradient tensor normalised on its root-mean-square value 𝑄/𝑄rms, evaluated at the particle position for increasing 𝑀𝑝. The normalised invariant 𝑄 is evaluated within a spherical shell with an e…

**Figure 12.** Figure 12: Joint probability density function (JPDF) of the normalised second invariant of the velocity gradient tensor 𝑄 (3.7a) evaluated around the particles, and the normalised Voronoi volume Vvoro/⟨Vvoro⟩ for (𝑎) 𝑀𝑝 = 0.1, (𝑏) 𝑀𝑝 = 0.2, (𝑐) 𝑀𝑝 = 0.4 and (𝑑) 𝑀𝑝 = 0.6. The mag…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

72 extracted references

[1]

Amsden, A. A. & Harlow, F. H.1970 A simplified MAC technique for incompressible fluid flow calculations. J. Comput. Phys.6(2), 322–325

1970
[2]

Ardekani, M. N. & Brandt, L.2019 Turbulence modulation in channel flow of finite-size spheroidal particles. J. Fluid Mech.859, 887–901

2019
[3]

I.1965 Sur la topologie des ´ecoulements stationnaires des fluides parfaits.C

Arnold, V. I.1965 Sur la topologie des ´ecoulements stationnaires des fluides parfaits.C. R. Acad. Sci. Paris 261, 17–20

1965
[4]

& Coletti, F.2017 Coherent clusters of inertial particles in homogeneous turbulence.J

Baker, L., Frankel, A., Mani, A. & Coletti, F.2017 Coherent clusters of inertial particles in homogeneous turbulence.J. Fluid Mech.833, 364–398

2017
[5]

& Eaton, J

Balachandar, S. & Eaton, J. K.2010 Turbulent dispersed multiphase flow.Ann. Rev. Fluid Mech.42(Vol. 42, 2010), 111–133

2010
[6]

& Lakhote, M.2019 Self-induced velocity correction for improved drag estimation in Euler–Lagrange point-particle simulations.J

Balachandar, S., Liu, K. & Lakhote, M.2019 Self-induced velocity correction for improved drag estimation in Euler–Lagrange point-particle simulations.J. Computat. Phys.376, 160–185

2019
[7]

& Succi, S.1993 Extended self-similarity in turbulent flows.Phys

Benzi, R., Ciliberto, S., Tripiccione, R., Baudet, C., Massaioli, F. & Succi, S.1993 Extended self-similarity in turbulent flows.Phys. Rev. E48, R29–R32

1993
[8]

Fluid Mech

Betchov, R.1956 An inequality concerning the production of vorticity in isotropic turbulence.J. Fluid Mech. 1(5), 497–504

1956

Show all 72 references

[9]

Birnstiel, T.2024 Dust growth and evolution in protoplanetary disks.Ann. Rev. Astron. Astrophys.62(Vol. 62, 2024), 157–202

2024
[10]

& Tuval, I.2019 Turbulence induces clustering and segregation of non-motile, buoyancy-regulating phytoplankton.J

Borgnino, M., Arrieta, J., Boffetta, G., De Lillo, F. & Tuval, I.2019 Turbulence induces clustering and segregation of non-motile, buoyancy-regulating phytoplankton.J. Royal Soc. Interface16(159), 20190324

2019
[11]

& Coletti, F.2022 Particle-laden turbulence: Progress and perspectives.Ann

Brandt, L. & Coletti, F.2022 Particle-laden turbulence: Progress and perspectives.Ann. Rev. Fluid Mech. 54(Vol. 54, 2022), 159–189

2022
[12]

Breugem, W.-P.2012 A second-order accurate immersed boundary method for fully resolved simulations of particle-laden flows.J. Comput. Phys.231(13), 4469–4498

2012
[13]

& Rosti, M

Cannon, I., Olivieri, S. & Rosti, M. E.2024 Spheres and fibers in turbulent flows at various Reynolds numbers.Phys. Rev. Fluids9, 064301

2024
[14]

J.1992 Exact solution of a restricted Euler equation for the velocity gradient tensor.Phys

Cantwell, B. J.1992 Exact solution of a restricted Euler equation for the velocity gradient tensor.Phys. Fluids A: Fluid Dyn.4(4), 782–793

1992
[15]

& Rosti, M

Chiarini, A., Cannon, I. & Rosti, M. E.2024 Anisotropic mean flow enhancement and anomalous transport of finite-size spherical particles in turbulent flows.Phys. Rev. Lett.132, 054005

2024
[16]

& Rosti, M

Chiarini, A. & Rosti, M. E.2024 Finite-size inertial spherical particles in turbulence.J. Fluid Mech.988, A17

2024
[17]

& Rosti, M

Chiarini, A., Tandurella, S. & Rosti, M. E.2025 Kolmogorov-size particles in homogeneous and isotropic turbulence.J. Fluid Mech.1007, A81

2025
[18]

& Uhlmann, M.2019 On the influence of forced homogeneous-isotropic turbulence on the settling and clustering of finite-size particles.Acta Mech.230(2), 387–412

Chouippe, A. & Uhlmann, M.2019 On the influence of forced homogeneous-isotropic turbulence on the settling and clustering of finite-size particles.Acta Mech.230(2), 387–412

2019
[19]

A., Meyer, D

Codispoti, L. A., Meyer, D. W. & Jenny, P.2025 Dissecting inertial clustering and sling dynamics in high- Reynolds number particle-laden turbulence.Phys. Fluids37(1), 015161. 0X0-26 Journal of Fluid Mechanics

2025
[20]

J., Westerweel, J

Costa, P., Boersma, B. J., Westerweel, J. & Breugem, W.-P.2015 Collision model for fully resolved simulations of flows laden with finite-size particles.Phys. Rev. E92, 053012

2015
[21]

& Picano, F.2021 Near-wall turbulence modulation by small inertial particles.J

Costa, P., Brandt, L. & Picano, F.2021 Near-wall turbulence modulation by small inertial particles.J. Fluid Mech.922, A9

2021
[22]

& Bourgoin, M.2012 Clustering of finite-size particles in turbulence.Phys

Fiabane, L., Zimmermann, R., Volk, R., Pinton, J.-F. & Bourgoin, M.2012 Clustering of finite-size particles in turbulence.Phys. Rev. E86, 035301

2012
[23]

B.2014 Numerical simulation of flowing blood cells.Ann

Freund, J. B.2014 Numerical simulation of flowing blood cells.Ann. Rev. Fluid Mech.46(Volume 46, 2014), 67–95

2014
[24]

& Richter, D

Gao, W., Samtaney, R. & Richter, D. H.2023 Direct numerical simulation of particle-laden flow in an open channel at𝑟𝑒 𝜏 =5186.J. Fluid Mech.957, A3

2023
[25]

& Crowe, C.T.1989 Effect of particle size on modulating turbulent intensity.Int

Gore, R.A. & Crowe, C.T.1989 Effect of particle size on modulating turbulent intensity.Int. J. Multiph. Flow 15(2), 279–285

1989
[26]

& Vassilicos, J

Goto, S. & Vassilicos, J. C.2006 Self-similar clustering of inertial particles and zero-acceleration points in fully developed two-dimensional turbulence.Phys. Fluids18(11), 115103

2006
[27]

& Vassilicos, J

Goto, S. & Vassilicos, J. C.2008 Sweep-stick mechanism of heavy particle clustering in fluid turbulence. Phys. Rev. Lett.100, 054503

2008
[28]

& Casciola, C.M.2023 Effect of Stokes number and particle-to- fluid density ratio on turbulence modification in channel flows.J

Gualtieri, P., Battista, F., Salvadore, F. & Casciola, C.M.2023 Effect of Stokes number and particle-to- fluid density ratio on turbulence modification in channel flows.J. Fluid Mech.974, A26

2023
[29]

Hassaini, R., Petersen, A. J. & Coletti, F.2023 Effect of two-way coupling on clustering and settling of heavy particles in homogeneous turbulence.J. Fluid Mech.976, A12

2023
[30]

Hori, N., Rosti, M. E. & Takagi, S.2022 An Eulerian-based immersed boundary method for particle suspensions with implicit lubrication model.Comput. Fluids236, 105278

2022
[31]

& Eaton, J

Hwang, W. & Eaton, J. K.2006 Homogeneous and isotropic turbulence modulation by small heavy (𝑆𝑡∼50) particles.J. Fluid Mech.564, 361–393

2006
[32]

& Beishuizen, N.2012 Modeling of turbulent dilute spray combustion.Prog

Jenny, P., Roekaerts, D. & Beishuizen, N.2012 Modeling of turbulent dilute spray combustion.Prog. Energy Combust. Sci.38(6), 846–887

2012
[33]

& Miyake, Y.2001 Turbulence structure of particle-laden flow in a vertical plane channel due to vortex shedding.JSME Int

Kajishima, T., Takiguchi, S., Hamasaki, H. & Miyake, Y.2001 Turbulence structure of particle-laden flow in a vertical plane channel due to vortex shedding.JSME Int. J. Ser. B Fluids Therm. Eng.44(4), 526–535

2001
[34]

& Fr¨ohlich, J.2012 An improved immersed boundary method with direct forcing for the simulation of particle laden flows.J

Kempe, T. & Fr¨ohlich, J.2012 An improved immersed boundary method with direct forcing for the simulation of particle laden flows.J. Comput. Phys.231(9), 3663–3684

2012
[35]

& Moser, R.1987 Turbulence statistics in fully developed channel flow at low Reynolds number.J

Kim, J., Moin, P. & Moser, R.1987 Turbulence statistics in fully developed channel flow at low Reynolds number.J. Fluid Mech.177, 133–166

1987
[36]

Raj Mohan, B., Byeon, S.-H., Lim, K.-S

Lee, B.-K., B. Raj Mohan, B., Byeon, S.-H., Lim, K.-S. & Hong, E.-P.2013 Evaluating the performance of a turbulent wet scrubber for scrubbing particulate matter.J. Air Waste Manag. Assoc.63(5), 499–506

2013
[37]

Lu, J., Lazar, E. A. & Rycroft, C. H.2023 An extension to Voro++ for multithreaded computation of Voronoi cells.Comput. Phys. Commun.291, 108832

2023
[38]

& Elghobashi, S.2010 Modulation of isotropic turbulence by particles of taylor length-scale size.J

Lucci, F., Ferrante, A. & Elghobashi, S.2010 Modulation of isotropic turbulence by particles of taylor length-scale size.J. Fluid Mech.650, 5–55

2010
[39]

& Fan, J.2017 Fully resolved simulations of turbulence modulation by high-inertia particles in an isotropic turbulent flow.Phys

Luo, K., Wang, Z., Li, D., Tan, J. & Fan, J.2017 Fully resolved simulations of turbulence modulation by high-inertia particles in an isotropic turbulent flow.Phys. Fluids29(11), 113301

2017
[40]

O., Magnaudet, J., Reeks, M., Simonin, O., Sommerfeld, M., Toschi, F., Wang, L.-P

Marchioli, C., Bourgoin, M., Coletti, F., Fox, R. O., Magnaudet, J., Reeks, M., Simonin, O., Sommerfeld, M., Toschi, F., Wang, L.-P. & Balachandar, S.2025 Particle-laden flows.Int. J. Multiph. Flow191, 105291

2025
[41]

Marchioli, C., Rosti, M. E. & Verhille, G.2026 Flexible fibers in turbulence.Ann. Rev. Fluid Mech. 58(Volume 58, 2026), 167–192

2026
[42]

& Schneider, K.2024 Heavy particle clustering in inertial subrange of high– Reynolds number turbulence.Phys

Matsuda, K., Yoshimatsu, K. & Schneider, K.2024 Heavy particle clustering in inertial subrange of high– Reynolds number turbulence.Phys. Rev. Lett.132, 234001

2024
[43]

R.1987 The gravitational settling of aerosol particles in homogeneous turbulence and random flow fields.J

Maxey, M. R.1987 The gravitational settling of aerosol particles in homogeneous turbulence and random flow fields.J. Fluid Mech.174, 441–465

1987
[44]

Maxey, M. R. & Riley, J. J.1983 Equation of motion for a small rigid sphere in a nonuniform flow.Phys. Fluids26(4), 883–889

1983
[45]

P.2017 Cloud-top entrainment in stratocumulus clouds.Ann

Mellado, J. P.2017 Cloud-top entrainment in stratocumulus clouds.Ann. Rev. Fluid Mech.49(Volume 49, 2017), 145–169

2017
[46]

Meneveau, C.2011 Lagrangian dynamics and models of the velocity gradient tensor in turbulent flows.Ann. Rev. Fluid Mech.43(Vol. 43, 2011), 219–245

2011
[47]

Phys.14(9), 095013

Monchaux, R.2012 Measuring concentration with Voronoi diagrams: the study of possible biases.New J. Phys.14(9), 095013. 0X0-27 L. Foss`a and M. E. Rosti

2012
[48]

& Cartellier, A.2010 Preferential concentration of heavy particles: A Vorono¨ı analysis.Phys

Monchaux, R., Bourgoin, M. & Cartellier, A.2010 Preferential concentration of heavy particles: A Vorono¨ı analysis.Phys. Fluids22(10), 103304

2010
[49]

Monti, A., Rathee, V., Shen, A. Q. & Rosti, M. E.2021 A fast and efficient tool to study the rheology of dense suspensions.Phys. Fluids33(10), 103314

2021
[50]

Mortimer, L. F. & Fairweather, M.2020 Density ratio effects on the topology of coherent turbulent structures in two-way coupled particle-laden channel flows.Phys. Fluids32(10), 103302

2020
[51]

& Bourgoin, M.2014 Preferential concentration of heavy particles in turbulence.J

Obligado, M., Teitelbaum, T., Cartellier, A., Mininni, P. & Bourgoin, M.2014 Preferential concentration of heavy particles in turbulence.J. Turbul.15(5), 293–310

2014
[52]

& Goto, S.2022 Attenuation of turbulence in a periodic cube by finite-size spherical solid particles

Oka, S. & Goto, S.2022 Attenuation of turbulence in a periodic cube by finite-size spherical solid particles. J. Fluid Mech.949, A45

2022
[53]

& Ooms, G.2007 Particle–fluid interactions in grid-generated turbulence.J

Poelma, C., Westerweel, J. & Ooms, G.2007 Particle–fluid interactions in grid-generated turbulence.J. Fluid Mech.589, 315–351

2007
[54]

Fluid Mech.791, 174–190

Sun, C.2016 Energy spectra in turbulent bubbly flows.J. Fluid Mech.791, 174–190

2016
[55]

& Risso, F.2024 Spectral analysis of dispersed multiphase flows in the presence of fluid interfaces.Int

Ramirez, G., Burlot, A., Zamansky, R., Bois, G. & Risso, F.2024 Spectral analysis of dispersed multiphase flows in the presence of fluid interfaces.Int. J. Multiph. Flow177, 104860

2024
[56]

Rose, W. I. & Durant, A. J.2009 Fine ash content of explosive eruptions.J. Volcanol. Geotherm. Res.186(1), 32–39

2009
[57]

& Kheirkhah, S.2024 Three-dimensional clustering characteristics of large-Stokes-number dilute sprays interacting with turbulent swirling co-flows.J

Rostami, A., Li, R. & Kheirkhah, S.2024 Three-dimensional clustering characteristics of large-Stokes-number dilute sprays interacting with turbulent swirling co-flows.J. Fluid Mech.999, A73

2024
[58]

Rosti, M. E. & Brandt, L.2020 Increase of turbulent drag by polymers in particle suspensions.Phys. Rev. Fluids5, 041301

2020
[59]

& Schr ¨oder, W.2017 Direct particle–fluid simulation of Kolmogorov-length- scale size particles in decaying isotropic turbulence.J

Schneiders, L., Meinke, M. & Schr ¨oder, W.2017 Direct particle–fluid simulation of Kolmogorov-length- scale size particles in decaying isotropic turbulence.J. Fluid Mech.819, 188–227

2017
[60]

& Wang, L.-P.2024 The influence of particle density and diameter on the interactions between the finite-size particles and the turbulent channel flow.Int

Shen, J., Peng, C., Lu, Z. & Wang, L.-P.2024 The influence of particle density and diameter on the interactions between the finite-size particles and the turbulent channel flow.Int. J. Multiph. Flow170, 104659

2024
[61]

Squires, K. D. & Eaton, J. K.1990 Particle response and turbulence modification in isotropic turbulence. Phys. Fluids A: Fluid Dyn.2(7), 1191–1203

1990
[62]

M., Capuano, T., Brandt, L., Iudicone, D

Sugathapala, T. M., Capuano, T., Brandt, L., Iudicone, D. & Sardina, G.2025 Vertical transport of buoyant microplastic particles in the ocean: The role of turbulence and biofouling.Environ. Pollut.369, 125819

2025
[63]

& Bourgoin, M.2017 Preferential concentration of inertial sub-Kolmogorov particles: The roles of mass loading of particles, Stokes numbers, and Reynolds numbers

Sumbekova, S., Cartellier, A., Aliseda, A. & Bourgoin, M.2017 Preferential concentration of inertial sub-Kolmogorov particles: The roles of mass loading of particles, Stokes numbers, and Reynolds numbers. Phys. Rev. Fluids2, 024302

2017
[64]

M., Prakash, Vivek N., Calzavarini, E., Sun, C

Tagawa, Y., Mercado, J. M., Prakash, Vivek N., Calzavarini, E., Sun, C. & Lohse, D.2012 Three- dimensional Lagrangian Vorono ¨ı analysis for clustering of particles and bubbles in turbulence.J. Fluid Mech.693, 201–215

2012
[65]

& Eaton, J

Tanaka, T. & Eaton, J. K.2010 Sub-Kolmogorov resolution partical image velocimetry measurements of particle-laden forced turbulence.J. Fluid Mech.643, 177–206

2010
[66]

Tavanashad, V., Passalacqua, A., Fox, R. O. & Subramaniam, S.2019 Effect of density ratio on velocity fluctuations in dispersed multiphase flow from simulations of finite-size particles.Acta Mech.230(2), 469–484. Tom´e, M. F., Duffy, B. & McKee, S.1996 A numerical technique fo...

2019
[67]

& Tanaka, T.1993 Discrete particle simulation of two-dimensional fluidized bed

Tsuji, Y., Kawaguchi, T. & Tanaka, T.1993 Discrete particle simulation of two-dimensional fluidized bed. Powder Tech.77(1), 79–87

1993
[68]

Uhlmann, M.2005 An immersed boundary method with direct forcing for the simulation of particulate flows. J. Comput. Phys.209(2), 448–476

2005
[69]

Fluids20(5), 053305

Uhlmann, M.2008 Interface-resolved direct numerical simulation of vertical particulate channel flow in the turbulent regime.Phys. Fluids20(5), 053305

2008
[70]

& Chouippe, A.2017 Clustering and preferential concentration of finite-size particles in forced homogeneous-isotropic turbulence.J

Uhlmann, M. & Chouippe, A.2017 Clustering and preferential concentration of finite-size particles in forced homogeneous-isotropic turbulence.J. Fluid Mech.812, 991–1023

2017
[71]

& Goto, S.2007 Self-similar clustering of inertial particles in homogeneous turbulence.J

Yoshimoto, H. & Goto, S.2007 Self-similar clustering of inertial particles in homogeneous turbulence.J. Fluid Mech.577, 275–286. 0X0-28 Journal of Fluid Mechanics

2007
[72]

& Wang, L.-P.2017 Effects of particle-fluid density ratio on the interactions between the turbulent channel flow and finite-size particles.Phys

Yu, Z., Lin, Z., Shao, X. & Wang, L.-P.2017 Effects of particle-fluid density ratio on the interactions between the turbulent channel flow and finite-size particles.Phys. Rev. E96, 033102. 0X0-29

2017

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.

[1] [1]

Amsden, A. A. & Harlow, F. H.1970 A simplified MAC technique for incompressible fluid flow calculations. J. Comput. Phys.6(2), 322–325

1970

[2] [2]

Ardekani, M. N. & Brandt, L.2019 Turbulence modulation in channel flow of finite-size spheroidal particles. J. Fluid Mech.859, 887–901

2019

[3] [3]

I.1965 Sur la topologie des ´ecoulements stationnaires des fluides parfaits.C

Arnold, V. I.1965 Sur la topologie des ´ecoulements stationnaires des fluides parfaits.C. R. Acad. Sci. Paris 261, 17–20

1965

[4] [4]

& Coletti, F.2017 Coherent clusters of inertial particles in homogeneous turbulence.J

Baker, L., Frankel, A., Mani, A. & Coletti, F.2017 Coherent clusters of inertial particles in homogeneous turbulence.J. Fluid Mech.833, 364–398

2017

[5] [5]

& Eaton, J

Balachandar, S. & Eaton, J. K.2010 Turbulent dispersed multiphase flow.Ann. Rev. Fluid Mech.42(Vol. 42, 2010), 111–133

2010

[6] [6]

& Lakhote, M.2019 Self-induced velocity correction for improved drag estimation in Euler–Lagrange point-particle simulations.J

Balachandar, S., Liu, K. & Lakhote, M.2019 Self-induced velocity correction for improved drag estimation in Euler–Lagrange point-particle simulations.J. Computat. Phys.376, 160–185

2019

[7] [7]

& Succi, S.1993 Extended self-similarity in turbulent flows.Phys

Benzi, R., Ciliberto, S., Tripiccione, R., Baudet, C., Massaioli, F. & Succi, S.1993 Extended self-similarity in turbulent flows.Phys. Rev. E48, R29–R32

1993

[8] [8]

Fluid Mech

Betchov, R.1956 An inequality concerning the production of vorticity in isotropic turbulence.J. Fluid Mech. 1(5), 497–504

1956