EconPapers    
Economics at your fingertips  
 

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

Juan Cruz Gonzalez Sembla, Camille Rambert, Fred Feudel and Laurette S. Tuckerman ()
Additional contact information
Juan Cruz Gonzalez Sembla: PMMH, CNRS, ESPCI Paris, Université PSL, Sorbonne Université, Université de Paris, 75005 Paris, France
Camille Rambert: PMMH, CNRS, ESPCI Paris, Université PSL, Sorbonne Université, Université de Paris, 75005 Paris, France
Fred Feudel: Institut für Physik und Astronomie, Universität Potsdam, 14476 Potsdam, Germany
Laurette S. Tuckerman: PMMH, CNRS, ESPCI Paris, Université PSL, Sorbonne Université, Université de Paris, 75005 Paris, France

Mathematics, 2025, vol. 13, issue 13, 1-26

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 Ek , 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.

Keywords: rotating fluids; Krylov methods; preconditioning; bifurcation; rotating waves (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2025
References: Add references at CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/13/13/2113/pdf (application/pdf)
https://www.mdpi.com/2227-7390/13/13/2113/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:13:y:2025:i:13:p:2113-:d:1689528

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().

 
Page updated 2025-06-28
Handle: RePEc:gam:jmathe:v:13:y:2025:i:13:p:2113-:d:1689528