The goal of this work is to conceive a model and algorithms for the realistic simulation of granular materials that are able capture the complex dynamics created by the interactions of millions of individual grains. Such algorithms will have make use of the full potential of the modern massively parallel architectures.