Non-axisymmetric magnetostatic equilibrium

This paper considers the structural properties of a sunspot-like magnetic flux tube which lacks perfect axisymmetry. The flux tube is taken to be in static equilibrium with an atmosphere in a uniform gravity. Assuming the departure from axisymmetry to be slight, the equations for the first order non-axisymmetric part of the equilibrium are derived in cylindrical coordinates. These first order equations reduce to a linear second order hyperbolic partial differential equation in the r-z plane. Whereas Cauchy type boundary conditions are appropriate for hyperbolic equations, physical considerations dictate the specification of boundary conditions on a closed surve for our problem of interest. The construction of solutions to this boundary value problem is illustrated with three analytically soluble cases, where the zero-order axisymmetric equilibria are chosen to have magnetic field geometry of different complexity. A physical discussion of the results is given.     

On non-axisymmetric magnetic equilibria in stars

In previous work, stable approximately axisymmetric equilibrium configurations for magnetic stars were found by numerical simulation. Here, I investigate the conditions under which more complex, non-axisymmetric configurations can form. I present numerical simulations of the formation of stable equilibria from turbulent initial conditions and demonstrate the existence of non-axisymmetric equilibria consisting of twisted flux tubes lying horizontally below the surface of the star, meandering around the star in random patterns. Whether such a non-axisymmetric equilibrium or a simple axisymmetric equilibrium forms depends on the radial profile of the strength of the initial magnetic field. The results could explain observations of non-dipolar fields on stars such as the B0.2 main-sequence star τ Sco or the pulsar 1E 1207.4-5209. The secular evolution of these equilibria due to Ohmic and buoyancy processes is also examined.     

Propagation in 3D spiral-arm cosmic-ray source distribution models and secondary particle production using Picard

We study the impact of possible spiral-arm distributions of Galactic cosmic-ray sources on the flux of various cosmic-ray nuclei throughout our Galaxy. We investigate model cosmic-ray spectra at the nominal position of the sun and at different positions within the Galaxy. The modelling is performed using the recently introduced numerical cosmic ray propagation code Picard. Assuming non-axisymmetric cosmic-ray source distributions yields new insights on the behaviour of primary versus secondary nuclei.We find that primary cosmic rays are more strongly confined to the vicinity of the sources, while the distribution of secondary cosmic rays is much more homogeneous compared to the primaries. This leads to stronger spatial variation in secondary to primary ratios when compared to axisymmetric source distribution models. A good fit to the cosmic-ray data at Earth can be accomplished in different spiral-arm models, although leading to decisively different spatial distributions of the cosmic-ray flux. These lead to different cosmic ray anisotropies, where even reproducing the data becomes possible. Consequently, we advocate directions to seek best fit propagation parameters that take into account the higher complexity introduced by the spiral-arm structure on the cosmic-ray distribution. We specifically investigate whether the flux at Earth is representative for a large fraction of the Galaxy. The variance among possible spiral-arm models allows us to quantify the spatial variation of the cosmic-ray flux within the Galaxy in presence of non-axisymmetric source distributions.     

Exact MHD solutions for rotating,magnetic stars

Simple exact solutions of the magnetohydrodynamic equations are found for rotating, magnetic stars. The velocity and magnetic field are axisymmetric and purely toroidal, and the magnetic energy density equals the kinetic energy density. For constant mass density, the solution reduces to that of Chandrasekhar (1956), which is stable even against non-axisymmetric perturbations. For an ideal gas equation of state, the condition for radiative thermal equilibrium is solved to lowest order in the non-spherical perturbation. The velocity, magnetic field and non-spherical pressure and temperature perturbations all vanish within cones centered around the rotation axis, |cos |>x _i a zero of a Legendre polynomial. Low-order, long-period stellar oscillations may be excited by MHD instabilities near the equatorial region which become damped near the axis.     

Self-similar propagation of axisymmetric radiative magnetogasdynamic shocks with increasing energy,II

