Thesis (Ph. D.)--University of Hawaii at Manoa, 2004. / Includes bibliographical references (leaves 134-139). / Also available by subscription via World Wide Web / xvi, 139 leaves, bound ill. 29 cm

Thesis (Ph. D.)--University of Texas at Austin, 1998. / Vita. Includes bibliographical references (leaves 124-131). Available also in a digital version from Dissertation Abstracts.

Thesis (Ph. D.)--University of Hawaii at Manoa, 2004. / Includes bibliographical references (leaves 134-139).

Flux is continually emerging on the Sun, making its way from the solar interior up into the atmosphere. Emergence occurs on small-scales in the quiet Sun where magnetic fragments emerge, interact and cancel and on large-scales in active regions where magnetic fields emerge and concentrate to form sunspots. This thesis has been concerned with the large-scale emergence process and in particular the results from previous solar flux emergence modelling endeavours. Modelling uses numerical methods to evolve a domain representing simplified layers of the Sun’s atmosphere, within which the subsurface layer contains magnetic flux. The flux is initialised such that it will rises towards the surface at the start of the simulation. Once the flux reaches the solar surface, it can only emerge into the atmosphere if a magnetic buoyancy instability occurs, after which it expands rapidly both vertically and horizontally. The aim of this thesis is to test the robustness of these general findings from simulations to date upon the seed magnetic field. More explicitly, we have used three-dimensional numerical simulations to investigate how variations in the subsurface magnetic field modify the emergence process and the resulting atmospheric field. We initially consider a simple constant twist flux tube for the seed field and vary the tube’s magnetic field strength and degree of twist. Additionally, we have examined the effects of using non-constant twist flux tubes as the seed field by choosing two different profiles for the twist that are functions of the tube’s radius. Finally, we have investigated the effects of increasing the complexity of the seed field by positioning two flux tubes below the solar surface and testing two different configurations for the tubes. In both cases, the magnetic fields of the two tubes are such that, once the tubes come into contact with each other, reconnection occurs and a combined flux system is formed. From our investigations, we conclude that the general emergence results given by previous simulations are robust. However, for constant twist tubes with low field strength and twist, the buoyancy instability fails to be launched when the tubes reach the photosphere and they remain trapped in the low atmosphere. Similarly, when the non-constant twist profile results in a low tension force throughout the tube, we find that the buoyancy instability is not initialised.

The new CRISP filter at the Swedish 1-m Solar Telescope provides opportunities for observing the solar atmosphere with unprecedented spatial resolution and cadence. In order to benefit from the high quality of observational data from this instrument, we have developed methods for calibrating and restoring polarized Stokes images, obtained at optical and near infrared wavelengths, taking into account field-of-view variations of the filter properties. In order to facilitate velocity measurements, a time series from a 3D hydrodynamical granulation simulation is used to compute quiet Sun spectral line profiles at different heliocentric angles. The synthetic line profiles, with their convective blueshifts, can be used as absolute references for line-of-sight velocities. Observations of the Ca II 8542 Å line are used to study magnetic fields in chromospheric fibrils. The line wings show the granulation pattern at mid-photospheric heights whereas the overlying chromosphere is seen in the core of the line. Using full Stokes data, we have attempted to observationally verify the alignment of chromospheric fibrils with the magnetic field. Our results suggest that in most cases fibrils are aligned along the magnetic field direction, but we also find examples where this is not the case. Detailed interpretation of Stokes data from spectral lines formed in the chromospheric data can be made using non-LTE inversion codes. For the first time, we use a realistic 3D MHD chromospheric simulation of the quiet Sun to assess how well NLTE inversions recover physical quantities from spectropolarimetric observations of Ca II 8542 Å. We demonstrate that inversions provide realistic estimates of depth-averaged quantities in the chromosphere, although high spectral resolution and high sensitivity are needed to measure quiet Sun chromospheric magnetic fields. / At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 1: Submitted. Paper 2: Manuscript. Paper 3: Submitted. Paper 4: Submitted.

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

Magnetic flux continually emerges from the Sun, rising through the solar interior, emerging at the photosphere in the form of sunspots and expanding into the atmosphere. Observations of sunspot rotations have been reported for over a century and are often accompanied by solar eruptions and flaring activity. In this thesis, we present 3D numerical simulations of the emergence of twisted flux tubes from the uppermost layers of the solar interior, examining the rotational movements of sunspots in the photospheric plane. The basic experiment introduces the mechanism and characteristics of sunspot rotation by a clear calculation of rotation angle, vorticity, magnetic helicity and energy, whereby we find an untwisting of the interior portion of the field, accompanied by an injection of twist into the atmospheric field. We extend this model by altering the initial field strength and twist of the sub-photospheric tube. This comparison reveals the rotation angle, helicity and current show a direct dependence on field strength. An increase in field strength increases the rotation angle, the length of fieldlines extending into the atmosphere, and the magnetic energy transported to the atmosphere. The fieldline length is crucial as we predict the twist per unit length equilibrates to a lower value on longer fieldlines, and hence possesses a larger rotation angle. No such direct dependence is found when varying the twist but there is a clear ordering in rotation angle, helicity, and energy, with more highly twisted tubes undergoing larger rotation angles. We believe the final angle of rotation is reached when the system achieves a constant degree of twist along the length of fieldlines. By extrapolating the size of the modelled active region, we find rotation angles and rates comparable with those observed. In addition, we explore sunspot rotation caused by sub-photospheric velocities twisting the footpoints of flux tubes.

