Search

The stability of magnetized flows of a nonuniformly rotating layer of an electrically conducting nanofluid is investigated with regard to the effects of Brownian diffusion and thermophoresis. In the absence of a temperature gradient, new types of magnetorotational instability (MRI) are considered in axial, azimuthal, and helical magnetic fields in thin layers of a nanofluid. The increments and development regions of these instabilities are obtained depending on the angular velocity profile (Rossby number Ro) and the radial wavenumber k. In the presence of temperature and nanoparticle concentration gradients, stationary regimes of nonuniformly rotating convection in axial and helical magnetic fields are studied. Expressions are obtained for the critical Rayleigh numbers Rast, and neutral stability curves are plotted depending on the angular velocity profile, the profile of the external azimuthal magnetic field (the magnetic Rossby number Rb), and the radial wavenumber k. Conditions for stabilization and destabilization of stationary convection in axial and helical magnetic fields are found.

Convection has been discussed in the field of accretion discs for several decades, both as a means of angular momentum transport and also because of its role in controlling discs' vertical structure via heat transport. If the gas is sufficiently ionized and threaded by a weak magnetic field, convection might interact in non-trivial ways with the magnetorotational instability (MRI). Recently, vertically stratified local simulations of the MRI have reported considerable variation in the angular momentum transport, as measured by the stress to thermal pressure ratio $\alpha$, when convection is thought to be present. Although MRI turbulence can act as a heat source for convection, it is not clear how the instabilities will interact dynamically. Here we aim to investigate the interplay between the two instabilities in controlled numerical experiments, and thus isolate the generic features of their interaction. We perform vertically stratified, 3D MHD shearing box simulations with a perfect gas equation of state with the conservative, finite-volume code PLUTO. We find two characteristic outcomes of the interaction between the two instabilities: straight MRI and MRI/convective cycles, with the latter exhibiting alternating phases of convection-dominated flow (during which the turbulent transport is weak) and MRI-dominated flow. During the latter phase we find that $\alpha$ is enhanced by nearly an order of magnitude, reaching peak values of $\sim 0.08$. In addition, we find that convection in the non-linear phase takes the form of large-scale and oscillatory convective cells. Convection can also help the MRI persist to lower Rm than it would otherwise do. Finally we discuss how our results help interpret simulations of Dwarf Novae., Comment: Accepted for publication in MNRAS (23 pages, 14 figures, 1 table)

Instability modes of viscoelastic Taylor–Couette flow in the Keplerian regime are investigated using both linear stability analysis and experimental detection of critical states. A generalised Rayl...

We investigate whether the regular Galilean satellites could have formed in the dead zone of a circumplanetary disc. A dead zone is a region of weak turbulence in which the magnetorotational instability (MRI) is suppressed, potentially an ideal environment for satellite formation. With the grid-based hydrodynamic code, FARGO3D, we examine the evolution of a circumplanetary disc model with a dead zone. Material accumulates in the dead zone of the disc leading to a higher total mass and but a similar temperature profile compared to a fully turbulent disc model. The tidal torque increases the rate of mass transport through the dead zone leading to a steady state disc with a dead zone that does not undergo accretion outbursts. We explore a range of disc, dead zone and mass inflow parameters and find that the maximum mass of the disc is around 0.001 MJ . Since the total solid mass of such a disc is much lower, we find that there is not sufficient material in the disc for in situ formation of the Galilean satellites and that external supplement is required., Comment: 10 pages, 6 figures

We explore the pulsationally driven orbital mass ejection mechanism for Be star disc formation using isothermal, 3D magnetohydrodynamic (MHD) and hydrodynamic simulations. Non-radial pulsations are added to a star rotating at 95\% of critical as an inner boundary condition that feeds gas into the domain. In MHD, the initial magnetic field within the star is weak. The hydrodynamics simulation has limited angular momentum transport, resulting in repeating cycles of mass accumulation into a rotationally-supported disc at small radii followed by fall-back onto the star. The MHD simulation, conversely, has efficient (Maxwell $\alpha_{\rm M}$ $\sim$ 0.04) angular momentum transport provided by both of turbulent and coherent magnetic fields; a slowly decreting midplane driven by the magnetorotational instability and a supersonic wind on the surface of the disc driven by global magnetic torques. The angle and time-averaged properties near the midplane agree reasonably well with a 1D viscous decretion disc model with a modified $\tilde\alpha=0.5$, in which the gas transitions from a subsonic thin disc to a supersonic spherical wind at the critical point. 1D models, however, cannot capture the multi-phase decretion/angular structure seen in our simulations. Our results demonstrate that, at least under certain conditions, non-radial pulsations on the surface of a rapidly rotating, weakly magnetized star can drive a Keplerian disc with the basic properties of the viscous decretion disc paradigm, albeit coupled to a laminar wind away from the midplane. Future modeling of Be star discs should consider the possible existence of such a surface wind., Comment: Accepted by MNRAS. 16 pages, 10 figures

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two and three dimensional numerical modelling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical ($z$) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal ($rz$) plane, driving magnetic reconnection when the Lundquist number $S > 3 \times 10^2$, making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at $S = 5 \times 10^3$, self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When non-axisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets further unstable, and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation., 12 pages, 6 figures, Accepted for publication to ApJ Letters

In the present study, we addressed the question of how the growth rate of the magnetorotational instability is modified when the radial component of the stellar dipole magnetic field is taken into account in addition to the vertical component. Considering a fiducial radius in the disc where diamagnetic currents are pronounced, we carried out a linear stability analysis to obtain the growth rates of the magnetorotational instability for various parameters such as the ratio of the radial-to-vertical component and the gradient of the magnetic field, the Alfvenic Mach number and the diamagnetization parameter. Our results show that the interaction between the diamagnetic current and the radial component of the magnetic field increases the growth rate of the magnetorotational instability and generates a force perpendicular to the disc plane which may induce a torque. It is also shown that considering the radial component of the magnetic field and taking into account a radial gradient in the vertical component of the magnetic field causes an increase in the magnitudes of the growth rates of both the axisymmetric ($m=0$) and the non-axisymmetric ($m=1$) modes., Comment: 9 pages, 5 figures, published in MNRAS, 2020, vol. 491, pages 5481-5488

