Effects of magnetic field on stellar differential rotation are studied by comparing magnetohydrodynamic (MHD) models and their hydrodynamic (HD) counterparts in the broad range of rotation rate and in varying initial rotation profile. Fully-compressible MHD simulations of rotating penetrative convection are performed in a full-spherical shell geometry. Critical conditions for the transition of the differential rotation between faster equator (solar-type) and slower equator (anti-solar type) are explored with focusing on the "Rossby number (${\rm Ro}$)" and the "convective Rossby number (${\rm Ro}_{\rm conv}$)". It is confirmed that the transition is more gradual and the critical value for it is higher in the MHD model than the HD model in the view of the ${\rm Ro}_{\rm conv}$-dependence. The rotation profile shows, as observed in earlier studies, the bistability near the transition in the HD model, while it disappears when allowing the growth of magnetic fields except for the model with taking anti-solar type solution as the initial condition. We find that the transition occurs at ${\rm Ro} \simeq 1$ both in the MHD and HD models independently of the hysteresis. Not only the critical value, the sharpness of the transition is also similar between the two models in the view of the ${\rm Ro}$-dependence. The influences of the dynamo-generated magnetic field and/or the hysteresis on convective motion are reflected in the ${\rm Ro}$. This would be the reason why the transition is unified in the view of the ${\rm Ro}$-dependence. We finally discuss the ${\rm Ro}$-dependence of magnetic dynamo activities with emphasis on its possible relation to the kinetic helicity profile.

Spherical solar dynamo simulations are performed. Self-consistent, fully compressible magnetohydrodynamic system with a stably stratified layer below the convective envelope is numerically solved with a newly developed simulation code based on the Yin-Yang grid. The effects of penetrative convection are studied by comparing two models with and without the stable layer. The differential rotation profile in both models is reasonably solar-like with equatorial acceleration. When considering the penetrative convection, a tachocline-like shear layer is developed and maintained beneath the convection zone without assuming any forcing. While turbulent magnetic field becomes predominant in the region where the convective motion is vigorous, mean-field component is preferentially organized in the region where the convective motion is less vigorous. Especially in the stable layer, the strong large-scale field with a dipole symmetry is spontaneously built up. The polarity reversal of the mean-field component takes place globally and synchronously throughout the system regardless the presence of the stable layer. Our results suggest that the stably stratified layer is a key component for organizing the large-scale strong magnetic field, but is not essential for the polarity reversal.

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.

Although the magnetorotational instability (MRI) has been widely accepted as a powerful accretion mechanism in magnetized accretion disks, it has not been realized in the laboratory. The possibility of studying MRI in a rotating liquid-metal annulus (Couette flow) is explored by local and global stability analysis and magnetohydrodynamic (MHD) simulations. Stability diagrams are drawn in dimensionless parameters, and also in terms of the angular velocities at the inner and outer cylinders. It is shown that MRI can be triggered in a moderately rapidly rotating table-top apparatus, using easy-to-handle metals such as gallium. Practical issues of this proposed experiment are discussed.