Azimuthal velocity profiles in Rayleigh-stable Taylor-Couette flow and implied axial angular momentum transport

We present azimuthal velocity profiles measured in a Taylor-Couette apparatus, which has been used as a model of stellar and planetary accretion disks. The apparatus has a cylinder radius ratio of $\eta = 0.716$, an aspect-ratio of $\Gamma = 11.74$, and the plates closing the cylinders in the axial direction are attached to the outer cylinder. We investigate angular momentum transport and Ekman pumping in the Rayleigh-stable regime. The regime is linearly stable and is characterized by radially increasing specific angular momentum. We present several Rayleigh-stable profiles for shear Reynolds numbers $Re_S \sim O(10^5) \,$, both for $\Omega_i > \Omega_o > 0$ (quasi-Keplerian regime) and $\Omega_o > \Omega_i > 0$ (sub-rotating regime) where $\Omega_{i,o}$ is the inner/outer cylinder rotation rate. None of the velocity profiles matches the non-vortical laminar Taylor-Couette profile. The deviation from that profile increased as solid-body rotation is approached at fixed $Re_S$. Flow super-rotation, an angular velocity greater than that of both cylinders, is observed in the sub-rotating regime. The velocity profiles give lower bounds for the torques required to rotate the inner cylinder that were larger than the torques for the case of laminar Taylor-Couette flow. The quasi-Keplerian profiles are composed of a well mixed inner region, having approximately constant angular momentum, connected to an outer region in solid-body rotation with the outer cylinder and attached axial boundaries. These regions suggest that the angular momentum is transported axially to the axial boundaries. Therefore, Taylor-Couette flow with closing plates attached to the outer cylinder is an imperfect model for accretion disk flows, especially with regard to their stability.Comment: 22 pages, 10 figures, 2 tables, under consideration for publication in Journal of Fluid Mechanics (JFM

On the Phenomenology of Hydrodynamic Shear Turbulence

The question of a purely hydrodynamic origin of turbulence in accretion disks is reexamined, on the basis of a large body of experimental and numerical evidence on various subcritical (i.e., linearly stable) hydrodynamic flows. One of the main points of this paper is that the length scale and velocity fluctuation amplitude which are characteristic of turbulent transport in these flows scale like $Re_m^{-1/2}$, where $Re_m$ is the minimal Reynolds number for the onset of fully developed turbulence. From this scaling, a simple explanation of the dependence of $Re_m$ with relative gap width in subcritical Couette-Taylor flows is developed. It is also argued that flows in the shearing sheet limit should be turbulent, and that the lack of turbulence in all such simulations performed to date is most likely due to a lack of resolution, as a consequence of the effect of the Coriolis force on the large scale fluctuations of turbulent flows. These results imply that accretion flows should be turbulent through hydrodynamic processes. If this is the case, the Shakura-Sunyaev $\alpha$ parameter is constrained to lie in the range $10^{-3}-10^{-1}$ in accretion disks, depending on unknown features of the mechanism which sustains turbulence. Whether the hydrodynamic source of turbulence is more efficient than the MHD one where present is an open question.Comment: 31 pages, 3 figures. Accepted for publication in Ap

Magnetorotational Instability in Liquid Metal Couette Flow

Despite the importance of the magnetorotational instability (MRI) as a fundamental mechanism for angular momentum transport in magnetized accretion disks, it has yet to be demonstrated in the laboratory. A liquid sodium alpha-omega dynamo experiment at the New Mexico Institute of Mining and Technology provides an ideal environment to study the MRI in a rotating metal annulus (Couette flow). A local stability analysis is performed as a function of shear, magnetic field strength, magnetic Reynolds number, and turbulent Prandtl number. The later takes into account the minimum turbulence induced by the formation of an Ekman layer against the rigidly rotating end walls of a cylindrical vessel. Stability conditions are presented and unstable conditions for the sodium experiment are compared with another proposed MRI experiment with liquid gallium. Due to the relatively large magnetic Reynolds number achievable in the sodium experiment, it should be possible to observe the excitation of the MRI for a wide range of wavenumbers and further to observe the transition to the turbulent state.Comment: 12 pages, 22 figures, 1 table. To appear in the Astrophysical Journa

Numerical and Experimental Investigation of Circulation in Short Cylinders

In preparation for an experimental study of magnetorotational instability (MRI) in liquid metal, we explore Couette flows having height comparable to the gap between cylinders, centrifugally stable rotation, and high Reynolds number. Experiments in water are compared with numerical simulations. Simulations show that endcaps corotating with the outer cylinder drive a strong poloidal circulation that redistributes angular momentum. Predicted azimuthal flow profiles agree well with experimental measurements. Spin-down times scale with Reynolds number as expected for laminar Ekman circulation; extrapolation from two-dimensional simulations at $Re\le 3200$ agrees remarkably well with experiment at $Re\sim 10^6$. This suggests that turbulence does not dominate the effective viscosity. Further detailed numerical studies reveal a strong radially inward flow near both endcaps. After turning vertically along the inner cylinder, these flows converge at the midplane and depart the boundary in a radial jet. To minimize this circulation in the MRI experiment, endcaps consisting of multiple, differentially rotating rings are proposed. Simulations predict that an adequate approximation to the ideal Couette profile can be obtained with a few rings

An hydrodynamic shear instability in stratified disks

We discuss the possibility that astrophysical accretion disks are dynamically unstable to non-axisymmetric disturbances with characteristic scales much smaller than the vertical scale height. The instability is studied using three methods: one based on the energy integral, which allows the determination of a sufficient condition of stability, one using a WKB approach, which allows the determination of the necessary and sufficient condition for instability and a last one by numerical solution. This linear instability occurs in any inviscid stably stratified differential rotating fluid for rigid, stress-free or periodic boundary conditions, provided the angular velocity $\Omega$ decreases outwards with radius $r$. At not too small stratification, its growth rate is a fraction of $\Omega$. The influence of viscous dissipation and thermal diffusivity on the instability is studied numerically, with emphasis on the case when $d \ln \Omega / d \ln r =-3/2$ (Keplerian case). Strong stratification and large diffusivity are found to have a stabilizing effect. The corresponding critical stratification and Reynolds number for the onset of the instability in a typical disk are derived. We propose that the spontaneous generation of these linear modes is the source of turbulence in disks, especially in weakly ionized disks.Comment: 19 pages, 13 figures, to appear in A&

Analog of Astrophysical Magnetorotational Instability in a Couette-Taylor Flow of Polymer Fluids

We report experimental observation of an instability in a Couette-Taylor flow of a polymer fluid in a thin gap between two coaxially rotating cylinders in a regime where their angular velocity decreases with the radius while the specific angular momentum increases with the radius. In the considered regime, neither the inertial Rayleigh instability nor the purely elastic instability are possible. We propose that the observed "elasto-rotational" instability is an analog of the magnetorotational instability which plays a fundamental role in astrophysical Keplerian accretion disks.Comment: 4 pages, 1 figur