No Arabic abstract
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.
Magnetorotational instability (MRI) is one of the fundamental processes in astrophysics, driving angular momentum transport and mass accretion in a wide variety of cosmic objects. Despite much theoretical/numerical and experimental efforts over the last decades, its saturation mechanism and amplitude, which sets the angular momentum transport rate, remains not well understood, especially in the limit of high resistivity, or small magnetic Prandtl numbers typical to interiors (dead zones) of protoplanetary disks, liquid cores of planets and liquid metals in laboratory. Using direct numerical simulations, in this paper we investigate the nonlinear development and saturation properties of the helical magnetorotational instability (HMRI) -- a relative of the standard MRI -- in a magnetized Taylor-Couette flow at very low magnetic Prandtl number (correspondingly at low magnetic Reynolds number) relevant to liquid metals. For simplicity, the ratio of azimuthal field to axial field is kept fixed. From the linear theory of HMRI, it is known that the Elsasser number, or interaction parameter determines its growth rate and plays a special role in the dynamics. We show that this parameter is also important in the nonlinear problem. By increasing its value, a sudden transition from weakly nonlinear, where the system is slightly above the linear stability threshold, to strongly nonlinear, or turbulent regime occurs. We calculate the azimuthal and axial energy spectra corresponding to these two regimes and show that they differ qualitatively. Remarkably, the nonlinear state remains in all cases nearly axisymmetric suggesting that HMRI turbulence is quasi two-dimensional in nature. Although the contribution of non-axisymmetric modes increases moderately with the Elsasser number, their total energy remains much smaller than that of the axisymmetric ones.
The helical magnetorotational instability is known to work for resistive rotational flows with comparably steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear are continuously connected when some axial electrical current is allowed to flow through the rotating fluid. Using a local approximation we demonstrate that the magnetohydrodynamic behavior of this dissipation-induced instability is intimately connected with the nonmodal growth and the pseudospectrum of the underlying purely hydrodynamic problem.
We consider the effect of stratification on systematic, large-scale flows generated in anelastic convection. We present results from three-dimensional numerical simulations of convection in a rotating plane layer in which the angle between the axis of rotation and gravity is allowed to vary. This model is representative of different latitudes of a spherical body. We consider two distinct parameter regimes: (i) weakly rotating and (ii) rapidly rotating. In each case, we examine the effect of stratification on the flow structure and heat transport properties focussing on the difference between Boussinesq and anelastic convection. Furthermore, we show that regimes (i) and (ii) generate very different large-scale flows and we investigate the role stratification has in modifying these flows. The stratified flows possess a net helicity not present in the Boussinesq cases which we suggest, when combined with the self-generated shear flows, could be important for dynamo action.
Origin of hydrodynamic turbulence in rotating shear flows is investigated. The particular emphasis is the flows whose angular velocity decreases but specific angular momentum increases with increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to explain observed data. Such a mismatch between the linear theory and observations/experiments is more severe when any hydromagnetic/magnetohydrodynamic instability and then the corresponding turbulence therein is ruled out. The present work explores the effect of stochastic noise on such hydrodynamic flows. We essentially concentrate on a small section of such a flow which is nothing but a plane shear flow supplemented by the Coriolis effect. This also mimics a small section of an astrophysical accretion disk. It is found that such stochastically driven flows exhibit large temporal and spatial correlations of perturbation velocities, and hence large energy dissipations of perturbation, which presumably generate instability. A range of angular velocity (Omega) profiles of background flow, starting from that of constant specific angular momentum (lambda = Omega r^2 ; r being the radial coordinate) to that of constant circular velocity (v_phi = Omega r), is explored. However, all the background angular velocities exhibit identical growth and roughness exponents of perturbations, revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear flows. This work, to the best of our knowledge, is the first attempt to understand origin of instability and turbulence in the three-dimensional Rayleigh stable rotating shear flows by introducing additive noise to the underlying linearized governing equations. This has important implications to resolve the turbulence problem in astrophysical hydrodynamic flows such as accretion disks.
We investigate the asymptotic properties of axisymmetric inertial modes propagating in a spherical shell when viscosity tends to zero. We identify three kinds of eigenmodes whose eigenvalues follow very different laws as the Ekman number $E$ becomes very small. First are modes associated with attractors of characteristics that are made of thin shear layers closely following the periodic orbit traced by the characteristic attractor. Second are modes made of shear layers that connect the critical latitude singularities of the two hemispheres of the inner boundary of the spherical shell. Third are quasi-regular modes associated with the frequency of neutral periodic orbits of characteristics. We thoroughly analyse a subset of attractor modes for which numerical solutions point to an asymptotic law governing the eigenvalues. We show that three length scales proportional to $E^{1/6}$, $E^{1/4}$ and $E^{1/3}$ control the shape of the shear layers that are associated with these modes. These scales point out the key role of the small parameter $E^{1/12}$ in these oscillatory flows. With a simplified model of the viscous Poincare equation, we can give an approximate analytical formula that reproduces the velocity field in such shear layers. Finally, we also present an analysis of the quasi-regular modes whose frequencies are close to $sin(pi/4)$ and explain why a fluid inside a spherical shell cannot respond to any periodic forcing at this frequency when viscosity vanishes.