MCNP : the Godfather

bjr, je viens d'installer le code mcnp, et je dois calculer le débit de dose avec le code; sachant que je né jamas travailler avec je suis arriver à faire cela
c     Created on: Thursday, December 19, 2013 at 12:35
    1     1       1 -1
    2     0         1
 
1            so 20
 
 mode p e
 
Sdef pos=-450 0 -100  par=2  erg=1.25
 f8: p, e
nps 1000

Bonsoir,
Comment dire, à part une "fatal" error" je ne vois pas bien ce que peu donner le code ....
Par exemple vous avez limité l'espace  à une sphère de 20 cm de rayon et émettez des particules en dehors, le f8 ne fait pas référence à une cellule, le matériau n'est pas décrit ....


PS pour les moderateurs peut on déplacer le post dans la bonne section :

https://www.forum-rpcirkus.com/t4128-mcnp-the-godfather

Bonsoir,

Merci pour votre réponse, je vien d'avoir le code il ya peu de temps, j'ai fait quelque modifications mé ca marche tjrs pas
c     Created on: Friday, December 20, 2013 at 18:53
  
  1     1 -0.00129 -1 
  
  2     0         1 

 
  
 1        so 20


 

m1    7000.02p      -0.755268  $l'air 
    
 8000.02p      -0.231781 18000.02p     -0.012827 6000.02p      -0.000124
mode p e
imp:p 1 0
imp:e 1 0
sdef pos 1 1 1  par=2  erg=1.25
f6:p,e  1
nps 100000

Bonsoir,
Cela devrait mieux marcher mais il faut respecter les "lignes vides" 1 seule entre les cellules et la géométrie et une entre la géométrie et les cards. Neanmoins même si mcnp donne un résultat je ne vois pas ce que vous calculez ...

Merci encore pour votre réponse, 

Au fait, je viens d'avoir le code par mon encadrant, lui qui connait rien sur le mcnp( tout comme moi), me demande de calculer l'énergie déposé par une source de cobalt 60 ponctuelle et isotrope dans une sphère de rayon 20, et comparer avec le calcul expérimental; 

Merci beaucoup

Re,
Etonnant d'encadrer quelqu'un sur un sujet non maîtrisé !
Je ne sais pas insérer un fichier texte donc je fais un copier/coller (attention faut faire du ascii avec mcnp). Voici un fichier qui fonctionne (testé à l'instant) et un peu plus "propre"

C fichier cirkus
1 1 -0.00129 -1
2 0 1

1 so 20

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

Bonsoir, merci beaucoup, ça marche très très bien merci beaucoup, j'ai deux autres questions si vous permettez  ( erg=d1 ); 


quelle est la valeurs de l'énergie dans votre fichier de sortie


j’espère que je ne vous dérange pas

Re,
Erg =d1, signifie que la donnée energie se trouve dans le drapeau 1. Ici si1. Si1 L signifie que les énergies sont données sous forme discrete : ici 2 énergies de 1, 17 et 1, 33 MeV (énergie du Co60). Enfin sp1 donne les intensités d'émission ici 1 et 1 : donc les énergies de 1, 17 et 1, 33 MeV seront émis à 100%.
Restera juste à bien normaliser le résultat.

lail a écrit:

quelle est la valeurs de l'énergie dans votre fichier de sortie

Je n'avais pas vu la question : je trouve 2, 095E-5 MeV/g/photon soit 4,19E-5 MeV/g/Bq

Citation :

je vous remercie infiniment pour votre aide

you're welcome

GA_mcnp (x) a écrit:

je trouve 2, 095E-5 MeV/g/photon soit 4,19E-5 MeV/g/Bq

Bien entendu il faut lire

2, 095E-5 MeV/g/photon soit 4,19E-5 MeV/g/Bq/s. Faut que j'arrête de poster si tard.

Bsr, je viens de voir votre message, merci infiniment 

Bonsoir, existe t-il une formation de courte durée sur le code mcnp en France
Merci

Bonjour,
Le LANL ayant bien bordé son affaire les seules formations mcnp autorisées sont celles dispensées par le LANL. Par chance, de temps en temps il y en a en Europe, et même en France : en 2014 elle a lieu du 31 mars au 4 avril 2014. L'inscription est ici : https://mcnpx.lanl.gov/ puis cliquer sur "classes"

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