IDEAS home Printed from https://ideas.repec.org/a/gam/jmathe/v13y2025i13p2113-d1689528.html
   My bibliography  Save this article

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

Author

Listed:
  • 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)

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.

Suggested Citation

  • Juan Cruz Gonzalez Sembla & Camille Rambert & Fred Feudel & Laurette S. Tuckerman, 2025. "Convection in a Rapidly Rotating Spherical Shell: Newton’s Method Using Implicit Coriolis Integration," Mathematics, MDPI, vol. 13(13), pages 1-26, June.
  • Handle: RePEc:gam:jmathe:v:13:y:2025:i:13:p:2113-:d:1689528
    as

    Download full text from publisher

    File URL: https://www.mdpi.com/2227-7390/13/13/2113/pdf
    Download Restriction: no

    File URL: https://www.mdpi.com/2227-7390/13/13/2113/
    Download Restriction: no
    ---><---

    More about this item

    Keywords

    ;
    ;
    ;
    ;
    ;

    Statistics

    Access and download statistics

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:gam:jmathe:v:13:y:2025:i:13:p:2113-:d:1689528. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    We have no bibliographic references for this item. You can help adding them by using this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: MDPI Indexing Manager (email available below). General contact details of provider: https://www.mdpi.com .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.