Suppressing wall modes in confined rotating turbulent convection
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel 2026-07-08 08:40 UTCglm-5.2pith:F7F7UWWOrecord.jsonopen to challenge →
The pith
Ring barriers suppress wall modes 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
The central discovery is that the suppression of wall modes by sidewall barriers is controlled by a scale-matching condition: the barrier width must exceed the inner Stewartson layer thickness (scaling as Ekman number to the one-third power) but remain below the outer wall-mode extent (scaling as Ekman number to the one-quarter power). When this condition is met, a single barrier eliminates the dominant azimuthal mode, suppresses near-wall velocity peaks, removes the boundary zonal flow, and recovers heat-transport values close to those of laterally periodic domains. The barriers also induce a secondary baroclinic flow caused by isotherm-isobar misalignment at the barrier, which is a direct,
What carries the argument
wall modes
If this is right
- Experiments in rapidly rotating convection can now isolate bulk dynamics from sidewall contamination by installing properly scaled ring barriers, enabling closer comparison with idealized periodic-domain simulations.
- The scale-matching condition (barrier width between Ek^{1/3} and Ek^{1/4} times the cell height) provides a concrete design rule for future experimental facilities targeting geostrophic convection regimes.
- The baroclinic flow induced by thermal-conductivity mismatch suggests that barrier material choice is a critical design parameter; high-conductivity barriers reduce secondary flows but may introduce other thermal artifacts.
- At extreme Rayleigh numbers, single barriers become insufficient and multiple barriers or stronger rotational confinement are needed, indicating that suppression strategies must be tailored to the specific region of parameter space.
Where Pith is reading between the lines
- If wall modes are topologically protected (as prior work suggests), the barriers likely do not destroy them but rather prevent their nonlinear saturation and spatial organization, effectively trapping them below threshold — this distinction between suppression and elimination is not fully resolved in the paper.
- The baroclinic flow asymmetry between DNS and experiments, attributed to centrifugal buoyancy, implies that even small physical effects normally neglected in simulations can qualitatively change local flow topology near obstacles, raising questions about the fidelity of volume-penalization methods for other barrier geometries.
- The finding that barrier effectiveness depends on the ratio of barrier width to wall-mode scale suggests a potential scaling law for the minimum number of barriers needed as a function of Rayleigh and Ekman numbers, though the paper does not derive such a law explicitly.
Load-bearing premise
The direct numerical simulations assume the volume penalization method with a fixed penalization strength accurately represents the fluid-barrier interaction, and that neglecting centrifugal buoyancy and assuming the barrier has the same thermal conductivity as the fluid are sufficient simplifications — yet the DNS-experiment mismatch in the baroclinic flow asymmetry shows these second-order effects qualitatively change the local dynamics near the barrier.
What would settle it
If barriers with width satisfying the stated scale-matching condition fail to reduce near-wall vertical velocity peaks or heat-transport bias in an independent experimental apparatus at the same parameters, the suppression mechanism would not be robust.
Figures
read the original abstract
In confined turbulent rotating convection, the largest vertical velocities are found near the sidewalls in the form of wave-like structures known as wall modes. These structures persist deep into the turbulent regime, bias heat transport, and disrupt bulk flow organisation through radial jets. Controlling or suppressing wall modes is, therefore, essential for accessing bulk dynamics free from wall-induced effects. Here, we combine experiments and direct numerical simulations to investigate wall modes control in cylindrical cells equipped with ring-shaped sidewall barriers. Barriers suppress vertical-velocity maxima near the sidewall and disrupt the characteristic wave-like pattern. Simulations further show that the barriers reduce the wall-mode-induced enhancement of heat transport, shifting it towards values characteristic of laterally periodic domains. The suppression efficiency is governed by the ratio of the barrier width to the wall-mode scale and is enhanced by the addition of a second barrier. In the horizontal plane, radial jet ejections are attenuated, while the time-averaged flow reveals suppression of the boundary zonal flow (BZF), a ring-shaped region of positive azimuthal velocity near the sidewall, provided measurements are taken away from the immediate vicinity of the barriers. In this region, isotherms bend toward the poorly conducting barrier, creating a local misalignment with the isobars and inducing a baroclinic flow adjacent to the barrier faces. This effect weakens with increasing barrier conductivity or smoother geometry. These results demonstrate that sidewall barriers provide a robust route for suppressing wall modes signatures in experimental turbulent rotating convection, while locally inducing secondary baroclinic flows near the barriers. Their use enables access to extreme rotating-convection regimes with reduced sidewall influence.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript investigates the suppression of wall modes in confined rotating Rayleigh–Bénard convection using ring-shaped sidewall barriers, combining PIV experiments and DNS. The central claim is that barriers suppress the characteristic signatures of wall modes—near-wall vertical velocity peaks, heat-transport bias, radial jet ejections, and the boundary zonal flow—when the barrier width exceeds the Stewartson inner-layer scale δ_m ~ Ek^{1/3}. The paper also characterizes a secondary baroclinic flow induced near the barriers due to thermal conductivity mismatch and proposes mitigation strategies. The experimental and DNS results show good quantitative agreement for the vertical, radial, and azimuthal velocity components, while the heat-transport results are DNS-only.
Significance. The problem addressed is well-motivated: wall modes contaminate heat-transport measurements and bulk flow organization in rotating convection experiments, and a practical suppression strategy is needed. The combination of laboratory experiments and DNS across multiple velocity components is a strength, as is the identification and partial mitigation of the secondary baroclinic flow. The scaling criterion d > δ_m = Ek^{1/3} provides a falsifiable, parameter-free design rule grounded in established Stewartson-layer theory. The extension to extreme parameter regimes (Appendix A) and the systematic variation of barrier width (Fig. 4, Fig. 5d) add value. The paper builds directly on the numerical proposal of Terrien et al. (2023) and provides the first experimental validation, which is a meaningful contribution to the field.
major comments (2)
- §3.2, Fig. 5: The heat-transport claim (pillar 2 of the central result) is supported by DNS alone; no experimental Nu measurements are presented. The DNS treats the barrier thermal diffusivity as identical to the fluid's (Eq. 2.4, §2.2), while the experimental Plexiglas barriers have κ approximately κ_water/3. The paper itself demonstrates in §3.5 that this conductivity mismatch substantially alters local barrier dynamics (isotherm bending, baroclinic flow intensity changing with aluminium vs. Plexiglas barriers, Fig. 9). Since the heat-transport modification depends on how the barrier redistributes convective transport in the near-wall region—precisely where the thermal-conductivity simplification matters most—the Nu results in Fig. 5 carry an unquantified systematic uncertainty. The authors should either (i) provide an estimate of the sensitivity of Nu to the barrier conductivity (e.g.
- a DNS comparison with a contrasting conductivity, even if approximate), or (ii) explicitly state this limitation in §3.2 and temper the claim that Nu shifts toward 'values characteristic of laterally periodic domains.' As written, the abstract and conclusion state the heat-transport shift as an established result without flagging that it rests on an idealized barrier model whose simplifications the paper itself shows to be consequential for local dynamics.
minor comments (8)
- §2.2: The penalization strength ζ=100 is stated as adequate based on a previous study (Martínez-Ortíz et al. 2026). A brief note on whether convergence with respect to ζ was checked for the present barrier geometry would strengthen confidence, since the barrier–fluid interaction is the central object of study.
- Fig. 3(b): The experimental and DNS profiles are at slightly different Ek (1.4×10⁻⁶ vs. 1.13×10⁻⁶) and different Ra. The agreement is described as 'excellent' but the parameter mismatch should be noted explicitly.
- §3.5, Fig. 8(a): The scaling δ_b ~ Ek^{0.20} is based on only four data points. The exponent should be labeled as preliminary, and the range of Ek over which it was measured should be stated. It is unclear whether this is a fitted exponent or merely empirical.
- Fig. 5(d): The barrier width axis uses d/D, while the scaling criterion is expressed in terms of d/δ_m or d/H. Adding a secondary axis in units of d/δ_m would make the figure self-consistent with the text and more interpretable.
- Appendix B, Fig. 11(f): The caption states 'Note the difference in Ek in bot cases' — the DNS and experimental Ek values used for this comparison should be stated explicitly in the caption or text, since the comparison involves centrifugal buoyancy effects that depend on Fr.
- §3.1: The statement 'predicting an exact relation for the number and characteristics of the barriers is not straightforward' is honest but could be sharpened. The paper does provide a useful criterion (δ_m < d < δ); clarifying whether this is intended as a necessary condition, a sufficient condition, or a heuristic would help the reader.
- Table 1: The grid resolution column uses notation like '193×385×[577,769]²' which is ambiguous. Clarify whether the bracket notation indicates two grid levels or a range.
- The abstract states 'shifting it towards values characteristic of laterally periodic domains' without qualification. Given that this is a DNS-only result with the conductivity simplification noted above, a qualifier such as 'in DNS' or 'under the idealized barrier model' would be appropriate.
Simulated Author's Rebuttal
We thank the referee for a careful and constructive report. The referee's single major comment concerns the heat-transport results in §3.2, which are DNS-only and rely on a barrier thermal diffusivity matched to the fluid rather than to the experimental Plexiglas. We agree this limitation should be made explicit and are happy to revise the manuscript accordingly.
read point-by-point responses
-
Referee: §3.2, Fig. 5: The heat-transport claim (pillar 2 of the central result) is supported by DNS alone; no experimental Nu measurements are presented. The DNS treats the barrier thermal diffusivity as identical to the fluid's (Eq. 2.4, §2.2), while the experimental Plexiglas barriers have κ approximately κ_water/3. The paper itself demonstrates in §3.5 that this conductivity mismatch substantially alters local barrier dynamics (isotherm bending, baroclinic flow intensity changing with aluminium vs. Plexiglas barriers, Fig. 9). Since the heat-transport modification depends on how the barrier redistributes convective transport in the near-wall region—precisely where the thermal-conductivity simplification matters most—the Nu results in Fig. 5 carry an unquantified systematic uncertainty. The authors should either (i) provide an estimate of the sensitivity of Nu to the barrier conductivity (e.g.
Authors: We thank the referee for raising this important point, which we fully acknowledge as a genuine limitation of the current study. The referee is correct on both counts: the Nu results in Fig. 5 are DNS-only, and the DNS assumes the barrier thermal diffusivity is identical to that of the fluid (Eq. 2.4), whereas the experimental Plexiglas barriers have a thermal conductivity roughly one-third that of water. We also agree that §3.5 demonstrates that the conductivity mismatch has consequential local effects on the isotherm bending and baroclinic flow intensity near the barrier. We address the referee's two suggested remedies in turn. Regarding option (i) — a DNS comparison with contrasting conductivity: we have in fact already begun such a comparison. In Appendix B, we include DNS with centrifugal buoyancy included (varying Froude number), which was necessary to recover the experimentally observed asymmetry in the baroclinic flow above and below the barrier. In that appendix, we explicitly note that differences in thermal diffusivity between the barrier and the fluid are an additional mechanism not captured in the main-text DNS and likely contribute to the observed symmetry breaking. However, a full DNS parameter study varying the barrier thermal diffusivity independently is not feasible within the revision timeframe, as each additional simulation at the relevant parameters (Ra ~ 10^9–10^10, Ek ~ 10^-6) requires substantial computational resources, and the volume penalisation method as currently implemented does not straightforwardly support a spatially varying thermal diffusivity field within the penalised region. We therefore adopt option (ii). We will revise §3.2 to explicitly state that the Nu results rest on an idealised barrier model in which the thermal diffusivity is revision: no
Circularity Check
No significant circularity: the paper's central claims are grounded in independent experimental measurements and DNS, with the scaling relation derived from established Stewartson layer theory rather than postulated to make results work.
full rationale
The paper's central claim — that sidewall barriers suppress wall modes — is supported by independent experimental PIV measurements (Figs. 2, 3b, 6, 7) and DNS (Figs. 3a, 4, 5), with quantitative agreement between the two. The scaling relation d > δ_m = Ek^{1/3} is derived from established Stewartson layer theory (Stewartson 1957, 1966; Kunnen et al. 2011), not defined in terms of the suppression result. The heat-transport shift toward periodic-domain Nu (Fig. 5) is DNS-only, but this is a limitation of experimental coverage, not circularity: the DNS solves the full Navier-Stokes equations (Eq. 2.2-2.4) with independently specified boundary conditions, and the Nu is computed from the resulting fields, not fitted to a target. The baroclinic flow mechanism (§3.5) is explained physically via isotherm-isobar misalignment (∇p × ∇θ), not defined circularly. The one self-citation (Martínez-Ortíz et al. 2026) merely justifies the penalization strength ζ=100, a numerical parameter choice that does not load-bear the central physical claims. The DNS-experiment mismatch in baroclinic flow asymmetry (§3.5, Appendix B) is openly acknowledged and investigated, not hidden. No step in the derivation chain reduces to its inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (3)
- Barrier width d =
0.0125 to 0.2 (d/D)
- Number of barriers =
1 or 2
- Penalisation strength ζ =
100
axioms (3)
- domain assumption Boussinesq approximation
- standard math Stewartson layer scaling (δ_m/H = Ek^(1/3), δ/H = Ek^(1/4))
- domain assumption Volume penalization method accurately represents solid barriers
Reference graph
Works this paper leans on
-
[1]
J. M. Aurnou and M. A. Calkins and J. S. Cheng and K. Julien and E. M. King and D. Nieves and K. M. Soderlund and S. Stellmach , title =. Phys. Earth Planet. Inter. , year =
- [2]
- [3]
-
[4]
R. P. J. Kunnen , title =. J. Turbul. , year =
- [5]
-
[6]
H. T. Rossby , title =. J. Fluid Mech. , year =
- [7]
- [8]
- [9]
-
[10]
A. E. Gill , title =
-
[11]
S. Stellmach and M. Lischper and K. Julien and G. Vasil and J. S. Cheng and A. Ribeiro and E. M. King and J. M. Aurnou , title =. Phys. Rev. Lett. , year =
-
[12]
A. M. Rubio and K. Julien and E. Knobloch and J. B. Weiss , title =. Phys. Rev. Lett. , year =
-
[13]
R. P. J. Kunnen and R. Ostilla-M\'onico and E. P. van der Poel and R. Verzicco and D. Lohse , title =. J. Fluid Mech. , year =
-
[14]
B. Favier and L. J. Silvers and M. R. E. Proctor , title =. Phys. Fluids , year =
-
[15]
Guervilly, C. and Hughes, D. W. and Jones, C. A. , title =. Journal of Fluid Mechanics , year =
- [16]
-
[17]
A. J. Competition between. Phys. Rev. Lett. , year =
-
[18]
A. van Kan and K. Julien and B. Miquel and E. Knobloch , title =. J. Fluid Mech. , year =
- [19]
-
[20]
R. P. J. Kunnen and R. J. A. M. Stevens and J. Overkamp and C. Sun and G. J. F. van Heijst and H. J. H. Clercx , title =. J. Fluid Mech. , year =
- [21]
-
[22]
Y. B. Hu and Y. C. Xie and K. Q. Xia , title =. J. Fluid Mech. , year =
-
[23]
Convection in a rotating cylinder
Goldstein, HF and Knobloch, E and Mercader, I and Net, M , journal=. Convection in a rotating cylinder. 1993 , publisher=
work page 1993
- [24]
-
[25]
Asymmetric modes and the transition to vortex structures in rotating
Zhong, Fang and Ecke, Robert and Steinberg, Victor , journal=. Asymmetric modes and the transition to vortex structures in rotating. 1991 , publisher=
work page 1991
-
[26]
R. E. Ecke and F. Zhong and E. Knobloch , title =. Europhys. Lett. , year =
-
[27]
X. Zhang and P. Reiter and O. Shishkina and R. E. Ecke , title =. Phys. Rev. Fluids , year =
- [28]
-
[29]
X. Zhang and D. P. M. van Gils and S. Horn and M. Wedi and L. Zwirner and G. Ahlers and R. E. Ecke and S. Weiss and E. Bodenschatz and O. Shishkina , title =. Phys. Rev. Lett. , year =
-
[30]
X. M. de Wit and A. J. Aguirre Guzm\'an and M. Madonia and J. S. Cheng and H. J. H. Clercx and R. P. J. Kunnen , title =. Phys. Rev. Fluids , year =
- [31]
-
[32]
M. Wedi and V. M. Moturi and D. Funfschilling and S. Weiss , title =. J. Fluid Mech. , year =
- [33]
- [34]
- [35]
-
[36]
X. M. de Wit and W. J. M. Boot and M. Madonia and A. J. Aguirre Guzm\'an and R. P. J. Kunnen , title =. Phys. Rev. Fluids , year =
-
[37]
Velocimetry in rapidly rotating convection: spatial correlations, flow structures and length scales , author=. Europhysics Letters , volume=. 2021 , publisher=
work page 2021
- [38]
-
[39]
L. Terrien and B. Favier and E. Knobloch , title =. Phys. Rev. Lett. , year =
-
[40]
R. E. Ecke and O. Shishkina , title =. Annu. Rev. Fluid Mech. , year =
- [41]
-
[42]
M. Raffel and C. E. Willert and F. Scarano and C. J. K. Particle Image Velocimetry: A Practical Guide , publisher =
- [43]
- [44]
-
[45]
Journal of Computational Physics , volume=
A volume penalization method for incompressible flows and scalar advection--diffusion with moving obstacles , author=. Journal of Computational Physics , volume=. 2012 , publisher=
work page 2012
-
[46]
T. Engels and D. Kolomenskiy and K. Schneider and J. Sesterhenn , title =. J. Comput. Phys. , year =
-
[47]
L. Mart. Persistent clock of turbulent thermal convection , journal =. 2026 , volume =
work page 2026
- [48]
-
[49]
M. Madonia and A. J. Reynolds number scaling and energy spectra in geostrophic convection , journal =. 2023 , volume =
work page 2023
-
[50]
J. S. Cheng and M. Madonia and A. J. Laboratory exploration of heat transfer regimes in rapidly rotating turbulent convection , journal =. 2020 , volume =
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.