1 / 44

Application de Simulation par MCNP (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. Salle ELDORDO 78 . source. y. x. Béton. 7 cm. Salle ELDORDO 78 . source. 4m. y. (0,0,0). x. 10m x 4m x 3m . 10m.

milla
Download Presentation

Application de Simulation par MCNP (vérification de la loi de l'inverse carrée de la distance)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 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

  2. Salle ELDORDO 78 source y x

  3. Béton 7 cm Salle ELDORDO 78 source 4m y (0,0,0) x 10m x 4m x 3m 10m

  4. Salle ELDORDO 78 10m Béton (0,0,0) 7 cm 3m source Z 10m x 4m x 3m x

  5. 1- open MCNP 3- change the title 2- clique input

  6. Béton 7 cm

  7. 1. Création des surfaces

  8. 4- clique surface 5- click Wizard

  9. 6- Click Macrobodies 7- Suivant 8- Click RPP 9- Suivant

  10. 10- remplir tableau 12- valeurs tableau 13- Suivant 14- comment15- Terminer

  11. 16- input Béton 7 cm 17- cube 2

  12. 2. Création des cellules

  13. Cel 3 (béton) 18- cel 1 ( inter cube2 ) (Air)1 0 -2 Cel 1 (Air) 19- cel 2 ( ext cube1 ) (Vide)2 0 1 Cel 2 (vide) Cube 2 20- cel 3 ( cunbe 2 intersection cube1 ) (Béton)1 0 -1 2 Cube 1 21- Input

  14. 3. Définition des Matériaux

  15. 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 indicateweight fractions

  16. 4. Création de source

  17. 22- modemode p e • 23- source (60Co)-Position: (x,y,z) ? • Type de particule: photons (2) • Energie: 1,25 MeV • Type de source: ponctuelleet isotrope vide

  18. Position de source 60Co: (x,y,z) ? y= 0 source 4m y (0,0,0) x X= 4,5 m = 450 cm 10m x 4m x 3m 10m

  19. Position de source 60Co: (x,y,z) ? Position : (x,y,z) = (-450, 0, -100) (0,0,0) 3m Z= 1m = 100 cm source Z 10m x 4m x 3m x

  20. 24- modemode p e • 25- source (60Co)-Position: (x,y,z) = (-450, 0, -100) • Type de particule: photons (2) • Energie: 1,25 MeV • Type de source: ponctuelleet isotrope vide

  21. 5. Création des détecteurs

  22. Lesdétecteurs utilisés sont des sphères (cellules)

  23. R= 1 cm -330 cm -350 cm -370 cm -390 cm (0,0,0) y X= -4,5 m = -450 cm x 60 80 100 10m x 4m x 3m 120

  24. (0,0,0) Z= 1m = 100 cm Z 10m x 4m x 3m x

  25. Création des surfaces & cellulespour les détecteurs

  26. 26- surfaces (détecteurs) • 27- cellules (détecteurs)

  27. 6.Importance

  28. 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

  29. 7.Tallies

  30. TalliesLes tallies sont les observables définies dans MCNP.Ces tallies se présentent de la manière suivante : Fkn:XC1C2 … 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)

  31. 8. NPS (nombre d’histoire)

  32. NPS

  33. 9.Run

  34. 10.Output

  35. Input même place avec xsdir.

  36. Click sur le fichieroutput

  37. Date Titre Input Mass/volume

  38. Unité Tally Energie Erreur

  39. Calcul de dose • *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.

More Related