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 : Précédent  1, 2, 3  Suivant
AuteurMessage
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/
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: Re: MCNP : the Godfather   Sam 21 Déc 2013 - 14:35

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



MessageSujet: Re: MCNP : the Godfather   Dim 22 Déc 2013 - 0:13

Bsr, je viens de voir votre message, merci infiniment Very Happy
Revenir en haut Aller en bas
Invité
Invité



MessageSujet: Re: MCNP : the Godfather   Sam 28 Déc 2013 - 22:25

Bonsoir, existe t-il une formation de courte durée sur le code mcnp en France
Merci
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: Re: MCNP : the Godfather   Dim 29 Déc 2013 - 14:38

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



MessageSujet: Re: MCNP : the Godfather   Jeu 2 Jan 2014 - 11:30

Bonjour,
pour répondre à lail voir : http://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 ( http://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)
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Gamma.Lambda
Trapéziste
Trapéziste



MessageSujet: Re: MCNP : the Godfather   Sam 4 Jan 2014 - 15:26

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



MessageSujet: Re: MCNP : the Godfather   Sam 4 Jan 2014 - 15:32

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



MessageSujet: Re: MCNP : the Godfather   Sam 4 Jan 2014 - 15:55

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



MessageSujet: Re: MCNP : the Godfather   Sam 4 Jan 2014 - 16:00

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



MessageSujet: Re: MCNP : the Godfather   Sam 4 Jan 2014 - 17:06

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



MessageSujet: Re: MCNP : the Godfather   Sam 4 Jan 2014 - 17:41

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



MessageSujet: Re: MCNP : the Godfather   Sam 4 Jan 2014 - 20:55

Encore merci,
finalement les cylindres fonctionnent aussi sous cette version, mais pas les cubes.
Bonne soirée.
Revenir en haut Aller en bas
Invité
Invité



MessageSujet: Re: MCNP : the Godfather   Dim 5 Jan 2014 - 20:50

Bonsoir, vous m'aidez bcp, je ne sais pas comment vous remercier,
 j'ai deux petite question par rapport à:
cut:e j .01 3j
ctme 15
Revenir en haut Aller en bas
Invité
Invité



MessageSujet: Re: MCNP : the Godfather   Lun 6 Jan 2014 - 19:57

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



MessageSujet: Re: MCNP : the Godfather   Lun 6 Jan 2014 - 20:45

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



MessageSujet: transport de particule   Mar 7 Jan 2014 - 13:41

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];
Revenir en haut Aller en bas
Invité
Invité



MessageSujet: Re: MCNP : the Godfather   Mar 7 Jan 2014 - 13:46

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



MessageSujet: Re: MCNP : the Godfather   Mar 7 Jan 2014 - 14:02

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



MessageSujet: Re: MCNP : the Godfather   Mar 7 Jan 2014 - 21:00

Bonsoir,
Je vais essayer de répondre à lail  sur le post http://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
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: Re: MCNP : the Godfather   Mar 7 Jan 2014 - 21:02

Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Invité
Invité



MessageSujet: Re: MCNP : the Godfather   Jeu 9 Jan 2014 - 13:25

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



MessageSujet: Re: MCNP : the Godfather   Jeu 9 Jan 2014 - 13:26

FICHIER DE SORTIE
1-       c fichier cirkus 2                                                             
    2-       1 1 -0.00129 -1 2                                                              
    3-       2 1 -0.00129 -2                                                                
    4-       999 0 1                                                                         
    5-                                                                                      
    6-       1 so 150                                                                       
    7-       2 sx 100 .5                                                                     
    8-                                                                                      
    9-       M1 7014 .781 8016 .212 7015 .0028 $ air                                        
   10-       cut:e j .01 3j                                                                 
   11-       mode p e                                                                       
   12-       imp:p 1  1 0                                                                    
   13-       imp:e 0 1 0                                                                    
   14-       sdef pos 0 0 0 par=2 erg=d1                                                    
   15-       si1 L 1.17 1.33                                                                 
   16-       Sp1 1 1                                                                        
   17-       *f6:p,e 2                                                                      
 fatal error.  tally 6 is not available for electrons.
   18-       ctme 15                                                                        
 
 warning. zero electron importance within non-void cell    1.
1source                                                                                                 print table 10
 
 values of defaulted or explicitly defined source variables
 
    cel        0.0000E+00
    sur        0.0000E+00
    tme        0.0000E+00
    dir        isotropic
    pos        0.0000E+00   0.0000E+00   0.0000E+00
    x          0.0000E+00
    y          0.0000E+00
    z          0.0000E+00
    rad        0.0000E+00
    ext        0.0000E+00
    axs        0.0000E+00   0.0000E+00   0.0000E+00
    vec        0.0000E+00   0.0000E+00   0.0000E+00
    ccc        0.0000E+00
    nrm        1.0000E+00
    ara        0.0000E+00
    wgt        1.0000E+00
    eff        1.0000E-02
    par        2.0000E+00
    tr         0.0000E+00
 
 
 
 probability distribution   1 for source variable erg
 unbiased discrete  distribution
 
 source        source           cumulative          probability
  entry         value           probability           of bin
 
    1       1.17000E+00        5.000000E-01        5.000000E-01
    2       1.33000E+00        1.000000E+00        5.000000E-01
 
     the mean of source distribution   1 is  1.2500E+00
 
 
 order of sampling source variables.
 pos erg tme
 
