In this work we present an efficient implementation of Eulerian TVD methods. We apply parallelization strategies based entirely on GPU for the solution of the 2D transport equation in heterogeneous porous media. Additionally, a parallel strategy is proposed for the generation of exponentially correlated lognormally distributed permeability fields in GPU. The programs are developed using C++/CUDA. The implemented methods are used to solve advective dominant problems, in a context of Monte Carlo type simulations to numerically determine the longitudinal and transversal macrodispersion coefficients averaging over 100 simulations for permeability fields for a large range of variances. The following types of transport are considered for testing: pure advection, advection-diffusion and advection-dispersion. The performance in terms of the computation time of explicit and implicit methods are compared. We show that the implemented algorithms allow to efficiently solve problems in computational domains of up to 134.5 million cells in a single GPU. L.B. acknowledges the support of UTN FRCON through “Scholarships for the training of doctors to strengthen the areas of R+D+i, Res. No. 1460/15”. Peer reviewed