Numerical simulations of dissipationless disk accretion

Our goal is to study the regime of disk accretion in which almost all of the angular momentum and energy is carried away by the wind outflowing from the disk in numerical experiments. For this type of accretion the kinetic energy flux in the outflowing wind can exceed considerably the bolometric luminosity of the accretion disk, what is observed in the plasma flow from galactic nuclei in a number of cases. In this paper we consider the nonrelativistic case of an outflow from a cold Keplerian disk. All of the conclusions derived previously for such a system in the self-similar approximation are shown to be correct. The numerical results agree well with the analytical predictions. The inclination angle of the magnetic field lines in the disk is less than 60°, which ensures a free wind outflow from the disk, while the energy flux per wind particle is greater than the particle rotation energy in its Keplerian orbit by several orders of magnitude, provided that the ratio r _A/r ? 1, where r _A is the Alfvénic radius and r is the radius of the Keplerian orbit. In this case, the particle kinetic energy reaches half the maximum possible energy in the simulation region. The magnetic field collimates the outflowing wind near the rotation axis and decollimates appreciably the wind outflowing from the outer disk periphery.