MCNP : the Godfather

Bonjour,

on continue sur les topics uniques sur les codes de calcul. Il n'existe pas de forum français sur MCNP. Donc maintenant, si.

Comme les autres sujets, au besoin on étoffera le topic en plusieurs fil de discussion.

Si des codeurs veulent animer ce forum, ils sont les bienvenues.

A vous de jouer maintenant.

klax

Bonjour,
Klax m'a demandé de vous rejoindre pour "animer" ce post sur mcnp (x). Force est de constater que les forums mcnp francophones ne sont pas légion et qu'il nous semble intéressant que votre forum, favori puisse accueillir un "essayons d'approcher la bete". Tous ceux qui ont débuté sur la chose se sont forcement vus affublés d'un "fatal error" et seul devant l'adversité vous y avez passé des nuits. C'est pour vous que ce post est fait ! Que vous soyez débutant (on l'a tous été) ou expert de classe international, on vous y attend. Venez parler des vos fichiers recalcitrants, aidez les nouveaux à débuter, parler, des publications qui associent mcnp et la radio ...
A vos claviers ...

Autre precision, des calculs mcnp pourraient être envisagés, si, évidemment, ils ne sont pas à but commerciaux,
GA

Bonjour,

Bonne initiative. Car malgré ce que pourrait laisser entendre certains post, pour ma part MCNP c'est "full respect". The code

L'un des premiers codes scientifiques qui ait tourné sur les premiers ordinateurs (ENIAC). Conçus par les plus grands phycisiens de l'époque.

On peut noter que pour beaucoup (j'ai fait ma petite enquète), très souvent ce sont les résultats de MCNP qui valident les mesures, et pas l'inverse.

Par contre effectivement pour les débutants il est très difficile , voire pour moi impossible, de s'y former tout seul, malgré les nombreux manuels qui existent.

J'ai eu la chance de pouvoir me faire payer une formation perso de 2 jours sur la bête, avec un Paganini de la chose (bises).Après ça, j'ai pratiqué un certain temps un peu tous les jours.  J'étais à deux doigts de pouvoir enlever mes petites roulettes et d'avancer tout seul sur mon petit vélo. mais là, paf, d'autres occupations m'ont aspirée à 100 %. lorsque 2 ans après j'ai voulu recommencer, j'étais incapable de seulement lancer un calcul. De désespoir j'ai définitivement renoncé.

Donc un conseil d'un 'ex-éternel débutant : si vous commencez, n'arrétez pas, et si vous en avez l'occasion, commencez jeune. C'est comme pour le piano ou le Mandarin, c'est plus fastoch de commencer petit

Gluonmou, ever old

Bonsoir
Vive les codes de calculs...  
KLOUG

Gluonmou a écrit:

 J'étais à deux doigts de pouvoir enlever mes petites roulettes et d'avancer tout seul sur mon petit vélo. mais là, paf, d'autres occupations m'ont aspirée à 100 %. lorsque 2 ans après j'ai voulu recommencer, j'étais incapable de seulement lancer un calcul. De désespoir j'ai définitivement renoncé

Bon résumé de ce que nous aimerions faire : aider à commencer puis vous vous lancez
GA

Bonjour,
dans ce nouveau post consacré à MCNP, je propose, afin de mieux connaitre ses capacités et de faire une "veille" sur les publications scientifiques utilisant ce code.
Je commence donc par le dernier numéro de "radioprotection" et la publication :
Monte Carlo simulation of a proton therapy beamline for intracranial treatments
http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=8995224&fulltextType=RA&fileId=S0033845112000543
Les auteurs ont simulé l'installation du centre de protontherapie d'Orsay, utilisant des protons de 180 MeV pour les traitements pédiatriques de tumeurs intracranials. Des calculs de la distribution de la dose dans un fantôme d'eau sont effectués et comparés à des mesures. De même les auteurs effectuent des calculs de la dose neutron sur le patient.
langue : anglais
difficulté : ++
mon point de vue : la description des paramètres pour une simulation de protons dans MCNPX est très bien decrite ce qui permet d'éviter certain piège de ce type de simulation.
GA

Bonjour,

Une question pour vous: ce logiciel a-t-il une version d'essai pour suivre éventuellement ce tuto/post? je pense que non mais bon sait-on jamais....

Bonjour,
Malheureusement, il n'existe pas de version test de ce code,
GA

mmm...en effet tous les deux ans (plus ou moins) une nouvelle version beta de MCNP/MCNPX est sortie. Lors de l’apparition de la nouvelle version beta, l'ancienne version du code est fournie.

Si vous êtes capables d'attendre le bon moment, la dernière version beta + l'avant-dernière version sont distribuées gratuitement pour une période limitée de l’année ici:

https://rsicc.ornl.gov/Default.aspx

Sinon le prix du code est d'environ 800 $.

Avec le code il est également distribué VisEd, logiciel qui permet de visualiser les géométries de MCNP.

Bonsoir,
 J'ai également oublié que l'on pouvait avoir le code gratuit en "beta testeur".  Tout est expliqué là : https://mcnpx.lanl.gov/ :

"MCNPX (source code, executables, data) is available from this WWW site to "beta testers" who have access to intermediate code versions. Beta versions of MCNPX are available from "The Code" tab at the left of this web site. Beta Testers are sponsors, collaborators, and those who take MCNPX workshops (see "Classes tab on the left.) For further information on the Beta Test program, contact mcnpx@lanl.gov."

