The kink instability of magnetohydrodynamics is believed to be fundamental to many aspects of the dynamic activity of the solar atmosphere, such as the initiation of flares and the heating of the solar corona. In this work, we investigate the importance of viscosity on the kink instability. In particular, we focus on two forms of viscosity; isotropic viscosity (independent of the magnetic field) and anisotropic viscosity (with a preferred direction following the magnetic field). Through the detailed analysis of magnetohydrodynamic simulations of the kink instability with both types of viscosity, we show that the form of viscosity has a significant effect on the nonlinear dynamics of the instability. The different viscosities allow for different flow and current structures to develop, thus affecting the behaviour of magnetic relaxation, the formation of secondary instabilities and the Ohmic and viscous heating produced. Our results have important consequences for the interpretation of solar observations of the kink instability.

Stellar radiative zones are typically assumed to be motionless in standard models of stellar structure but there is sound theoretical and observational evidence that this cannot be the case. We investigate by direct numerical simulations a three-dimensional and time-dependent model of stellar radiation zones consisting of an electrically-conductive and stably-stratified anelastic fluid confined to a rotating spherical shell and driven by a baroclinic torque. As the baroclinic driving is gradually increased a sequence of transitions from an axisymmetric and equatorially-symmetric time-independent flow to flows with a strong poloidal component and lesser symmetry are found. It is shown that all flow regimes characterised with significant non-axisymmetric components are capable of generating self-sustained magnetic field. As the value of the Prandtl number is decreased and the value of the Ekman number is decreased flows become strongly time-dependent with progressively complex spatial structure and dynamos can be generated at lower values of the magnetic Prandtl number.

The dynamics of stably stratified stellar radiative zones is of considerable interest due to the availability of increasingly detailed observations of Solar and stellar interiors. This article reports the first non-axisymmetric and time-dependent simulations of flows of anelastic fluids driven by baroclinic torques in stably stratified rotating spherical shells -- a system serving as an elemental model of a stellar radiative zone. With increasing baroclinicity a sequence of bifurcations from simpler to more complex flows is found in which some of the available symmetries of the problem are broken subsequently. The poloidal component of the flow grows relative to the dominant toroidal component with increasing baroclinicity. The possibility of magnetic field generation thus arises and this paper proceeds to provide some indications for self-sustained dynamo action in baroclinically-driven flows. We speculate that magnetic fields in stably stratified stellar interiors are thus not necessarily of fossil origin as it is often assumed.

Numerical MHD simulations play increasingly important role for understanding mechanisms of stellar magnetism. We present simulations of convection and dynamos in density-stratified rotating spherical fluid shells. We employ a new 3D simulation code for the solution of a physically consistent anelastic model of the process with a minimum number of parameters. The reported dynamo simulations extend into a "buoyancy-dominated" regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients and strong anti-solar differential rotation develops as a result. We find that the self-generated magnetic fields, despite being relatively weak, are able to reverse the direction of differential rotation from anti-solar to solar-like. We also find that convection flows in this regime are significantly stronger in the polar regions than in the equatorial region, leading to non-oscillatory dipole-dominated dynamo solutions, and to concentration of magnetic field in the polar regions. We observe that convection has different morphology in the inner and at the outer part of the convection zone simultaneously such that organized geostrophic convection columns are hidden below a near-surface layer of well-mixed highly-chaotic convection. While we focus the attention on the buoyancy-dominated regime, we also demonstrate that conical differential rotation profiles and persistent regular dynamo oscillations can be obtained in the parameter space of the rotation-dominated regime even within this minimal model.

Aims: To investigate the relationship between surges and magnetic reconnection during the emergence of small-scale active regions. In particular, to examine how the large-scale geometry of the magnetic field, shaped by different phases of reconnection, guides the flowing of surges. Methods: We present three flux emergence models. The first model, and the simplest, consists of a region emerging into a horizontal ambient field that is initially parallel to the top of the emerging region. The second model is the same as the first but with an extra smaller emerging region which perturbs the main region. This is added to create a more complex magnetic topology and to test how this complicates the development of surges compared to the first model. The last model has a non-uniform ambient magnetic field to model the effects of emergence near a sunspot field and impose asymmetry on the system through the ambient magnetic field. At each stage, we trace the magnetic topology to identify the locations of reconnection. This allows for field lines to be plotted from different topological regions, highlighting how their geometry affects the development of surges. Results: In the first model, we identify distinct phases of reconnection. Each phase is associated with a particular geometry for the magnetic field and this determines the paths of the surges. The second model follows a similar pattern to the first but with a more complex magnetic topology and extra eruptions. The third model highlights how an asymmetric ambient field can result in preferred locations for reconnection, subsequently guiding the direction of surges.

