Convection in a Rapidly Rotating Spherical Shell: Newton's Method Using Implicit Coriolis Integration
Fuente:
arXiv
Salvato in:
| Autori principali: | , , , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2025
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866911117218938880 |
|---|---|
| author | Sembla, Juan Cruz Gonzalez Rambert, Camille Feudel, Fred Tuckerman, Laurette S. |
| author_facet | Sembla, Juan Cruz Gonzalez Rambert, Camille Feudel, Fred Tuckerman, Laurette S. |
| contents | 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_16635 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Convection in a Rapidly Rotating Spherical Shell: Newton's Method Using Implicit Coriolis Integration Sembla, Juan Cruz Gonzalez Rambert, Camille Feudel, Fred Tuckerman, Laurette S. Geophysics Fluid Dynamics 35B32, 37M20, 65F08, 65M70, 65N35, 65P30, 76E06, 76U60 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. |
| title | Convection in a Rapidly Rotating Spherical Shell: Newton's Method Using Implicit Coriolis Integration |
| topic | Geophysics Fluid Dynamics 35B32, 37M20, 65F08, 65M70, 65N35, 65P30, 76E06, 76U60 |
| url | https://arxiv.org/abs/2508.16635 |