We investigate acceleration and propagation processes of high-energy particles inside hot accretion flows. The magnetorotational instability (MRI) creates turbulence inside accretion flows, which triggers magnetic reconnection and may produce non-thermal particles. They can be further accelerated stochastically by the turbulence. To probe the properties of such relativistic particles, we perform magnetohydrodynamic simulations to obtain the turbulent fields generated by the MRI, and calculate orbits of the high-energy particles using snapshot data of the MRI turbulence. We find that the particle acceleration is described by a diffusion phenomenon in energy space with a diffusion coefficient of the hard-sphere type: $D_\epsilon\propto \epsilon^2$, where $\epsilon$ is the particle energy. Eddies in the largest scale of the turbulence play a dominant role in the acceleration process. On the other hand, the stochastic behaviour in configuration space is not usual diffusion but superdiffusion: the radial displacement increases with time faster than that in the normal diffusion. Also, the magnetic field configuration in the hot accretion flow creates outward bulk motion of high-energy particles. This bulk motion is more effective than the diffusive motion for higher energy particles. Our results imply that typical active galactic nuclei that host hot accretion flows can accelerate CRs up to $\epsilon\sim 0.1-10$ PeV., Comment: 17 pages, 16 figures, Figure 4 is corrected, conclusions are unchanged, published in MNRAS

The Magnetorotational Instability (MRI) has long been considered a plausibly ubiquitous mechanism to destabilize otherwise stable Keplerian flows to support radially outward transport of angular momentum. Such an efficient transport process would allow fast accretion in astrophysical objects such as stars and black holes to release copious kinetic energy that powers many of the most luminous sources in the universe. But the standard MRI under a purely vertical magnetic field has heretofore never been directly measured despite numerous efforts over more than a decade. Here we report an unambiguous laboratory demonstration of the spring-mass analogue to the standard MRI by comparing motion of a spring-tethered ball within different rotating flows. The experiment corroborates the theory: efficient outward angular momentum transport manifests only for cases with a weak spring in quasi-Keperian flow. Our experimental method accomplishes this in a new way, thereby connecting solid and fluid mechanics to plasma astrophysics., Comment: 7 pages; 4 figures; accepted by Communications Physics

Many high-energy astrophysical systems contain magnetized collisionless plasmas with relativistic particles, in which turbulence can be driven by an arbitrary mixture of solenoidal and compressive motions. For example, turbulence in hot accretion flows may be driven solenoidally by the magnetorotational instability or compressively by spiral shock waves. It is important to understand the role of the driving mechanism on kinetic turbulence and the associated particle energization. In this work, we compare particle-in-cell simulations of solenoidally driven turbulence with similar simulations of compressively driven turbulence. We focus on plasma that has an initial beta of unity, relativistically hot electrons, and varying ion temperature. Apart from strong large-scale density fluctuations in the compressive case, the turbulence statistics are similar for both drives, and the bulk plasma is described reasonably well by an isothermal equation of state. We find that nonthermal particle acceleration is more efficient when turbulence is driven compressively. In the case of relativistically hot ions, both driving mechanisms ultimately lead to similar power-law particle energy distributions, but over a different duration. In the case of non-relativistic ions, there is significant nonthermal particle acceleration only for compressive driving. Additionally, we find that the electron-to-ion heating ratio is less than unity for both drives, but takes a smaller value for compressive driving. We demonstrate that this additional ion energization is associated with the collisionless damping of large-scale compressive modes via perpendicular electric fields., 29 pages, 28 figures, accepted for publication in ApJ

We study the convective and absolute forms of azimuthal magnetorotational instability (AMRI) in a Taylor-Couette (TC) flow with an imposed azimuthal magnetic field. We show that the domain of the convective AMRI is wider than that of the absolute AMRI. Actually, it is the absolute instability which is the most relevant and important for magnetic TC flow experiments. The absolute AMRI, unlike the convective one, stays in the device, displaying a sustained growth that can be experimentally detected. We also study the global AMRI in a TC flow of finite height using DNS and find that its emerging butterfly-type structure -- a spatio-temporal variation in the form of upward and downward traveling waves -- is in a very good agreement with the linear stability analysis, which indicates the presence of two dominant absolute AMRI modes in the flow giving rise to this global butterfly pattern., Comment: 11 pages, 6 figures, accepted for publication in J. Fluid Mech. Rapids

Extensive simulations of the Princeton Magnetorotational Instability (MRI) Experiment with the Spectral/Finite Element code for Maxwell and Navier-Stokes Equations (SFEMaNS) have been performed to map the MRI-unstable region as a function of inner cylinder angular velocity and applied vertical magnetic field. The angular velocities of the outer cylinder and the end-cap rings follow the inner cylinder in fixed ratios optimized for MRI. We first confirm the exponential growth of the MRI linear phase using idealized conducting vertical boundaries (end caps) rotating differentially with a Taylor-Couette profile. Subsequently, we run a multitude of simulations to scan the experimental parameter space and find that the normalized volume-averaged mean-square radial magnetic field, our main instability indicator, rises significantly where MRI is expected. At various locations, the local radial components of fluid velocity and generated magnetic field are well correlated with the volume-averaged indicator. Based on this correlation, a diagnostic system that will measure the radial magnetic field at several locations on the inner cylinder is proposed as the main comparison between simulation and experiment. A detailed analysis of poloidal mode structures in the SFEMaNS code indicates that MRI, rather than Ekman circulation or Rayleigh instability, dominates the fluid behavior in the region where MRI is expected.

In an earlier paper we showed that the combination of azimuthal magnetic fields and super-rotation in Taylor-Couette flows of conducting fluids can be unstable against non-axisymmetric perturbations if the magnetic Prandtl number of the fluid is $Pm\neq 1$. Here we demonstrate that the addition of a weak axial field component may allow axisymmetric perturbation patterns for $Pm$ of order unity depending on the boundary conditions. The axisymmetric modes only occur for magnetic Mach numbers (of the azimuthal field) of order unity, while higher values are necessary for non-axisymmetric modes. The typical growth time of the instability and the characteristic time scale of the axial migration of the axisymmetric mode are long compared with the rotation period, but short compared with the magnetic diffusion time. The modes travel in the positive or negative $z$-direction along the rotation axis depending on the sign of $B_\phi B_z$. We also demonstrate that the azimuthal components of flow and field perturbations travel in phase if $|B_\phi|\gg |B_z|$, independent of the form of the rotation law. Within a short-wave approximation for thin gaps it is also shown (in an Appendix) that for {\em ideal} fluids the considered helical magnetorotational instability (HMRI) only exists for rotation laws with negative shear., Comment: 18 pages, 10 Figures

