Transcript vérification de la loi de l`inverse carrée de la distance

Application de Simulation par MCNP
(vérification de la loi de l'inverse
carrée de la distance)
Cas de l’unité d’irradiation
ELDORDO 78
(CRNA)
Par: A. SIDI MOUSSA
y
x
7 cm
y
(0,0,0)
x
10m x 4m x 3m
10m
4m
10m
(0,0,0)
3m
Z
10m x 4m x 3m
x
7 cm
1- open MCNP
3- change the title
2- clique input
7 cm
1. Création des surfaces
4- clique surface
5- click Wizard
6- Click Macrobodies
7- Suivant
8- Click RPP
9- Suivant
10- remplir tableau
12- valeurs tableau
13- Suivant
14- comment
15- Terminer
16- input
17- cube 2
7 cm
2. Création des cellules
Cel 3 (béton)
18- cel 1 ( inter cube2 ) (Air)
1 0 -2
Cel 1 (Air)
Cube 2
Cel 2 (vide)
19- cel 2 ( ext cube1 ) (Vide)
2 0 1
Cube 1
20- cel 3 ( cunbe 2 intersection cube1 ) (Béton)
1 0 -1 2
21- Input
3. Définition des Matériaux
22- Air
d=1.29 mg/cm3 = 0.00129 g/cm3
Nitrogen=0.755268
Oxygen=0.231781
Argon=0.012827
Carbon=0.000124
The negative values indicate weight fractions
4. Création de source
22- mode
mode p e
23- source (60Co)
- Position: (x,y,z) ?
-Type de particule: photons (2)
-Energie: 1,25 MeV
-Type de source: ponctuelle
et isotrope
vide
Position de source 60Co: (x,y,z) ?
y= 0
y
(0,0,0)
x
X= 4,5 m = 450 cm
10m x 4m x 3m
10m
4m
Position de source 60Co: (x,y,z) ?
Position : (x,y,z) = (-450, 0, -100)
(0,0,0)
3m
Z= 1m = 100 cm
Z
10m x 4m x 3m
x
24- mode
mode p e
25- source (60Co)
- Position: (x,y,z) = (-450, 0, -100)
-Type de particule: photons (2)
-Energie: 1,25 MeV
-Type de source: ponctuelle
et isotrope
vide
5. Création des détecteurs
Les détecteurs utilisés sont des
sphères (cellules)
R= 1 cm
-330 cm
-350 cm
-370 cm
-390 cm
(0,0,0)
y
X= -4,5 m = -450 cm
x
10m x 4m x 3m
60
80
100
120
(0,0,0)
Z= 1m = 100 cm
60 80 100
Z
120
10m x 4m x 3m
x
Création des surfaces & cellules
pour les détecteurs
26- surfaces (détecteurs)
27- cellules (détecteurs)
6. Importance
Les particules propagées ont un poids, généralement égal à 1.
Ce poids intervient dans toutes les formules comme élément
de pondération : une particule de poids 0.1 influençant 10 fois
moins un tally qu’une particule de poids 1.
L’importance d’une cellule est définie par la carte IMP :X=imp
avec X la particule (N, P ou E).
Cette définition peut se faire sur la même ligne que la
déclaration de cellule ou tout à la fin du fichier MCNP.
imp:p 1 0 0 1 1 1 1
imp:e 1 0 0 1 1 1 1
7. Tallies
Tallies
Les tallies sont les observables définies dans MCNP.
Ces tallies se présentent de la manière suivante :
Fkn:X C1 C2 … Ci (tally sur une cellule)
avec : k un nombre entre 0 et 99 destiné à différencier les tallies de
même type
n un chiffre entre 1 et 8 destiné à indiqué le type de tally calculé
X le type de particule
Ci une cellule sur laquelle on veut calculer le tally
On a :
- F1 : courant surfacique
- F2 : flux surfacique
- F4 : flux volumique
- F6 : énergie déposée
- F8 : « pulse height tally »
On peut rajouter une * devant le F pour changer l’unité du résultat (se reporter à la documentation
MCNP)
8. NPS (nombre d’histoire)
NPS
9. Run
10. Output
• Input même place
avec xsdir.
Click sur
le fichier
output
Date
Titre
Input
Mass/volume
Unité
Tally
Energie
Erreur
Calcul de dose
∗
F8 . A . t
D=
1,6. 10−10
m
• *F8: énergie en (Mev).
• A: Activité en (Bq) Unité usuelle: 1[Ci] = 3.7E+10 [Bq].
• t: temps d’irradiation en (s).
• m: mass en (g).
• 1,6.10-10 : facteur de conversion (Mev/g) --˃ Gy.