Forum technique de RadioProtection Cirkus

Le portail de la RadioProtection pratique et opérationnelle - www.rpcirkus.org
 
RP CirkusRP Cirkus  AccueilAccueil  FAQFAQ  RechercherRechercher  S'enregistrerS'enregistrer  Connexion  
Partagez | 
 

 MCNP : the Godfather

Voir le sujet précédent Voir le sujet suivant Aller en bas 
Aller à la page : 1, 2, 3  Suivant
AuteurMessage
Klax
Homme-canon
Homme-canon



MessageSujet: 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...
Revenir en haut Aller en bas
http://www.rpcirkus.org
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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 ...
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Gluonmou
Contorsionniste
Contorsionniste



MessageSujet: 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
Revenir en haut Aller en bas
Invité
Invité



MessageSujet: Re: MCNP : the Godfather   Sam 19 Oct 2013 - 19:55

Bonsoir
Vive les codes de calculs... study 
KLOUG
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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=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
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
SheepeR
Trapéziste
Trapéziste



MessageSujet: 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....
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: Re: MCNP : the Godfather   Lun 21 Oct 2013 - 11:54

Bonjour,
Malheureusement, il n'existe pas de version test de ce code,
GA
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Invité
Invité



MessageSujet: 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.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.

Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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 Laughing 
GA
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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 : http://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
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Fred
Homme-canon
Homme-canon



MessageSujet: Re: MCNP : the Godfather   Mer 23 Oct 2013 - 20:34

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

_________________
Fred :pig:
Revenir en haut Aller en bas
Gluonmou
Contorsionniste
Contorsionniste



MessageSujet: 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?
Revenir en haut Aller en bas
Invité
Invité



MessageSujet: 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
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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 :

http://www.forum-rpcirkus.com/t4128-mcnp-the-godfather
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Invité
Invité



MessageSujet: 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
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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 ...
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Invité
Invité



MessageSujet: 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
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Invité
Invité



MessageSujet: 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
Revenir en haut Aller en bas
Invité
Invité



MessageSujet: 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
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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.
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Invité
Invité



MessageSujet: Re: MCNP : the Godfather   Ven 20 Déc 2013 - 22:47

 Very Happy je vous remercie infiniment pour votre aide
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: 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
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
 
MCNP : the Godfather
Voir le sujet précédent Voir le sujet suivant Revenir en haut 
Page 1 sur 3Aller à la page : 1, 2, 3  Suivant
 Sujets similaires
-
» Présentation The Godfather
» Accès MCNP et MERCURAD

Permission de ce forum:Vous ne pouvez pas répondre aux sujets dans ce forum
Forum technique de RadioProtection Cirkus :: Un peu de théorie :: Code de calcul et transport de particules-
Sauter vers: