MCNP : the Godfather

Bonjour,
Merci klax.
@lail : il vous faut prendre le fichier "fichier entr,.txt" le renommer en fichier de 7 lettres sans extension, par exemple "flaile"et lancer mcnp "mcnpx n=flaile" (en fonction de la version de mcnp la commande peut changer".
GA

Bonjour, je débute avec le code mcnp, je suis arrivé à faire des trucs avec, et là je dois mesurer l'efficacité d'un détecteur NaI, je bloque sur la partie tally car je dois en même temps tiré le spectre tracé, 

Merci d'avance pour vos réponses


Sujet déplacé : ce topic rassemble toutes les questions sur MCNP.

Bonsoir,
Il suffit de pondérer ton tally par E.
Par exemple E8 0 1024i 3.072 pour un 1024 canaux.

Bonjour a tous,

Enfin un Forum francais sur MCNP. Je vais tenter ma Chance!

Je simule actuellement un detecteur a scintillation NaI(Tl). Je souhaite obtenir le spectre du Cs137 pour commencer.
Ma Simulation est plutot complete: géometrie interne du detecteur, environnement...
Mes résultats sont jolis par leur forme. Mais le compton continuum est trop Faible par Rapport a mes essais tandis que le photopic est trop élevé. En jouant un peu Avec la densité de l'air ou du cristal, j'obtiens des résultats "satisfaisants". Mais voila, l'air a vu sa densité multipliée par 20.

Je suis tout ouie pour des pistes de reflections.

PS: je ne connais pas le vocabulaire technique... en francais!

Bonjour,
Le mieux est de mettre en pj votre fichier d'entrée afin de regarder si il n'a pas des évidences à corriger.

Bonjour,

Etes vous sur des dimensions de votre cristal?
Un plateau Compton plus faible et un pic plus grand signe un détecteur plus efficace que prévu en absorption totale, c'est à dire globalement plus volumineux ou avec un rapport Volume/surface plus grand à volume identique

Cordialement

Gluonmou

Après verification, l'activité de la soruce a une erreur maximale de 3%. Apres Simulation, seul le photopic est trop elevé, le continuum est tres satisfaisant.

Mais treve de blabla, voici la bete (les annotations sont en allemand)

c BLOCK 1: CELLS ***********************************************
   1 0          101          imp:p=0 $ Au�er des Problems
   2 15 -.0012 -101 119 121 111 129 127 130 131 132 133 134 135 136 137 &
     138 139 140 141 147 150 151 imp:p=1 $ Au�er des Detektors
   3 15 -.0012 (104 109 110 -122 123)-120 118 116 114:-115)-113 110): &
     (-112 110)-124 104)-125 104) imp:p=1 $ Au�er des Housing    
  11 1 -3.667  -102          imp:p=1 $ Probe aus NaJ (Dichte = 3,667 g/cm�)
  12 2 -2.7    -103 144      imp:p=1 $ Alu extrabeschichtung
  13 2 -2.7    -104 103      imp:p=1 $ Kristall Housing (AlMg3)
  14 3 -2.23   -105 106      imp:p=1 $ Borosilicate glass
  15 4 -8.7    -107 105      imp:p=1 $ Mu Metall
  16 5 -.93    -108 107      imp:p=1 $ PE Housing
  17 2 -2.7    -109 108      imp:p=1 $ Housing (AlMg3)
  18 15 -.0012 -106 126 143 142 imp:p=1 $ Luft in Szintillator
  19 5 -0.93  (-110 146)-110 -145)   imp:p=1 $ Plastik Sockel (PE)
  20 6 -8.47   -111 112 113  imp:p=1 $ Sockel (CuZn39Pl3)
  21 2 -2.7    -114 115      imp:p=1 $ Scheibe (Alu)
  22 7 -2.7    -116 117      imp:p=1 $ Tr�ger (Alu)
  23 2 -2.7    -118          imp:p=1 $ Scheibe (Alu)
  24 2 -2.7    -119 120      imp:p=1 $ Hintere Haube
  25 2 -2.7    -121 122      imp:p=1 $ Vornere Haube
  26 8 -1.38   -123 124 125  imp:p=1 $ Kappe (PET)
  27 9 -3.5    -126          imp:p=1 $ Dynode aus Stahl
  28 10 -1     -127          imp:p=1 $ Tisch aus Holz $ dichte 1 statt .75
  29 11 -2.25  -129:-151     imp:p=1 $ Wand und Boden aus Concrete MCNP
  30 2 -2.7 (-130:-131:-132) 140.1 141.1 imp:p=1 $ Unter Tr�ger
  31 2 -2.7   -133 140.1 141.1 121.1 imp:p=1 $ Tr�ger
  32 10 -1 (-134:-135) 136.1 imp:p=1 $ Drehtisch $ dichte 1 statt .75
  33 13 -8.07  -136          imp:p=1 $ Drehtisch Achse
  34 9 -5.5   -137:-138:-139 imp:p=1 $ Tisch Struktur
  35 9 -7.87   -140:-141     imp:p=1 $ Schrauben
  36 9 -7.87   -142          imp:p=1 $ Kollimator
  37 14 -4.5   -143          imp:p=1 $ 1. Dynode
  38 12 -3.58  -144 102      imp:p=1 $ Reflektor (MgO)
  39 9 -7      -146 145 -110 imp:p=1 $ Ring in Sockel
  40 9 -7      -147          imp:p=1 $ Buchse
  41 16 -1.8   -117          imp:p=1 $ Leiterblatt
  42 17 -1.873 -148          imp:p=1 $ Quelle
  43 10 -1     -149 148      imp:p=1 $ Papier um die Quelle
  44 5 -.93    -150 149      imp:p=1 $ Plastik um die Quelle
 
