The processes of planet formation and migration depend intimately on the interaction between planetesimals and the gaseous disks in which they form. The formation of gaps in the disk can severely limit the mass of the planet and its migration toward the protostar. We investigate the process of gap formation through magnetohydrodynamic simulations in which internal stress arises self-consistently from turbulence generated by the magnetorotational instability. The simulations investigate three different planetary masses and two disk temperatures to bracket the tidal (thermal) and viscous gap opening conditions. The results are in general qualitative agreement with previous simulations of gap formation, but show significant differences. In the presence of MHD turbulence, the gaps produced are shallower and asymmetrically wider than those produced with pure hydrodynamics. The rate of gap formation is also slowed, with accretion occurring across the developing gap. Viscous hydrodynamics does not adequately describe the evolution, however, because planets capable of producing gaps also may be capable of affecting the level MHD turbulence in different regions of the disk.

Chaotic flow is studied in a series of numerical magnetohydrodynamical simulations that use the shearing box formalism. This mimics important features of local accretion disk dynamics. The magnetorotational instability gives rise to flow turbulence, and quantitative chaos parameters, such as the largest Lyapunov exponent, can be measured. Linear growth rates appear in these exponents even when the flow is fully turbulent. The extreme sensitivity to initial conditions associated with chaotic flows has practical implications, the most important of which is that hundreds of orbital times are needed to extract a meaningful average for the stress. If the evolution time in a disk is less than this, the classical $\alpha$ formalism will break down.

We employ a variety of numerical simulations in the local shearing box system to investigate in greater depth the local hydrodynamic stability of Keplerian differential rotation. In particular we explore the relationship of Keplerian shear to the nonlinear instabilities known to exist in simple Cartesian shear. The Coriolis force is the source of linear stabilization in differential rotation. We exploit the formal equivalence of constant angular momentum flows and simple Cartesian shear to examine the transition from stability to nonlinear instability. The manifestation of nonlinear instability in shear flows is known to be sensitive to initial perturbation and to the amount of viscosity; marginally (linearly) stable differentially rotating flows exhibit this same sensitivity. Keplerian systems, however, are completely stable; the strength of the stabilizing Coriolis force easily overwhelms any destabilizing nonlinear effects. In fact, nonlinear effects speed the decay of applied turbulence by producing a rapid cascade of energy to high wavenumbers where dissipation occurs. Our conclusions are tested with grid resolution experiments and by comparison with results from a code that employs an alternative numerical algorithm. The properties of hydrodynamic differential rotation are contrasted with magnetohydrodynamic differential rotation. The kinetic stress couples to the vorticity which limits turbulence, while the magnetic stress couples to the shear which promotes turbulence. Thus magnetohydrodynamic turbulence is uniquely capable of acting as a turbulent angular momentum transport mechanism in disks.