pith. sign in

arxiv: 2605.22326 · v2 · pith:ZOP6QGKEnew · submitted 2026-05-21 · ⚛️ physics.flu-dyn

Lagrangian single-particle, multi-particle and topological analyses in turbulent Rayleigh-B\'enard convection

Matti Ettel , Roshan J. Samuel , Michael Chertkov , J\"org Schumacher This is my paper

Pith reviewed 2026-05-22 02:43 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords Lagrangian turbulenceRayleigh-Bénard convectionparticle dispersionheat transfer intermittencyvelocity gradient topologyplume dynamics
0
0 comments X

The pith

Lagrangian particle tracking reveals extreme intermittency and organized dispersion in Rayleigh-Bénard convection

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

This paper tracks individual fluid particles through direct numerical simulations of Rayleigh-Bénard convection across Rayleigh numbers from 10^5 to 10^10 at Prandtl number 0.7. Heat transport along single trajectories proves far more variable than global averages indicate, with some particles carrying fluxes hundreds of times the mean value. Velocity-gradient invariants along paths mark specific regions tied to vortex stretching and plume release. Pair separation does not settle into one power law; instead a brief rapid ejection phase gives way to sustained scaling, sequenced by buoyancy then shear. These trajectory-based diagnostics therefore isolate the actual transport mechanisms that Eulerian statistics average away.

Core claim

The central claim is that Lagrangian heat transfer remains extremely intermittent, with individual massless particles achieving convective fluxes up to 500 times the global Eulerian mean while higher-order moments relax toward Gaussian values at higher Ra; Q-R analysis of the velocity gradient along trajectories isolates a distinct topological footprint of dust-devil-like convective vortices in the Q>0, R<0 quadrant linked to stretching, detachment, and localized heat transfer; and conditioned pair-dispersion plus cloud PCA show temporally organized transport consisting of a short t^5-like ejection episode followed by sustained Richardson-like t^3 scaling rather than extended plateaus.

What carries the argument

Q-R invariant topology of the velocity gradient tensor sampled along material trajectories together with principal-component analysis of dense Lagrangian particle clouds

Load-bearing premise

The chosen aspect ratio of 4:4:1 and the numerical resolution at Rayleigh numbers up to 10^10 are sufficient to capture the full range of Lagrangian statistics without dominant boundary or discretization artifacts affecting the reported intermittency and dispersion regimes.

What would settle it

A higher-resolution simulation or laboratory experiment at the same or larger Rayleigh numbers that finds substantially lower peak Lagrangian heat fluxes or that eliminates the distinct Q>0, R<0 population would falsify the claimed mechanism-resolving power of the topological and cloud-geometry diagnostics.

Figures

Figures reproduced from arXiv: 2605.22326 by J\"org Schumacher, Matti Ettel, Michael Chertkov, Roshan J. Samuel.

