Rossby wave instability of thin accretion disks. III. Nonlinear simulations

We study the nonlinear evolution of the Rossby wave instability in thin dis ks using global two-dimensional hydrodynamic simulations. The detailed line ar theory of this nonaxisymmetric instability was developed earlier by Love lace et al. and Li et al., who found that the instability can be excited wh en there is an extremum in the radial profile of an entropy-modified versio n of potential vorticity. The key questions we are addressing in this paper are the following : (1) What happens when the instability becomes nonlinea r ? Specifically, does it lead to vortex formation? (2) What is the detaile d behavior of a vortex ? (3) Can the instability sustain itself and can the vortex last a long time? Among various initial equilibria that we have exa mined, we generally find that there are three stages of the disk evolution : (1) The exponential growth of the initial small amplitude perturbations. This is in excellent agreement with the linear theory; (2) The production o f large-scale vortices and their interactions with the background flow, inc luding shocks. Significant accretion is observed owing to these vortices. ( 3) The coupling of Rossby waves/vortices with global spiral waves, which fa cilitates further accretion throughout the whole disk. Even after more than 20 revolutions at the radius of vortices, we find that the disk maintains a state that is populated with vortices, shocks, spiral waves/shocks, all o f which transport angular momentum outward. We elucidate the physics at eac h stage and show that there is an efficient outward angular momentum transp ort in stages (2) and (3) over most parts of the disk, with an equivalent S hakura-Sunyaev angular momentum transport parameter alpha in the range from 10(-4) to 10(-2). By carefully analyzing the flow structure around a vorte x, we show why such vortices prove to be almost ideal "units" in transporti ng angular momentum outward, namely by positively correlating the radial an d azimuthal velocity components. In converting the gravitational energy to the internal energy, we find some special cases in which entropy can remain the same while angular momentum is transported. This is different from the classical a-disk model, which results in the maximum dissipation (or entro py production). The dependence of the transport efficiency on various physi cal parameters are examined and effects of radiative cooling are briefly di scussed as well. We conclude that Rossby wave/vortex instability is an effi cient, purely hydrodynamic mechanism for angular momentum transport in thin disks, and may find important applications in many astrophysical systems.