High-frequency quasi-periodic oscillations (HFQPOs) observed in the emission of black-hole X-ray binary systems promise insight into strongly curved spacetime. `Discoseismic' oscillations with frequencies set by the intrinsic properties of the central black hole, in particular `trapped inertial waves' (r-modes), offer an attractive explanation for HFQPOs. To produce an observable signature, however, such oscillations must be excited to sufficiently large amplitudes. Turbulence driven by the magnetorotational instability (MRI) fails to provide the necessary amplification, but r-modes may still be excited via interaction with accretion disc warps or eccentricities. We present 3D global hydrodynamic simulations of relativistic accretion discs, which demonstrate for the first time the excitation of trapped inertial waves by an imposed eccentricity in the flow. While the r-modes' saturated state depends on the vertical boundary conditions used in our unstratified, cylindrical framework, their excitation is unambiguous in all runs with eccentricity >0.005 near the ISCO. These simulations provide a proof of concept, demonstrating the robustness of trapped inertial wave excitation in a non-magnetized context. In a companion paper, we explore the competition between this excitation, and damping by magnetohydrodynamic turbulence., 16 pages, 15 figures, accepted in MNRAS. Updated with grammatical and typographical corrections

We consider a differentially rotating flow of an incompressible electrically conducting and viscous fluid subject to an external axial magnetic field and to an azimuthal magnetic field that is allowed to be generated by a combination of an axial electric current external to the fluid and electrical currents in the fluid itself. In this setting we derive an extended version of the celebrated Hain-Lüst differential equation for the radial Lagrangian displacement that incorporates the effects of the axial and azimuthal magnetic fields, differential rotation, viscosity, and electrical resistivity. We apply the Wentzel-Kramers-Brillouin method to the extended Hain-Lüst equation and derive a comprehensive dispersion relation for the local stability analysis of the flow to three-dimensional disturbances. We confirm that in the limit of low magnetic Prandtl numbers, in which the ratio of the viscosity to the magnetic diffusivity is vanishing, the rotating flows with radial distributions of the angular velocity beyond the Liu limit, become unstable subject to a wide variety of the azimuthal magnetic fields, and so is the Keplerian flow. In the analysis of the dispersion relation we find evidence of a new long-wavelength instability which is caught also by the numerical solution of the boundary value problem for a magnetized Taylor-Couette flow.

Previous axisymmetric investigations in two dimensions (2D) have shown that rings and gaps develop naturally in non-ideal magnetohydrodynamic (MHD) disk-wind systems, especially in the presence of ambipolar diffusion (AD). Here we extend these 2D simulations to three dimensions (3D) and find that rings and gaps are still formed in the presence of a moderately strong ambipolar diffusion. The rings and gaps form from the same basic mechanism that was identified in the 2D simulations, namely, the redistribution of the poloidal magnetic flux relative to the disk material as a result of the reconnection of a sharply pinched poloidal magnetic field lines. Thus, the less dense gaps are more strongly magnetized with a large poloidal magnetic field compared to the less magnetized (dense) rings. The rings and gaps start out rather smoothly in 3D simulations that have axisymmetric initial conditions. Non-axisymmetric variations arise spontaneously at later times, but they do not grow to such an extent as to disrupt the rings and gaps. These disk substructures persist to the end of the simulations, lasting up to 3000 orbital periods at the inner edge of the simulated disk. The longevity of the perturbed yet still coherent rings make them attractive sites for trapping large grains that would otherwise be lost to rapid radial migration due to gas drag. As the ambipolar diffusivity decreases, both the disk and the wind become increasingly turbulent, driven by the magnetorotational instability, with tightly-wound spiral arms becoming more prominent in the disk., Comment: 19 pages, 13 figures, accepted by MNRAS

Magnetic fields of planets, stars and galaxies are generated by self-excitation in moving electrically conducting fluids. Once produced, magnetic fields can play an active role in cosmic structure formation by destabilising rotational flows that would be otherwise hydro-dynamically stable. For a long time, both hydromagnetic dynamo action as well as magnetically triggered flow instabilities had been the subject of purely theoretical research. Meanwhile, however, the dynamo effect has been observed in large-scale liquid sodium experiments in Riga, Karlsruhe and Cadarache. In this paper, we summarise the results of liquid metal experiments devoted to the dynamo effect and various magnetic instabilities such as the helical and the azimuthal magnetorotational instability and the Tayler instability. We discuss in detail our plans for a precession-driven dynamo experiment and a large-scale Tayler–Couette experiment using liquid sodium, and on the prospects to observe magnetically triggered instabilities of flows with positive shear.

The nonaxisymmetric azimuthal magnetorotational instability is studied for hydromagnetic Taylor-Couette flows between cylinders of finite electrical conductivity. We find that the magnetic Prandtl number Pm determines whether perfectly conducting or insulating boundary conditions lead to lower Hartmann numbers for the onset of instability. Regardless of the imposed rotation profile, for small Pm the solutions for perfectly conducting cylinders become unstable for weaker magnetic fields than the solutions for insulating cylinders. The critical Hartmann and Reynolds numbers form monotonic functions of the ratio of the electrical conductivities of the cylinders and the fluid, such that a ratio of about 10 provides a very good approximation to perfectly conducting cylinders, and a ratio of about 0.1 a very good approximation to insulating cylinders. These results are of particular relevance for the super-rotating case where the outer cylinder rotates faster than the inner one; in this case the critical onset values are substantially different for perfectly conducting versus insulating boundary conditions. An experimental realization of the super-rotating instability, with liquid sodium as the fluid and cylinders made of copper, would require an electric current of at least 33.5 kAmp running along the central axis., Comment: 18 pages with 15 figures

