Global gyrokinetic simulations have revealed an important nonlinear flow generation process due to the residual stress produced by electrostatic turbulence of ion temperature gradient (ITG) modes and trapped electron modes (TEMs). In collisionless TEM (CTEM) turbulence, nonlinear residual stress generation by both the fluctuation intensity and the intensity gradient in the presence of broken symmetry in the parallel wavenumber spectrum is identified for the first time. Concerning the origin of the symmetry breaking, turbulence self-generated low frequency zonal flow shear has been identified to be a key, universal mechanism in various turbulence regimes. Simulations reported here also indicate the existence of other mechanisms beyond E×B shear. The ITG turbulence driven "intrinsic" torque associated with residual stress is shown to increase close to linearly with the ion temperature gradient, in qualitative agreement with experimental observations in various devices. In CTEM dominated regimes, a net toroidal rotation is driven in the cocurrent direction by intrinsic torque, consistent with the experimental trend of observed intrinsic rotation. The finding of a "flow pinch" in CTEM turbulence may offer an interesting new insight into the underlying dynamics governing the radial penetration of modulated flows in perturbation experiments. Finally, simulations also reveal highly distinct phase space structures between CTEM and ITG turbulence driven momentum, energy, and particle fluxes, elucidating the roles of resonant and non-resonant particles. © 2010 American Institute of Physics.