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.

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.

Radial substructures in circumstellar disks are now routinely observed by ALMA. There is also growing evidence that disk winds drive accretion in such disks. We show through 2D (axisymmetric) simulations that rings and gaps develop naturally in magnetically coupled disk-wind systems on the scale of tens of au, where ambipolar diffusion (AD) is the dominant non-ideal MHD effect. In simulations where the magnetic field and matter are moderately coupled, the disk remains relatively laminar with the radial electric current steepened by AD into a thin layer near the midplane. The toroidal magnetic field sharply reverses polarity in this layer, generating a large magnetic torque that drives fast accretion, which drags the poloidal field into a highly pinched radial configuration. The reconnection of this pinched field creates magnetic loops where the net poloidal magnetic flux (and thus the accretion rate) is reduced, yielding dense rings. Neighbouring regions with stronger poloidal magnetic fields accrete faster, forming gaps. In better magnetically coupled simulations, the so-called `avalanche accretion streams' develop continuously near the disk surface, rendering the disk-wind system more chaotic. Nevertheless, prominent rings and gaps are still produced, at least in part, by reconnection, which again enables the segregation of the poloidal field and the disk material similar to the more diffusive disks. However, the reconnection is now driven by the non-linear growth of MRI channel flows. The formation of rings and gaps in rapidly accreting yet laminar disks has interesting implications for dust settling and trapping, grain growth, and planet formation.

Rings and gaps are being observed in an increasing number of disks around young stellar objects. We illustrate the formation of such radial structures through idealized, 2D (axisymmetric) resistive MHD simulations of coupled disk-wind systems threaded by a relatively weak poloidal magnetic field (plasma-$\beta \sim 10^3$). We find two distinct modes of accretion depending on the resistivity and field strength. A small resistivity or high field strength promotes the development of rapidly infalling `avalanche accretion streams' in a vertically extended disk envelope that dominates the dynamics of the system, especially the mass accretion. The streams are suppressed in simulations with larger resistivities or lower field strengths, where most of the accretion instead occurs through a laminar disk. In these simulations, the disk accretion is driven mainly by a slow wind that is typically accelerated by the pressure gradient from a predominantly toroidal magnetic field. Both wind-dominated and stream-dominated modes of accretion create prominent features in the surface density distribution of the disk, including rings and gaps, with a strong spatial variation of the magnetic flux relative to the mass. Regions with low mass-to-flux ratios accrete quickly, leading to the development of gaps, whereas regions with higher mass-to-flux ratios tend to accrete more slowly, allowing matter to accumulate and form dense rings. In some cases, avalanche accretion streams are observed to produce dense rings directly through continuous feeding. We discuss the implications of ring and gap formation driven by winds and streams on grain growth and planet formation.