Geometrically thick tori with constant specific angular momentum have been widely used in the last decades to construct numerical models of accretion flows onto black holes. Such discs are prone to a global non-axisymmetric hydrodynamic instability, known as Papaloizou-Pringle instability (PPI), which can redistribute angular momentum and also lead to an emission of gravitational waves. It is, however, not clear yet how the development of the PPI is affected by the presence of a magnetic field and by the concurrent development of the magnetorotational instability (MRI). We present a numerical analysis using three-dimensional GRMHD simulations of the interplay between the PPI and the MRI considering, for the first time, an analytical magnetized equilibrium solution as initial condition. In the purely hydrodynamic case, the PPI selects as expected the large-scale $m=1$ azimuthal mode as the fastest growing and non-linearly dominant mode. However, when the torus is threaded by a weak toroidal magnetic field, the development of the MRI leads to the suppression of large-scale modes and redistributes power across smaller scales. If the system starts with a significantly excited $m=1$ mode, the PPI can be dominant in a transient phase, before being ultimately quenched by the MRI. Such dynamics may well be important in compact star mergers and tidal disruption events., Comment: 15 pages, 13 figures, submitted to MNRAS

Helical and azimuthal magnetorotational instabilities operate in rotating magnetized flows with relatively steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear, which determine the threshold of modal growth of these instabilities, are continuously connected when some axial electrical current is allowed to pass through the rotating fluid. We investigate the nonmodal dynamics of these instabilities arising from the non-normality of shear flow in the local approximation, generalizing the results of the modal approach. It is demonstrated that moderate transient/nonmodal amplification of both types of magnetorotational instability occurs within the Liu limits, where the system is stable according to modal analysis. We show that for the helical magnetorotational instability this magnetohydrodynamic behavior is closely connected with the nonmodal growth of the underlying purely hydrodynamic problem.

The evolution of a large-scale poloidal magnetic field in an accretion disc is an important problem because it determines the launching of winds and the feasibility of the magnetorotational instability to generate turbulence or channel flows. Recent studies, both semi-analytical calculations and numerical simulations, have highlighted the crucial role non-ideal MHD effects (Ohmic resistivity, Hall drift and ambipolar diffusion), relevant in the protoplanetary disc context, might play in magnetic flux evolution in the disc. We investigate the flux transport in discs through the use of two one-dimensional semi-analytic models in the vertical direction, exploring regimes where different physical source terms and effects dominate. The governing equations for both models are derived by performing an asymptotic expansion in the limit of a thin disc, with the different regimes isolated through setting the relative order of the leading terms between variables. Flux transport rates and vertical structure profiles are calculated for a range of diffusivities and disc magnetisations. We found that Ohmic and ambipolar diffusivities drive radially outward flux transport with an outwardly inclined field. A wind outflow drives inward flux transport, which is significantly enhanced in the presence of Hall drift in the positive polarity case, $\eta_H (\boldsymbol{B}_z \cdot \boldsymbol{\Omega}) > 0$, an effect which has only been briefly noted before. Coupled only with outward inclination, the Hall effect reduces the flux transport given by a background Ohmic and/or ambipolar diffusivity, but drives no flux transport when it is the only non-ideal effect present., Comment: In press. Accepted by MNRAS

We develop a framework for magnetohydrodynamical (MHD) simulations in a local cylindrical shearing box by extending the formulation of the Cartesian shearing box. We construct shearing-periodic conditions at the radial boundaries of a simulation box from the conservation relations of the basic MHD equations, taking into account the explicit radial dependence of physical quantities. We demonstrate quasi-steady mass accretion, which cannot be handled by the standard Cartesian shearing box model, with an ideal MHD simulation in a vertically unstratified cylindrical shearing box up to 200 rotations. In this demonstrative run we set up (i) net vertical magnetic flux, (ii) a locally isothermal equation of state, and (iii) a sub-Keplerian equilibrium rotation, whereas the sound velocity and the initial Alfven velocity have the same radial dependence as that of the Keplerian velocity. Inward mass accretion is induced to balance with the outward angular momentum flux of the MHD turbulence triggered by the magnetorotational instability in a self-consistent manner. We discuss detailed physical properties of the saturated magnetic field, in comparison to the results of a Cartesian shearing box simulation., 27 pages, 10 figures, accepted by PASJ, movie is available at http://ea.c.u-tokyo.ac.jp/astro/Members/stakeru/research/cylshbx/

The gas temperature in protoplanetary disks (PPDs) is determined by a combination of irradiation heating and accretion heating, with the latter conventionally attributed to turbulent dissipation. However, recent studies have suggested that the inner disk (a few AU) is largely laminar, with accretion primarily driven by magnetized disk winds, as a result of nonideal magnetohydrodynamic (MHD) effects from weakly ionized gas, suggesting an alternative heating mechanism by Joule dissipation. We perform local stratified MHD simulations including all three nonideal MHD effects (ohmic, Hall, and ambipolar diffusion) and investigate the role of Joule heating and the resulting disk vertical temperature profiles. We find that in the inner disk, as ohmic and ambipolar diffusion strongly suppress electrical current around the midplane, Joule heating primarily occurs at several scale heights above the midplane, making the midplane temperature much lower than that with the conventional viscous heating model. Including the Hall effect, Joule heating is enhanced/reduced when the magnetic fields threading the disks are aligned/anti-aligned with the disk rotation, but it is overall ineffective. Our results further suggest that the midplane temperature in the inner PPDs is almost entirely determined by irradiation heating, unless viscous heating can trigger thermal ionization in the disk innermost region to self-sustain magnetorotational instability turbulence., Comment: 23 pages, 11 figures; matches the published version

The magnetorotational instability (MRI) occurs when a weak magnetic field destabilises a rotating, electrically conducting fluid with inwardly increasing angular velocity. The MRI is essential to astrophysical disk theory where the shear is typically Keplerian. Internal shear layers in stars may also be MRI unstable, and they take a wide range of profiles, including near-critical. We show that the fastest growing modes of an ideal magnetofluid are three-dimensional provided the shear rate, $S$, is near the two-dimensional onset value, $S_c$. For a Keplerian shear, three-dimensional modes are unstable above $S\approx0.10S_c$, and dominate the two-dimensional modes until $S\approx2.05S_{c}$. These three-dimensional modes dominate for shear profiles relevant to stars and at magnetic Prandtl numbers relevant to liquid-metal laboratory experiments. Significant numbers of rapidly growing three-dimensional modes remain well past $2.05S_{c}$. These finding are significant in three ways. First, weakly nonlinear theory suggests that the MRI saturates by pushing the shear rate to its critical value. This can happen for systems, like stars and laboratory experiments, that can rearrange their angular velocity profiles. Second, the non-normal character and large transient growth of MRI modes should be important whenever three-dimensionality exists. Finally, three-dimensional growth suggests direct dynamo action driven from the linear instability., Comment: 15 pages, 5 figures. Proceedings of the Royal Society A, accepted