c BLOCK 2: SURFACE CARDS ***************************************
 101 RPP -100 100 -100 200 -100 100    $ Problem Zylinder
 102 RCC -2.54 0 0 5.08 0 0 2.54       $ Kristall Zylinder
 103 RCC -2.7737 0 0 5.3137 0 0 2.7737 $ Extra beschichtung (0,2337cm dick)
 104 RCC -2.8245 0 0 5.3645 0 0 2.8245 $ Housing Zylinder (0,0508cm dick)
 105 RCC 2.54 0 0 13.0251 0 0 2.54     $ Au�er Glass
 106 RCC 2.64 0 0 12.8251 0 0 2.44     $ Inner Glass (1mm dick)
 107 RCC 2.54 0 0 13.0251 0 0 2.5908   $ Mu Metall (0,0508cm)
 108 RCC 2.54 0 0 13.0251 0 0 2.8702   $ PE Housing
 109 RCC 2.54 0 0 13.0251 0 0 2.921    $ Alu Housing (0,0508cm)
 110 RCC 15.5651 0 0 4 0 0 2.8         $ Plastik Sockel
 111 RCC 18.5651 0 0 1.9 0 0 3.25      $ CuZn39Pb3 Sockel
 112 RCC 18.5651 0 0 1.9 0 0 1.7       $ Inner CuZn39Pb3 Sockel
 113 RCC 18.5651 0 0 1.5 0 0 2.85      $ Inner CuZn39Pb3 Sockel
 114 RCC 20.4651 0 0 0.4 0 0 2.95      $ Alu Scheibe
 115 RCC 20.4651 0 0 0.4 0 0 1.5       $ Alu Scheibe
 116 RPP 20.8651 29.0951 -2.8 2.8 -0.75 0.75 $ Alu Tr�ger
 117 RPP 20.8651 29.0951 -2.3 2.3 -0.75 0.75 $ Inner Alu Tr�ger
 118 RCC 29.0951 -0 0 0.4 0 0 2.95     $ Alu Scheibe
 119 RCC 20.4651 0 0 10.05 0 0 3.25    $ Hintere Haube
 120 RCC 20.4651 0 0 9.75 0 0 3        $ Inner hintere Haube
 121 RCC -3.2745 0 0 21.8396 0 0 3.25  $ Vordere Haube             
 122 RCC -3.1745 0 0 21.7396 0 0 3     $ Inner vordere Haube
 123 RCC -3.1745 0 0 1 0 0 3           $ PE Kappe  3,5mm
 124 RCC -3.1745 0 0 1 0 0 2.55        $ Inner PE Kappe
 125 RCC -2.8245 0 0 0.65 0 0 2.8245   $ Inner PE Kappe
 126 RPP 6.2 10.7 1 2 -1 1             $ Dynode
 127 RPP -100 100 -25.5 -22.5 -50 50   $ Tisch
