Pith. sign in

REVIEW 2 major objections 5 minor 48 references

Influence of Boundary Conditions and Heating Modes on the Onset of Columnar Convection in Rotating Spherical Shells

T0 review · 2 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash

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

Pith's one-line read No-slip boundaries trigger convection before stress-free ones in thick rotating shells or at high Prandtl numbers.

desk verdict A careful, well-validated linear-stability scan mapping how boundary-condition preference depends on geometry and heating; the continuation-restricted m-search deserves a robustness check near transitions. read the letter →

arxiv 2509.06632 v1 pith:XBP2ANEY submitted 2025-09-08 physics.flu-dyn astro-ph.EPastro-ph.SRphysics.geo-ph

William Seeley , Francesca Coke , Radostin Simitev , Robert Teed This is my paper · ORCID
classification physics.flu-dynastro-ph.EPastro-ph.SRphysics.geo-ph
keywords rotatingsphericalshellsthermalconvectionboundaryconditionsEkmanlayerinternalheatingcriticalRayleighnumberPrandtlcolumnar
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 asks which mechanical boundary conditions and which heating mechanism make thermal convection easiest to start in a rotating spherical shell, the canonical model for fluid motion in planetary cores and stellar interiors. The central finding is that the preferred boundary condition is not a fixed modelling choice: for thick shells or large Prandtl numbers, the Ekman boundary layer at the outer surface becomes destabilising, so no-slip walls give a lower critical Rayleigh number than stress-free walls, while the reverse holds elsewhere. It also finds that internal heating, compared with differential heating, raises the critical Rayleigh number, selects larger azimuthal wavenumbers and oscillation frequencies, and moves the first convective columns away from the tangent cylinder. If these results hold, numerical models of planetary and stellar convection need to justify their boundary-condition and heating choices rather than treating either as a neutral convention.

What carries the argument

The argument is carried by the linear stability problem of the rotating-shell equations, solved by a spectral expansion in Chebyshev polynomials and spherical harmonics and by a continuation in parameter space that tracks the marginal mode with the largest growth rate. The physical mechanism invoked is the Ekman boundary layer, whose thickness scales as $E^{1/2}$: it restrains the flow at low Prandtl number but becomes destabilising for thick shells or large Prandtl numbers, reversing which boundary condition lowers the critical Rayleigh number. A second structural element is the critical cylinder radius, located by the maximum of the azimuthally averaged radial velocity in the equatorial plane; internal heating detaches the onset from the tangent cylinder and moves that radius with Prandtl number and Ekman number.

What would settle it

Re-run the linear stability search at parameter values straddling the reported transitions, for example with Prandtl number near 1.7 in differential heating and radius ratio near 0.4 or inside the boundary-condition reversal region in internal heating at Ekman number $3\times10^{-5}$, scanning all azimuthal wavenumbers from several independent starting guesses instead of a small neighbourhood of the previous critical mode; any mismatch would show the reported critical Rayleigh number is not the global value.

Watch

Extended reading notes

Core claim

At linear onset in a Boussinesq rotating spherical shell at moderate Ekman numbers ($10^{-4}$ to $10^{-5}$), the critical Rayleigh number, azimuthal wavenumber, and oscillation frequency are computed for Prandtl numbers 0.1 to 100, radius ratios 0.2 to 0.8, differential and internal heating, and no-slip, stress-free, and mixed boundary conditions. The paper's main discovery is that the outer Ekman boundary layer can be either stabilising or destabilising: under differential heating, no-slip boundaries give a lower critical Rayleigh number than stress-free boundaries once the Prandtl number exceeds about 1.7, and under internal heating the transition occurs near a Prandtl number of about 0.9 and also appears as a function of radius ratio, with the relative difference between the two critical values changing sign in an Ekman-number-dependent way. Internal heating raises the critical Rayleigh number by up to an order of magnitude, increases the wavenumber and frequency at onset, and places the critical convection column at a radius that departs from the previously predicted asymptotic value when the Prandtl number differs from 1. Mixed boundaries with no-slip on the inner sphere and stress-free outside closely reproduce the purely stress-free results, because the convection columns interact mainly with the outer surface.

Load-bearing premise

The calculation assumes that the most unstable azimuthal wavenumber changes only gradually as parameters move in small steps; if the true minimum jumps abruptly between steps, the reported critical Rayleigh number could be a local rather than global minimum.

Editorial extensions