In this paper we studied the weakly nonlinear stage of stationary convective instability in a nonuniformly rotating layer of an electrically conductive fluid in an axial uniform magnetic field under the influence of: a) temperature modulation of the layer boundaries; b) gravitational modulation; c) modulation of the magnetic field; d) modulation of the angular velocity of rotation. As a result of applying the method of perturbation theory for the small parameter of supercriticality of the stationary Rayleigh number nonlinear non-autonomous Ginzburg-Landau equations for the above types of modulation were obtaned. By utilizing the solution of the Ginzburg-Landau equation, we determined the dynamics of unsteady heat transfer for various types of modulation of external fields and for different profiles of the angular velocity of the rotation of electrically conductive fluid., Comment: 44 pages, 28 figures

We analyze the sustaining mechanism of nonlinear perturbations/turbulence in spectrally stable smooth shear flows. The essence of the sustenance is a subtle interplay of linear transient growth of Fourier harmonics and nonlinear processes. In spectrally stable shear flows, the transient growth of perturbations is strongly anisotropic in spectral (k-)space. This, in turn, leads to anisotropy of nonlinear processes ink-space and, as a result, the main (new) nonlinear process appears to be not a direct/inverse, but rather a transverse/angular redistribution of harmonics in Fourier space referred to as the nonlinear transverse cascade. It is demonstrated that nonlinear state is sustained owing to the interplay of the linear nonmodal growth and the transverse cascade. The possibility of such course of events has been described ink-space byG. Chagelishvili, J.-P. Zahn, A. Tevzadze and J. Lominadze, A&A, 402, 401 (2003)that reliably exemplifies the well-known bypass scenario of subcritical turbulence in spectrally stable shear flows. We present selected results of the simulations performed in different (HD and MHD; 2D and 3D; plane and Keplerian) shear flows to demonstrate the transverse cascade in action.

In this paper the stability of the non-uniformly rotating cylindrical plasma in the axial uniform magnetic field with the vertical temperature gradient is investigated. In the approximation of geometrical optics a dispersion equation for small axisymmetric perturbations is obtained with the effects of viscosity, ohmic and heat conductive dissipation taken into account. The stability criteria for azimuthal plasma flows are obtained in the presence of the vertical temperature gradient and the constant magnetic field. The Rayleigh-Benard problem for stationary convection in the non-uniformly rotating layer of the electrically conducting fluid in the axial uniform magnetic field is considered. In the linear theory of stationary convection the critical value of the Rayleigh number subject to the profile of the inhomogeneous rotation (Rossby number) is obtained. It is shown that the negative values of the Rossby number have a destabilizing effect, since the critical Rayleigh number becomes smaller, than in the case of the uniform rotation , or of the rotation with positive Rossby numbers . To describe the nonlinear convective phenomena the local Cartesian coordinate system was used, where the inhomogeneous rotation of the fluid layer was represented as the rotation with a constant angular velocity and azimuthal shear with linear dependence on the coordinate. As a result of applying the method of perturbation theory for the small parameter of supercriticality of the stationary Rayleigh number a nonlinear Ginzburg-Landau equation was obtaned. This equation describes the evolution of the finite amplitude of perturbations by utilizing the solution of the Ginzburg-Landau equation. It is shown that the weakly nonlinear convection based on the equations of the six-mode Lorentz model transforms into the identical Ginzburg-Landau equation. By utilizing the solution of the Ginzburg-Landau equation, we determined the dynamics of unsteady heat transfer for various profiles of the angular velocity of the rotation of electrically conductive fluid.

We solve the two-dimensional hydrodynamic equations of hot accretion flow in the presence of the thermal conduction. The flow is assumed to be in steady-state and axisymmetric, and self-similar approximation is adopted in the radial direction. In this hydrodynamic study, we consider the viscous stress tensor to mimic the effects of the magnetorotational instability for driving angular momentum. We impose the physical boundary conditions at both the rotation axis and the equatorial plane and obtain the solutions in the full $ r-\theta $ space. We have found that thermal conduction is indispensable term for investigating the inflow-wind structure of the hot accretion flows with very low mass accretion rates. One of the most interesting results here is that the disc is convectively stable in hot accretion mode and in the presence of the thermal conduction. Furthermore, the properties of wind and also its driving mechanisms are studied. Our analytical results are consistent with previous numerical simulations of hot accretion flow., Comment: 20 pages, 10 figures, accepted for publication in ApJ

We study the dynamo generation (exponential growth) of large scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal ($x$-$y$) averaging, we demonstrate the presence of large scale fields when either horizontal or vertical ($y$-$z$) averaging is employed. By computing planar averaged fields and power spectra, we find large scale dynamo action in the early MRI growth phase---a previously unidentified feature. Fast growing horizontal low modes and fiducial vertical modes over a narrow range of wave numbers amplify these planar averaged fields in the MRI growth phase, before turbulence sets in. The large scale field growth requires linear fluctuations but not nonlinear turbulence (as defined by mode-mode coupling) and grows as a direct global mode of the MRI. Only by vertical averaging, can it be shown that the growth of horizontal low wavenumber MRI modes directly feed-back to the initial vertical field providing a clue as to why the large scale vertical field sustains against turbulent diffusion in the saturation regime. We compute the terms in the planar averaged mean field equations to identify the individual contributions to large scale field growth for both vertical and horizontal averaging. The large scale fields obtained from such vertical averaging are found to compare well with global cylindrical simulations and quasilinear analytical analysis from a previous study by Ebrahimi \& Blackman. We discuss the potential implications of these new results for understanding large scale MRI dynamo saturation and turbulence., 12 pages, 9 figues, submitted to MNRAS, comments are welcome

The dynamics of solid bodies in protoplanetary disks are subject to the properties of any underlying gas turbulence. Turbulence driven by disk self-gravity shows features distinct from those driven by the magnetorotational instability (MRI). We study the dynamics of solids in gravito-turbulent disks with two-dimensional (in the disk plane), hybrid (particle and gas) simulations. Gravito-turbulent disks can exhibit stronger gravitational stirring than MRI-active disks, resulting in greater radial diffusion and larger eccentricities and relative speeds for large particles (those with dimensionless stopping times $t_{stop} \Omega > 1$, where $\Omega$ is the orbital frequency). The agglomeration of large particles into planetesimals by pairwise collisions is therefore disfavored in gravito-turbulent disks. However, the relative speeds of intermediate-size particles $t_{stop} \Omega \sim 1$ are significantly reduced as such particles are collected by gas drag and gas gravity into coherent filament-like structures with densities high enough to trigger gravitational collapse. First-generation planetesimals may form via gravitational instability of dust in marginally gravitationally unstable gas disks., Comment: MNRAS accepted

