Large-scale, fully three-dimensional particle simulations have been carried out for the collisionless drift instabilities in a cylindrical geometry. It is found that nonlinear excitation of convective cells (w = 0,& !( = 0) due to drift instabilites gives rise to enhanced turbulence and anomalous plasma diffusion. The observed power spectra have resemblance with the recent measurements in toroidal devices.Here we report our recent results of numerical simulations on collisionless drift instabilities using large-scale fully three-dimensional cylindrical plasma models developed earlier. 1 Theoretical interpretations of the observed results on plasma turbulence and diffusion are also given. Drift-wave turbulence has been of current interest in toroidal confinement devices such as tokamaks, 2 stellarators, 3 and internal ring devices 45 because of the strong correlation between the observed density fluctuations and anomalous plasma transport.The simulation model used is a straight cylinder in a uniform external magnetic field B 0 along the z direction with its length L^/p i = 640, p t = (Tjm i ) 1/2 /^2 i being the ion gyroradius. In the cross section, a 64x64 (L 2 ) spatial grid is used for numerical computation with its physical length L/pi = 32. The plasma is periodic in z while it is surrounded by a conducting wall at the boundary of the cross section. 1 Initially both ions and electrons have Maxwellian velocity distributions with T e /T i = 4. Initial plasma density is taken to be n e (r) = n i (r) = w 0 exp(-4r 2 /a 2 ) with the average density given by Q, e /o) Pe = 5, m { /m. e -100 9 a = L/2; and there is no initial temperature gradient. Seven Fourier modes w = 0, ±1, ±2, and ±3 are kept in the z direction with k z~ 2m/L z .Linear theory predicts that the collisionless drift instability (universal mode) is strongly unstable for n = ±l and becomes stable rapidly with increasing \n\ because of the onset of ion Landau damping. This situation is commonly observed in Q devices. In terms of physical size, the simulation plasma is certainly smaller than a tokamak plasma but may be larger than a conventional Q device; and one can study quite a few problems with the model using realistic plasma parameters.Let us first look at the gross behavior of the instability. Figure 1 shows the particle diffusion, heat transfer, electron velocity distribution, and radial-mode structure associated with the insta-