MCNP : the Godfather

Bonjour,
pour répondre à lail voir : https://www.forum-rpcirkus.com/t4259-l-energie-depose#48413
voici le fichier MCNP* permettant le calcul du kerma à un mètre d'une source de 1,17 et 1,33 MeV émettant un photon. Pour le calcul de lail il faut normaliser tout ça. A noter le résultat est en jerk ( https://www.youtube.com/watch?v=Xyzi28ckyU8 **) soit en GJ/g

c fichier cirkus 2
1 1 -0.00129 -1 2
2 1 -0.00129 -2
999 0 1

1 so 150
2 sx 100 .5

M1 7014 .781 8016 .212 7015 .0028 $ air
cut:e j .01 3j
mode p e
imp:p 1  1 0
imp:e 0 1 0
sdef pos 0 0 0 par=2 erg=d1
si1 L 1.17 1.33
Sp1 1 1
*f6:p,e 2
ctme 15

n'hésitez pas si vous avez des questions.

nota * : c'est luxueux d'utiliser MCNP pour ce type de calcul mais bon ...
nota ** : c'est pour que le post soit moins austère (et non pas pour rivaliser avec l'humour avec G...u)

Bonjour,
j'essaye de m'initier à mcnp.
Je veux juste dans un premier temps faire  un calcul simple de débit de dose à 10cm d'une source de Cs137.
Aprés avoir lu 2 ou 3 truc, j'ai créer un fichier input mais qui me sort des fatal error...Pourriez-vous m'aiguiller sur ce qui ne va pas dans mon raisonnement:
Voici ce que je pensais faire: Je me suis définit une aire de jeu cubique de 1m³ dedans on calcul dehors on ne calcul rien, l'interieur est rempli d'air.
Dans cette aire au centre il y a une source de Cs137
Pour calculer le débit de dose à 10 cm j'ai positionné un cylindre de 0,5cm de rayon et 1cm de hauteur selon l'axe OX. Comme je voudrai un débit de dose j'ai mis de l'eau assimilable au tissus mou en première approche. Je ne sais déjà pas si ceci est la bonne méthode pour définir un débit de dose efficace.
S'agissant d'un débit de dose j'ai choisis un tally en F6, et je pulse 10^5 photons pour commencer bref voici ce que je propose en fichier input:

C TEST CS137
10 0 100 IMP:P=0
11 1 -0.0013 -100 101 IMP:P=1
12 2 -1 -100 -101 IMP:P=1

100 RPP -50.00 50.00 -50.00 50.00 -50.00 50.00
101 RCC 10.00 0.00 0.00 1.00 0.00 0.00 0.5

MODE P
M1 7014 -0.80 8016 -0.2
M2 1000 -0.11 8016 -0.89
SDEF POS 0.0 0.0 0.0 PAR=2 ERG=D1
SI1 L 0.662
SP1 D 0.852
F6:P 101
NPS 100000

En lançant mcnp j'ai:
fatal error.  f card cell   101 of bin   1 tally     6 not found.

Je ne trouve pas dans le fichier de sortie de résultat, je ne vois d'ailleurs pas trop ou cela se situe...vu la quantité de données qui sont transmises.
Pourriez-vous me mettre sur la voie de mon erreur s'il vous plait (version utilisée 4c3)
Merci

Bonjour,
La fatal error vient du f6:p 101 car 101 n'est pas une cellule mais un élément de géométrie.  Il eut fallut (pas sûr du t à la fin de fallu) mettre 12. Je regarde le reste voir si il y a d'autres erreurs moins évidentes.
GA

Bonjour,
Merci Ga de cette réponse ultra rapide,
logique on calcul l'énergie déposée dans une cellule et non un volume! Trop fort ce GA.
Pour le "t" je suis super nul en orthographe alors avec ou sans pour moi c'est pareil...
Encore merci.

Autre remarque, je ne suis pas sur que la version 4c3 "accepte" les macrobodies (RPP et RCC), contrairement aux versions X et 5. Si ce n'est pas le cas il vous faudra construire vos cubes à partir dintersection de plan type PX PY PZ.

Re-bonjour GA,
En effet cela ne fonctionnais toujours pas j'ai fait un test avec des sphères et cela fonctionne ceci est donc certainement lié aux macrobodies.
Merci.

Gamma.Lambda a écrit:

j'ai fait un test avec des sphères

Solution de facilité.  En 4c3 il eût fallu faire (désolé pour le flood klax mais comme on ne peut pas insérer de txt ...) :
C TEST CS137
10 0 -100:110:-120:130:-140:150 IMP:P=0
11 1 -0.0013 100 -110 120 -130 140 -150 #12 IMP:P=1
12 2 -1  -200 210 -220 IMP:P=1 vol .7854

100 PX -50
110 PX 50
120 PY -50
130 PY 50
140 PZ -50
150 PZ 50
200 CX .5
210 PX 10
220 PX 11

MODE P
M1 7014 -0.80 8016 -0.2
M2 1000 -0.11 8016 -0.89
SDEF POS 0 0 0 PAR=2 ERG=.662
c SI1 L .662
c SP1 .852
*F6:P 12
ctme 20

Encore merci,
finalement les cylindres fonctionnent aussi sous cette version, mais pas les cubes.
Bonne soirée.

Bonjour, mejuhrci bcp pour votre aide, ce que j'ai compris, c'est que cela permet de calculer le débit de kerma à un métre, dans une sphère de rayon 100,5 (2 sx 100 .5); 
jé pas compris tt ca
999 0 1???
cut:e j .01 3j??
ctme 15??
mode p e??