In this paper self-similar solutions have been investigated for the propagation of axisymmetric radiative gasdynamic shocks caused by an explosion into an inhomogeneous ideal gas permeated by a current free azimuthal magnetic field. The effects of radiation flux and magnetic field together have been seen in the region of interest on the other flow variables. The total energy of the flow between the inner expanding surface and the shock is taken to be dependent on shock radius obeying a power law. The radiative pressure and energy have been neglected.     

Investigations of spherical kinematic mean-field dynamo models

The paper supplements an earlier one on the mean-field approach to spherical kinematic dynamo models (Rädler 1980a) by results of numerical investigations. A number of dynamo models working on the basis of the α²-mechanism are considered. Cases of pure α²-mechanism are studied, which includes only the simplest form of α-effect and no other induction effect, as well as cases with several additional effects due to fluctuating or mean motions. By the pure α²-mechanism axisymmetric and non-axisymmetric fields, can be excited and maintained with nearly equal ease. Part of the additional induction effects, however, clearly favour axisymmetric fields, and others non-axisymmetric fields. The non-axisymmetric fields are waves which travel in azimuthal direction, eastward or westward, depending on the models. For special dynamo models the transition from α² to αω-mechanism and properties of the latter are investigated. The results support the presumption that the αω-mechanism is able to maintain only axisymmetric but never non-axisymmetric fields. Conditions for the occurrence of non-oscillatory or oscillatory fields are discussed, and again the influence of additional induction effects is studied. There are further presented a model with βω-mechanism maintaining an axisymmetric non-oscillatory field, and models with two kinds of δω-mechanisms allowing axisymmetric non-oscillatory and oscillatory fields. Some ideas concerning dynamo models for the Earth, the Sun and magnetic stars are discussed. It seems possible to construct dynamo models for the Earth, on the basis of the α²-mechanism which explain not only the presence of a magnetic field with a strong dipole part but also the inclination of the dipole axis against the axis of rotation, the occurrence of higher multipoles and the westward drift of the non-axisymmetric parts. Models with αω, βω or δω-mechanism, which have to be considered in the case of a strong differential rotation inside the core, provide an explanation at first only of the axisymmetric parts of the field, and the non-axisymmetric parts have then to be interpreted, for example, as MAC-waves. As far as dynamo models for the Sun are concerned, in addition to the possibility of an αω-mechanism also that of a βω or δω-mechanism is discussed, which, however, does not look not very promising. In the models developed so far, which work with the αω-mechanism, only a few of the induction effects of fluctuating motions have been included; it seems necessary to investigate also influences of other effects. The sectorial structure of the solar magnetic field can hardly be understood in terms of the traditional mean-field concept. The magnetic stars possess fields which strongly deviate from symmetry with respect to the axis of rotation. The occurrence of such fields seems understandable only if there is no noticeable differential rotation. They can be maintained by an α²-mechanism but hardly by αω, βω or δω-mechanisms.     

A dynamo model for axisymmetric and non-axisymmetric solar magnetic fields

More and more observations are showing a relatively weak, but persistent, non-axisymmetric magnetic field co-existing with the dominant axisymmetric field on the Sun. Its existence indicates that the non-axisymmetric magnetic field plays an important role in the origin of solar activity. A linear non-axisymmetric  α²– Ω  dynamo model is derived to explore the characteristics of the axisymmetric  ( m = 0)  and the first non-axisymmetric  ( m = 1)  modes and to provide a theoretical basis with which to explain the 'active longitude', 'flip-flop' and other non-axisymmetric phenomena. The model consists of an updated solar internal differential rotation, a turbulent diffusivity varying with depth, and an α-effect working at the tachocline in a rotating spherical system. The difference between the  α²–Ω  and the  α–Ω  models and the conditions that favour the non-axisymmetric modes under solar-like parameters are also presented.     

Large-Scale Structure of the Sunspot Activity and the Background Magnetic Field Formation

