Anomalous momentum transport in a typical astrophysical return-current beam plasma system is studied by means of two-dimensional PIC code simulations. A forward going hot electron beam compensated by a cold return beam is considered. A linear dispersion analysis predicts the linerarly unstable wave modes. Our simulation reveals that the nonlinerly generated waves and the consequent wave-particle interactions cause the electron heating and the relaxation of the electron drifts. Both, the developments of electrostatic and electromagnetic waves are analyzed as well as the roles they play in energy conversion. In particular it is found that the relaxation of electron drifts is stronger if the electromagnetic turbulence is taken into account.