Numerical simulations of rotating axisymmetric sunspots

Gert Botha, F. H. Busse, Neal Hurlburt, Alistair Rucklidge

Research output: Contribution to journalArticlepeer-review

10 Citations (Scopus)
16 Downloads (Pure)


A numerical model of axisymmetric convection in the presence of a vertical magnetic flux bundle and rotation about the axis is presented. The model contains a compressible plasma described by the non-linear MHD equations, with density and temperature gradients simulating the upper layer of the Sun’s convection zone. The solutions exhibit a central magnetic flux tube in a cylindrical numerical domain, with convection cells forming collar flows around the tube. When the numerical domain is rotated with a constant angular velocity, the plasma forms a Rankine vortex, with the plasma rotating as a rigid body where the magnetic field is strong, as in the flux tube, while experiencing sheared azimuthal flow in the surrounding convection cells, forming a free vortex. As a result, the azimuthal velocity component has its maximum value close to the outer edge of the flux tube. The azimuthal flow inside the magnetic flux tube and the vortex flow is prograde relative to the rotating cylindrical reference frame. A retrograde flow appears at the outer wall. The most significant convection cell outside the flux tube is the location for the maximum value of the azimuthal magnetic field component. The azimuthal flow and magnetic structure are not generated spontaneously, but decay exponentially in the absence of any imposed rotation of the cylindrical domain.
Original languageEnglish
Pages (from-to)1445-1462
JournalMonthly Notices of the Royal Astronomical Society
Issue number4
Publication statusPublished - 2008


Dive into the research topics of 'Numerical simulations of rotating axisymmetric sunspots'. Together they form a unique fingerprint.

Cite this