Sharp transitions in rotating turbulent convection: Lagrangian acceleration statistics reveal a second critical Rossby number
Pith reviewed 2026-05-25 09:07 UTC · model grok-4.3
The pith
Lagrangian acceleration statistics reveal a second critical Rossby number for flow structure change in rotating convection
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In rotating Rayleigh-Benard convection, the dominant flow structure changes from a domain-filling large-scale circulation to a collection of rotation-aligned plumes at a critical Rossby number Ro_c2 ≈ 2.25. This is revealed by the sudden increase in the root-mean-square values and the kurtosis of the horizontal acceleration of Lagrangian tracers near the top plate. This structural transition precedes the sharp transition in the Nusselt number, which occurs at Ro_c1 ≈ 2.7 for the same parameter settings.
What carries the argument
Lagrangian acceleration statistics of tracers near the top plate, with focus on the rms and kurtosis of their horizontal components, which capture small-scale flow properties.
If this is right
- The flow structure transition occurs at Ro ≈ 2.25, earlier than the heat transfer transition at Ro ≈ 2.7.
- The large-scale circulation is replaced by aligned plumes at the lower critical value.
- The boundary layer change from Prandtl-Blasius to Ekman type follows the structural reorganization.
- Two distinct critical Rossby numbers govern the dynamics at these parameters.
Where Pith is reading between the lines
- The increase in heat transfer may be a consequence of the plume alignment rather than the driver of the transition.
- Similar Lagrangian diagnostics could uncover multiple thresholds in other parameter regimes of rotating convection.
- Experiments measuring particle accelerations near boundaries could confirm the separation of the two critical values.
Load-bearing premise
The abrupt change in tracer acceleration statistics at Rossby number 2.25 corresponds to an independent shift in the dominant flow structures, separate from the Nusselt number transition, and is not due to limitations in the simulation resolution or sampling.
What would settle it
A higher-resolution simulation or laboratory experiment showing continuous variation in horizontal acceleration rms and kurtosis through Rossby number 2.25, or an exact coincidence with the Nusselt transition, would disprove the existence of a distinct second critical Rossby number.
Figures
read the original abstract
In RB convection for fluids with Prandtl number $Pr\gtrsim 1$, rotation beyond a critical (small) rotation rate is known to cause a sudden enhancement of heat transfer which can be explained by a change in the character of the BL dynamics near the top and bottom plates of the convection cell. Namely, with increasing rotation rate, the BL signature suddenly changes from Prandtl--Blasius type to Ekman type. The transition from a constant heat transfer to an almost linearly increasing heat transfer with increasing rotation rate is known to be sharp and the critical Rossby number $Ro_{c}$ occurs typically in the range $2.3\lesssim Ro_{c}\lesssim 2.9$ (for Rayleigh number $Ra=1.3\times 10^9$, $Pr=6.7$, and a convection cell with aspect ratio $\Gamma=\frac{D}{H}=1$, with $D$ the diameter and $H$ the height of the cell). The explanation of the sharp transition in the heat transfer points to the change in the dominant flow structure. At $1/Ro\lesssim 1/Ro_c$ (slow rotation), the well-known LSC is found: a single domain-filling convection roll made up of many individual thermal plumes. At $1/Ro\gtrsim 1/Ro_c$ (rapid rotation), the LSC vanishes and is replaced with a collection of swirling plumes that align with the rotation axis. In this paper, by numerically studying Lagrangian acceleration statistics, related to the small-scale properties of the flow structures, we reveal that this transition between these different dominant flow structures happens at a second critical Rossby number, $Ro_{c_2}\approx 2.25$ (different from $Ro_{c_1}\approx 2.7$ for the sharp transition in the Nusselt number $Nu$; both values for the parameter settings of our present numerical study). When statistical data of Lagrangian tracers near the top plate are collected, it is found that the root-mean-square (rms) values and the kurtosis of the horizontal acceleration of these tracers show a sudden increase at $Ro_{c_2}$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports direct numerical simulations of rotating Rayleigh-Bénard convection (Ra=1.3×10^9, Pr=6.7, Γ=1) and uses Lagrangian tracer statistics to argue for a second critical Rossby number Ro_c2≈2.25. At this value the rms and kurtosis of horizontal acceleration for tracers near the top plate increase abruptly, which the authors interpret as the transition from a domain-filling large-scale circulation to rotation-aligned plumes; this is stated to be distinct from the known Nusselt-number transition at Ro_c1≈2.7.
Significance. If the reported discontinuity is shown to be robust and independent of the heat-transport transition, the result would refine the picture of how rotation reorganizes the flow before the boundary-layer regime changes. The Lagrangian-acceleration diagnostic is a potentially useful probe of small-scale structure that is not directly accessible from Eulerian fields.
major comments (2)
- [Abstract / Numerical methods] Abstract and numerical-setup section: no grid resolution, tracer count, integration length, spatial definition of “near the top plate,” or convergence tests are supplied. Because the central claim rests on the sudden jump in rms and kurtosis being physical rather than a sampling or binning artifact, these quantities must be documented and shown to be insensitive to reasonable variations in tracer ensemble size and location.
- [Results (tracer statistics)] Results section describing the acceleration statistics: the manuscript must demonstrate that the location of the jump at Ro≈2.25 is statistically distinguishable from Ro_c1≈2.7 and that the same simulations reproduce the known Nu(Ro) transition at the higher value. Without error bars on the acceleration moments or a direct overlay of Nu(Ro) and acceleration(Ro) from identical runs, the claim of two distinct critical values remains under-supported.
minor comments (2)
- [Abstract] Abstract: the notation “Prandtl--Blasius” uses an en-dash; standard usage is “Prandtl–Blasius.”
- [Abstract] Abstract: the parenthetical remark “both values for the parameter settings of our present numerical study” should be expanded to state explicitly that Ro_c1 was also measured in the same runs.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. We will revise the paper to address the concerns regarding numerical details and the presentation of the two transitions. Below we respond point by point to the major comments.
read point-by-point responses
-
Referee: [Abstract / Numerical methods] Abstract and numerical-setup section: no grid resolution, tracer count, integration length, spatial definition of “near the top plate,” or convergence tests are supplied. Because the central claim rests on the sudden jump in rms and kurtosis being physical rather than a sampling or binning artifact, these quantities must be documented and shown to be insensitive to reasonable variations in tracer ensemble size and location.
Authors: We agree that the numerical parameters and convergence information are necessary to substantiate the claim. The revised manuscript will include these details in the methods section: grid resolution, tracer count, integration length, the spatial definition of the near-plate region, and convergence tests with respect to tracer number and binning. These additions will demonstrate that the jump at Ro_c2 is robust and not due to sampling artifacts. revision: yes
-
Referee: [Results (tracer statistics)] Results section describing the acceleration statistics: the manuscript must demonstrate that the location of the jump at Ro≈2.25 is statistically distinguishable from Ro_c1≈2.7 and that the same simulations reproduce the known Nu(Ro) transition at the higher value. Without error bars on the acceleration moments or a direct overlay of Nu(Ro) and acceleration(Ro) from identical runs, the claim of two distinct critical values remains under-supported.
Authors: The manuscript presents data from the same simulations for both the Nusselt number and the Lagrangian statistics. However, to better demonstrate the distinction between the two critical values, we will add error bars to the acceleration statistics and a direct comparison plot of Nu(Ro) and the acceleration moments in the revised version. This will show that the transitions occur at different Rossby numbers. revision: yes
Circularity Check
No circularity: critical values extracted directly from simulation data
full rationale
The paper reports an empirical observation from direct numerical simulations: a sudden increase in rms and kurtosis of horizontal Lagrangian acceleration at Ro_c2≈2.25, identified as a distinct transition point separate from the Nu transition at Ro_c1≈2.7. This is presented as a numerical finding from tracer statistics near the top plate, with no equations, fitted parameters, or self-citations that define the reported critical value in terms of itself or reduce the claim to a renaming or construction from prior inputs. The derivation chain consists of simulation output analysis rather than any self-referential structure.
Axiom & Free-Parameter Ledger
free parameters (3)
- Rayleigh number Ra =
1.3e9
- Prandtl number Pr =
6.7
- Aspect ratio Gamma =
1
axioms (2)
- standard math The flow obeys the incompressible Navier-Stokes equations under the Boussinesq approximation.
- domain assumption Lagrangian tracers faithfully sample the acceleration field of the underlying Eulerian flow near the top plate.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.lean (Jcost uniqueness)washburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
by numerically studying Lagrangian acceleration statistics... rms values and the kurtosis of the horizontal acceleration... sudden increase at Ro_c2≈2.25
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.
Reference graph
Works this paper leans on
-
[1]
Root-mean-square values of acceleration The rms values of the horizontal and vertical acceleration are computed asarms i = √ ⟨(ai−⟨ai⟩)2⟩, wherei =xy or i =z and the average is taken over time and over the top and center measurement volumes, respectively, as sketched in figure 1. For the horizontal component, arms xy , the average is taken over a statistic...
-
[2]
Kurtosis and skewness of acceleration The probability density functions (PDFs) of acceleration in turbulent flows are characterized by exponential tails [47, 51]. In (rotating) RBC, the coherent flow structures are influencing the tails of the acceleration PDFs [18]. The relative importance of the tails in such a PDF can be quantified by the kurtosis, Ki =⟨(a...
work page 2020
-
[3]
J.-C. Gascard, A.J. Watson, M.-J. Messias, K.A. Olsson, T. Johannessen, and K. Simonsen, Long-lived vortices as a mode of deep ventilation in the Greenland Sea, Nature 416, 525 (2002)
work page 2002
-
[4]
J. Marshall and F. Schott, Open-ocean convection: Observations, theory, and models, Rev. Geophys. 37, 1 (1999)
work page 1999
-
[5]
Hadley, Concerning the cause of the general trade-winds, Phil
G. Hadley, Concerning the cause of the general trade-winds, Phil. Trans. 39, 58 (1753)
-
[6]
K.A. Emanuel, J.D. Neelin, and C.S. Bretherton, On largescale circulations in convecting atmospheres, Q. J. Royal Meteorol. Soc. 120, 1111 (1994)
work page 1994
- [7]
-
[8]
D. Lohse and K.-Q. Xia, Small-scale properties of turbulent Rayleigh-B´ enard convection, Annu. Rev. Fluid Mech.42, 335 (2010)
work page 2010
-
[9]
F. Chilla and J. Schumacher, New perspectives in turbulent Rayleigh-B´ enard convection, Eur. J. Phys. E 35, 58 (2012)
work page 2012
-
[10]
Pope, Turbulent Flows (Cambridge University Press, Cambridge, 2000)
S.B. Pope, Turbulent Flows (Cambridge University Press, Cambridge, 2000)
work page 2000
-
[11]
Y. Liu and R.E. Ecke, Heat transport scaling in turbulent Rayleigh–B´ enard convection: Effects of rotation and Prandtl number, Phys. Rev. Lett. 79, 2257 (1997)
work page 1997
-
[12]
R.P.J. Kunnen, H.J.H. Clercx, and B.J. Geurts, Breakdown of large-scale circulation in turbulent rotating convection, Europhys. Lett. 84, 24001 (2008)
work page 2008
-
[13]
R.J.A.M. Stevens, J.-Q. Zhong, H.J.H. Clercx, G. Ahlers, and D. Lohse, Transitions between turbulent states in rotating Rayleigh–B´ enard convection, Phys. Rev. Lett.103, 024503 (2009)
work page 2009
-
[14]
J.-Q. Zhong, R.J.A.M. Stevens, H.J.H. Clercx, R. Verzicco, D. Lohse, and G. Ahlers, Prandtl-, Rayleigh-, and Rossby- number dependence of heat transport in turbulent rotating Rayleigh–B´ enard convection, Phys. Rev. Lett. 102, 044502 (2009). 15
work page 2009
-
[15]
R.P.J. Kunnen, B.J. Geurts, H.J.H. Clercx, Experimental and numerical investigation of turbulent convection in a rotating cylinder, J. Fluid Mech. 642, 445 (2010)
work page 2010
-
[16]
S. Weiss, R.J.A.M. Stevens, J.-Q. Zhong, H.J.H. Clercx, D. Lohse, and G. Ahlers, Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh–B´ enard convection, Phys. Rev. Lett.105, 224501 (2010)
work page 2010
-
[17]
J.-Q. Zhong and G. Ahlers, Heat transport and the large-scale circulation in rotating turbulent Rayleigh-B´ enard convection, J. Fluid Mech. 665, 300 (2010)
work page 2010
-
[18]
R.J.A.M. Stevens, J. Overkamp, D. Lohse, and H.J.H. Clercx, Effect of aspect ratio on vortex distribution and heat transfer in rotating Rayleigh–B´ enard convection, Phys. Rev. E84, 056313 (2011)
work page 2011
-
[19]
P. Wei, S. Weiss, and G. Ahlers, Multiple transitions in rotating turbulent Rayleigh–B´ enard convection, Phys. Rev. Lett. 114, 114506 (2015)
work page 2015
- [20]
-
[21]
K.L. Chong, Y. Yang, S.-D. Huang, J.-Q. Zhong, R.J.A.M. Stevens, R. Verzicco, D. Lohse, and K.-Q. Xia, Confined Rayleigh–B´ enard, rotating Rayleigh-B´ enard, and double diffusive convection: a unifying view on turbulent transport enhancement through coherent structure manipulation, Phys. Rev. Lett. 119, 064501 (2017)
work page 2017
-
[22]
G.S. Lewis and H.L. Swinney, Velocity structure functions, scaling, and transitions in high-Reynolds-number Couette– Taylor flow, Phys. Rev. E 59, 5457 (1999)
work page 1999
-
[23]
S. Grossmann, D. Lohse, and C. Sun, High-Reynolds number Taylor–Couette turbulence, Annu. Rev. Fluid Mech. 48, 53 (2016)
work page 2016
-
[24]
F. Ravelet, L. Mari´ e, A. Chiffaudel, and F. Daviaud, Multistability and memory effect in a highly turbulent flow: experi- mental evidence for a global bifurcation, Phys. Rev. Lett. 93, 164501 (2004)
work page 2004
-
[25]
F. Ravelet, A. Chiffaudel, and F. Daviaud, Supercritical transition to turbulence in an inertially driven von K´ arm´ an closed flow, J. Fluid Mech. 601, 339 (2008)
work page 2008
-
[26]
S.G. Huisman, R.C.A. van der Veen, C. Sun, and D. Lohse, Multiple states in highly turbulent Taylor–Couette flow, Nat. Commun. 5, 3820 (2014)
work page 2014
-
[27]
R.P.J. Kunnen, H.J.H. Clercx, and B.J. Geurts, Heat flux intensification by vortical flow localization in rotating convection, Phys. Rev. E 74, 056306 (2006)
work page 2006
-
[28]
Rossby, A study of B´ enard convection with and without rotation, J
H.T. Rossby, A study of B´ enard convection with and without rotation, J. Fluid Mech. 36, 309 (1969)
work page 1969
-
[29]
R.J.A.M. Stevens, H.J.H. Clercx, and D. Lohse, Optimal Prandtl number for heat transfer in rotating Rayleigh–B´ enard convection, New J. Phys. 12, 075005 (2010)
work page 2010
-
[30]
P. Vorobieff and R.E. Ecke, Vortex structure in rotating Rayleigh–B´ enard convection, Physica D123, 153 (1998)
work page 1998
-
[31]
R.P.J. Kunnen, R.J.A.M. Stevens, J. Overkamp, C. Sun, G.J.F. van Heijst, and H.J.H. Clercx, The role of Stewartson and Ekman layers in turbulent rotating Rayleigh–B´ enard convection, J. Fluid Mech.688, 422 (2011)
work page 2011
-
[32]
R.J.A.M. Stevens, H.J.H. Clercx, and D. Lohse, Boundary layers in rotating weakly turbulent Rayleigh–B´ enard convection, Phys. Fluids 22, 085103 (2010)
work page 2010
- [33]
- [34]
- [35]
-
[36]
S. Stellmach, M. Lischper, K. Julien, G. Vasil, J.S. Cheng, A. Ribeiro, E.M. King, and J.M. Aurnou, Approaching the asymptotic regime of rapidly rotating convection: boundary layers versus interior dynamics, Phys. Rev. Lett. 113, 254501 (2014)
work page 2014
-
[37]
R.J.A.M. Stevens, H.J.H. Clercx, and D. Lohse, Heat transport and flow structure in rotating Rayleigh–B´ enard convection, Eur. J. Mech. B/Fluids 40, 41-49 (2013)
work page 2013
-
[38]
O. Shishkina and C. Wagner, Analysis of sheet-like thermal plumes in turbulent Rayleigh–B´ enard convection, J. Fluid Mech. 599, 383 (2008)
work page 2008
-
[39]
E.M. King, S. Stellmach, and J.M. Aurnou, Heat transfer by rapidly rotating Rayleigh-B´ enard convection, J. Fluid Mech. 691, 568 (2012)
work page 2012
- [40]
-
[41]
S. Horn and O. Shishkina, Rotating non-Oberbeck Boussinesq Rayleigh-B´ enard convection in water, Phys. Fluids 26, 055111 (2014)
work page 2014
-
[42]
R.E. Ecke and J.J. Niemela, Heat transport in the geostrophic regime of rotating Rayleigh-B´ enard convection, Phys. Rev. Lett. 113, 114301 (2014)
work page 2014
-
[43]
H. Rajaei, K.M.J. Alards, R.P.J. Kunnen, and H.J.H. Clercx, Velocity and acceleration statistics in rapidly rotating Rayleigh-B´ enard convection, J. Fluid Mech.857, 374 (2018)
work page 2018
-
[44]
M.A. Green, C.W. Rowley, and G. Haller, Detection of Lagrangian coherent structures in three-dimensional turbulence, J. Fluid Mech. 572, 111 (2007)
work page 2007
-
[45]
Schumacher, Lagrangian studies in convective turbulence, Phys
J. Schumacher, Lagrangian studies in convective turbulence, Phys. Rev. E 79, 056301 (2009)
work page 2009
-
[46]
Y. Choi, Y. Park, and C. Lee, Helicity and geometric nature of particle trajectories in homogeneous isotropic turbulence, Int. J. Heat Fluid Flow 31, 482 (2010)
work page 2010
-
[47]
L. Biferale, G. Boffetta, A. Celani, B.J. Devenish, A. Lanotte, and F. Toschi, Multifractal statistics of Lagrangian velocity 16 and acceleration in turbulence, Phys. Rev. Lett. 93, 064502 (2004)
work page 2004
-
[48]
F. Moisy and J. Jim´ enez, Geometry and clustering of intense structures in isotropic turbulence, J. Fluid Mech. 513, 111 (2004)
work page 2004
-
[49]
G.A. Voth, A. La Porta, A.M. Crawford, J. Alexander, and E. Bodenschatz, Measurement of particle accelerations in fully developed turbulence, J. Fluid Mech. 469, 121 (2002)
work page 2002
-
[50]
L. Del Castello and H.J.H. Clercx, Lagrangian acceleration of passive tracers in statistically steady rotating turbulence, Phys. Rev. Lett. 107, 214502 (2011)
work page 2011
-
[51]
R. Verzicco and R. Camussi, Numerical experiments on strongly turbulent thermal convection in a slender cylindrical cell, J. Fluid Mech. 477, 19 (2003)
work page 2003
-
[52]
R. Verzicco and P. Orlandi, A finite-difference scheme for three-dimensional incompressible flows in cylindrical coordinates, J. Comput. Phys. 123, 402 (1996)
work page 1996
-
[53]
F. Toschi and E. Bodenschatz, Lagrangian properties of particles in turbulence, Annu. Rev. Fluid Mech. 41, 375 (2009)
work page 2009
- [54]
-
[55]
Q. Zhou, C. Sun, and K.-Q. Xia, Experimental investigation of homogeneity, isotropy, and circulation of the velocity field in buoyancy-driven turbulence, J. Fluid Mech. 598, 361 (2008)
work page 2008
-
[56]
O. Shishkina, S. Horn, S. Wagner, Falkner–Skan boundary layer approximation in Rayleigh–B´ enard convection, J. Fluid Mech. 730, 442 (2013)
work page 2013
-
[57]
S. Sterl, H.-M. Li, and J.-Q. Zhong, Dynamical and statistical phenomena of circulation and heat transfer in periodically forced rotating turbulent Rayleigh–B´ enard convection, Phys. Rev. Fluids1, 084401 (2016)
work page 2016
-
[58]
H. Schlichting and K. Gersten, Boundary-layer theory (Springer-Verlag Berlin, Heidelberg, 2000)
work page 2000
-
[59]
P.K. Kundu and I.M. Cohen, Fluid Mechanics (Academic Press, London, 2002)
work page 2002
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.