Large-scale distribution of the sunspot activity of the Sun has been analyzed by using a technique worked out previously (Erofeev, 1997) to study long-lived, non-axisymmetric magnetic structures with different periods of rotation. Results of the analysis have been compared with those obtained by analyzing both the solar large-scale magnetic field and large-scale magnetic field simulated by means of the well-known flux transport equation using the sunspot groups as a sole source of new magnetic flux in the photosphere. A 21-year period (1964–1985) has been examined.The rotation spectra calculated for the total time interval of two 11-year cycles indicate that sunspot activity consists of a series of discrete components (modes) with different periods of rotation. The largest-scale component of the sunspot activity reveals modes with 27-day and 28-day periods of rotation situated, correspondingly, in the northern and southern hemispheres of the Sun, and two modes with rotation periods of about 29.7 days situated in both hemispheres. Such a modal structure of the sunspot activity agrees well with that of the large-scale solar magnetic field. Moreover, the magnetic field distribution simulated with the flux transport equation also reveals the same modal structure. However, such an agreement between the large-scale solar magnetic field and both the sunspot activity and simulated magnetic field is unstable in time; so, it is absent in the northern hemisphere of the Sun during solar cycle No. 20. Thus the sources of magnetic flux responsible for formation of the large-scale, rigidly rotating magnetic patterns appear to be closely connected, but are not identical with the discrete modes of the sunspot activity.     

Magnetostatic configuration of a confined magnetic flux tube

In this paper an isolated magnetic flux tube confined in stratified atmosphere is studied for slender and axisymmetric model. The functions of the pressure, density, and temperature are expanded as a Taylor series of magnetic surface function . Several models of an isolated magnetic flux tube confined in a stratified atmosphere are constructed, and the external pressure of the stratified atmosphere decreases reasonably with increasing height. The distribution of thermal dynamic quantities and the magnetic pressure in the flux tube are also obtained.     

Sausage MHD Waves in Incompressible Flux Tubes with Twisted Magnetic Fields

Twisted magnetic flux tubes are of considerable interest because of their natural occurrence from the Sun’s interior, throughout the solar atmosphere and interplanetary space up to a wide range of applicabilities to astrophysical plasmas. The aim of the present work is to obtain analytically a dispersion equation of linear wave propagation in twisted incompressible cylindrical magnetic waveguides and find appropriate solutions for surface, body and pseudobody sausage modes (i.e. m = 0) of a twisted magnetic flux tube embedded in an incompressible but also magnetically twisted plasma. Asymptotic solutions are derived in long- and short-wavelength approximations. General solutions of the dispersion equation for intermediate wavelengths are obtained numerically. We found, that in case of a constant, but non-zero azimuthal component of the equilibrium magnetic field outside the flux tube the index ν of Bessel functions in the dispersion relation is not integer any more in general. This gives rise to a rich mode-structure of degenerated magneto-acoustic waves in solar flux tubes. In a particular case of a uniform magnetic twist the total pressure is found to be constant across the boundary of the flux tube. Finally, the effect of magnetic twist on oscillation periods is estimated under solar atmospheric conditions. It was found that a magnetic twist will increase, in general, the periods of waves approximately by a few percent when compared to their untwisted counterparts.     

Anisotropic Magnetic Confinement. Applications to Field-Reversed Configurations

Introducing an auxiliary function of the usual poloidal magnetic stream function, it is possible to obtain axisymmetric solutions of the ideal anisotropic magnetohydrodynamic equations for steady rotating plasmas, in terms of solutions of the Maschke and Perrin equation for isotropic plasmas, with temperature as a surface function. For vanishing rotation, the problem is reduced to the classical Grad–Schlüter–Shafranov equation for static equilibria. Some applications of the equilibrium models to the study of tilting stability and classical particle transport in field-reversed configurations are presented.     

On the origin of strong downdrafts associated with the birth of sunspots

During the emerging phase of sunspot pores strong downdrafts are observed in photospheric lines, whereas the upward flow is not detected within the observational accuracy. For the study of the origin of these downdrafts, we present the numerical solutions of non-steady hydrodynamic equations in one dimension along the rising magnetic flux tubes. In these solutions, it is assumed that the solar gas inside the tube is initially in the hydrostatic HSRA and then the tube rises with a given velocity and configuration. The results reproduce well the observed radial velocity distributions, hence it is concluded that the observed strong downdrafts originate from the sliding matter along the rising magnetic flux tubes.     