If this is right

  • Planetary and stellar convection models must treat no-slip versus stress-free boundaries as a physical parameter with a measurable effect, not a numerical convenience: at thick shells or high Prandtl numbers, no-slip onset is easier.
  • Models using no-slip inner and stress-free outer boundaries can, near onset, be approximated by purely stress-free conditions because the outer boundary dominates.
  • Internal heating stabilises the shell relative to differential heating, so planets or stars with strong volumetric heating require stronger thermal forcing to begin convection.
  • At Ekman numbers relevant to planetary interiors, the sign of the boundary-condition preference depends on shell thickness under internal heating, with no-slip favoured for small radius ratios and stress-free for large ones.
  • The reported critical values across the parameter grid provide a benchmark for nonlinear simulations, whose onset is expected to fall close to these linear critical parameters.

Reading between the lines

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

  • If the outer boundary dominates onset, then for stars modelled with mixed boundary conditions the stress-free approximation may also hold in weakly nonlinear regimes, though whether it survives into strongly nonlinear dynamo states is untested.
  • The Prandtl-number drift of the critical radius found here suggests the asymptotic global theory's fixed value is only the unit-Prandtl-number case; a Prandtl-dependent correction to that theory would be a natural next step.
  • Because the trend of the critical Rayleigh number with radius ratio reverses when the outer radius rather than the shell thickness sets the length scale, comparisons of convective ease between planets with different inner-core sizes should use the outer-radius scaling to avoid a purely geometric artefact.
  • A combined heating model with both differential and internal sources would likely interpolate between the two preference maps, and the boundary-condition transition could then occur at intermediate Prandtl numbers and radius ratios.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The manuscript studies the linear onset of thermal convection in rotating spherical shells using a spectral method in spherical harmonics and Chebyshev polynomials. It systematically varies mechanical boundary conditions (no-slip, stress-free, and mixed with no-slip inner/stress-free outer) and heating modes (differential versus internal) over Pr ∈ [0.1, 100], χ ∈ [0.2, 0.8], and E ∈ {10^-4, 3×10^-5, 10^-5}. The main claims are that the preferred mechanical boundary condition at onset is not universal: for sufficiently thick shells or large Pr the no-slip outer boundary gives a lower critical Rayleigh number than the stress-free boundary, attributed to a destabilising Ekman layer; internal heating generally raises Rac, shifts onset to larger azimuthal wavenumbers and frequencies, and moves the critical column away from the tangent cylinder; and mixed boundary conditions with no-slip inner/stress-free outer closely mimic the purely stress-free case. The paper also analyses how the choice of nondimensionalisation can reverse apparent trends of Rac with χ, and provides fits for scaling laws.

Significance. The study is potentially valuable because it maps, over a wide parameter range, where the conventional choice of mechanical boundary condition in rotating-shell convection models changes the onset threshold. If correct, the results imply that the preferred boundary condition is a function of Pr, χ, E, and heating type rather than a fixed modelling convention, with direct implications for nonlinear convection and dynamo simulations. The numerical work is quantitatively anchored by comparisons with several prior studies (Dormy et al. 2004, Al-Shamali et al. 2004, Barik et al. 2023, Fan et al. 2024) and by convergence tests showing Rac converged to within 0.003% at the most demanding parameters (Appendix B); the data are openly available. The principal weakness is methodological: the global minimisation over the azimuthal wavenumber m relies on a continuation search restricted to a small neighbourhood of the previous m_c, and this is not tested at the reported mode transitions, where the boundary-condition preference reversals are decided by small differences in Rac.

major comments (2)
  1. [Section 2.2] The global critical mode is defined as Rac = min_m Ra_c^(m), but the implementation restricts the minimisation over m to 'a small neighbourhood of the previously identified mc' during continuation. Near reported transitions—for example the kink at χ≈0.4 in the internally heated case (Section 3.2.2, Figure 9b) and the sign changes of 1−Ra_c^SF/Ra_c^NS in Figure 10—the true minimizer may jump azimuthal wavenumber by more than that neighbourhood, so the reported Rac could be a local rather than global minimum. Because the preferred-boundary-condition conclusion is determined by the sign of a small difference between Ra_c^NS and Ra_c^SF, a missed global mode at these points could flip the ordering. The paper does not report any test that restarts the m-search from an unconstrained set or otherwise verifies the global minimum at these parameter values; the authors should add such a test and state whether the preference conclusions survive it.
  2. [Abstract, Section 3.1.1, Section 4] The mechanistic claim that the Ekman boundary layer at the outer boundary 'becomes destabilising' for sufficiently thick shells or large Pr is inferred solely from the inequality Ra_c^NS < Ra_c^SF. No direct computation or scaling analysis of the Ekman layer's contribution (e.g., boundary-layer suction or viscous dissipation) is provided, so the mechanism is an interpretation rather than an established result. The observed ordering of critical Rayleigh numbers is a robust numerical finding, but the causal language in the abstract overstates what the calculations show; either add a supporting boundary-layer argument or rephrase the claim as 'consistent with a destabilising role of the Ekman layer' rather than asserting it.