This thesis studies the acceleration of particles to super-thermal energies in explosive solar events as well as the magnetic changes in connectivity that may be responsible for changes in the morphology of quiescent filaments. Firstly a review of some of the observations of solar flare dynamics is given, as well as an introduction to the competing theories attempting to explain both particle acceleration and filament formation. An explanation of the numerical FORTRAN code that is used to calculate the trajectories of particle distribution functions in prescribed electromagnetic fields is given. Examples of known fields are used to test the accuracy of the code and the simple example of the well-known Litvinenko current sheet field is investigated. The results of charged particle orbit calculations in prescribed electric and magnetic fields motivated by magnetic reconnection models are then presented. The electromagnetic fields are chosen to resemble a current sheet with a localised reconnection region. The dependence of the model on the important physical parameters is considered. An introduction to the mathematical formulation of a collapsing magnetic trap is given. The same numerical code is used to calculate single electron orbits in this more complicated time dependent electromagnetic field. Consideration of important previous work is given before describing the best attempts to model the movement of flare loops in a realistic fashion. Finally the process of flux cancellation and filament formation is studied using a range of data including ground-based Hα and SoHO MDI magnetograms. It is found that the cancellation occurs at the ends of Hα sections of the filament and is accompanied by a noticeable increase in the Hα intensity and linkage of the sections. Measurements of the amount of flux cancelled at each site show it is in agreement with an estimate of the axial flux contained in the filament.

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.

In magnetised plasmas, magnetic reconnection is the process of magnetic field merging and recombination through which considerable amounts of magnetic energy may be converted into other forms of energy. Reconnection is a key mechanism for solar flares and coronal mass ejections in the solar atmosphere, it is believed to be an important source of heating of the solar corona, and it plays a major role in the acceleration of particles in the Earth's magnetotail. For reconnection to occur, the magnetic field must, in localised regions, be able to diffuse through the plasma. Ideal locations for diffusion to occur are electric current layers formed from rapidly changing magnetic fields in short space scales. In this thesis we consider the formation and nature of these current layers in magnetised plasmas. The study of current sheets and current layers in two, and more recently, three dimensions, has been a key field of research in the last decades. However, many of these studies do not take plasma pressure effects into consideration, and rather they consider models of current sheets where the magnetic forces sum to zero. More recently, others have started to consider models in which the plasma beta is non-zero, but they simply focus on the actual equilibrium state involving a current layer and do not consider how such an equilibrium may be achieved physically. In particular, they do not allow energy conversion between magnetic and internal energy of the plasma on their way to approaching the final equilibrium. In this thesis, we aim to describe the formation of equilibrium states involving current layers at both two and three dimensional magnetic null points, which are specific locations where the magnetic field vanishes. The different equilibria are obtained through the non-resistive dynamical evolution of perturbed hydromagnetic systems. The dynamic evolution relaxes via viscous damping, resulting in viscous heating. We have run a series of numerical experiments using LARE, a Lagrangian-remap code, that solves the full magnetohydrodynamic (MHD) equations with user controlled viscosity and resistivity. To allow strong current accumulations to be created in a static equilibrium, we set the resistivity to be zero and hence simply reach our equilibria by solving the ideal MHD equations. We first consider the relaxation of simple homogeneous straight magnetic fields embedded in a plasma, and determine the role of the coupling between magnetic and plasma forces, both analytically and numerically. Then, we study the formation of current accumulations at 2D magnetic X-points and at 3D magnetic nulls with spine-aligned and fan-aligned current. At both 2D X-points and 3D nulls with fan-aligned current, the current density becomes singular at the location of the null. It is impossible to be precisely achieve an exact singularity, and instead, we find a gradual continuous increase of the peak current over time, and small, highly localised forces acting to form the singularity. In the 2D case, we give a qualitative description of the field around the magnetic null using a singular function, which is found to vary within the different topological regions of the field. Also, the final equilibrium depends exponentially on the initial plasma pressure. In the 3D spine-aligned experiments, in contrast, the current density is mainly accumulated along and about the spine, but not at the null. In this case, we find that the plasma pressure does not play an important role in the final equilibrium. Our results show that current sheet formation (and presumably reconnection) around magnetic nulls is held back by non-zero plasma betas, although the value of the plasma pressure appears to be much less important for torsional reconnection. In future studies, we may consider a broader family of 3D nulls, comparing the results with the analytical calculations in 2D, and the relaxation of more complex scenarios such as 3D magnetic separators.