c 128 RPP -140 -130 -100 205 -50 50     $ W�nde
 129 RPP -100 100 -100 -90 -100 100    $ Boden
 130 RPP 20 24 -6.5 -5 -13.25 13.25    $ Tr�ger Lineal
 131 RCC 22 -16.5 -11.25 0 10 0 2      $ Tr�ger Bein Links
 132 RCC 22 -16.5 11.25 0 10 0 2       $ Tr�ger Bein Recht
 133 RPP 20.5 23.5 -5 5 -7.25 7.25     $ Tr�ger
 134 RCC 0 -16.5 0 0 -3 0 40           $ �bere Drehtisch
 135 RCC 0 -19.5 0 0 -3 0 45           $ Untere Drehtisch
 136 RCC 0 -16.5 0 0 -6 0 3            $ Drehtisch Achse
 137 RPP -100 -98 -30.5 -25.5 -50 50   $ Tisch Struktur 1
 138 RPP 98 100 -30.5 -25.5 -50 50     $ Tisch Struktur 2
 139 RPP -98 98 -28.5 -25.5 -1 1       $ Tisch Struktur 3
 140 RCC 22 -8.5 -5 0 16 0 .75         $ Schrauben 1
 141 RCC 22 -8.5 5 0 16 0 .75          $ Schrauben 2
 142 RCC 6.04 0 0 .06 0 0 2            $ Kollimator
 143 RPP 6.2 7.2 -.75 .75 -1 1         $ 1. Dynode
 144 RCC -2.5585 0 0 5.0985 0 0 2.5585 $ Luft zwischen kristall und reflektor
 145 CX 2                              $ Ring in Sockel
 146 CX 2.5                            $ Ring in Sockel au�er
 147 RCC 30.5151 0 0 .86 0 0 .95       $ Buchse
 148 RCC -50.01 0 0 .02 0 0 3.51       $ Quelle
 149 RCC -50.025 0 0 .05 0 0 3.51      $ Papier um die Quelle
 150 RCC -50.0365 0 0 .075 0 0 4       $ Plastik um die Quelle
 151 RPP -100 100 195 200 -100 100     $ Dach
c BLOCK 3:DATA CARDS ******************************************
SDEF erg=d1 par=2 pos=-50 0 0 axs=1 0 0 ext=d3 rad=d2 WGT=4.31e8
SI1 L .66166
SP1 D 1.2e5
si2 0 3.51
sp2 -21 1
si3 0 0.05
sp3 -21 0
mode p
Phys:p 1 1 0 0 1 0
CUT:p j .01
nps  1e8
c
f18:p 11
e18 .06 861i .8
ft18 geb -.01 .07 -.08
c
m1   11000 0.45 $ Na
     53000 0.45 $ I
     81000 0.1  $ Tl (0,1% mole)
m2   13000 -0.97     $ Al (AlMg3)
     12000 -0.03     $ Mg
m3   5000 -0.040064  $ Glass Borosilicate
     8000 -0.539562
     11000 -0.02819
     13000 -0.011644
     14000 -0.377220
     19000 -0.00332
m4   28000 -0.77     $ Ni (Mu Metall - NiFe15Cu5Mo3)#(-111 112 113)#(-114 115)#(-116 117)
     26000 -0.15     $ Fe
     29000 -0.05     $ Cu
     42000 -0.03     $ Mo
