Nonlinear evolution of unstable solar inertial modes: The case of viscous modes on a differentially rotating sphere
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2026
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866908873938436096 |
|---|---|
| author | Mushtaq, Muneeb Fournier, Damien Ayoub, Rama Schmid, Peter J. Gizon, Laurent |
| author_facet | Mushtaq, Muneeb Fournier, Damien Ayoub, Rama Schmid, Peter J. Gizon, Laurent |
| contents | On the Sun, the inertial mode with the largest observed amplitude (rms velocity exceeding $10$ m/s) is the high-latitude mode with longitudinal wavenumber $m=1$. In two dimensions, on the sphere, linear theory predicts that this mode is unstable due to a shear instability associated with latitudinal differential rotation (fast equator, slower polar regions). We investigate the evolution of this instability numerically and theoretically.
The nonlinear vorticity equation is solved using direct numerical simulations in the time domain. The only control parameter is the Ekman number $E$. For $10^{-3}\lesssim E< E_c \approx 1.5\times10^{-3}$, only the high-latitude $m=1$ mode is unstable. We extract its saturation amplitude as a function of $E$ and compare the results with predictions from two perturbative approaches in nonlinear stability theory.
The simulations reveal a supercritical Hopf bifurcation. Near onset, the mode amplitude is well described by the Landau equation $d|A|/dt=σ_I |A|+β_I |A|^3$, with a positive linear growth rate $σ_I$ and a negative nonlinear coefficient $β_I$. The coefficient $β_I$ depends weakly on $E$, implying that the saturated amplitude scales approximately as $|A|\proptoσ_I^{1/2}$. The equilibrium mode contains the $m=1$ fundamental and harmonics $m=2$ and $m=3$, whose amplitudes scale as $σ_I^{m/2}$. Saturation results from Reynolds stresses that smooth the latitudinal differential rotation. For $E=4\times10^{-4}$, consistent with solar-like turbulent viscosity, the saturated velocity reaches $28$ m/s, comparable to solar observations.
These results should be interpreted cautiously, since in three dimensions the instability is baroclinic and involves different physics. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2603_08528 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Nonlinear evolution of unstable solar inertial modes: The case of viscous modes on a differentially rotating sphere Mushtaq, Muneeb Fournier, Damien Ayoub, Rama Schmid, Peter J. Gizon, Laurent Solar and Stellar Astrophysics Fluid Dynamics On the Sun, the inertial mode with the largest observed amplitude (rms velocity exceeding $10$ m/s) is the high-latitude mode with longitudinal wavenumber $m=1$. In two dimensions, on the sphere, linear theory predicts that this mode is unstable due to a shear instability associated with latitudinal differential rotation (fast equator, slower polar regions). We investigate the evolution of this instability numerically and theoretically. The nonlinear vorticity equation is solved using direct numerical simulations in the time domain. The only control parameter is the Ekman number $E$. For $10^{-3}\lesssim E< E_c \approx 1.5\times10^{-3}$, only the high-latitude $m=1$ mode is unstable. We extract its saturation amplitude as a function of $E$ and compare the results with predictions from two perturbative approaches in nonlinear stability theory. The simulations reveal a supercritical Hopf bifurcation. Near onset, the mode amplitude is well described by the Landau equation $d|A|/dt=σ_I |A|+β_I |A|^3$, with a positive linear growth rate $σ_I$ and a negative nonlinear coefficient $β_I$. The coefficient $β_I$ depends weakly on $E$, implying that the saturated amplitude scales approximately as $|A|\proptoσ_I^{1/2}$. The equilibrium mode contains the $m=1$ fundamental and harmonics $m=2$ and $m=3$, whose amplitudes scale as $σ_I^{m/2}$. Saturation results from Reynolds stresses that smooth the latitudinal differential rotation. For $E=4\times10^{-4}$, consistent with solar-like turbulent viscosity, the saturated velocity reaches $28$ m/s, comparable to solar observations. These results should be interpreted cautiously, since in three dimensions the instability is baroclinic and involves different physics. |
| title | Nonlinear evolution of unstable solar inertial modes: The case of viscous modes on a differentially rotating sphere |
| topic | Solar and Stellar Astrophysics Fluid Dynamics |
| url | https://arxiv.org/abs/2603.08528 |