1tally   6                                                                                              print table 30
           tally type 6*   track length estimate of heating.                             
           tally for  photons   electrons
1material composition                                                                                   print table 40
 
 the sum of the fractions of material    1 was 9.958000E-01
 
 material
  number     component nuclide, atom fraction
 
     1            7014, 7.84294E-01       8016, 2.12894E-01       7015, 2.81181E-03
 
 
 
 
 
 material
  number     component nuclide, mass fraction
 
     1            7014, 7.61094E-01       8016, 2.35983E-01       7015, 2.92291E-03
 
 warning.   1 materials had unnormalized fractions. print table 40.
1cell volumes and masses                                                                                print table 50
 
        cell     atom          gram         input       calculated                            reason volume
                density       density       volume        volume         mass       pieces    not calculated
 
     1     1  5.38356E-05   1.29000E-03   0.00000E+00   1.41372E+07   1.82369E+04      1                     
     2     2  5.38356E-05   1.29000E-03   0.00000E+00   5.23599E-01   6.75442E-04      1                      
     3   999  0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00      0      infinite       
1surface areas                                                                                          print table 50
 
         surface    input      calculated    reason area
                    area          area       not calculated
 
     1        1  0.00000E+00   2.82743E+05                  
     2        2  0.00000E+00   3.14159E+00                  
 
 warning. p. and e. importances are inconsistent in    1 cells.
1cells                                                                                                  print table 60
 
                      atom        gram                                     photon     electron                        
        cell  mat    density     density     volume       mass     pieces importance importance                       
 
     1     1    1  5.38356E-05 1.29000E-03 1.41372E+07 1.82369E+04     1  1.0000E+00 0.0000E+00                        
     2     2    1  5.38356E-05 1.29000E-03 5.23599E-01 6.75442E-04     1  1.0000E+00 1.0000E+00                       
     3   999    0  0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00     0  0.0000E+00 0.0000E+00                        
 
 total                                     1.41372E+07 1.82369E+04
1surfaces                                                                                               print table 70
 
         surface     trans  type   surface coefficients
 
     1      1                so    1.5000000E+02
     2      2                sx    1.0000000E+02   5.0000000E-01
 
    minimum source weight = 1.0000E+00    maximum source weight = 1.0000E+00
 
 ***************************************************
 * Random Number Generator  =                    1 *
 * Random Number Seed       =       19073486328125 *
 * Random Number Multiplier =       19073486328125 *
 * Random Number Adder      =                    0 *
 * Random Number Bits Used  =                   48 *
 * Random Number Stride     =               152917 *
 ***************************************************
 
 
         3 warning messages so far.
 
          1 fatal error so far.
1physical constants                                                                                     print table 98
 
 
   name                    value     description
 
   huge      1.0000000000000E+36     infinity
    pie      3.1415926535898E+00     pi
  euler      5.7721566490153E-01     euler constant
 avogad      6.0220434469282E+23     avogadro number (molecules/mole)
  aneut      1.0086649670000E+00     neutron mass (amu)
  avgdn      5.9703109000000E-01     avogadro number/neutron mass (1.e-24*molecules/mole/amu)
  slite      2.9979250000000E+02     speed of light (cm/shake)
 planck      4.1357320000000E-13     planck constant (mev shake)
  fscon      1.3703930000000E+02     inverse fine structure constant h*c/(2*pi*e**2)
 gpt(1)      9.3958000000000E+02     neutron mass (mev)
 gpt(3)      5.1100800000000E-01     electron mass (mev)
 
 
 fission q-values:     nuclide    q(mev)   nuclide    q(mev)
                         90232    171.91     91233    175.57
                         92233    180.84     92234    179.45
                         92235    180.88     92236    179.50
                         92237    180.40     92238    181.31
                         92239    180.40     92240    180.40
                         93237    183.67     94238    186.65
                         94239    189.44     94240    186.36
                         94241    188.99     94242    185.98
                         94243    187.48     95241    190.83
                         95242    190.54     95243    190.25
                         96242    190.49     96244    190.49
                         other    180.00
 
 
 the following compilation options were used:
 
     cheap
     dec
     plot
     mcplot
     xlib
     default datapath: C:\Program Files\LANL\MCNPDATA                                                  
                       C:\Program Files\LANL\MCNP5DATA                                                
 
 bad trouble in subroutine unique of xact                            
 
 cannot create teruntp                      
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: Re: MCNP : the Godfather   Jeu 9 Jan 2014 - 20:40

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



MessageSujet: Re: MCNP : the Godfather   Ven 10 Jan 2014 - 14:58

et voici les fichiers...
klax
Fichiers joints
Desktop.zip
Vous n'avez pas la permission de télécharger les fichiers joints.
(6 Ko) Téléchargé 12 fois

_________________
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 11 Jan 2014 - 8:42

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
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 2 sur 3Aller à la page : Précédent  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: