Thesis (Ph.D.) - University of St Andrews, January 2007.

The solar corona is a highly complex and active plasma environment, containing many exotic phenomena such as solar flares, coronal mass ejections, prominences, coronal loops, and bright points. The fundamental element giving coherence to all this apparent diversity is the strong coronal magnetic field, the dominant force shaping the plasma there. In this thesis, I model the 3D magnetic fields of various coronal features using the techniques of magnetic charge topology (MCT) in a potential field. Often the real coronal field has departures from its potential state, but these are so small that the potential field method is accurate enough to pick out the essential information about the structure and evolution of the magnetic field. First I perform a topological analysis of the magnetic breakout model for an eruptive solar flare. Breakout is represented by a topological bifurcation that allows initially enclosed flux from the newly emerging region in my MCT model of a delta sunspot to reconnect out to large distances. I produce bifurcation diagrams showing how this behaviour can be caused by changing the strength or position of the emerging flux source, or the force-free parameter α. I also apply MCT techniques to observational data of a coronal bright point, and compare the results to 3D numerical MHD simulations of the effects of rotating the sources that underlie the bright point. The separatrix surfaces that surround each rotating source are found to correspond to locations of high parallel electric field in the simulations, which is a signature of magnetic reconnection. The large-scale topological structure of the magnetic field is robust to changes in the method of deriving point magnetic sources from the magnetogram. Next, I use a Green’s function expression for the magnetic field to relax the standard topological assumption of a flat photosphere and extend the concept of MCT into a spherical geometry, enabling it to be applied to the entire global coronal magnetic field. I perform a comprehensive study of quadrupolar topologies in this new geometry, producing several detailed bifurcation diagrams. These results are compared to the equivalent study for a flat photosphere. A new topological state is found on the sphere which has no flat photosphere analogue; it is named the dual intersecting state because of its twin separators joining a pair of magnetic null points. The new spherical techniques are then applied to develop a simple six-source topological model of global magnetic field reversal during the solar cycle. The evolution of the large-scale global magnetic field is modelled through one complete eleven-year cycle, beginning at solar minimum. Several distinct topological stages are exhibited: active region flux connecting across the equator to produce transequatorial loops; the dominance of first the leading and then the following polarities of the active regions; the magnetic isolation of the poles; the reversal of the polar field; the new polar field connecting back to the active regions; the polar flux regaining its dominance; and the disappearance of the transequatorial loops.

The problem of how the Sun's corona is heated is of central importance to solar physics research. In this thesis we model three main areas. The first, annihilation, is a feature of non-ideal MHD and focusses on how magnetic field of opposite polarity meets at a null point and annihilates, after having been advected with plasma toward a stagnation point in the plasma flow. Generally, the null point of the field and the stagnation point of the flow are coincident at the origin, but in chapter 2 a simple extension is considered where an asymmetry in the boundary conditions of the field moves the null point away from the origin. Chapter 3 presents a model of reconnective annihilation in three dimensions. It represents flux being advected through the fan plane of a 3D null, and diffusing through a thin diffusion region before being annihilated at the spine line, and uses the method of matched asymptotic expansions to find the solution for small values of the resistivity. The second area of the thesis covers null collapse. This is when the magnetic field in close proximity to a null point is disturbed, causing the field to fold up on itself and collapse. This is a feature of ideal MHD, and causes a strong current to build up, allowing non-ideal effects to become important. When using linearised equations for the collapse problem, we are in fact looking at a linear instability. If this instability initiates a collapse, this is only a valid model until non-linear effects become important. By talking about collapse in chapters 4 and 5 (as it is talked about in the literature), we mean that the linear instability initiates collapse, which in principle, non-linear effects could later stop. Chapter 4 introduces a two-dimensional model for collapse, using the ideal, compressible, linearised MHD equations. It is a general solution in which all spatially linear nulls and their supporting plasma flows and pressure gradients can be checked for susceptibility to collapse under open boundary conditions. Chapter 5 uses the model introduced in chapter 4 to investigate the collapse of three-dimensional, potential nulls (again, spatially linear) for all possible supporting plasma flows and pressure gradients. Using this model, all nulls under consideration are found to collapse and produce large currents, except for a group of 2D O-type nulls supported by highly super-Alfvenic plasma flows. The third area of this thesis involves numerically simulating a model of heating by coronal tectonics (Priest et al, 2002). A simple magnetic field is created and the boundary is driven, also in a simple manner. Current sheets which scale with grid resolution are seen to build up on the quasi-separatrix layers, and there is some evidence of magnetic reconnection.

The dynamics of the solar corona are dominated by the magnetic field which creates its structure. The magnetic field in most of the corona is ‘frozen’ to the plasma very effectively. The exception is in small localised regions of intense current concentrations where the magnetic field can slip through the plasma and a restructuring of the magnetic field can occur. This process is known as magnetic reconnection and is believed to be responsible for a wide variety of phenomena in the corona, from the rapid energy release of solar flares to the heating of the high-temperature corona. The coronal field itself is three-dimensional (3D), but much of our understanding of reconnection has been developed through two-dimensional (2D) models. This thesis describes several models for fully 3D reconnection, with both kinematic and fully dynamic models presented. The reconnective behaviour is shown to be fundamentally different in many respects from the 2D case. In addition a numerical experiment is described which examines the reconnection process in coronal magnetic flux tubes whose photospheric footpoints are spun, one type of motion observed to occur on the Sun. The large-scale coronal field itself is thought to be generated by a magnetohydrodynamic dynamo operating in the solar interior. Although the dynamo effect itself is not usually associated with reconnection, since the essential element of the problem is to account for the presence of large-scale fields, reconnection is essential for the restructuring of the amplified small-scale flux. Here we examine some simple models of the solar-dynamo process, taking advantage of their simplicity to make a full exploration of their behaviour in a variety of parameter regimes. A wide variety of dynamic behaviour is found in each of the models, including aperiodic modulation of cyclic solutions and intermittency that strongly resembles the historic record of solar magnetic activity.