Models of magnetic field generation in partly stable planetary cores: Applications to Mercury and Saturn

A substantial part of Mercury's iron core may be stably stratified because the temperature gradient is subadiabatic. A dynamo would operate only in a deep sublayer. We show that such a situation arises for a wide range of values for the heat flow and the sulfur content in the core. In Saturn the upper part of the metallic hydrogen core could be stably stratified because of helium depletion. The magnetic field is unusually weak in the case of Mercury and unusually axisymmetric at Saturn. We study numerical dynamo models in rotating spherical shells with a stable outer region. The control parameters are chosen such that the magnetic Reynolds number is in the range of expected Mercury values. Because of its slow rotation, Mercury may be in a regime where the dipole contribution to the internal magnetic field is weak. Most of our models are in this regime, where the dynamo field consists mainly of rapidly varying higher multipole components. They can hardly pass the stable conducting layer because of the skin effect. The weak low-degree components vary more slowly and control the structure of the field outside the core, whose strength matches the observed field strength at Mercury. In some models the axial dipole dominates at the planet's surface and in others the axial quadrupole is dominant. Differential rotation in the stable layer, representing a thermal wind, is important for attenuating non-axisymmetric components in the exterior field. In some models that we relate to Saturn the axial dipole is intrinsically strong inside the dynamo. The surface field strength is much larger than in the other cases, but the stable layer eliminates non-axisymmetric modes. The Messenger and Bepi Colombo space missions can test our predictions that Mercury's field is large-scaled, fairly axisymmetric, and shows no secular variations on the decadal time scale.     

Non-Axisymmetric Oscillations of Thin Prominence Fibrils

We study non-axisymmetric oscillations of thin prominence fibrils. A fibril is modeled by a straight thin magnetic tube with the ends frozen in dense plasmas. The density inside and outside the tube varies only along the tube and it is discontinuous at the tube boundary. Making a viable assumption that the tube radius is much smaller than its length, we show that the squares of the frequencies of non-axisymmetric tube oscillations are given by the eigenvalues of the Sturm–Liouville problem for a second-order ordinary differential equation on a finite interval with the zero boundary conditions. For an equilibrium density that is constant outside the tube and piecewise constant inside we derived a simple dispersion equation determining the frequencies of non-axisymmetric oscillations. We carry out a parametric study of this equation both analytically and numerically, restricting our analysis to the first even mode and the first odd mode. In particular, we obtained a criterion that allows to find out if each of these modes is a normal or leaky mode.     

Three-dimensional model for generation of the mean solar magnetic field

A three-dimensional (non-axisymmetric) model for the solar mean magnetic field generation is studied. The sources of generation are the differential rotation and mean helicity in the convective shell. The system is described by two equations of the first order in time and the fourth order in space coordinates. The solution is sought for in the form of expansion over the spherical function Y_n^m. The modes of different m are separated. A finite-difference scheme similar to the Peaceman-Rachford scheme is constructed so to find coefficients of the expansion depending on the time and radial coordinates. It is shown that a mode with a smaller azimuthal number m is primarily excited. The axisymmetric mode m = o describes the 22 year solar cycle oscillations. The modes of m o have no such periodicity, the oscillate with a period of rotation of the low boundary of the solar convective shell, The solutions which are symmetric relative to the equator plane are excited more easily compared with the antisymmetrical ones. The results obtained are confronted to the observational picture of the non-axisymmetric large-scale solar magnetic fields.     

Synoptic Solar Cycle 24 in Corona,Chromosphere, and Photosphere Seen by the Solar Dynamics Observatory

