Secondary tearingmode in the nonlinear evolution of magnetorotational instability

Numerical investigation of the two‐dimensional magnetic reconnection is given in the context of the nonlinear evolution of the Magneto‐Rotational Instability (MRI). With a careful comparison to various theories using both one‐ and twodimensional analysis, it is found that a new stabilizing effect of the centrifugal force on tearing instability must be present in the specific geometry of the MRI.Magnetic reconnection might play a key role to the formation of the nonaxisymmetric structures observed in MRI experiments. The results may also be useful for the estimate of the accretion rate in various astrophysical objects. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Shear‐driven magnetic buoyancy oscillations

The effects of uniform horizontal shear on a stably stratified layer of gas is studied. The system is initially destabilized by a magnetically buoyant flux tube pointing in the cross‐stream direction. The shear amplifies the initial field to Lundquist numbers of about 200–400, but then its value drops to about 100–300, depending on the value of the sub‐adiabatic gradient. The larger values correspond to cases where the stratification is strongly stable and nearly isothermal. At the end of the runs the magnetic field is nearly axisymmetric, i.e. uniform in the streamwise direction. In view of Cowling's theorem the sustainment of the field remains a puzzle and may be due to subtle numerical effects that have not yet been identified in detail. In the final state the strength of the magnetic field decreases with height in such a way that the field is expected to be unstable. Low amplitude oscillations are seen in the vertical velocity even at late times, suggesting that they might be persistent (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solar dynamo models with α ‐effect and turbulent pumping from local 3D convection calculations

Results from kinematic solar dynamo models employing α ‐effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2–3 m s^–1, having scaled the local quantities with the convective velocity at the bottom of the convection zone from a solar mixing‐length model. Rotation profile of the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface shear layer on the dynamo solutions. With these choices of the small‐ and large‐scale velocity fields, we obtain estimate of the ratio of the two induction effects, C _α /C _Ω ≈ 10^–3, which we keep fixed in all models. We also include a one‐cell meridional circulation pattern having a magnitude of 10–20 m s^–1 near the surface and 1–2 m s^–1 at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the α ‐effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of latitude. Turbulent pumping of the mean fields is predominantly down‐ and equatorward. The anisotropies in the turbulent diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B _r B _ϕ . Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the α ‐effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5° … 60°). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar‐type solutions. The role of the surface shear layer is found to be important only in the case where the α ‐effect has an appreciable magnitude near the surface. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Reduction of boundary effects in the spiralMRI experiment PROMISE

Magnetorotational instability (MRI) is one of the most important and most common instabilities in astrophysics. It is widely accepted that it serves as a source of turbulent viscosity in accretion disks – the most energy efficient objects in the Universe. However it is very difficult to bring this process down on earth and model it in a laboratory experiment. Several different approaches have been proposed, one of the most recent is PROMISE (Potsdam‐ROssendorf Magnetorotational InStability Experiment). It consists of a flow of a liquid metal between two rotating cylinders under applied current‐free spiral magnetic field. The cylinders must be covered with plates which introduce additional end‐effects which alter the flow and make it more difficult to clearly distinguish between MRI stable and unstable state. In this paper we propose simple and inexpensive improvement to the PROMISE experiment which would reduce those undesirable effects. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Collapse of weakly ionized rotating turbulent Cloud Cores

In view of the Turbulent Cooling Flows scenario we carry out several 3D axisymmetric calculations to follow the evolution of magnetically subcritical weakly ionized and rotating turbulent cloud cores. Turbulent Cooling Flows appear to pronounce the effects of ambipolar diffusion considerably, inducing thereby a runaway collapse of the core already on a diluted free-fall time scale. Ambipolar diffusion significantly weakens the efficiency of magnetic braking. This implies that most of the rotational energy is trapped into the dynamically collapsing core and that initiation of outflows is prevented at least in the early isothermal phases. The trapped rotational energy is found to enhance the formation of rings that may afterwards fragment. It is shown that the central region of a strongly ionized magnetically subcritical core is principally overdense, with central density up to one order of magnitude larger than the surroundings. These results confirm that large scale magnetic fields threading a cloud core relax the supersonic random motions on an Alfvén wave crossing time. Moreover, ambipolar diffusion enhances dissipation of supersonic turbulence even more.     

Whistler wave instabilities in collisionless current sheet

Instability of whistler wave in collisionless current sheet is studied with numerical solution of the general dispersion relation obtained in Ref.[4] for the physical model A. As revealed by the results, the whistler wave can be directly absorbed by collisionless current sheets. On the neutral sheet (z/d_i = 0) oblique whistler waves over a rather wide range of wave numbers can propagate, while they are basically stable. In the ionic inertial region (z/d_i < 1), the obliquely propagating whistler wave is unstable. On the edge of the ionic inertial region (z/d_i = 1), the whistler wave is still unstable, with an increase in the growth rate, and in the frequency of the unstable wave. The growth rate is larger for the whistler wave propagating towards the neutral sheet (k_zd_i < 0) than away from the neutral sheet (k_zd_i > 0).     

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)     

Stability of the line preserving flows

We examine the equations that are used to describe flows which preserve field lines. We study what happens if we introduce perturbations to the governing equations. The stability of the line preserving flows in the case of the magneto-fluids permeated by magnetic fields is strictly connected to the non-null magnetic reconnection processes. In most of our study we use the Euler potential representation of the external magnetic field. We provide general expressions for the perturbations of the Euler potentials that describe the magnetic field. Similarly, we provide expressions for the case of steady flow as well as we obtain certain conditions required for the stability of the flow. In addition, for steady flows we formulate conditions under which the perturbations of the external field are negligible and the field may be described by its initial unperturbed form. Then we consider the flow equation that transforms quantities from the laboratory coordinate system to the related external field coordinate system. We introduce perturbations to the equation and obtain its simplified versions for the case of a steady flow. For a given system, use of this method allows us to simplify the considerations provided that some part of the system may be described as a perturbation. Next, to study regions favourable for the magnetic reconnection to occur we introduce a deviation vector to the basic line preserving flows condition equation. We provide expressions of the vector for some simplifying cases. This method allows us to examine if given perturbations either stabilise the system or induce magnetic reconnection. To illustrate some of our results we study two examples, namely a simple laboratory plasma flow and a simple planetary magnetosphere model.