Codes Nucléaires - ChainSolver - PDF

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 2

Université Ferhat Abbas Sétif-1- Année universitaire : 2020 / 2021

Faculté des Sciences Master II : Physique des Rayonnements


Département de Physique Module : Codes Nucléaires

Le code de calcul de transmutation « ChainSolver »

Le code de calcul ChainSolver :


Ce code, développé par E.G.Romanov, est destiné à la réalisation de calculs concrets pour la transmutation
des isotopes. Il prend en considération : la structure initiale des produits à transmuter, les flux de neutrons et
la période d’irradiation. Il permet de résoudre les équations de Bateman, équations différentielles ordinaires
(ODE), qui décrivent la chaîne d’évolution d’un nucléide sous irradiation :
𝑑𝑁𝑖
= − 𝜎𝑖𝑎𝑏𝑠 𝜑 𝑁𝑖 + ∑ 𝜎𝑗→𝑖 𝜑 𝑁𝑗 − 𝜆𝑖 𝑁𝑖 + ∑ 𝜆𝑗 →𝑖 𝑁𝑗
𝑑𝑡
𝑖 ≠𝑗 𝑗

Avec : Ni : le nombre de noyaux i ; 𝜎𝑖𝑎𝑏𝑠 𝜑 : le taux de réaction d’absorption du noyau i ; 𝜎𝑗→𝑖 𝜑 : le taux de
réaction du noyau j produisant i ; λi : la constante de décroissance totale du noyau i ; λj → i : la constante de
décroissance du noyau j vers le noyau i.

L’interface du programme ChainSolver est construite d’un composant visuel Delphi, qui comporte cinq
pages complémentaires : Chain, Conditions, Time, Calc, et Answers.

1- La page « Chain » :
C’est la première page qui permet d’éditer les chaînes de transmutation, voir la figure 1. La chaîne peut être
construite par l’introduction des nucléides dans la zone d'édition de "States", en vérifiant l’instruction
"Build" pour choisir le type de liens : "Decay" désintégration radioactive, "Capture" capture neutronique,
"Fission" fission induite, "Threshold" réaction de seuil. Les données doivent être initialement chargées dans
le programme par le bouton "Load DB".

Fig.1. La page « Chain ».

2- La Page « Conditions » :
Elle permet d’introduire les caractéristiques de l’échantillon à irradier, voir la figure 2. Le tableau "Masses,
g" est rempli par les poids initiaux des nucléides présents dans la chaîne. La case "masse totale, g" qui est
remplie par la masse totale de l’échantillon.

Page - 1 - Préparé par : Dr. Mounira HOUAS


Fig.2. La page « Conditions ».
3- La page « Time » :
Cette page (figure 3) permet d’introduire les conditions d’irradiation : la période d’irradiation (jours), la
puissance du réacteur (MW) et les flux neutroniques (thermique, de résonance et rapide). Le programme
permet aussi d’effectuer des irradiations multiples avec différentes puissances et positions d’irradiation.

Fig.3. La page « Time ».


4- La page « Calc » :
Elle donne l'interface aux quatre sous-programmes utilisés pour la résolution du système des équations de
Bateman : VODE, LSODA, RADAU et MEBDF, voir la figure 4. Le tableau ensuite affiche les masses
calculées des isotopes de la chaîne pour chaque intervalle de temps.

Fig.4. La page « Calc ».


5- La page « Answers » :
C’est la dernière page (figure 5) qui présente les résultats de calcul, comme l’évolution des masses des
radionucléides présents dans la chaîne, les activités, …. On sauvegarde ces résultats dans un fichier Excel.

Fig.5. La page « Answers ».

Page - 2 - Préparé par : Dr. Mounira HOUAS

Vous aimerez peut-être aussi