In the past few years we have been involved in the development and implementation of parallel finite element and pseudo-spectral methods for porous media flow. In this paper we present a brief review of our methods and show several examples of highly nonlinear simulations in two and three-dimensions, on different current high performance computers. New issues on finite element techniques and pseudo-spectral methods are also addressed