The Solar Dynamics Observatory provides multiwavelength imagery from extreme ultraviolet (EUV) to visible light as well as magnetic-field measurements. These data enable us to study the nature of solar activity in different regions of the Sun, from the interior to the corona. For solar-cycle studies, synoptic maps provide a useful way to represent global activity and evolution by extracting a central meridian band from sequences of full-disk images over a full solar Carrington rotation (≈?27.3 days). We present the global evolution during Solar Cycle 24 from 20 May 2010 to 31 August 2013 (CR?2097?–?CR?2140), using synoptic maps constructed from full-disk, line-of-sight magnetic-field imagery and EUV imagery (171 Å, 193 Å, 211 Å, 304 Å, and 335 Å). The synoptic maps have a resolution of 0.1 degree in longitude and steps of 0.001 in sine of latitude. We studied the axisymmetric and non-axisymmetric structures of solar activity using these synoptic maps. To visualize the axisymmetric development of Cycle 24, we generated time–latitude (also called butterfly) images of the solar cycle in all of the wavelengths, by averaging each synoptic map over all longitudes, thus compressing it to a single vertical strip, and then assembling these strips in time order. From these time–latitude images we observe that during the ascending phase of Cycle 24 there is a very good relationship between the integrated magnetic flux and the EUV intensity inside the zone of sunspot activities. We observe a North–South asymmetry of the EUV intensity in high-latitudes. The North–South asymmetry of the emerging magnetic flux developed and resulted in a consequential asymmetry in the timing of the polar magnetic-field reversals.     

Instability of a streaming jet under the influence of vacuum magnetic fields

The magnetohydrodynamic stability of a streaming liquid cylinder subject to surface tension and pervaded by a magnetic vacuum field has been elaborated for all axisymmetric and non-axisymmetric disturbances. The dispersion relation is obtained and studied analytically and numerically. the streaming has always a destabilizing effect. The axial magnetic fields inside and outside the jet have always stabilizing effects for all perturbations. The transverse magnetic field has a destabilizing effect. However, if the axial field intensity is so high and paramount over that the transverse field, the destabilizing character of the model is suppressed. The latter is satisfied if the Alfvén wave velocity is greater than the equilibrium liquid velocity.     

Analytical solutions for the Newtonian gravitational field induced by matter within axisymmetric boundaries

An analytical method originally applied to the problem of the actuator disc in fluid mechanics has been applied to the closely analogous problem of constructing the classical Newtonian potential and attractions. The method can treat axisymmetric problems and also non-axisymmetric cases where matter is confined within axisymmetric boundaries. The potential and attractions for the generalized thin finite disc can be given in closed form in terms of elliptic integrals and elementary functions. For the general case of matter within an axisymmetric boundary, the potentials and attractions can be evaluated as one-dimensional integrals of albeit complex analytical expressions. These expressions represent the fields induced by matter in an extended region as a distribution of gravitating discs. For certain special cases, such as matter bounded by a circular cylinder and also for matter distributed in a spherical region, closed-form solutions can be given that appear to be new. Some non-axisymmetric results are also given for the thin disc of infinite radial extent.     

Numerical simulations of rotating axisymmetric sunspots

A numerical model of axisymmetric convection in the presence of a vertical magnetic flux bundle and rotation about the axis is presented. The model contains a compressible plasma described by the non-linear MHD equations, with density and temperature gradients simulating the upper layer of the Sun's convection zone. The solutions exhibit a central magnetic flux tube in a cylindrical numerical domain, with convection cells forming collar flows around the tube. When the numerical domain is rotated with a constant angular velocity, the plasma forms a Rankine vortex, with the plasma rotating as a rigid body where the magnetic field is strong, as in the flux tube, while experiencing sheared azimuthal flow in the surrounding convection cells, forming a free vortex. As a result, the azimuthal velocity component has its maximum value close to the outer edge of the flux tube. The azimuthal flow inside the magnetic flux tube and the vortex flow is prograde relative to the rotating cylindrical reference frame. A retrograde flow appears at the outer wall. The most significant convection cell outside the flux tube is the location for the maximum value of the azimuthal magnetic field component. The azimuthal flow and magnetic structure are not generated spontaneously, but decay exponentially in the absence of any imposed rotation of the cylindrical domain.