Les alliages utilisés dans les applications nucléaires subissent une irradiation aux neutrons, laquelle introduit un grand nombre de défauts lacunaires et interstitiels. Au cours du temps, ces défauts migrent, se recombinent et s’agglomèrent pour former des amas. Ce phénomène physique affecte les propriétés mécaniques des aciers et conduit à sa fragilisation. Dans ce contexte, il est important de pouvoir simuler l’évolution de la microstructure à l’aide de la méthode de dynamique d’amas. Malheureusement, cette méthode devient inefficace lorsque plusieurs éléments d’alliage doivent être pris en compte. La difficulté provient du nombre trop élevé de variables de simulation à gérer. Le projet a pour objectif d’optimiser l’efficacité du code sur une architecture parallèle distribuée en faisant appel à des fonctions dédiées, vectorielles et matricielles, de la bibliothèque SUNDIALS. Cette librairie est utilisée pour intégrer l’équation différentielle ordinaire décrivant les réactions entre amas. Un autre aspect du travail, plus théorique, consistera à reformuler le problème non-linéaire de recherche de zéros du schéma d’intégration en tirant profit de la réversibilité des réactions chimiques. Cette propriété doit permettre la mise en oeuvre de solveurs directs et itératifs pour matrices creuses, symétriques et définies positives. Un axe de recheche explorera la combinaison des approches directs et itératives, en utilisant une méthode de factorisation multi-frontale de type Cholesky pour préconditionner des itérations de gradient conjugué.