III. 2. Présentation du code MCNP
Une présentation exhaustive du code MCNP sort du cadre
de ce manuscrit, on se contentera plutôt de la description et la
présentation du fichier d'entrée "Input". C'est le fichier
à travers lequel on arrive à introduire nos commandes et nos
directives pour obtenir la simulation requise.
III. 2. 1. Structure du fichier d'entrée
Tout utilisateur du code doit tout d'abord écrire un
fichier d'entrée (input file) qui sera lu par le
simulateur MCNP. Le contenu du fichier d'entrée MCNP
nécessite plusieurs
Chapitre III Code Simulation Monte-Carlo
24
paramètres décrits à l'aide de cartes
réparties en bloc. Sa structure est constituée de trois grands
blocs séparés par une ligne vide : les cartes
des cellules, les cartes des surfaces et les cartes des
données dont y existent :
- Les données pour la source,
- Les données pour le matériau,
- Les données pour la grandeur à calculer ou
tally.
La structure du fichier d'entrée est la suivante :
Figure 12 : Schématisation du fichier
input du code MCNP.
Une autre caractéristique du code MCNP5,
c'est que celui-ci permet la préparation des fichiers input
moyennant une interface graphique appelée « Vised ». Ceci
permet de visualiser instantanément le modèle introduit dans
l'input (Figure 12).
Chapitre III Code Simulation Monte-Carlo
25
Figure 13 : interface Vised.
III. 2. 1. a) Définition des cellules [23]
Ce sont des boîtes constituant la
géométrie à simuler. Elles se constituent des
différents paramètres que sont : les numéros de cellule et
de matériau attribués par l'utilisateur, la densité de
matériau, les numéros des surfaces désignant la cellule
avec les signes adéquats et les importances pour chaque particule. Un
numéro de surface muni d'un signe (-) indique l'intérieur de la
cellule, sinon son extérieur. Un zéro à la place du
numéro de matériau et de la densité indique une cellule
vide. La cellule peut être définie par deux opérateurs
d'intersection et d'union. Ils sont représentés respectivement
par un espace et le caractère (:) entre deux numéros de surface
dans la carte cellule. Il est aussi possible d'exclure une cellule d'un espace
par l'opérateur symbolisé par (#). Les cellules
sont définies verticalement et se déclarent de la manière
suivant :
j m d geom imp: x= params
Avec :
j : est un nombre compris entre 1
et 99999 (colonne 1-5) et définit le
numéro de la cellule.
m : désigne la matière constituant
cette cellule. Si on choisit m=0, la cellule est
déclarée vide ;
si m?0, ce paramètre indique le
matériau constituant la cellule.
d : densité du matériau de la
cellule :
o
Chapitre III Code Simulation Monte-Carlo
26
Pas valeur de densité pour le vide donc égal
zéro,
o Valeur positive = densité atomique (atoms
/cm3),
o Valeur négative = densité de masse
(g/cm3),
geom : spécifie la
géométrie de la cellule par les surfaces qui la
délimitent. imp:x : importance de la cellule pour la
particule x. params : paramètres optionnels de la
cellule.
Par exemple
Commentaire
C Cell card
4 1 1.234 E-3 -1 2 -3 imp : n=1
Numéro de la cellule densité surface
paramètres (importance de la
particule dans le calcul)
Numéro du matériau de la
cellule
|