Bonsoir,

lail a écrit:

j'ai compris, c'est que cela permet de calculer le débit de kerma à un métre, dans une sphère de rayon 100,5

Non en fait sx 100 0.5 est une sphère sur l'axe ox à la position 100, 0, 0 qui a un rayon de 0, 5 cm. Donc je calcule le kerma à 100 cm dans une sphère de 0, 5 cm.

Citation :

999 0 1???

Signifie que la cellule 999 est constituée de vide et à l'extérieur de la géométrie 1 (sphère centrée en 0, 0, 0 et de rayon 150 cm)

Citation :

cut:e j .01 3j??

Permet de ne pas transporter les electrons de moins de 10 keV (le transport des électrons est très gourmand en temps et empêche la convergence du code

Citation :

ctme 15??

Je fais tourner le code 15 minutes (au lieu d'utiliser un nps qui transporte un nombre de particules)

Citation :

mode p e??

On transporte les photons et les électrons

Voilou

Bonjour,j'ai un petit problème que je n'arrive pas à saisir,  que signifie un mode de transport de particule; est ce qu'il est lié à la désintégration ? ou au milieu interagissant avec le rayonnement?

pourquoi le mode de transport pour le césium 137 est mode p [photon], et pour le cobalt c'est mode pe[photon électron];

Bonjour
S'agit-il d'une discussion en référence au code MCNP ?
Si oui les modérateurs pourraient déplacer le sujet.
Parce que "mode de transport" ne me parle pas d'emblée.
KLOUG

Sujet déplacé par Scully   

Bonjour, justement c'est un problème que j'ai rencontré en essayons de faire des modélisations simple avec le code " pour faire rappeler que je suis débutante à travailler avec"; 
Merci

Bonsoir,
Je vais essayer de répondre à lail  sur le post https://www.forum-rpcirkus.com/t4270-transport-de-particule#48472

lail a écrit:

Bonjour,j'ai un petit problème que je n'arrive pas à saisir,  que signifie un mode de transport de particule; est ce qu'il est lié à la désintégration ? ou au milieu interagissant avec le rayonnement? pourquoi le mode de transport pour le césium 137 est mode p [photon], et pour le cobalt c'est mode pe[photon électron];

Tout d'abord mcnp (pour faire très simple) est ce que l'on appelle un code de transport (au sens propre, du terme), c'est à, dire qu'il transporte des particules d'un point A à un point B et si les particules rencontrent des matériaux,  elles interagissent avec celui-ci. Je suppose que vous avez des bases solides de physique (notamment en IRM -sinon l'utilisation de ce type de code est déconseillé par la médecine légale).  En interagissant ces particules primaires engendrent d'autres particules (selon les processus classique de l'IRM) qui elles mêmes sont transportées et peuvent à leur tour interagir et engendrer des particules ... Le processus s' arrête lorsque les particules sortent du domaine (dans notre cas la cellule 999). Ainsi un mode p va transporter des photons et lors d'une interaction de type photoelectrique il considere l'electron à l'endroit de l'interaction (il peut meme créer le photon de freinage associé) mais ne transporte pas l'électron.  Par contre un mode p e va faire que l'electron va etre transporté dans le milieu (avec toute la physique lourde du transport des électrons et les temps de calcul associés). Ainsi bien qu'un mode p e soit physiquement plus "propre" il est beaucoup plus gourmand en temps qu'un mode p qui, certes moins pointu, reste suffisant pour les calculs.
Pour les forts en thème je suis sur qu'ils ont déjà intuité qu'avec cette subtilité on peut soit calculer un kerma soit une dose.
En esperant ne pas avoir été trop abscons je vais me refaire la saga starwars.
GA

Bonsoir,
J'ai essayé de répondre ici : https://www.forum-rpcirkus.com/t4128p30-mcnp-the-godfather#48480

Bonjour, 
étant le seul forum  en français sur MCNP( comme je suis super nulle en anglais), j pense  que je vais continuer à posé de questions jusqu'à ce que vous arrêtez de me répondre ; 
alors voilà, j'ai comme même des erreurs dans mon fichier de sortie[ calcul du kerma à un  mètre ]

 fatal error.  tally 6 is not available for electrons. 


Merci

Bonjour lail,
je suis bien embêté car votre fichier fonctionne très bien sur mon MCNP ! peut être il y a t'il des caractères non ascii dans votre fichier ASCII (je ne peux pas le vérifier sans votre fichier d'entrée) ou votre MCNP est mal installé (quelle version utilsez vous ?)
Je vais essayer d'envoyer mes fichiers d'entrée et de sortie à Klax voir si vous pouvez les télécharger (sinon cela fait des post interminables et trop de flood agace klax)
En attendant essayer :
c fichier cirkus 2
1 1 -0.00129 -1 2
2 1 -0.00129 -2
999 0 1

1 so 150
2 sx 100 .5

M1 7014 .781 8016 .212 7015 .0028 $ air
mode p
imp:p 1 1 0
sdef pos 0 0 0 par=2 erg=d1
si1 L 1.17 1.33
Sp1 1 1
*f6:p 2
ctme 15

J'ai juste supprimer le transport des électrons, ce qui ne change rien dans votre cas puisque vous êtes à l'équilibre électronique et K=D.

et voici les fichiers...
klax

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