Turbulence in the protoplanetary disks induces collisions between dust grains, and thus facilitates grain growth. We investigate the two fundamental assumptions of the turbulence in obtaining grain collisional velocities -- the kinetic energy spectrum and the turbulence autocorrelation time -- in the context of the turbulence generated by the magneto-rotational instability (MRI). We carry out numerical simulations of the MRI as well as driven turbulence, for a range of physical and numerical parameters. We find that the convergence of the turbulence $\alpha$-parameter does not necessarily imply the convergence of the energy spectrum. The MRI turbulence is largely solenoidal, for which we observe a persistent kinetic energy spectrum of $k^{-4/3}$. The same is obtained for solenoidal driven turbulence with and without magnetic field, over more than 1 dex near the dissipation scale. This power-law slope appears to be converged in terms of numerical resolution, and to be due to the bottleneck effect. The kinetic energy in the MRI turbulence peaks at the fastest growing mode of the MRI. In contrast, the magnetic energy peaks at the dissipation scale. The magnetic energy spectrum in the MRI turbulence does not show a clear power-law range, and is almost constant over approximately 1 dex near the dissipation scale. The turbulence autocorrelation time is nearly constant at large scales, limited by the shearing timescale, and shows a power-law drop close to $k^{-1}$ at small scales, with a slope steeper than that of the eddy crossing time. The deviation from the standard picture of the Kolmogorov turbulence with the injection scale at the disk scale height can potentially have a significant impact on the grain collisional velocities., Comment: Accepted by ApJ

Stability and nonlinear evolution of rotating magnetohydrodynamic flows in the Princeton magnetorotational instability (MRI) experiment are examined using three-dimensional non-axisymmetric simulations. In particular, the effect of axial boundary conductivity on a free Stewartson-Shercliff layer (SSL) is numerically investigated using the spectral finite-element Maxwell and Navier Stokes (SFEMaNS) code. The free SSL is established by a sufficiently strong magnetic field imposed axially across the differentially rotating fluid with two rotating rings enforcing the boundary conditions. Numerical simulations show that the response of the bulk fluid flow is vastly different in the two different cases of insulating and conducting end caps. We find that, for the insulating end caps, there is a transition from stability to instability of a Kelvin-Helmholtz-like mode that saturates at an azimuthal mode number m=1, whereas for the conducting end caps, the reinforced coupling between the magnetic field and the bulk fluid generates a strong radially localized shear in the azimuthal velocity resulting in axisymmetric Rayleigh-like modes even at reduced thresholds for the axial magnetic field. For reference, three-dimensional nonaxisymmetric simulations have also been performed in the MRI unstable regime to compare the modal structures.

We study local instabilities of a differentially rotating viscous flow of electrically conducting incompressible fluid subject to an external azimuthal magnetic field. A hydrodynamically stable flow can be destabilized by the magnetic field both in an ideal and a viscous and resistive system giving rise to the azimuthal magnetorotational instability. A special solution to the equations of ideal magnetohydrodynamics characterized by the constant total pressure, the fluid velocity parallel to the direction of the magnetic field, and by the magnetic and kinetic energies that are finite and equal—the Chandrasekhar equipartition solution—is marginally stable in the absence of viscosity and resistivity. Performing a local stability analysis, we find the conditions under which the azimuthal magnetorotational instability can be interpreted as a dissipation-induced instability of the Chandrasekhar equipartition solution.

We investigate numerically the flow of an electrically conducting fluid in a rapidly rotating spherical shell where the inner boundary spins slightly faster than the outer one. The magnetic field evolves self-consistently from an initial dipolar configuration of weak amplitude, and a toroidal field is produced by winding this poloidal field through the internal differential rotation. First, we characterize the axisymmetric field solutions obtained at long times when the Lorentz force is negligible and the flow follows the steady, purely hydrodynamical solution. We then examine the stability of these solutions, focusing on the regime of large magnetic Reynolds numbers where the field is dominantly toroidal. When the ratio of the azimuthal Alfven frequency to the rotation frequency exceeds a certain value, a nonaxisymmetric instability develops. We show that the instability properties are compatible with those expected for the magnetorotational instability. Finally, we compare the instability properties with predictions obtained from a local linear stability analysis. The linear analysis agrees well with the numerical simulation results, except in a number of cases where the discrepancies are attributed to shearing effects on the unstable modes.

We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu limit, in that it works instead for a wide range of the positive shear when ${\rm (i)}$ a combination of axial/poloidal and azimuthal/toroidal magnetic fields is applied and ${\rm (ii)}$ the magnetic Prandtl number is not too close to unity. We study this instability first with radially local WKB analysis and then confirm its existence using a global stability analysis of the magnetized flow between two rotating cylinders with conducting or insulating boundaries. From an experimental point of view, we also demonstrate the presence of the new instability in a magnetized viscous and resistive Taylor-Couette flow with positive shear for such values of the flow parameters, which can be realized in upcoming experiments at the DRESDYN facility. Finally, this instability might have implications for the dynamics of the equatorial parts of the solar tachocline and dynamo action there, since the above two necessary conditions for the instability to take place are satisfied in this region. Our global stability calculations for the tachocline-like configuration, representing a thin rotating cylindrical layer with the appropriate boundary conditions -- conducting inner and insulating outer cylinders -- and the values of the flow parameters, indicate that it can indeed arise in this case with a characteristic growth time comparable to the solar cycle period., 15 pages, 13 figures, accepted for publication in Phys. Rev. Fluids