minor comments (5)
  1. [Section 3.2, Scaling laws and Table 3] Equation (11) defines C1 = χ^{-1}Rac, but the caption of Table 3 states C1 = Rac; the notation should be made consistent.
  2. [Section 3.2, Flow patterns] 'radius radio' should read 'radius ratio'.
  3. [Section 3.1.1] 'whether or nor such a situation' contains a typo ('nor' should be 'not').
  4. [Section 4, Remarks on nondimsionalisation] The heading and text use 'nondimsionalisation' and 'nondimensionalistaion'; these should be corrected to 'nondimensionalisation'.
  5. [Table 1 and Section 2.3] The text notes that frequency rescaling is needed when the timescale changes, but no explicit conversion formula for ω is given; adding the conversion to Table 1 would improve usability.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central onset computations are self-contained, independently benchmarked, and the only self-cited scaling laws are explicitly used as fits, not as predictions.

full rationale

The paper's central claims are the numerical determination of critical Rayleigh numbers, azimuthal wavenumbers, and frequencies from the linearized eigenvalue problem, Equations (3) and (5), under different boundary conditions and heating modes. These quantities are computed directly from the discretized governing equations, and the preference between no-slip and stress-free boundaries is read off from the computed ratio 1 - RaSF_c / RaNS_c, not from any fitted or assumed formula. The code is validated against independent prior studies including Dormy et al., Al-Shamali et al., Barik et al., Fan et al., and Christensen and Aubert, with conversions described in Section 2.3 and selected comparisons in Figures 2 and 7. The scaling laws in Equations (10) and (11) are explicitly introduced as interpolations of the paper's own best-fit curves, with the text saying 'To interpolate the best-fit curves of Figure 2, we use functions based on the scaling laws found in [7]'; they are not used to generate the critical values or to infer the boundary-condition preference, so they are not fitted inputs disguised as predictions. The only self-citation to the authors' earlier work appears in this curve-fitting context and in general references to prior pattern studies; it is not load-bearing for the physical conclusions. The continuation-restricted minimization over azimuthal wavenumber m noted in Section 2.2 is a possible numerical robustness issue near mode jumps, but it is a methodological risk rather than a circularity: it does not make any derived quantity equal to an input by construction. Overall, the derivation chain from the governing equations to the reported onset parameters is self-contained, and the externally benchmarked numerical results provide independent content for the paper's central claims.

Assumptions & free parameters 2 free parameters · 6 assumptions · 0 invented entities

The central numerical comparisons are computed from the discretized linear equations with no fitting to the claimed outcome. The listed free parameters are interpolation coefficients used only to draw smooth curves through computed points. The main physical assumptions are the Boussinesq model, the uniform volumetric heating idealization, and accepted asymptotic scalings from prior literature. No new entities are postulated.

free parameters (2)
  • Equation (10) fit coefficients a_i, b_i for Rac, mc, and |omega_c| versus Pr = Table 2: for E=1e-4, a1=1.556e5, b1=0.522, a2=6.938, b2=0.073, a3=1.365e2, b3=0.364; values vary with E
    Fitted to the computed critical values solely to draw best-fit curves in Figure 2; they do not enter the physical comparisons.
  • Equation (11) fit coefficients a_i, b_i, c_i for Rac, mc, and |omega_c| versus chi = Table 3: for E=1e-4, a1=3.407e5, b1=-2.355, c1=0.123, a2=4.172e2, b2=10.147, c2=-4.273, a3=9.751e4, b3=13.165…
    Fitted to the computed critical values to draw best-fit curves in Figure 7; they are not used in the central onset comparisons.
