La thèse comporte deux parties. La première traite de la convergence vers l'équilibre de flots de gradients discrets ou plus généralement de discrétisations d'un système autonome admettant une fonction de Lyapunov. En se plaçant dans une cadre pour lequel les olutions du problème continu converge vers un état stationnaire en temps grand, il est démontré sous des hypothèses générale que le système discret a la même propriété. Ce résultat conduit à des conclusions nouvelles sur le comportement en temps grand de schémas numériques anciens. La seconde partie concerne la simulation numérique de particules en suspension dans un fluide visqueux. Il est montré que les méthodes utilisées actuellement pour simuler l'interaction hydrodynamique entre particules perdent de leur précision quand de grandes forces non-hydrodynamiques sont en jeu et que au moins deux particules sont proches l'une de l'autre. Ce cas survient, dans le contexte de l'ingénierie biomédicale, lors de la conception de nano-robots capables de nager. Cette perte de précision est due au caractère singulier de l'écoulement de Stokes dans les zones de presque contact. Une nouvelle méthode est introduite ici. Des expérimentations numériques sont effectuées pour mettre en évidence sa grande précision.