| MCNP : the Godfather | |
|
+4SheepeR Gluonmou GA_mcnp (x) Klax 8 participants |
|
Auteur | Message |
---|
Klax Homme-canon
![Homme-canon Homme-canon](https://i.servimg.com/u/f49/16/29/75/03/1110.gif)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: MCNP : the Godfather Sam 19 Oct 2013 - 18:17 | |
| 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 _________________ J'aime pas les matins... ![MCNP : the Godfather Zzzzz-9758](https://2img.net/h/smileys.sur-la-toile.com/repository/Endormi/zzzzz-9758.gif) |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Sam 19 Oct 2013 - 18:51 | |
| 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 ...
|
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Sam 19 Oct 2013 - 19:08 | |
| Autre precision, des calculs mcnp pourraient être envisagés, si, évidemment, ils ne sont pas à but commerciaux, GA |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Gluonmou Contorsionniste
![Contorsionniste Contorsionniste](https://i.servimg.com/u/f49/16/29/75/03/1010.gif)
![Gluonmou](https://2img.net/u/2811/10/03/89/avatars/1761-26.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Sam 19 Oct 2013 - 19:27 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Sam 19 Oct 2013 - 19:55 | |
| Bonsoir Vive les codes de calculs... KLOUG |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Sam 19 Oct 2013 - 20:46 | |
| - 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: publications utilisant MCNP : protonthérapie Dim 20 Oct 2013 - 9:12 | |
| 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=S0033845112000543Les 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
SheepeR Trapéziste
![Trapéziste Trapéziste](https://i.servimg.com/u/f49/16/29/75/03/810.gif)
![SheepeR](https://2img.net/u/2811/10/03/89/avatars/1211-51.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Lun 21 Oct 2013 - 9:26 | |
| 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.... |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Lun 21 Oct 2013 - 11:54 | |
| Bonjour, Malheureusement, il n'existe pas de version test de ce code, GA |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Mar 22 Oct 2013 - 15:16 | |
| 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.aspxSinon 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. |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Mar 22 Oct 2013 - 19:20 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Mer 23 Oct 2013 - 20:03 | |
| 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-mammographiej'ai fait quelques calculs MCNP. on trouve alors pour une HT de 50 kV : anode W calcul à 1m (dans le vide) ![MCNP : the Godfather Gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=10~mSv.min^{^{-1}}](http://latex.codecogs.com/gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=10~mSv.min^{^{-1}}.mA^{-1}) anode Mo calcul à 1m (dans le vide) ![MCNP : the Godfather Gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=13~mSv.min^{^{-1}}](http://latex.codecogs.com/gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=13~mSv.min^{^{-1}}.mA^{-1}) anode Rh calcul à 1m (dans le vide) ![MCNP : the Godfather Gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=12~mSv.min^{^{-1}}](http://latex.codecogs.com/gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=12~mSv.min^{^{-1}}.mA^{-1}) anode Cu calcul à 1m (dans le vide) ![MCNP : the Godfather Gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=4,6~mSv.min^{^{-1}}](http://latex.codecogs.com/gif.latex?\dot{H}^{*}\left&space;(&space;10&space;\right&space;)=4,6~mSv.min^{^{-1}}.mA^{-1}) 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Fred Homme-canon
![Homme-canon Homme-canon](https://i.servimg.com/u/f49/16/29/75/03/1110.gif)
![Fred](https://2img.net/u/2811/10/03/89/avatars/72-35.gif)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Mer 23 Oct 2013 - 20:34 | |
| Bonjour, Effectivement la HT est de l'ordre de 28 / 30 kV. _________________ Fred I IRSN |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Gluonmou Contorsionniste
![Contorsionniste Contorsionniste](https://i.servimg.com/u/f49/16/29/75/03/1010.gif)
![Gluonmou](https://2img.net/u/2811/10/03/89/avatars/1761-26.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Dim 10 Nov 2013 - 11:53 | |
| 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? |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: calcul de 'énérgie déposé en MeV par mcnp Ven 20 Déc 2013 - 16:19 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 19:37 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 20:08 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 20:21 | |
| 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 ... |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 20:33 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 21:39 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 22:13 | |
| 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 22:19 | |
| j'ai pas compris ca c'est un intervalle d'énergie? erg=d1 si1 L 1.17 1.33 Sp1 1 1
????
Merci |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 22:33 | |
| 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. |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
Invité Invité
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | |
![Aller en bas](https://2img.net/i/empty.gif) | |
GA_mcnp (x) Ventriloque
![Ventriloque Ventriloque](https://i.servimg.com/u/f49/16/29/75/03/510.gif)
![GA_mcnp (x)](https://2img.net/u/2811/10/03/89/avatars/2097-89.jpg)
![MCNP : the Godfather Empty](https://2img.net/i/empty.gif) | Sujet: Re: MCNP : the Godfather Ven 20 Déc 2013 - 23:09 | |
| - 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 |
|
![Aller en bas](https://2img.net/i/empty.gif) | |
| MCNP : the Godfather | |
|