assumptions (6)
  • domain assumption Boussinesq approximation with constant material properties except buoyancy
    Adopted in Section 2.1; its validity for thick shells and internal heating in planetary and stellar contexts is assumed.
  • domain assumption Uniform volumetric heat source S and self-consistently determined boundary temperatures for internal heating
    Section 2.1 justifies this by convective mixing homogenizing heat sources, but it is a modeling idealization rather than a derived result.
  • standard math Busse symmetry and Schmidt-normalized spherical harmonic trial expansions
    Appendix A.2 relies on prior derivations that onset modes satisfy this symmetry, following Busse and Riahi.
  • domain assumption Asymptotic scalings Rac ~ E^-4/3, mc ~ E^-1/3, |omega_c| ~ E^-2/3 and the fitting forms of Equations (10) and (11)
    Used as expected trends and as templates for fitting curves; these scalings come from prior asymptotic theory rather than being derived here.
  • standard math Sturm-Liouville eigenvalue ordering and inverse power method convergence
    Section 2.2 and Appendix A inherit this numerical framework from Jones et al.; no proof is given in this paper.
  • domain assumption Asymptotic critical radii sM and sL from global and local theories
    Sections 3.1.2 and 3.2.2 compare computed critical radii with sM from Jones et al. and Dormy et al. and with sL from local theory; these prior results are taken as benchmarks.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Influence of Boundary Conditions and Heating Modes on the Onset of Columnar Convection in Rotating Spherical Shells." pith.science (2026). https://pith.science/paper/XBP2ANEY

@misc{pith2026250906632,
  author       = {Pith},
  title        = {Pith review of: Influence of Boundary Conditions and Heating Modes on the Onset of Columnar Convection in Rotating Spherical Shells},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XBP2ANEY}},
  note         = {Machine review of arXiv:2509.06632}
}
abstract

We investigate the linear onset of thermal convection in rotating spherical shells with a focus on the influence of mechanical boundary conditions and thermal driving modes. Using a spectral method, we determine critical Rayleigh numbers, azimuthal wavenumbers, and oscillation frequencies over a wide range of Prandtl numbers and shell aspect ratios at moderate Ekman numbers. We show that the preferred boundary condition for convective onset depends systematically on both aspect ratio and Prandtl number: for sufficiently thick shells or for large $\text{Pr}$, the Ekman boundary layer at the outer boundary becomes destabilising, so that no-slip boundaries yield a lower $\text{Ra}_c$ than stress-free boundaries. Comparing differential and internal heating, we find that internal heating generally raises $\text{Ra}_c$, shifts the onset to larger wavenumbers and frequencies, and relocates the critical column away from the tangent cylinder. Mixed boundary conditions with no-slip on the inner boundary behave similarly to purely stress-free boundaries, confirming the dominant influence of the outer surface. These results demonstrate that boundary conditions and heating mechanisms play a central role in controlling the onset of convection and should be carefully considered in models of planetary and stellar interiors.

Figures

Figures reproduced from arXiv: 2509.06632 by the authors.

