Cette thèse comporte trois parties. <br />Dans la première, nous présentons une méthode de simulation d'écoulements fluide/particules. Nous montrons que la pénalisation du tenseur des contraintes, associée à une discrétisation en temps par la méthode des caractéristiques, conduit à une formulation variationnelle de type Stokes généralisée. Des tests numériques sont effectués sous FreeFem++ afin d'étudier la convergence. Nous en présentons également trois exemples d'utilisation. <br />Dans la seconde partie nous proposons un modèle permettant de prendre en compte les forces de lubrification dans les simulations directes d'écoulements fluide/particules. Nous présentons d'abord un modèle de contact visqueux dans le cas particule/plan, obtenu comme limite, à viscosité nulle, du modèle de lubrification. Nous décrivons ensuite un algorithme reposant sur une étape de projection des vitesses, à chaque instant, sur un espace dit de vitesses admissibles. On montre alors la convergence du schéma et on généralise l'algorithme au cas multi-particules. Nous en présentons également un exemple de programmation orientée objet. <br />Dans la dernière partie, nous considérons un système discret de sphères (boulier en 1D) qui interagissent à travers la force de lubrification. Le modèle microscopique repose sur le développement de cette force à courte distance. Nous proposons une équation constitutive marcoscopique, de type Newtonien, reposant sur une viscosité linéique proportionnelle à l'inverse de la fraction locale de fluide. Nous établissons la convergence du modèle microscopique vers le modèle macroscopique proposé.