Convection in a Rapidly Rotating Spherical Shell: Newton's Method Using Implicit Coriolis Integration

Abstract

Geophysical flows are characterized by rapid rotation. Simulating these flows requires small timesteps to achieve stability and accuracy. Numerical stability can be greatly improved by the implicit integration of the terms that are most responsible for destabilizing the numerical scheme. We have implemented an implicit treatment of the Coriolis force in a rotating spherical shell driven by a radial thermal gradient. We modified the resulting timestepping code to carry out steady-state solving via Newton's method, which has no timestepping error. The implicit terms have the effect of preconditioning the linear systems, which can then be rapidly solved by a matrix-free Krylov method. We computed the branches of rotating waves with azimuthal wavenumbers ranging from 4 to 12. As the Ekman number (the non-dimensionalized inverse rotation rate) decreases, the flows are increasingly axially independent and localized near the inner cylinder, in keeping with well-known theoretical predictions and previous experimental and numerical results. The advantage of the implicit over the explicit treatment also increases dramatically with decreasing Ekman numbers, reducing the cost of computation by as much as a factor of 20 for Ekman numbers of order of $10^{-5}$. We carried out continuation for both the Rayleigh and Ekman numbers and obtained interesting branches in which the drift velocity remained unchanged between pairs of saddle-node bifurcations.

References in corpus (17)

• High Reynolds number Taylor-Couette turbulence
• The Significance of Simple Invariant Solutions in Turbulent Flows
• Turbulent geodynamo simulations: a leap towards Earth's core
• Efficient Spherical Harmonic Transforms aimed at pseudo-spectral numerical simulations
• Equilibrium and traveling-wave solutions of plane Couette flow
• Multiple states in highly turbulent Taylor-Couette flow
• Relative periodic orbits in transitional pipe flow
• Patterns of convection in rotating spherical shells
• Extreme multiplicity in cylindrical Rayleigh-Benard convection: II. Bifurcation diagram and symmetry classification
• Extreme multiplicity in cylindrical Rayleigh-Bénard convection: I. Time-dependence and oscillations
• Stability and Decay Rates of Non-Isotropic Attractive Bose-Einstein Condensates
• A blue sky catastrophe in double-diffusive convection
• Exact relations between Rayleigh-Bénard and rotating plane Couette flow in 2D
• Bridging the Rossby number gap in rapidly rotating thermal convection
• Natural convection in a vertical channel. Part 1. Wavenumber interaction and Eckhaus instability in a narrow domain
• Natural convection in a vertical channel. Part 2. Oblique solutions and global bifurcations in a spanwise-extended domain
• Seven decades of exploring planetary interiors with rotating convection experiments