Cette thèse s'inscrit dans le cadre de la modélisation en bio-mathématiques et dans celui de l'analyse numérique et du calcul scientifique. Le modèle bidomaine décrit l'activité électrique du coeur. Cette activité est complexe : elle relève à l'échelle cellulaire de processus biochimiques et à l'échelle macroscopique de la structure anisotrope des tissus cardiaques, des caractéristiques du thorax. Une application fondamentale du modèle est la simulation d'électrocardiogrammes. Des méthodes de calcul type volumes finis sont développées pour la résolution du modèle. Dans un premier temps, la stabilité et la convergence de schémas volumes finis classiques est établie, en théorie et numériquement, pour une version simplifiée du modèle bidomaine. Pour faire face à des difficultés conceptuelles et pratiques du modèle complet (anisotropie des tissus, conditions limites, maillages non structurés distordus), une seconde classe de schémas 2D-3D, cell-vertex centered, est mise au point et testée.

L'analyse numérique se développe en un outil puissant dans l'étude des équations aux dérivées partielles (EDPs), permettant d'illustrer des théorèmes existants et de trouver des conjectures. En utilisant des techniques sophistiquées, des questions apparaissant inaccessibles avant, comme des oscillations rapides ou un blow-up des solutions, peuvent être étudiées. Des oscillations rapides dans les solutions sont observées dans des EDPs dispersives sans dissipation ou les solutions des EDPs correspondantes sans dispersion ont des chocs. Pour résoudre numériquement ces oscillations, l'application de méthodes efficaces introduisant peu de dissipation numérique artificielle est impérative, en particulier pour l'étude d' EDPs en plusieurs dimensions. Comme les EDPs étudiées dans ce contexte sont typiquement raides, l'intégration efficace dans le temps représente le principal problème. Une analyse des intégrants exponentiels et symplectiques a permis de déterminer les méthodes les plus efficaces pour chaque EDP étudiée. L'apprentissage et l'utilisation de techniques de parallélisation de codes numériques permet de nos jours de grandes avancées, plus précisément dans ce travail d'étudier numériquement la stabilité des solutions et l'apparition de blow-up dans l'équation de Davey-Stewartson.