Sedimenting suspensions exist in a varity of natural phenomena and industrial applications. It is already observed in experiments that the dilute fibre suspensions experience a concentration instability under gravity at low Reynolds numbers. Initially well-mixed suspensions become inhomogeneous and anisotropic due to this instability. This project is focused on the development and validation of numerical models to understand the instability in a dilute fibre suspension by means of the mixture model and the point-particle model. For periodic boundary condition, we use a linear stability analysis to show that inertia and hydro dynamic translational diffusion damp perturbations at long wavelengths and short wavelengths, respectively, leading to a wavenumber selection. However, numerical simulations indicate a weak wavenumber selection even at zero Reynolds number. Numerical simulations also show that the induced flow may either die or saturate on a finite amplitude. The characterof this long time behaviour is dictated by the wavenumber, the presence or absence of the translational diusivity, rotational diffusivity, and the fluid inertia on particle motions. Moreover, the most unstable wavenumber decreases with time and the maximum amplitude increases. The smallest wavenumber obtains the largest amplitude at steady state. For a vessel bounded by sidewalls, the near-wall convection is an upward back flow in the very beginning, due to the combined effects of the steric-depleted layer and a hydrodynamiclly-depleted region near the wall. However, the evolution of the near-wall convection at later times depends on the aspect ratio of the bres, the translational diffusivity and the initial perturbations. The steric-depleted layer in the mixture model can be neglected for large widths. Multiple streamers are obtained due to the sidewalls, implying that the sidewalls can generate a wavelength which is smaller than the channel width. The suspension ends up with a single streamer on one side of the container, consistent with the results of the cases with periodic boundary condition but different from the experimental results. This might be due to the absence of the botton wall in the mixture model. Moreover, the global structure evolution of a suspension is dependent on the width of the vessel and the amplitude ofthe initial perturbations. / <p>QC 20140207</p>