The onset of convection in a rotating cylindrical annulus with parallel ends filled with a compressible fluid is studied in the anelastic approximation. Thermal Rossby waves propagating in the azimuthal direction are found as solutions. The analogy to the case of Boussinesq convection in the presence of conical end surfaces of the annular region is emphasized. As in the latter case the results can be applied as an approximation for the description of the onset of anelastic convection in rotating spherical fluid shells. Reasonable agreement with three-dimensional numerical results published by Jones et al. (J. Fluid Mech., vol. 634, 2009, pp. 291-319) for the latter problem is found. As in those results the location of the onset of convection shifts outward from the tangent cylinder with increasing number $N_\rho$ of density scale heights until it reaches the equatorial boundary. A new result is that at a much higher number $N_\rho$ the onset location returns to the interior of the fluid shell.

Aims. Linear magnetohydrostatic (MHS) models of solar magnetic fields balance plasma pressure gradients, gravity and Lorentz forces where the current density is composed of a linear force-free component and a cross-field component that depends on gravitational stratification. In this paper, we investigate an efficient numerical procedure for calculating such equilibria. Methods. The MHS equations are reduced to two scalar elliptic equations - one on the lower boundary and the other within the interior of the computational domain. The normal component of the magnetic field is prescribed on the lower boundary and a multigrid method is applied on both this boundary and within the domain to find the poloidal scalar potential. Once solved to a desired accuracy, the magnetic field, plasma pressure and density are found using a finite difference method. Results. We investigate the effects of the cross-field currents on the linear MHS equilibria. Force-free and non-force-free examples are given to demonstrate the numerical scheme and an analysis of speed-up due to parallelization on a graphics processing unit (GPU) is presented. It is shown that speed-ups of x30 are readily achievable.

This document presents an accurate record of the scientific programme of the 2013 UK National Conference on Geophysical, Astrophysical and Industrial Magnetohydrodynamics held on 23rd and 24th of May 2013 at the School of Mathematics and Statistics of the University of Glasgow. The abstracts of presented contributions are listed in chronological order and a full list of participants and their affiliations is included.

A physically consistent model of magnetic field generation by convection in a rotating spherical shell with a minimum of parameters is applied to the Sun. Despite its unrealistic features the model exhibits a number of properties resembling those observed on the Sun. The model suggests that the large scale solar dynamo is dominated by a non-axisymmetric $m=1$ component of the magnetic field.

Simple models of magnetic field generation by convection in rotating spherical shells exhibit properties resembling those observed on the sun. The {assumption of the Boussinesq approximation made in these models} prevents a realistic description of the solar cycle, but through a physically motivated change in the boundary condition for the differential rotation the propagation of dynamo waves towards higher latitudes can be reversed at least at low latitudes.

Possibilities and difficulties of applying the theory of magnetic field generation by convection flows in rotating spherical fluid shells to the Giant Planets are outlined. Recent progress in the understanding of the distribution of electrical conductivity in the Giant Planets suggests that the dynamo process occurs predominantly in regions of semiconductivity. In contrast to the geodynamo the magnetic field generation in the Giant Planets is thus characterized by strong radial conductivity variations. The importance of the constraint on the Ohmic dissipation provided by the planetary luminosity is emphasized. Planetary dynamos are likely to be of an oscillatory type, although these oscillations may not be evident from the exterior of the planets.

A convection-driven MHD dynamo in a rotating spherical shell, with clearly defined structural elements in the flow and magnetic field, is simulated numerically. Such dynamos can be called deterministic, in contrast to those explicitly dependent on the assumed properties of turbulence. The cases most interesting from the standpoint of studying the nature of stellar magnetism demonstrate the following features. On a global scale, the convective flows can maintain a ``general'' magnetic field with a sign-alternating dipolar component. Local (in many cases, bipolar) magnetic structures are associated with convection cells. Disintegrating local structures change into background fields, which drift toward the poles. From time to time, reversals of the magnetic fields in the polar regions occur, as ``new'' background fields expel the ``old'' fields.