Figure 1
Figure 1. Figure 1: Lagrangian particle acceleration statistics. Probability density functions (PDFs) of the vertical ( [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Evaluation of the Heisenberg-Yaglom constant [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Lagrangian Nusselt number statistics. Probability density functions (PDFs) of the total Lagrangian [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Parameter plane spanned by the invariants [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Analysis of dust-devil-like convective vortices. (a) Projected trajectories ( [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Evolution of particle pair and single-particle dispersion. Relative particle pair dispersion [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Scale-dependent Lagrangian eddy viscosity [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Geometric overview of particle clouds at [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Geometric evolution of a representative particle cloud seeded in the mixing zone ( [PITH_FULL_IMAGE:figures/full_fig_p024_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Ensemble statistics of the temporal evolution of Lagrangian particle clouds. Ensemble averages [PITH_FULL_IMAGE:figures/full_fig_p025_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Cubic smoothing splines for Lagrangian particle accelerations at [PITH_FULL_IMAGE:figures/full_fig_p029_11.png] view at source ↗
read the original abstract

We present three-dimensional direct numerical simulations of turbulent Rayleigh-B\'enard convection (RBC) in the Lagrangian frame of reference for Rayleigh numbers $10^5 \leq Ra \leq 10^{10}$ and a Prandtl number $Pr=0.7$ in a plane layer at an aspect ratio $L:L:H=4:4:1$ with a horizontal length $L$ and height $H$. We use particle accelerations, Lagrangian heat transfer, $Q$-$R$ invariant topology, Lagrangian particle pair dispersion, scale-dependent Lagrangian eddy viscosity, and principal component analysis (PCA) of dense particle clouds to characterise convective transport along material trajectories. By computing particle accelerations at the integration time step and controlling spectral element method signatures, we obtain robust acceleration statistics and recover Heisenberg-Yaglom behaviour. Lagrangian heat transfer is extremely intermittent: individual massless Lagrangian particles can carry convective heat fluxes up to $500$ times the global Eulerian mean, although higher-order heat flux moments decrease toward Gaussian values with increasing $Ra$. The analysis of velocity gradient invariants in the $Q$-$R$ plane along trajectories identifies a distinct topological footprint of dust-devil-like convective vortices in the quadrant of $Q>0$, $R<0$, associated with vortex stretching, plume detachment, and intense localised heat transfer. Global unconditioned pair dispersion exhibits neither extended Richardson nor Bolgiano-Obukhov scaling plateaus. Rather, scale-dependent eddy viscosity and conditioned PCA of dense particle clouds reveal that buoyancy- and shear-driven dispersion are temporally organised: rapid plume-driven ejection produces a short $t^5$-like episode, followed by sustained Richardson-like $t^3$-scaling. Thus, Lagrangian topology and cloud geometry provide mechanism-resolving diagnostics for active-scalar turbulence beyond RBC-specific global scaling laws.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper presents 3D DNS of turbulent Rayleigh-Bénard convection in the Lagrangian frame for 10^5 ≤ Ra ≤ 10^10 at Pr=0.7 in a 4:4:1 aspect-ratio layer. Lagrangian particle tracking is used to compute accelerations, heat fluxes, Q-R velocity-gradient invariants, pair dispersion, scale-dependent eddy viscosity, and PCA of dense particle clouds. Key results include recovery of Heisenberg-Yaglom acceleration scaling, extreme intermittency with individual particles carrying up to 500 times the mean heat flux, a distinct Q>0, R<0 vortex-stretching footprint linked to plume detachment, and temporally sequenced dispersion (short t^5-like ejection followed by sustained t^3 scaling) revealed only after conditioning. The central claim is that Lagrangian topology and cloud geometry supply mechanism-resolving diagnostics for active-scalar turbulence that go beyond RBC global scaling laws.

Significance. If the reported dispersion regimes, intermittency levels, and topological signatures prove robust to domain size and resolution, the work supplies concrete Lagrangian diagnostics that can distinguish buoyancy-driven ejection from shear-driven spreading in convective turbulence. Explicit recovery of Heisenberg-Yaglom scaling and the Q-R footprint associated with dust-devil-like vortices are strengths that could be useful for model validation and for connecting single-particle statistics to multi-particle cloud evolution.

major comments (2)
  1. [Abstract and simulation-setup description] Abstract and simulation-setup description: the central claim that the short t^5-like ejection followed by sustained t^3 dispersion and the extreme heat-flux intermittency are intrinsic mechanism-resolving features requires that these statistics are insensitive to the chosen aspect ratio 4:4:1 and to the grid resolution at Ra=10^10. In a 4:4:1 domain the large-scale circulation spans a distance comparable to the horizontal box size, which can organize plume ejections and particle-cloud geometry; at Ra=10^10 the Kolmogorov scale is small enough that modest under-resolution of the dissipation range can alter acceleration PDFs and pair-separation curves. No grid-convergence or aspect-ratio sensitivity data are supplied for the conditioned PCA or scale-dependent eddy-viscosity results, leaving open the possibility that the reported temporal organization is partly an artifact of domain or grid
  2. [Dispersion and PCA results] Dispersion and PCA results: the manuscript states that global unconditioned pair dispersion shows neither extended Richardson nor Bolgiano-Obukhov plateaus, yet conditioned PCA of dense particle clouds recovers the t^5-then-t^3 sequence. Because the conditioning is performed post-hoc on plume events identified from the same trajectories, it is unclear whether the temporal organization is an independent diagnostic or is partly selected by the conditioning procedure itself; an a-priori, trajectory-independent definition of the conditioning criterion would be needed to establish that the t^5/t^3 sequence is a genuine dynamical feature rather than a consequence of the analysis choice.
minor comments (2)
  1. [Abstract] The phrase 'controlling spectral-element-method signatures' used to justify robust acceleration statistics is not defined quantitatively; a brief description of the filtering or interpolation procedure applied at the integration time step would clarify how the Heisenberg-Yaglom recovery is achieved.
  2. [Figure captions] Figure captions and axis labels for the Q-R plane plots and the conditioned dispersion curves should explicitly state the number of particles, the time window used for conditioning, and the precise definition of the scale-dependent eddy viscosity to allow direct reproduction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment below with clarifications and indicate planned revisions. The responses focus on substance and maintain the integrity of the reported results.

read point-by-point responses

  1. Referee: Abstract and simulation-setup description: the central claim that the short t^5-like ejection followed by sustained t^3 dispersion and the extreme heat-flux intermittency are intrinsic mechanism-resolving features requires that these statistics are insensitive to the chosen aspect ratio 4:4:1 and to the grid resolution at Ra=10^10. In a 4:4:1 domain the large-scale circulation spans a distance comparable to the horizontal box size, which can organize plume ejections and particle-cloud geometry; at Ra=10^10 the Kolmogorov scale is small enough that modest under-resolution of the dissipation range can alter acceleration PDFs and pair-separation curves. No grid-convergence or aspect-ratio sensitivity data are supplied for the conditioned PCA or scale-dependent eddy-viscosity results, leaving open the possibility that the reported temporal organization is partly an artifact of domain or grid

    Authors: We agree that explicit grid-convergence tests and aspect-ratio sensitivity studies for the conditioned PCA and scale-dependent eddy-viscosity statistics are not included in the manuscript. The recovery of Heisenberg-Yaglom scaling in the acceleration PDFs serves as an indirect resolution check, and the 4:4:1 aspect ratio follows standard practice in RBC DNS to accommodate the large-scale circulation. We will revise the manuscript to add a dedicated paragraph discussing resolution criteria (based on Kolmogorov and Batchelor scales) and the consistency of key statistics across the simulated Ra range. We will also explicitly note the absence of aspect-ratio variations as a limitation. New simulations for sensitivity tests are computationally demanding and beyond the scope of the current revision but will be flagged for future work. revision: partial

  2. Referee: Dispersion and PCA results: the manuscript states that global unconditioned pair dispersion shows neither extended Richardson nor Bolgiano-Obukhov plateaus, yet conditioned PCA of dense particle clouds recovers the t^5-then-t^3 sequence. Because the conditioning is performed post-hoc on plume events identified from the same trajectories, it is unclear whether the temporal organization is an independent diagnostic or is partly selected by the conditioning procedure itself; an a-priori, trajectory-independent definition of the conditioning criterion would be needed to establish that the t^5/t^3 sequence is a genuine dynamical feature rather than a consequence of the analysis choice.

    Authors: The conditioning is anchored in the Q-R topological signature (Q>0, R<0 quadrant) and associated intense heat-flux events that we identify as characteristic of plume detachment and dust-devil-like vortices. This criterion is physically motivated by the Eulerian-Lagrangian connection developed in the paper rather than an arbitrary post-selection. Nevertheless, we acknowledge the referee's point that greater transparency is needed. In the revision we will provide an explicit, a-priori definition of the conditioning threshold (including the precise Q-R and heat-flux criteria) and will demonstrate that the t^5-like to t^3 transition remains visible when the same criterion is applied to an independent subset of trajectories not used for topology identification. This will strengthen the claim that the temporal staging is a dynamical feature. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results follow from DNS trajectory post-processing

full rationale

The paper's claims rest on direct numerical simulation of RBC at specified Ra and aspect ratio, followed by post-processing to extract particle accelerations, heat fluxes, Q-R invariants, pair dispersion, scale-dependent eddy viscosity, and PCA of clouds. No derivation step reduces a claimed prediction or diagnostic to a fitted parameter or self-citation by construction; the reported t^5-like ejection, t^3 dispersion, vortex topology, and intermittency emerge as observed statistics from the trajectories. The scale-dependent eddy viscosity is introduced as an analysis tool rather than a self-consistent closure that forces the dispersion result. The work is self-contained against external benchmarks of Lagrangian turbulence statistics and does not rely on load-bearing self-citations or ansatzes imported from prior author work.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The study rests on standard incompressible fluid equations and numerical methods for DNS; no new free parameters, axioms, or postulated entities are introduced beyond the conventional setup of Rayleigh-Bénard convection.

axioms (2)
  • standard math Incompressible Navier-Stokes equations govern the flow at the simulated Prandtl and Rayleigh numbers
    Implicit foundation of all direct numerical simulations of RBC
  • domain assumption Lagrangian particles faithfully represent material trajectories without significant numerical diffusion or interpolation errors
    Required for all reported acceleration and heat-flux statistics

pith-pipeline@v0.9.0 · 5879 in / 1474 out tokens · 44320 ms · 2026-05-22T02:43:10.039641+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    Global unconditioned pair dispersion exhibits neither extended Richardson nor Bolgiano-Obukhov scaling plateaus. Rather, scale-dependent eddy viscosity and conditioned PCA of dense particle clouds reveal that buoyancy- and shear-driven dispersion are temporally organised: rapid plume-driven ejection produces a short t^5-like episode, followed by sustained Richardson-like t^3-scaling.

  • IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    The analysis of velocity gradient invariants in the Q–R plane along trajectories identifies a distinct topological footprint of dust-devil-like convective vortices in the quadrant of Q>0, R<0

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.