The prevalence and consequences of convection perpendicular to the plane of accretion discs have been discussed for several decades. Recent simulations combining convection and the magnetorotational instability have given fresh impetus to the debate, as the interplay of the two processes can enhance angular momentum transport, at least in the optically thick outburst stage of dwarf novae. In this paper we seek to isolate and understand the most generic features of disc convection, and so undertake its study in purely hydrodynamical models. First, we investigate the linear phase of the instability, obtaining estimates of the growth rates both semi-analytically, using one-dimensional spectral computations, as well as analytically, using WKBJ methods. Next we perform three-dimensional, vertically stratified, shearing box simulations with the conservative, finite-volume code PLUTO, both with and without explicit diffusion coefficients. We find that hydrodynamic convection can, in general, drive outward angular momentum transport, a result that we confirm with ATHENA, an alternative finite-volume code. Moreover, we establish that the sign of the angular momentum flux is sensitive to the diffusivity of the numerical scheme. Finally, in sustained convection, whereby the system is continuously forced to an unstable state, we observe the formation of various coherent structures, including large- scale and oscillatory convective cells, zonal flows, and small vortices., Accepted for publication in MNRAS (20 pages, 16 figures, 4 tables)

Recent numerical simulations of magnetized accretion discs show that the radial-azimuthal component of the stress tensor due to the magnetorotational instability (MRI) is well represented by a power-law function of the gas pressure rather than a linear relation which has been used in most of the accretion disc studies. The exponent of this power-law function which depends on the net flux of the imposed magnetic field is reported in the range between zero and unity. However, the physical consequences of this power-law stress-pressure relation within the framework of the standard disc model have not been explored so far. In this study, the structure of an accretion disc with a power-law stress-pressure relation is studied using analytical solutions in the steady-state and time-dependent cases. The derived solutions are applicable to different accreting systems, and as an illustrative example, we explore structure of protoplanetary discs using these solutions. We show that the slopes of the radial surface density and temperature distributions become steeper with decreasing the stress exponent. However, if the disc opacity is dominated by icy grains and value of the stress exponent is less than about $0.5$, the surface density and temperature profiles become so steep that make them unreliable. We also obtain analytical solutions for the protoplanetary discs which are irradiated by the host star. Using these solutions, we find that the effect of the irradiation becomes more significant with decreasing the stress exponent., Accepted for publication in MNRAS

We studied dynamical balances in magnetorotational instability (MRI) turbulence with a net vertical field in the shearing box model of disks. Analyzing the turbulence dynamics in Fourier (${\bf k}$-)space, we identified three types of active modes that define turbulence characteristics. These modes have lengths similar to the box size, i.e., lie in the small wavenumber region in Fourier space labeled the vital area and are: (i) the channel mode - uniform in the disk plane with the smallest vertical wavenumber,(ii) the zonal flow mode - azimuthally and vertically uniform with the smallest radial wavenumber and (iii) the rest modes. The rest modes comprise those harmonics in the vital area whose energies reach more than $50 \%$ of the maximum spectral energy. The rest modes individually are not so significant compared to the channel and zonal flow modes, however, the combined action of their multitude is dominant over these two modes. These three mode types are governed by interplay of the linear and nonlinear processes, leading to their interdependent dynamics. The linear processes consist in disk flow nonmodality-modified classical MRI with a net vertical field. The main nonlinear process is transfer of modes over wavevector angles in Fourier space - the transverse cascade. The channel mode exhibits episodic bursts supplied by linear MRI growth, while the nonlinear processes mostly oppose this, draining the channel energy and redistributing it to the rest modes. As for the zonal flow, it does not have a linear source and is fed by nonlinear interactions of the rest modes., 28 pages, 16 figures, published in ApJ

The streaming instability and Kelvin--Helmholtz instability are considered the two major sources causing clumping of dust particles and turbulence in the dust layer of a protoplanetary disk as long as we consider the dead zone where the magneto-rotational instability does not grow. Extensive numerical simulations have been carried out in order to elucidate the condition for the development of particle clumping caused by the streaming instability. In this paper, a set of two parameters suitable for classifying the numerical results is proposed. One is the Stokes number that has been employed in previous works and the other is the dust particle column density that is nondimensionalized using the gas density in the midplane, Keplerian angular velocity, and difference between the Keplerian and gaseous orbital velocities. The magnitude of dust clumping is a measure of the behavior of the dust layer. Using three-dimensional numerical simulations of dust particles and gas based on Athena code v. 4.2, it is confirmed that the magnitude of dust clumping for two disk models are similar if the corresponding sets of values of the two parameters are identical to each other, even if the values of the metallicity (i.e., the ratio of the columns density of the dust particles to that of the gas) are different.

Magnetorotational instability (MRI) is an instability that is responsible for accretion, the phenomenon observed in astrophysical disks, for example around black holes. MRI can be modeled using the equations of MHD. A typical linear analysis in the past made use of the small Prandtl number approximation, which results in dropping one term from these equations. Previously, one of the authors (Herron) showed that the small magnetic Prandtl number approximation suppresses MRI in axisymmetric viscous resistive magnetized Taylor–Couette flow when one has no-slip velocity boundary conditions, and insulating magnetic boundary conditions. We follow up here with a proof that MRI is still suppressed with perfectly conducting magnetic boundary conditions on the cylinders.

We present an approach to turbulence closure based on mixing length theory with three-dimensional fluctuations against a two-dimensional background. This model is intended to be rapidly computable for implementation in stellar evolution software and to capture a wide range of relevant phenomena with just a single free parameter, namely the mixing length. We incorporate magnetic, rotational, baroclinic and buoyancy effects exactly within the formalism of linear growth theories with nonlinear decay. We treat differential rotation effects perturbatively in the corotating frame using a novel controlled approximation which matches the time evolution of the reference frame to arbitrary order. We then implement this model in an efficient open source code and discuss the resulting turbulent stresses and transport coefficients. We demonstrate that this model exhibits convective, baroclinic and shear instabilities as well as the magnetorotational instability (MRI). It also exhibits non-linear saturation behaviour, and we use this to extract the asymptotic scaling of various transport coefficients in physically interesting limits., Comment: 17 pages, 20 figures, published in MNRAS. Figures 9-19 and C1 and associated text have been updated relative to MNRAS version to reflect minor changes due to corrections to the numerical code. Typographic errors have been corrected in equations 8, 10 and 11 and Figure 1 caption. Figure 3 legend updated to be clearer, and colors have been coordinated between Figures 1-3

