The stability of MHD Taylor‐Couette flow with current‐free spiral magnetic fields between conducting cylinders

We study the magnetorotational instability in cylindrical Taylor‐Couette flow, with the (vertically unbounded) cylinders taken to be perfect conductors, and with externally imposed spiral magnetic fields. The azimuthal component of this field is generated by an axial current inside the inner cylinder, and may be slightly stronger than the axial field. We obtain an instability beyond the Rayleigh line, for Reynolds numbers of order 1000 and Hartmann numbers of order 10, and independent of the (small) magnetic Prandtl number. For experiments with R_out = 2R_in = 10 cm and Ω_out = 0.27 Ω_in, the instability appears for liquid sodium for axial fields of ∼20 Gauss and axial currents of ∼1200 A. For gallium the numbers are ∼50 Gauss and ∼3200 A. The vertical cell size is about twice the cell size known for nonmagnetic experiments. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Energy distribution in nonaxisymmetric magnetic Taylor-Couette flow

Azimuthal magnetorotational instability is a mechanism that generates nonaxisymmetric field pattern. Nonlinear simulations in an infinite Taylor-Couette system with current-free external field show, that not only the linearly unstable mode m = 1 appears, but also an inverse cascade transporting energy into the axisymmetric field is possible. By varying the Reynolds number of the flow and the Hartmann number for the magnetic field, we find that the ratio between axisymmetric (m = 0) and dominating nonaxisymmetric mode (m = 1) can be nearly free chosen. On the surface of the outer cylinder this mode distribution appears similarly, but with weaker axisymmetric fields.We do not find significant differences in the case that a constant current within the flow is added. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Helical magnetorotational instability of Taylor‐Couette flows in the Rayleigh limit and for quasi‐Kepler rotation

The magnetorotational instability (MRI) of differential rotation under the simultaneous presence of axial and azimuthal components of the (current‐free) magnetic field is considered. For rotation with uniform specific angular momentum the MHD equations for axisymmetric perturbations are solved in a local short‐wave approximation. All the solutions are overstable for B_z · B_ϕ ≠ 0 with eigenfrequencies approaching the viscous frequency. For more flat rotation laws the results of the local approximation do not comply with the results of a global calculation of the MHD instability of Taylor‐Couette flows between rotating cylinders. – With B_ϕ and B_z of the same order the traveling‐mode solutions are also prefered for flat rotation laws such as the quasi‐Kepler rotation. For magnetic Prandtl number Pm → 0 they scale with the Reynolds number of rotation rather than with the magnetic Reynolds number (as for standard MRI) so that they can easily be realized in MHD laboratory experiments. – Regarding the nonaxisymmetric modes one finds a remarkable influence of the ratio B_ϕ/B_z only for the extrema. For B_ϕ ≫ B_z and for not too small Pm the nonaxisymmetric modes dominate the traveling axisymmetric modes. For standard MRI with B_z ≫ B_ϕ, however, the critical Reynolds numbers of the nonaxisymmetric modes exceed the values for the axisymmetric modes by many orders so that they are never prefered. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonlinear simulations of magnetic Taylor‐Couette flow with currentfree helical magnetic fields

Themagnetorotational instability (MRI) in cylindrical Taylor‐Couette flow with external helical magnetic field is simulated for infinite and finite aspect ratios. We solve the MHD equations in their small Prandtl number limit and confirm with timedependent nonlinear simulations that the additional toroidal component of the magnetic field reduces the critical Reynolds number from O (10⁶) (axial field only) to O (10³) for liquid metals with their small magnetic Prandtl number. Computing the saturated state we obtain velocity amplitudes which help designing proper experimental setups. Experiments with liquid gallium require axial field ∼50 Gauss and axial current ∼4 kA for the toroidal field. It is sufficient that the vertical velocity u_z of the flow can be measured with a precision of 0.1 mm/s.We also show that the endplates enclosing the cylinders do not destroy the traveling wave instability which can be observed as presented in earlier studies. For TC containers without and with endplates the angular momentum transport of the MRI instability is shown as to be outwards. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Tayler instability with Hall effect in young neutron stars

