Les méthodes de simulations numériques ont suscité une attention particulière et un développement considérable ces dernières années [1-4] grâce à la croissance explosive de la puissance de l'outil informatique (ordinateur). Ces méthodes sont considérées comme étant une révolution scientifique. En effet, elles constituent une troisième branche complémentaire aux deux premières traditionnelles qui sont les branches théorique et expérimentale. Souvent, les informations qui proviennent des théories analytiques ne sont exactes que dans des cas rares, car pour les autres cas des approximations incontrôlées sont souvent nécessaires. Quant aux informations issues de l'expérience, le fait qu'elles ne soient pas généralement caractérisées avec précision est lié: (i) aux phénomènes expérimentaux observés qui sont parfois intrinsèques et/ou (ii) à des effets d'impuretés inconnus (exemple, la composition chimique du matériau est connue approximativement). De ce fait, la comparaison des résultats expérimentaux et théoriques ne conduit pas forcément à des réponses conclusives d'où l'utilité de la simulation numérique. Cette dernière qui se situe à mi-chemin entre la théorie et l'expérience, sans pour autant aller se substituer à l'une ou à l'autre, constitue une approche scientifique valable pour comprendre certaines lois de la nature et apporter des réponses aux questions soulevées par les approches théoriques et expérimentales. Il existe ainsi un rapport d'interdépendance entre la théorique, la simulation et l'expérimental. Généralement, la connaissance des différentes configurations atomiques locales et par conséquent la microstructure, est nécessaire pour prévoir le comportement à long terme des matériaux intervenant dans plusieurs branches de l'industrie. Les méthodes d'interactions hyperfines telles que la Résonance Magnétique Nucléaire (RMN) et la spectrométrie Mِssbauer (SM) qui sont directement sensibles aux voisinages immédiats de l'atome résonnant, se sont avérées être un cutil d'investigatior: important pour la détermination de l'ordre local, par l'étude de la distribution de champ hyperfin dont la forme est liée à la nature des atomes solutés. En pratique, les distributions de champs hyperfins sont déduites de l'ajustement des spectres Mِssbauer par des programmes basés sur la méthode des moindres carrées non linéaires [5-8]. Toutefois, l'interprétation des résultats obtenus est rendue difficile: (i) à cause de l'existence de plusieurs éléments d'alliage dont les effets sont additifs et/ou antagonistes; et/ou (ii) par le manque d'information, dans la littérature, concernant les systèmes étudiés, d'où le recours parfois à des méthodes de simulations numériques. Parmi ces dernières,