Vous pouvez donc essayer d'envoyer un courriel à mcnpx@lanl.gov. Mais je doute que ce soit très simple pour entrer dans le beta test.
Noter que si vous trouvez un bugg dans mcnpx vous avez le  droit à 20 $ voir https://mcnpx.lanl.gov/opendocs/misc/Cashbugs.doc . 100% des gagnants ont tenté leur chance  
GA

Bonsoir ,
pour répondre au post concernant les doses dans le cas d'un appareil de mammographie : https://www.forum-rpcirkus.com/t4088-valeur-du-rendement-du-tube-dans-le-cas-du-examen-mammographie
j'ai fait quelques calculs MCNP.
on trouve alors pour une HT de 50 kV :

anode W calcul à 1m (dans le vide)


anode Mo calcul à 1m (dans le vide)


anode Rh calcul à 1m (dans le vide)


anode Cu calcul à 1m (dans le vide)


Pour le calcul j'ai considéré une pente d'anode à 45°. Je n'ai pas considéré de filtration (si vous avez des idées sur le sujet n'hésitez pas).
A noter également qu'il me semble que les appareils de mammographie ont des HT très faibles de l'ordre de 25 kV, dans ce cas le calcul peut être complexe car le spectre s'étend alors en deçà de 10 keV, région où l'ICRU/ICRP ne s'aventure pas !
Seule solution recalculer les facteurs de dose à ces faibles énergies.

Notez que vous avez appris ou lu que la dose dû au freinage varie linéairement avec le Z du matériau. Or on a pour le W Z=74, le Rh Z=45, le Mo Z=42, le Cu Z=29.
Etonnant non ? En fait cela fonctionne avec le cuivre et le tungstène et non avec Rh et Mo. Le soucis avec le deux derniers sont les raies caractéristiques autour de 20 keV pour le Rh et 17 et 19 keV qui perturbent cette régle.
GA

Bonjour,
Effectivement la HT est de l'ordre de 28 / 30 kV.

Cher Ga-mcnp(x)


Par curiosité curiosante, pourriez-vous svp calculer le rendement en H*(10) avec du Mo à 100 keV et 2 mm d'alu?

Prenons le pari que l'on devrait alors retrouver la règle en Z .

Et effectivement si quelqu'un a une idée de la filtration utilisée en mammographie. Il doit bien y en avoir une, pour filtrer les très basses énergies et éviter de décoller la peau

Merci encore

Gluonmou


PS : Super. Gamma.lambda a rejoint sa base polaire secrète pour coder les changements d'anodes. Mais où s’arrêtera-t-il?

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