m5   1000 -0.143716  $ PE ohne Bore
     6000 -0.856284
m6   29000 -0.58 $ CuZn39Pb3
     30000 -0.39
     82000 -0.03
m7   13000 -0.92 $ AlMgSi
     14000 -0.04
     12000 -0.04
m8   1000 -0.041960 $ PET (1,38)
     6000 -0.625016
     8000 -0.333024
m9   6000 -.005 $ Stahl
     26000 -.995
m10  1000 -.059642 $ Holz 0,75
     6000 -.497018
     7000 -.00497
     8000 -.427435
     12000 -.001988
     16000 -.004970
     19000 -.001988
     20000 -.001988
m11  1000 -.00453 $ Concrete MCNP
     8000 -.5126
     11000 -.01527
     13000 -.03555
     14000 -.36036
     18000 -.125
     20000 -.05791
     26000 -.01378
m12  8000  -.39696 $ Reflektor aus MgO
     12000 -.60304
m13  26000 -.000868 $ Messing
     29000 -.665381
     30000 -.325697
     50000 -.002672
     82000 -.005377
m14  29000 -1      $ Kupfer
m15  6000 -.000124 $ Luft
     7000 -.755268
     8000 -.231781
     18000 -.012827
m16  1000 -.25 $ Epoxy
     6000 -.625
     8000 -.125
m17  1000 -.037764
     6000 -.5625144
     8000 -.2997216
     55000 -.1 $ Cesium

J'ai trouvé!!!! Trop de Tallium dans mon cristal!
Merci tout de meme pour les pistes de reflections

xln10 a écrit:

J'ai trouvé!!!! Trop de Tallium dans mon cristal!
Merci tout de meme pour les pistes de reflections

Bonjour,
Cela doit jouer (perso je ne ne sais pas combien il y a de Tl dans un NaI).
Néanmoins,  le plus étonnant c'est que vous ayez simulé un spectro de 862 canaux fonctionnant entre 60 keV et 800 keV. En principe (sauf à simuler un spectro de votre invention ;-)) vous auriez dû simuler un 1024 canaux entre 0 et 3072 keV. Cela a une importance colossale car vos bandes d'énergie simulées sont différentes de celles mesurées donc vous n'avez pas les mêmes résultats.

Bonjour,
dans mon cas, il y a 0,1 %/mol de Tallium, ma Simulation en comportait 10%!
Pour ce qui est de mes canaux, j'ai d'abord calibré mes mesures Avec differentes sources. Jai ensuite déduis une loi linéaire reliant chaque canal à une énergie et je me suis focalisé sur l'intervale qui m'interesse : 60keV-800keV (bien plus pratique pour exporter mes données Sous Excel)
Par contre je ne comprend pas votre 3072 keV.