Magnetic fields of planets, stars, and galaxies are generated by self-excitation in moving electrically conducting fluids. Once produced, magnetic fields can play an active role in cosmic structure formation by destabilizing rotational flows that would be otherwise hydrodynamically stable. For a long time, both hydromagnetic dynamo action and magnetically triggered flow instabilities had been the subject of purely theoretical research. Meanwhile, however, the dynamo effect has been observed in large-scale liquid sodium experiments in Riga, Karlsruhe, and Cadarache. In this chapter, we summarize the results of some smaller liquid metal experiments devoted to various magnetic instabilities, such as the helical and the azimuthal magnetorotational instability, the Tayler instability, and the different instabilities that appear in a magnetized spherical Couette flow. We conclude with an outlook on a large scale Tayler-Couette experiment using liquid sodium, and on the prospects to observe magnetically triggered instabilities of flows with positive shear.

Space and astrophysical plasmas are frequently found in the regime of differential rotation, where the presence of a magnetic field can result in the magnetorotational instability, directly responsible for important phenomena such as turbulent angular momentum transport in accretion disks. In the absence of an imposed magnetic field, a nonlinear dynamo is necessary for this transport mechanism to take place. In protoplanetary disks there are regions with high density and very low temperatures, which are two necessary conditions for the Hall effect to operate, affecting the development of the dynamo and the associated turbulence. In this work we perform local magnetohydrodynamic (MHD) simulations to study transition to weak turbulence in Keplerian shear flows with Hall effect. The Hall effect is shown to lead the system to long-lived turbulent transients whose decay time follows an exponential dependence on the magnetic Reynolds number and the Hall parameter., Comment: 6 pages, 7 figures Submitted to EPL

A laboratory experiment on studying magnetorotational instability (MRI) has been numerically simulated. The experiment involves driving liquid sodium placed in a vertical magnetic field in rotation by passing a radial electric current between the surfaces of coaxial cylinders. The numerical simulation yields the parameters of the future experimental device making it possible to develop MRI. These parameters include the size of the working chamber, the value of the radial current, and the strength of the initial magnetic field.

The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars., 22 pages, 15 figures, accepted for publication in MNRAS

We investigate the evolution of hydromagnetic perturbations in a small section of accretion disks. It is known that molecular viscosity is negligible in accretion disks. Hence, it has been argued that Magnetorotational Instability (MRI) is responsible for transporting matter in the presence of weak magnetic field. However, there are some shortcomings, which question effectiveness of MRI. Now the question arises, whether other hydromagnetic effects, e.g. transient growth (TG), can play an important role to bring nonlinearity in the system, even at weak magnetic fields. Otherwise, whether MRI or TG, which is primarily responsible to reveal nonlinearity to make the flow turbulent? Our results prove explicitly that the flows with high Reynolds number (Re), which is the case of realistic astrophysical accretion disks, exhibit nonlinearity by best TG of perturbation modes faster than that by best modes producing MRI. For a fixed wavevector, MRI dominates over transient effects, only at low Re, lower than its value expected to be in astrophysical accretion disks, and low magnetic fields. This seriously questions (overall) persuasiveness of MRI in astrophysical accretion disks., 6 pages including 1 figure; to appear in the Proceedings of 14th Marcel Grossman Meeting (MG14), Rome, Italy, 12-18 July 2015; based on the talk given in the meeting in the session AC2

We examine the properties of strongly magnetized accretion discs in a global framework, with particular focus on the evolution of magnetohydrodynamic instabilities such as the magnetorotational instability (MRI). Work by Pessah and Psaltis showed that MRI is stabilized beyond a critical toroidal field in compressible, differentially rotating flows and, also, reported the appearance of two new instabilities beyond this field. Their results stemmed from considering geometric curvature effects due to the suprathermal background toroidal field, which had been previously ignored in weak-field studies. However, their calculations were performed under the local approximation, which poses the danger of introducing spurious behavior due to the introduction of global geometric terms in an otherwise local framework. In order to avoid this, we perform a global eigenvalue analysis of the linearized MHD equations in cylindrical geometry. We confirm that MRI indeed tends to be highly suppressed when the background toroidal field attains the Pessah-Psaltis limit. We also observe the appearance of two new instabilities that emerge in the presence of highly suprathermal toroidal fields. These results were additionally verified using numerical simulations in PLUTO. There are, however, certain differences between the the local and global results, especially in the vertical wavenumber occupancies of the various instabilities, which we discuss in detail. We also study the global eigenfunctions of the most unstable modes in the suprathermal regime, which are inaccessible in the local analysis. Overall, our findings emphasize the necessity of a global treatment for accurately modeling strongly magnetized accretion discs., 23 pages, 7 Figures and 1 Table. Accepted for publication in MNRAS

The project A2 of the LIMTECH Alliance aimed at a better understanding of those magnetohydrodynamic instabilities that are relevant for the generation and the action of cosmic magnetic fields. These comprise the hydromagnetic dynamo effect and various magnetically triggered flow instabilities, such as the magnetorotational instability and the Tayler instability. The project was intended to support the experimental capabilities to become available in the framework of the DREsden Sodium facility for DYNamo and thermohydraulic studies (DRESDYN). An associated starting grant was focused on the dimensioning of a liquid metal experiment on the newly found magnetic destabilization of rotating flows with positive shear. In this paper, the main results of these two projects are summarized., 17 pages, 10 figures

Magnetic dynamo action caused by the magnetorotational instability is studied in the shearing-box approximation with no imposed net magnetic flux. Consistent with recent studies, the dynamo action is found to be sensitive to the aspect ratio of the box: it is much easier to obtain in tall boxes (stretched in the direction normal to the disk plane) than in long boxes (stretched in the radial direction). Our direct numerical simulations indicate that the dynamo is possible in both cases, given a large enough magnetic Reynolds number. To explain the relatively larger effort required to obtain the dynamo action in a long box, we propose that the turbulent eddies caused by the instability most efficiently fold and mix the magnetic field lines in the radial direction. As a result, in the long box the scale of the generated strong azimuthal (stream-wise directed) magnetic field is always comparable to the scale of the turbulent eddies. In contrast, in the tall box the azimuthal magnetic flux spreads in the vertical direction over a distance exceeding the scale of the turbulent eddies. As a result, different vertical sections of the tall box are permeated by large-scale nonzero azimuthal magnetic fluxes, facilitating the instability. In agreement with this picture, the cases when the dynamo is efficient are characterized by a strong intermittency of the local azimuthal magnetic fluxes., Comment: 7 pages, 9 figures, 1 table