Collapse calculations indicate that the hot young neutron stars rotate differentially so that strong toroidal magnetic field components should exist in the outer shell where also the Hall effect appears to be important when the Hall parameter = ω_Bτ exceeds unity. The amplitudes of the induced toroidal magnetic fields are limited by the current‐induced Tayler instability. An important characteristics of the Hall effect is its distinct dependence on the sign of the magnetic field. We find for fast rotation that positive (negative) Hall parameters essentially reduce (increase) the stability domain. It is thus concluded that the toroidal field belts in young neutron stars induced by their differential rotation should have different amplitudes in both hemispheres which later are frozen in. Due to the effect of magnetic suppression of the heat conductivity also the brightness of the two hemispheres should be different. As a possible example for our scenario the isolated neutron star RBS 1223 is considered which has been found to exhibit different X‐ray brightness at both hemispheres (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Axisymmetry vs. nonaxisymmetry of a Taylor‐Couette flow with azimuthal magnetic fields

The instability of a supercritical Taylor‐Couette flow of a conducting fluid with resting outer cylinder under the influence of a uniform axial electric current is investigated for magnetic Prandtl number Pm = 1. In the linear theory the critical Reynolds number for axisymmetric perturbations is not influenced by the current‐induced axisymmetric magnetic field but all axisymmetric magnetic perturbations decay. The nonaxisymmetric perturbations with m = 1 are excited even without rotation for large enough Hartmann numbers (“Tayler instability”). For slow rotation their growth rates scale with the Alfvén frequency of the magnetic field but for fast rotation they scale with the rotation rate of the inner cylinder. In the nonlinear regime the ratio of the energy of the magnetic m = 1 modes and the toroidal background field is very low for the non‐rotating Tayler instability but it strongly grows if differential rotation is present. For super‐Alfv´enic rotation the energies in the m = 1 modes of flow and field do not depend on the molecular viscosity, they are almost in equipartition and contain only 1.5 % of the centrifugal energy of the inner cylinder. The geometry of the excited magnetic field pattern is strictly nonaxisymmetric for slow rotation but it is of the mixed‐mode type for fast rotation – contrary to the situation which has been observed at the surface of Ap stars. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Subcritical dynamos in shear flows

Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non‐rotating wall‐bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non‐axisymmetric instabilities of shearinduced axisymmetric toroidal velocity or magnetic fields, such as Kelvin‐Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high‐Reynolds number shear flows could act as a large‐scale driving mechanism of turbulent flows that would in turn generate an independent small‐scale dynamo. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Tayler instability of toroidal magnetic fields in a columnar gallium experiment

The nonaxisymmetric Tayler instability of toroidal magnetic fields due to axial electric currents is studied for conducting incompressible fluids between two coaxial cylinders without endplates. The inner cylinder is considered as so thin that the limit of R_in → 0 can be computed. The magnetic Prandtl number is varied over many orders of magnitudes but the azimuthal mode number of the perturbations is fixed to m = 1. In the linear approximation the critical magnetic field amplitudes and the growth rates of the instability are determined for both resting and rotating cylinders. Without rotation the critical Hartmann numbers do not depend on the magnetic Prandtl number but this is not true for the corresponding growth rates. For given product of viscosity and magnetic diffusivity the growth rates for small and large magnetic Prandtl number are much smaller than those for Pm = 1. For gallium under the influence of a magnetic field at the outer cylinder of 1 kG the resulting growth time is 5 s. The minimum electric current through a container of 10 cm diameter to excite the instability is 3.20 kA. For a rotating container both the critical magnetic field and the related growth times are larger than for the resting column (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The equilibrium and the stability of the Jeans spheroids with toroidal magnetic fields

We examine the effect of a toroidal magnetic field on the equilibrium and stability of homogeneous masses distorted by the tidal effects of a secondary (of massM at a distanceR). It is shown that if the toroidal magnetic field is assumed to be axisymmetric about the direction of the line joining the centres of mass of the primary and the secondary, then the equilibrium configuration is a prolate spheroid. Also determined are the characteristic frequencies of the various modes of oscillations belonging to the second harmonics. It is found that the magnetic field shifts these frequencies to higher values than the ones which prevail in the absence of a magnetic field.     

Properties of the negative effective magnetic pressure instability

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large‐scale instability. In this study, hydromagnetic mean‐field modelling is performed for an isothermally stratified layer in the presence of a horizontal magnetic field. The negative effective magnetic pressure instability (NEMPI) is comprehensively investigated. It is shown that, if the effect of turbulence on the mean magnetic tension force vanishes, which is consistent with results from direct numerical simulations of forced turbulence, the fastest growing eigenmodes of NEMPI are two‐dimensional. The growth rate is found to depend on a parameter β_* characterizing the turbulent contribution of the effective mean magnetic pressure for moderately strong mean magnetic fields. A fit formula is proposed that gives the growth rate as a function of turbulent kinematic viscosity, turbulent magnetic diffusivity, the density scale height, and the parameter β_*. The strength of the imposed magnetic field does not explicitly enter provided the location of the vertical boundaries are chosen such that the maximum of the eigenmode of NEMPI fits into the domain. The formation of sunspots and solar active regions is discussed as possible applications of NEMPI (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The role of toroidal magnetic field in the stability of accretion disks with alpha viscosity and other viscosities

The standard thin accretion disk model predicts that the inner regions of alpha model disks, where radiation pressure is dominant, are thermally and viscously unstable. However, observations show that the bright X-ray binaries and AGN accretion disks, corresponding to radiation-pressure thin disks, are stable. In this paper, we reconsider the linear and local instability of accretion disks in the presence of a toroidal magnetic field. In the basic equations, we consider physical quantities such as advection, thermal conduction, arbitrary viscosity, and an arbitrary cooling function also. A fifth order diffusion equation is obtained and is solved numerically. The solutions are compared to non-magnetic cases. The results show that the toroidal magnetic field can make the thermal instability in radiation pressure-dominated slim disks disappear if ? _m ≥0.3. However, it causes a more thermal instability in radiation pressure alpha disks without advection. Also, we consider the thermal instability in accretion disks with other values of the viscosity and obtain a general criterion for thermal instability in the long-wavelength limit and in the presence of a toroidal magnetic field.     

Non-radial oscillations and convective instability of a polytrope with a toroidal magnetic field

We examine the non-radial modes of oscillation, belonging to spherical harmonics of ordersl=1 andl=3, of a gaseous polytrope with a toroidal magnetic field. We find that a toroidal magnetic field increases the growth rate of convective instability for deformations belonging to the spherical harmonicl=1 whereas it decreases the growth rate of convective instability for deformations belonging to the harmonicsl=2 andl=3. The frequencies of the ‘acoustic’ mode and the ‘Kelvin’ mode are decreased by the presence of the toroidal magnetic field.     

Wound magnetostatics and flaring

We discuss the winding of a force-free axisymmetric magnetic field rooted on a heavy conductor onz=0. In quadrupolar symmetry the field expands in the half-spacez>0 and the toroidal flux concentrates on a conical surface. After a mean twist of 208°, the conical layer hosts large toroidal current loops with reversal of the magnetic flux on either side. The evolution of the field structure is described by scale-free static solutionsBr ^–(p+2), withp taking values between 0 and 2. The large expansion factor of the field structure is suggestive of flaring originating on the solar photosphere.     

A numerical study of the rotational behavior of the white dwarf in DQ Herculis

It is well-known that the optical pulsations in DQ Her are due to emission from the magnetic poles of the white dwarf. As the white dwarf spins on its axis, the magnetic poles sweep into and out of the line of sight due to the fact that the magnetic axis and the spin axis are not aligned, that is, the DQ Her white dwarf is an `oblique rotator'. So, a central question is if an initially axisymmetric model simulating the DQ Her white dwarf before its `turn-over' (where the term `turn-over' describes the process by which the magnetic axis gets inclining relative to the spin axis at a progressively increasing angle, the so-called `turn-over angle') is indeed susceptible to turn-over. For the puprose of resolving this problem, we compute several axisymmetric models of the DQ Her white dwarf. Our results show that, for both the rotation periods proposed on the basis of the observational evidence regarding the optical pulsations of DQ Her (i.e.,71 s or 142 s), the moment of inertia along the rotation axis is less than the corresponding moment of inertia along the remaining two principal axes of the axisymmetric configuration, I ₃₃ > I ₁₁(=I ₂₂). This is because toroidal magnetic field (tending to derive prolate equidensity surfaces) dominates over rotation (tending, in turn, to derive oblate equidensity surfaces), mainly in the interior of the star. The situation I ₁₁ < I ₃₃ is known as `dynamical asymmetry', and can cause a turn-over of the magnetic symmetry axis with respect to the rotation axis, eventually deriving a nonaxisymmetric configuration corresponding to the so-called `perpendicular rotator' with turn-over angle almost equal to 90^°. In this view, our results explain why the DQ Her white dwarf is now an oblique rotator. This revised version was published online in July 2006 with corrections to the Cover Date.     

Three-Dimensional Particle Acceleration in Electromagnetic Cylinder and Torus

Particle acceleration via Poynting vector with toroidal magnetic field is studied in 3D PIC simulation of electron-positron plasma. We choose two different initial magnetic field configurations to compare how the particle acceleration is affected by the expansion of electromagnetic wave. In the cylindrical case, the electromagnetic field strength decays as (ct)⁻², and particles are accelerated in the radial direction as well as the axial direction. Rayleigh-Taylor instability is also observed at the center of the cylinder. In the torus case, the field strength decays as (ct)⁻³, making the acceleration less efficient. Particles accelerated in the axial direction by E × B force creates strong charge separation.     

Stability of toroidal magnetic fields in the solar tachocline and beneath

Stability of toroidal magnetic field in a stellar radiation zone is considered for the cases of uniform and differential rotation. In the rigidly rotating radiative core shortly below the tachocline, the critical magnetic field for instability is about 600 G. The unstable disturbances for slightly supercritical fields have short radial scales ∼1 Mm. Radial mixing produced by the instability is estimated to conclude that the internal field of the sun can exceed the critical value of 600 G only marginally. Otherwise, the mixing is too strong and not compatible with the observed lithium abundance. Analysis of joint instability of differential rotation and toroidal field leads to the conclusion that axisymmetric models of the laminar solar tachocline are stable to nonaxisymmetric disturbances. The question of whether sun-like stars can posses tachoclines is addressed with positive answer for stars with rotation periods shorter than about two months. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetorotational supernovae. Magnetorotational instability. Jet formation

2D numerical simulations of magnetorotational (MR) supernova mechanism are described. It is shown that magnetic field is amplified due to the differential rotation after core collapse. When magnetic pressure reaches some level, a compression wave starts to move outwards. Moving along steeply decreasing density profile the compression wave transforms quickly into fast MHD shock. The magnetorotational instability (MRI) was found in our simulations. MRI leads to the exponential growth of the components of the magnetic field. The MRI significantly reduces MR supernova explosion time. Configuration of the initial magnetic field qualitatively defines the shape of MR supernova explosion. For the quadrupole-like initial poloidal field the MR supernova explosion develops mainly along equatorial plane, the dipole-like initial field results in MR supernova developing as mildly collimated jet along axis of rotation. The explosion energy of MR supernova found in our simulations is ∼0.5–0.6×10⁵¹ erg.     

Magnetohydrodynamical non-radiative accretion flows in two dimensions

We present the results of axisymmetric, time-dependent magnetohydrodynamic simulations of accretion flows around black holes. The calculations begin from a rotationally supported thick torus which contains a weak poloidal field. Accretion is produced by growth and saturation of the magnetorotational instability (MRI) provided that the wavelength of the fastest growing mode is less than the thickness of the torus. Using a computational grid that spans more than two decades in radius, we compare the time-averaged properties of the flow with previous hydrodynamical simulations. The net mass accretion rate is small compared with the mass inflow and outflow rates at large radii associated with turbulent eddies. Turbulence is driven by the MRI rather than convection. The two-dimensional structure of the time-averaged flow is significantly different compared with the hydrodynamical case. We discuss the limitations imposed on our results by the assumption of axisymmetry and the relatively small radial domain.     

Stability of relative equilibria in arbitrary axisymmetric gravitational and magnetic fields

Relative equilibria occur in a wide variety of physical applications, including celestial mechanics, particle accelerators, plasma physics, and atomic physics. We derive sufficient conditions for Lyapunov stability of circular orbits in arbitrary axisymmetric gravitational (electrostatic) and magnetic fields, including the effects of local mass (charge) and current density. Particularly simple stability conditions are derived for source‐free regions, where the gravitational field is harmonic (∇²U = 0) or the magnetic field irrotational (∇ × B = 0). In either case the resulting stability conditions can be expressed geometrically (coordinate‐free) in terms of dimensionless stability indices. Stability bounds are calculated for several examples, including the problem of two fixed centers, the J₂ planetary model, galactic disks, and a toroidal quadrupole magnetic field. This revised version was published online in July 2006 with corrections to the Cover Date.