Call number: LD2668 .T4 1966 S562 / Master of Science

The magnetic field in the solar atmosphere is highly inhomogeneous. In the photosphere, the field is concentrated into intense flux tubes and the coronal magnetic field consists of many loops and regions of open field. This thesis investigates some of the basic properties of inhomogeneous solar magnetic fields. First of all, the equilibrium properties of untwisted flux tubes, confined by a spatially varying external pressure distribution, are investigated. The behaviour of thick flux tubes, including the effects of a transverse field component and a variation in the field across the tube, is compared with slender flux tube theory. It is shown that slender tube theory is accurate for tubes which are approximately slender, but that completely misleading results can be obtained by applying slender tube theory if the pressure distribution is not slowly varying. Twisted flux tubes are then studied, with the aim of finding how twisting affects a tube confined by an inhomogeneous pressure distribution. It is shown that, in general, a tube expands as it is twisted; this is illustrated both by extensions to slender tube theory and by some exact analytical solutions. A family of linear solutions is used to model the evolution of a finite tube confined by a falling external pressure. It is shown that, if the confining pressure falls too low, the tube may burst, with some dynamic process ensuing. The equilibrium properties of a flux tube with a curved axis are then investigated, with the main aim of modelling coronal loops. Previous theory for the equilibrium of a curved slender flux tube in a gravitationally stratified atmosphere, with a balance between magnetic buoyancy and tension forces, is extended to take into account an external field and the effects of twist. Increasing the magnitude of the external field tends to lower the summit height of the tube. It is found that non-equilibrium sets in if the footpoints are separated more than a certain critical width, which does not depend on the magnitude of the external field. It is found that two possible equilibrium heights can exist for a twisted tube; however, if the tube is twisted too far, or if the footpoints are moved apart, non-equilibrium can set in. The critical width at which non-equilibrium occurs is lower for a twisted tube than for an untwisted one. This is suggested as an explanation for the rise of a filament prior to a two ribbon flare, and as a mechanism for coronal transients. An alternative description of the coronal magnetic field is given, using a perturbation expansion for an almost potential field, with small pressure gradients. The field is assumed to be line-tied at the photospheric base. Then the equilibrium properties of the global magnetic field of a star are investigated. A linear and non-linear family of solutions to the magnetostatic equilibrium equation are found. The linear solutions are used to investigate the twisting up of force-free dipolar and quadrupolar fields, including in a simple manner the effects of a stellar wind. In both cases, it was found that the field becomes physically unreasonable if it is twisted too far, with field lines detached from the star being formed, which would be pulled out by the stellar wind. Thus, if the field is twisted more than a critical amount, non-equilibrium sets in and some catastrophic behaviour takes place. This is suggested as a possible mechanism for stellar flares. Similar results are found in a study of the effects of increasing the pressure gradients at the stellar surface of a magnetostatic dipole-like field. The linear solutions are also used to study the equilibrium of a finite magnetosphere, and multiple equilibria are found. Finally, one aspect of the propagation of waves in an inhomogeneous magnetic field is studied, with particular reference to the problem of heating the solar corona. The mechanism of phase-mixing, which provides a means of dissipating shear Alfven waves that propagate in an inhomogeneous magnetic field, is investigated. The onset of Kelvin-Helmholtz instability, which could disrupt the wave and thus enhance the dissipation, is studied. First, the dispersion relation of the instability is calculated for the case of fully developed phase-mixing. Then, the onset of the instability is investigated, to find out whether the instability can grow before the phase-mixing is fully developed. It is found that instability can set in after only a very few wave periods. It is suggested that the instability triggers off a turbulent cascade which dissipates the wave energy. The heating rates that could be produced by such a process are calculated, and are found to be more than adequate for coronal heating.

Magnetic fields play a key role in a wide variety of phenomena found on the Sun. One such phenomena is the Coronal Mass Ejection (CME) where a large amount of material is ejected from the Sun. CME’s may directly affect the earth, therefore understanding their origin is of key importance for space weather and the near-Earth environment. In this thesis, the nature and evolution of solar magnetic fields is considered through a combination of Magnetic Flux Transport Simulations and Potential Field Source Surface Models. The Magnetic Flux Transport Simulations produce a realistic description of the evolution and distribution of the radial magnetic field at the level of the solar photosphere. This is then applied as a lower boundary condition for the Potential Field Source Surface Models which prescribe a coronal magnetic field. Using these two techniques, the location and variation of coronal null points, a key element in the Magnetic Breakout Model of CMEs, are determined. Results show that the number of coronal null points follow a cyclic variation in phase with the solar cycle. In addition, they preferentially form at lower latitudes as a result of the complex active latitude field. Although a significant number of coronal nulls may exist at any one time (≈ 17), it is shown that only half may satisfy the necessary condition for breakout. From this it is concluded that while the Magnetic Breakout Model of CMEs is an important model in understanding the origin of the CMEs, other processes must occur in order to explain the observed number of CMEs. Finally, the Magnetic Flux Transport Simulations are applied to stellar magnetic fields and in particular to the fast rotating star HD171488. From this speculative study it is shown that the Magnetic Flux Transport Simulations constructed for the Sun may be applied in very different stellar circumstances and that for HD171488 a significantly higher rate of meridional flow (1200-1400 ms⁻¹) is required to match observed magnetic field distributions.

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 solar corona is a highly conductive plasma which is dominated by the coronal magnetic field. Observations show that important solar phenomena like flares or the heating of the corona are driven by magnetic energy, probably through the process of magnetic reconnection. The release of magnetic energy by reconnection requires that non-ideal processes take place in contradiction to the high conductivity of the corona. One possibility to overcome this problem is to generate strong electrical currents in strongly localised regions. In this thesis we investigate how such localised currents can be formed by slow ideal evolution of topologically simple magnetic fields. To this purpose numerical simulations are carried out using an Eulerian and a Lagrangian MHD relaxation code. We first use a simple example (twisting of a uniform field) to investigate the advantages and disadvantages of both codes and to discover possible limitations for their application. We show that for the problems addressed in this thesis the Lagrangian code is more suited because it can resolve the localised current densities much better than the Eulerian code. We then focus in particular on magnetic fields containing a so-called Hyperbolic Flux Tube (HPT). A recently proposed analytical theory predicts that HFT's are sites where under certain conditions strong current build-up can be expected. We use our code to carry out a systematic parametric study of the dependence of current growth for a typical HFT configuration. We have also developed a completely new version of the analytical theory which is directly based on the set-up of our numerical simulations. We find that the simulations agree with the analytical prediction in a quantitative way but that the analytical theory underestimates the current growth quite substantially, probably by not taking into account the non-linear character of the full problem.

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.