Bonjour,
J'ai pris une électronique classique de spectro, c'est à dire un 1024 canaux de 0 à 3072 keV. Si vous voulez comparer vos résultats de  mesure par rapport à la simulation vous devez impérativement avoir la même largeur de bande en énergie.  Sinon il faut faire une comparaison en imp/keV (vous divisez votre nombre d'évènements par la largeur du canal.
Bonne continuation

Bonjour.

Après avoir parcouru le forum à de nombreuses reprises pour m'aider à m'initier à MCNP (merci pour votre aide !), je m'inscris enfin suite à un problème dont je n'arrive pas à sortir !

Je cherche à modéliser une source de Pu américiée afin d'estimer les débits de dose neutron et gamma en un point. Lorsque je fais tourner MCNP sur la partie photon de ma problématique, je score très peu, mes photons étant essentiellement stoppé dans ma source de Pu (un orthocylindre).

Est-ce quelqu'un aurait une idée sur une méthode me permettant d'améliorer ma convergence sans fausser le résultat de mon débit de dose ? Je me demandais par exemple si je pourrais sourcer uniquement en surface de ma source (quelle épaisseur serait alors maximale ?) pour n'avoir que des photons produits par MCNP capables de sortir de la source avec une probabilité suffisante.

Merci d'avance pour votre aide et vos réponses.

Bien cordialement,

Bonjour,

pour la dose générée par le 59 keV de l'Am241, le libre parcours moyen (lpm) est égal à 66 µm (seulement)
Prenez une couche surfacique de 5 lpm, soit 330 µm, ce sera suffisant pour estimer la dose

Au delà de 5 lpm, les photons qui tirent droit vers la sortie ont moins de 1 %  de chances de sortir. Donc inutile de les générer

Si votre source fait seulement 1 cm sur l'axe de calcul, vous gagnez un facteur 10000/330=30 de temps de calcul avec cette méthode

Gluonmou

Merci beaucoup pour cette réponse. Je vais essayer tout de suite et vous tient au courant !

Je viens de recevoir un PC avec MCNPX pour un mois. Le temps de faire un calcul.
Mon premier problème est l'installation. 
J'avais jusqu'à présent utilisé des ordinateurs avec MCNP déjà installé.
J'ai essayé de lire et suivre les consignes du manuel mais c'est pas gagné.

Déjà apparemment je dois télécharger un compilateur FORTRAN et C.
Ensuite je dois dézziper le dossier v25f.tar.gz mais je n'ai pas ce fichier. J'ai le : MCNPX_v250_EXE.tar.gz à l'interrieur se trouve le fichier MCNPX_v250_EXE.tar et à l'intérieur de celui-ci se trouve le dossier v250 à l'intérieur se trouve : Data, Docs, exe_dir et autres.

Puis je dois lancer cette commande :
CD MCNPX_BLD\MCNPX\configure

Mais je n'ai pas de dossier MCNPX_BLD et pas de fichier configure.

Pouvez-vous me renseigner ? 
Merci beaucoup !

Bonjour Samael22,

Bizarre, votre problème me fait penser penser à chante-sloubi :
http://kaamelott.co/livre-3/perceval-chante-sloubi/


Il y a dans votre post comme qui dirait du comique involontaire.
Ou pas?

On a ouvert un thème là -dessus récemment. Coïncidence?

Sinon si c'est pour de vrai ben bon courage quand même.

moi aussi il y a quelque années je me suis tapé le mode d'emploi de MCNP
Pour des figues...

Du coup avec Gamma-lambda on a pondu Dosimex.

heureusement que H*(10) n'et plus là pour voir ça; peuchère.

Bon c'est pas tout ça, j'ai un flan de courgettes à attaquer

En tout cas se regarder chante-sloubi à 7h10 avant de commencer à travailler et galérer sur MCNP ça fait plaisir.

J'ai un élément de réponse. Je n'avais pas tous les fichiers.
Il me manque plus que les droits admin et je pourrai l'installer.

ça, les droits admin, c'est encore autre chose, un autre monde..
Re-bon courage

Bonjour.

Nouvelle année, nouveau problème !

Donc tout d'abord, meilleurs voeux à tous !

Ensuite, je poste sur le forum pour un problème de géométrie dans MCNP.

J'ai besoin de modéliser 2 réseaux imbriqués pour modéliser une petite matrice d'objets qui se répètent, petite matrice qui se répète elle-même dans une plus grande matrice.

J'ai modélisé mon environnement et la géométrie "fonctionne" (aucun cooky cutter, la géométrie correspond à ce dont j'ai besoin). Par contre, lorsque j'essaie de distribuer ma source dans toutes les mailles de l'imbrication des matrices : "zero lattice element hit" ! Je suis reparti dans ma matrice de niveau 2, elle fonctionne. C'est lorsque je la mets dans la matrice de niveau 1 que les particules se perdent...

Quelqu'un pourrait-il m'aider sur les types de problèmes qui pourraient être à l'origine de mon blocage ?

Merci d'avance.

Bonne journée.

Bonjour,

je suis en Master 2 Radioprotection et j'aimerais savoir si vous auriez par hasard des TD et/ou des cours sur l'initiation au calcul Monte- Carlo. 

Merci,

Bien à vous.