Wednesday 2 November 2011

The next step: Modelling a realistic input based on observations

For realistic solar cycle simulations, corresponding inputs have to be modeled and fed to the Surface flux transport code.

The modelling was done with reference to the study carried out by Jie Jiang et al. (http://arxiv.org/abs/1102.1266v1) and van Ballegooijen et al. 1998.

The input parameters needed are:
(1) Time of BMR eruption
(2) Latitude of eruption
(3) Tilt angle (with respect to the solar equator)
(4) Longitude of eruption
(5) Radius of individual spots in a BMR
(6) Peak magnetic field (Bmax)
(7) Separation between the centers of the individual spots.


One solar cycle of 11 years was divided in 120 phases. And the sunspots were placed at their respective locations after each phase. The following figure shows the latitude of eruption(in degrees) vs phase that was fed into the surface flux transport code:







An initial polar field of +-4.5 Gauss was placed within 23 degrees of the poles, and the surface flux transport code was given a run for 12 sunspot cycles(132 years). This is what was the output:





This is again a butterfly diagram, but you can really see butterfly like structures in it. The polar magnetic field reversal is clearly evident from the diagram.

But, the problem here is, the magnitude of the peak polar field is not constant. It varies from cycle to cycle.   

No. of Cycle     Peak Magnetic field in the northern hemisphere
        1                                             -4.5
        2                                            +3.8
        3                                            -2.8
        4                                            +5.0
        5                                            -3.8
        6                                            +3.4
        7                                            -3.9
        8                                            +2.7
        9                                            -3.9
       10                                           +3.9
       11                                           -4.9
       12                                           +2.3
       13                                           -5.0

Work up till now: Polar magnetic field reversal

An initial polar field of +-10 Gauss was placed within 20 degrees of the poles, and a BMR was placed in each hemisphere, such that the trailing spot has magnetic field of a polarity opposite to the polar magnetic field in that hemisphere. This automatically set the leading spots of the BMR to be opposite in polarity of magnetic flux.

Such BMRs were placed at regular intervals till the polar magnetic field was totally cancelled, and then became opposite.

The following figure is a butterfly diagram of the simulation. There is nothing like a butterfly in it, it is named like that for a different reason.

On the X-axis is the no. of time steps,
on the Y-axis is the latitude of the surface of the Sun,
and the colors indicate the value of magnetic field on a particular latitude at a particular time(averaged over all longitudes).


Work up till now: Building up of poloidal field

In this video, you can see,
When 2 BMRs in different hemispheres are such that, their leading spots have magnetic flux of the opposite sign, then their flux gets cancelled over the equator. While, the flux of the trailing sunspots is carried to the poles by the meridional flow.


But when the tilt of the BMRs is such that the leading spots have magnetic flux of the same sign, there is no cancellation over the equator, and both the poles end up getting the same polarity of magnetic field.

Tuesday 1 November 2011

Work up till now

I have written the Surface flux transport code in FORTRAN77.






It seems to be working fine when I load an input with a single Bipolar Magnetic Region(BMR) placed near the equator. The diffusion seems to blow up the sunspots in size, and the meridional flow(11m/s) carries the poloidal component of the flux to the poles.

The above figures are contour plots of the radial magnetic field on the surface of the Sun. The X-axis is the azimuthal angle(phi) and the Y-axis is the polar angle(theta).

Purpose

I was introduced to research in astrophysics in the summer of 2011 when I went to IISER Kolkata for a summer project in Solar Physics, under the guidance of Dr. Dibyendu Nandi. He is one of the nicest person I have known.

My project was to develop a Surface Magnetic Flux Transport code that would simulate the evolution of radial magnetic field on the surface of the Sun.


Dr. Nandi, has himself worked on Solar Dynamo code. On combining the Dynamo code with the surface flux transport code, we will have a better understanding of the star nearest to us. And in the future, we can hope to extrapolate these fields to and simulate the interplanetary magnetic field in the Solar System.

The purpose of this blog is to post the recent advances in the development of the Surface flux transport code.