On s'intéresse dans cette thèse à l'étude de méthodes numériques pour les systèmes de réaction-diffusion. Tout d'abord, on étudie le schéma par régularisation du résidu et ses extrapolations; ce schéma introduit un préconditionneur en espace lors de la discrétisation en temps. On prouve la stabilité en norme usuelle et la convergence en norme d'énergie de cette méthode et on l'applique au préconditionnement de méthodes spectrales par des méthodes d'éléments finis. Cette application nécessite le calcul d'asymptotiques précises des polynômes de Legendre et de leurs extrema. On prouve aussi la convergence et l'ordre deux d'une méthode de splitting semi-discrétisée en temps pour les systèmes de réaction-diffusion, l'approximation de Peaceman-Rachford. Enfin, on applique ces méthodes à la simulation d'une équation parabolique non linéaire pour modéliser la croissance de grains et à une équation parabolique non locale venant de la mécanique statistique et modélisant les systèmes autogravitants de fermions.