Although the solar corona is one of the most studied areas in solar physics, its activity, such as flares, prominence eruptions and CMEs, is far from understood. Since the solar corona is a low-ß plasma, its structure and dynamics are driven by the magnetic field. The aim of this PhD thesis to study the magnetic field in the solar corona. Unfortunately, high quality direct measurements of the coronal magnetic field are not available and theoretical extrapolation using the observed photospheric magnetic field is required. The thesis is mainly divided in two parts. The first part deals with the comparison between theoretical models of magnetic fields and observed structures in the corona. For any theoretical model, a quantitative method to fit magnetic field lines to observed coronal loops is introduced. This method provides a quantity C that measures how closely a theoretical model can reproduce the observed coronal structures. Using linear force-free field extrapolation, the above field line fitting method is used to study the evolution of an active region. The method is also illustrated when the theoretical magnetic field depends on more than one parameter. The second part of the thesis focuses on the linear force-free field assumption using two different geometric configurations. Firstly a vertical rigid magnetic flux tube is considered. The analytical expression of the magnetic field is obtained as an expansion in terms of Bessel functions. The main properties of this system are discussed and compared with two cylindrically symmetric twist profiles. For the second system, the photosphere is assumed to be an infinite plane. Using translational geometry, the analytical expression of the linear force-free magnetic field that matches a prescribed line of sight magnetic field component is obtained. This solution is compared with the non-linear solution obtained by Roumeliotis (1993).

Interest In the vertical fault problem for electromagnetic fields has been recently revived by the papers of I. d'Erceville and G. Kunetz (1962) and D. Rankin (1962). In the derivation of his equations Rankin used d'Erceville1s theory which contains some fallacious assumptions. These have been pointed out by J.T. Weaver (1962) and also in this thesis. This thesis follows the lines of mathematical attack first employed by d'Erceville and Kunetz, and later developed by Weaver, in applying the theory of integral transforms to the partial differential equations satisfied by land and sea conductors. The problem of both a vertical fault and also a sloping fault, i.e. 0 < α < 90° where α is the angle of dip of the fault are considered. The results in the general case are Inconclusive, no solution has been found and no solution is suggested. The case of α = 90° has proved to be equally indeterminate, but a solution has been suggested, which, although it has not been proved rigourously, does not appear to violate any physical principles and also seems to represent the field equations on the surface of the land and the sea. / Science, Faculty of / Earth, Ocean and Atmospheric Sciences, Department of / Graduate

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.

Magnetic fields are fundamental to the structure and dynamics of the Sun’s corona. Observations show them to be locally complex, with highly sheared and twisted fields visible in solar filaments/prominences. The free magnetic energy contained in such fields is the primary source of energy for coronal mass ejections, which are important—but still poorly understood drivers of space weather in the near-Earth environment. In this thesis, a new model is developed for the evolution of the large-scale magnetic field in the global solar corona. The model is based on observations of the radial magnetic field on the solar photosphere (visible surface). New active regions emerge, and their transport and dispersal by surface motions are simulated accurately with a surface flux transport model. The 3D coronal magnetic field is evolved in response to these photospheric motions using a magneto-frictional technique. The resulting sequence of nonlinear force-free equilibria traces the build-up of magnetic helicity and free energy over many months. The global model is applied to study two phenomena: filaments and coronal mass ejections. The magnetic field directions in a large sample of observed filaments are compared with a 6-month simulation. Depending on the twist of newly-emerging active regions, the correct chirality is simulated for up to 96% of filaments tested. On the basis of these simulations, an explanation for the observed hemispheric pattern of filament chirality is put forward, including why exceptions occur for filaments in certain locations. Twisted magnetic flux ropes develop in the simulations, often losing equilibrium and lifting off, removing helicity. The physical basis for such losses of equilibrium is demonstrated through 2D analytical models. In the 3D global simulations, the twist of emerging regions is a key parameter controlling the number of lift-offs, which may explain around a third of observed coronal mass ejections.