Figure 1
Figure 1. Schematic diagram of the problem. The inner spherical surface at r = ri is held at temperature T = Ti and the outer spherical surface at r = ro is held at temperature T = To < Ti . The system rotates about the vertical (z-direction) with rotation rate, Ω. The gap width is d = ro − ri and gravity acts radially, g = −grˆr. Thermal driving. To model convection driven by thermal gradients, we adopt the Boussinesq approx… view at source ↗
Figure 2
Figure 2. Dependence of critical values on the Prandtl number, Pr, for several values of the Ekman number, E, in the case of differential heating with purely no-slip boundary conditions and radius ratio χ = 0.35. (a) Critical Rayleigh number, (b) critical frequency, (c) critical azimuthal wavenumber. Critical values found by Christensen and Aubert (☀) [12] and Fan et al (▲) [10] are also shown. The best-fit curves are given b… view at source ↗
Figure 3
Figure 3. The radial component of the velocity field at onset plotted in an equatorial slice for purely no-slip boundary conditions, χ = 0.35, differential heating (top half), internal heating (bottom half), E = 10−4 (left half of each panel), E = 10−5 (right half of each panel), and for (a) Pr = 0.1, (b) Pr = 1.0, and (c) Pr = 100.0. Red and blue represent positive and negative values, respectively. 3.1.1. Effect of Mechanic… view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: Comparison of critical parameter values as functions of the Prandtl number, Pr, for several values of the Ekman number, E, and different boundary conditions with radius ratio, χ = 0.35. Results are shown for each heating configuration: (a) differential heating (hollow …
Figure 5
Figure 5. Figure 5: The relative difference between the critical Rayleigh numbers of purely no-slip and purely stress-free boundary conditions. The quantity 1 − RaSF c /RaNS c is plotted against the Prandtl number, Pr, for different values of the Ekman number, E, with χ = 0.35 and the cas…
Figure 6
Figure 6. Figure 6: The critical radius, sc, as a function of the Prandtl number, Pr, for two values of the Ekman number, E, with χ = 0.35, purely no-slip boundary conditions, and for the case of internal heating. 3.2. Dependence on Radius Ratio All results presented in this subsection ar…
Figure 7
Figure 7. Figure 7: c shows that the critical wavenumber increases with increasing χ. This is partly due to the fact that, for differential heating, convection onsets in the form of columns at the tangent cylinder to the inner core, and as this inner core gets larger, the convection colum…
Figure 8
Figure 8. Figure 8: The radial component of the velocity field at onset plotted in an equatorial slice for purely no-slip boundary conditions, Pr = 1, differential heating (top half), internal heating (bottom half), E = 10−4 (left half of each panel), E = 10−5 (right half of each panel), …
Figure 9
Figure 9. Figure 9: Comparison of critical parameter values as functions of radius ratio, χ, for several values of the Ekman number, E, and different boundary conditions with Prandtl number, Pr = 1. Results are shown for each heating configuration: (a) differential heating (hollow symbols…
Figure 10
Figure 10. Figure 10: The relative difference between the critical Rayleigh numbers of purely no-slip and purely stress-free boundary conditions. The quantity 1 − RaSF c /RaNS c is plotted against radius ratio χ for different values of the Ekman number E with Prandtl number Pr = 1 and the …
Figure 11
Figure 11. Figure 11: The critical radius, sc, as a function of χ for two values of the Ekman number, E, with Pr = 1, purely no-slip boundary conditions, and for the case of internal heating (+) and differential heating (◇). The grey dashed line is given by sc/ro = χ, indicating the positi…
Figure 12
Figure 12. Figure 12: Dependence of the critical Rayleigh number on radius ratio for various definitions of Rac based on nondimensional set-ups discussed in Section 2.3. All results shown are for E = 10−5 and Pr = 1 with purely no-slip boundary conditions. Results are displayed for AN-1 (a…

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

48 extracted references · 36 canonical work pages

  1. [1]

    Planetary Dynamos

    Busse, F.H.; Simitev, R.D. Planetary Dynamos. In Treatise on Geophysics; Elsevier: Amsterdam, The Nether- lands, 2015; pp. 239–254. https://doi.org/10.1016/b978-0-444-53802-4.00172-x

    work page doi:10.1016/b978-0-444-53802-4.00172-x 2015

  2. [2]

    Rapidly Rotating Magnetohydrodynamics and the Geodynamo

    Dormy, E. Rapidly Rotating Magnetohydrodynamics and the Geodynamo. Annu. Rev. Fluid Mech. 2025, 57, 335–362. https://doi.org/10.1146/annurev-fluid-031224-121649

    work page doi:10.1146/annurev-fluid-031224-121649 2025

  3. [3]

    Solar Dynamo Theory

    Charbonneau, P . Solar Dynamo Theory. Annu. Rev. Astron. Astrophys. 2014, 52, 251–290. https: //doi.org/10.1146/annurev-astro-081913-040012

    doi:10.1146/annurev-astro-081913-040012 2014

  4. [4]

    Planetary Magnetic Fields and Fluid Dynamos

    Jones, C.A. Planetary Magnetic Fields and Fluid Dynamos. Annu. Rev. Fluid Mech. 2011, 43, 583–614. https://doi.org/10.1146/annurev-fluid-122109-160727

    doi:10.1146/annurev-fluid-122109-160727 2011

  5. [5]

    Dynamos in the Inner Solar System

    Tikoo, S.M.; Evans, A.J. Dynamos in the Inner Solar System. Annu. Rev. Earth Planet. Sci. 2022, 50, 99–122. https://doi.org/10.1146/annurev-earth-032320-102418

    work page doi:10.1146/annurev-earth-032320-102418 2022

  6. [6]

    Finite-amplitude convection in rotating spherical fluid shells

    Tilgner, A.; Busse, F.H. Finite-amplitude convection in rotating spherical fluid shells. J. Fluid Mech. 1997, 332, 359–376. https://doi.org/10.1017/S0022112096004259

    work page doi:10.1017/s0022112096004259 1997

  7. [7]

    Patterns of convection in rotating spherical shells

    Simitev, R.D.; Busse, F.H. Patterns of convection in rotating spherical shells. New J. Phys. 2003, 5, 97. https://doi.org/10.1088/1367-2630/5/1/397

    doi:10.1088/1367-2630/5/1/397 2003

  8. [8]

    Zonal flow driven by strongly supercritical convection in rotating spherical shells

    Christensen, U.R. Zonal flow driven by strongly supercritical convection in rotating spherical shells. J. Fluid Mech. 2002, 470, 115–133. https://doi.org/10.1017/S0022112002002008

    doi:10.1017/s0022112002002008 2002

Show all 48 references
  1. [9]

    Turbulent convection in liquid metal with and without rotation

    King, E.M.; Aurnou, J.M. Turbulent convection in liquid metal with and without rotation. Proc. Natl. Acad. Sci. USA 2013, 110, 6688–6693. https://doi.org/10.1073/pnas.1217553110

    2013 doi

  2. [10]

    Scaling behaviour of rotating convection in a spherical shell with different Prandtl numbers

    Fan, W.; Wang, Q.; Lin, Y. Scaling behaviour of rotating convection in a spherical shell with different Prandtl numbers. J. Fluid Mech. 2024, 998, A20. https://doi.org/10.1017/jfm.2024.893

    2024 doi

  3. [11]

    Dynamos driven by convection in rotating spherical shells

    Busse, F.H.; Simitev, R.D. Dynamos driven by convection in rotating spherical shells. Astron. Nachrichten Astron. Notes 2005, 326, 231–240. https://doi.org/10.1002/asna.200410382

    2005 doi

  4. [12]

    Scaling properties of convection-driven dynamos in rotating spherical shells and application to planetary magnetic fields

    Christensen, U.R.; Aubert, J. Scaling properties of convection-driven dynamos in rotating spherical shells and application to planetary magnetic fields. Geophys. J. Int. 2006, 166, 97–114. https://doi.org/10.1111/j.13 65-246X.2006.03009.x

    2006 arXiv

  5. [13]

    Double-diffusive convection in a rotating cylindrical annulus with conical caps

    Simitev, R.D. Double-diffusive convection in a rotating cylindrical annulus with conical caps. Phys. Earth Planet. Inter. 2011, 186, 183–190. https://doi.org/10.1016/j.pepi.2011.04.007

    2011 doi

  6. [14]

    On the thermal instability of a rotating-fluid sphere containing heat sources.Philos

    Roberts, P .H. On the thermal instability of a rotating-fluid sphere containing heat sources.Philos. Trans. R. Soc. London. Ser. A Math. Phys. Sci. 1968, 263, 93–117. https://doi.org/10.1098/rsta.1968.0007

    1968

  7. [15]

    Thermal instabilities in rapidly rotating systems

    Busse, F.H. Thermal instabilities in rapidly rotating systems. J. Fluid Mech. 1970, 44, 441–460. https: //doi.org/10.1017/S0022112070001921

    1970 doi

  8. [16]

    Hydrodynamic and Hydromagnetic Stability; Clarendon: Oxford, UK, 1961

    Chandrasekhar, S. Hydrodynamic and Hydromagnetic Stability; Clarendon: Oxford, UK, 1961

    1961

  9. [17]

    On the Finite amplitude thermal instability of a rapidly rotating fluid sphere

    Soward, A.M. On the Finite amplitude thermal instability of a rapidly rotating fluid sphere. Geophys. Astrophys. Fluid Dyn. 1977, 9, 19–74. https://doi.org/10.1080/03091927708242315

    1977 doi

  10. [18]

    The onset of thermal convection in a rapidly rotating sphere

    Jones, C.A.; Soward, A.M.; Mussa, A.I. The onset of thermal convection in a rapidly rotating sphere. J. Fluid Mech. 2000, 405, 157–179. https://doi.org/10.1017/S0022112099007235. 27 of 28

    2000 doi

  11. [19]

    The onset of thermal convection in rotating spherical shells

    Dormy, E.; Soward, A.M.; Jones, C.A.; Jault, D.; Cardin, P . The onset of thermal convection in rotating spherical shells. J. Fluid Mech. 2004, 501, 43–70. https://doi.org/10.1017/S0022112003007316

    2004 doi

  12. [20]

    On coupling between the Poincaré equation and the heat equation.J

    Zhang, K. On coupling between the Poincaré equation and the heat equation.J. Fluid Mech. 1994, 268, 211–229. https://doi.org/10.1017/s0022112094001321

    1994 doi

  13. [21]

    Thermal convection in rotating spherical shells

    Ardes, M.; Busse, F.; Wicht, J. Thermal convection in rotating spherical shells. Phys. Earth Planet. Inter. 1997, 99, 55–67. https://doi.org/10.1016/s0031-9201(96)03200-1

    1997 doi

  14. [22]

    Inertial convection in rotating fluid spheres

    Busse, F.H.; Simitev, R.D. Inertial convection in rotating fluid spheres. J. Fluid Mech. 2004, 498, 23–30. https://doi.org/10.1017/S0022112003006943

    2004 doi

  15. [23]

    On inertial waves in a rotating fluid sphere.J

    Zhang, K.; Earnshaw, P .; Liao, X.; Busse, F.H. On inertial waves in a rotating fluid sphere.J. Fluid Mech. 2001, 437, 103–119. https://doi.org/10.1017/s0022112001004049

    2001 doi

  16. [24]

    A new asymptotic method for the analysis of convection in a rapidly rotating sphere

    Zhang, K.; Liao, X. A new asymptotic method for the analysis of convection in a rapidly rotating sphere. J. Fluid Mech. 2004, 518, 319–346. https://doi.org/10.1017/s0022112004001120

    2004 doi

  17. [25]

    Asymptotic solutions of convection in rapidly rotating non-slip spheres

    Zhang, K.; Liao, X.; Busse, F.H. Asymptotic solutions of convection in rapidly rotating non-slip spheres. J. Fluid Mech. 2007, 578, 371–380. https://doi.org/10.1017/s002211200700523x

    2007 doi

  18. [26]

    Equations governing convection in earth’s core and the geodynamo.Geophys

    Braginsky, S.I.; Roberts, P .H. Equations governing convection in earth’s core and the geodynamo.Geophys. Astrophys. Fluid Dyn. 1995, 79, 1–97. https://doi.org/10.1080/03091929508228992

    1995 doi

  19. [27]

    Puzzles in Planetary Dynamos: Implications for Planetary Interiors

    Soderlund, K.M.; Stanley, S.; Cao, H.; Calkins, M.A.; Browning, M.K. Puzzles in Planetary Dynamos: Implications for Planetary Interiors. Annu. Rev. Earth Planet. Sci. 2025, 53, 305–37. https://doi.org/10.1146/ annurev-earth-111523-081635

    2025

  20. [28]

    Turbulence, Magnetism, and Shear in Stellar Interiors

    Miesch, M.S.; Toomre, J. Turbulence, Magnetism, and Shear in Stellar Interiors. Annu. Rev. Fluid Mech. 2009, 41, 317–345. https://doi.org/10.1146/annurev.fluid.010908.165215

    2009

  21. [29]

    Onset of low Prandtl number thermal convection in thin spherical shells

    Garcia, F.; Chambers, F.R.N.; Watts, A.L. Onset of low Prandtl number thermal convection in thin spherical shells. Phys. Rev. Fluids 2018, 3. https://doi.org/10.1103/physrevfluids.3.024801

    2018 doi

  22. [30]

    The onset of thermo-compositional convection in rotating spherical shells

    Silva, L.; Mather, J.F.; Simitev, R.D. The onset of thermo-compositional convection in rotating spherical shells. Geophys. Astrophys. Fluid Dyn. 2019, 113, 377–404. https://doi.org/10.1080/03091929.2019.1640875

    2019

  23. [31]

    Onset of convection in rotating spherical shells: Variations with radius ratio

    Barik, A.; Triana, S.; Calkins, M.; Stanley, S.; Aurnou, J. Onset of convection in rotating spherical shells: Variations with radius ratio. Earth Space Sci. 2023, 10, e2022EA002606. https://doi.org/10.1029/2022EA002 606

    2023 doi

  24. [32]

    Onset of thermal convection in a solid spherical shell with melting at either or both boundaries

    Morison, A.; Labrosse, S.; Deguen, R.; Alboussière, T. Onset of thermal convection in a solid spherical shell with melting at either or both boundaries. Geophys. J. Int. 2024, 238, 1121–1136. https://doi.org/10.1093/ gji/ggae208

    2024

  25. [33]

    Parameter dependences of convection-driven dynamos in rotating spherical fluid shells

    Busse, F.H.; Simitev, R.D. Parameter dependences of convection-driven dynamos in rotating spherical fluid shells. Geophys. Astrophys. Fluid Dyn. 2006, 100, 341–361. https://doi.org/10.1080/03091920600784873

    2006 doi

  26. [34]

    Varying the spherical shell geometry in rotating thermal convection

    Al-Shamali, F.; Heimpel, M.; Aurnou, J. Varying the spherical shell geometry in rotating thermal convection. Geophys. Astrophys. Fluid Dyn. 2004, 98, 153–169. https://doi.org/10.1080/03091920410001659281

    2004 doi

  27. [35]

    Effects of Shell Thickness on Cross-Helicity Generation in Convection-Driven Spherical Dynamos

    Silva, L.; Gupta, P .; MacTaggart, D.; Simitev, R.D. Effects of Shell Thickness on Cross-Helicity Generation in Convection-Driven Spherical Dynamos. Fluids 2020, 5, 245. https://doi.org/10.3390/fluids5040245

    2020 doi

  28. [36]

    A buoyancy profile for the Earth’ s core

    Davies, C.; Gubbins, D. A buoyancy profile for the Earth’ s core. Geophys. J. Int. 2011, 187, 549–563. https://doi.org/10.1111/j.1365-246X.2011.05144.x

    2011 arXiv

  29. [37]

    Linear theory of compressible convection in rapidly rotating spherical shells, using the anelastic approximation

    Jones, C.A.; Kuzanyan, K.M.; Mitchell, R.H. Linear theory of compressible convection in rapidly rotating spherical shells, using the anelastic approximation. J. Fluid Mech. 2009, 634, 291–319. https://doi.org/10.101 7/S0022112009007253

    2009

  30. [38]

    Convection in rotating spherical fluid shells with inhomogeneous heat flux at the outer boundary

    Gibbons, S.; Gubbins, D.; Zhang, K. Convection in rotating spherical fluid shells with inhomogeneous heat flux at the outer boundary. Geophys. Astrophys. Fluid Dyn. 2007, 101, 347–370. https://doi.org/10.1080/0309 1920701472550

    2007 doi

  31. [39]

    Solenoidal force balances in numerical dynamos

    Teed, R.J.; Dormy, E. Solenoidal force balances in numerical dynamos. J. Fluid Mech. 2023, 964, A26. https://doi.org/10.1017/jfm.2023.332

    2023 doi

  32. [40]

    On the onset of low-Prandtl-number convection in rotating spherical shells: Non-slip boundary conditions

    Net, M.; Garcia, F.; Sánchez, J. On the onset of low-Prandtl-number convection in rotating spherical shells: Non-slip boundary conditions. J. Fluid Mech. 2008, 601, 317–337. https://doi.org/10.1017/S0022112008000 61X

    2008 doi

  33. [41]

    Convection in a rapidly rotating spherical shell at infinite Prandtl number: Steadily drifting rolls

    Zhang, K. Convection in a rapidly rotating spherical shell at infinite Prandtl number: Steadily drifting rolls. Phys. Earth Planet. Inter. 1991, 68, 156–169. https://doi.org/10.1016/0031-9201(91)90015-A

    1991 doi

  34. [42]

    The influence of Ekman boundary layers on rotating convection

    Zhang, K.; Jones, C.A. The influence of Ekman boundary layers on rotating convection. Geophys. Astrophys. Fluid Dyn. 1993, 71, 145–162. https://doi.org/10.1080/03091929308203600. 28 of 28

    1993 doi

  35. [43]

    The ekman layer

    Cushman-Roisin, B.; Beckers, J.M. The ekman layer. In International Geophysics; Elsevier: Amsterdam, The Netherlands, 2011; Volume 101, pp. 239–270. https://doi.org/10.1016/B978-0-12-088759-0.00008-0

    2011 doi

  36. [44]

    Convection in rapidly rotating spherical fluid shells.Geophys

    Busse, F.H.; Cuong, P . Convection in rapidly rotating spherical fluid shells.Geophys. Astrophys. Fluid Dyn. 1977, 8, 17–41. https://doi.org/10.1080/03091927708240369

    1977 doi

  37. [45]

    Convection in rotating annulus: An asymptotic theory and numerical solutions

    Zhang, K.; Greed, G. Convection in rotating annulus: An asymptotic theory and numerical solutions. Phys. Fluids 1998, 10, 2396–2404. https://doi.org/10.1063/1.869756

    1998 doi

  38. [46]

    Influence of Boundary Conditions and Heating Modes on Convective Onset in Rotating Spherical Shells, [Data Collection], 2025

    Seeley, W.; Coke, F.; Simitev, R.D.; Teed, R.J. Influence of Boundary Conditions and Heating Modes on Convective Onset in Rotating Spherical Shells, [Data Collection], 2025. https://doi.org/10.5525/gla. researchdata.2015

    2025 doi

  39. [47]

    InPotential Theory in Gravity and Magnetic Applications; Cambridge University Press: Cambridge, UK, 1995; pp

    Blakely, R.J., Spherical Harmonic Analysis. InPotential Theory in Gravity and Magnetic Applications; Cambridge University Press: Cambridge, UK, 1995; pp. 100–127. https://doi.org/10.1017/CBO9780511549816.007

    1995 doi

  40. [48]

    Patterns of convection in spherical shells

    Busse, F.H.; Riahi, N. Patterns of convection in spherical shells. Part 2. J. Fluid Mech. 1982, 123, 283–301. https://doi.org/10.1017/S0022112082003061. Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individ...

    1982 doi

Pith tools

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