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, 4
AuteurMessage
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/
lail
Saltimbanque
Saltimbanque



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
lail
Saltimbanque
Saltimbanque



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/
serge
Curieux
Curieux



MessageSujet: Simulation par mcnp   Mer 5 Fév 2014 - 2:03

Bonjour, je débute avec le code mcnp, je suis arrivé à faire des trucs avec, et là je dois mesurer l'efficacité d'un détecteur NaI, je bloque sur la partie tally car je dois en même temps tiré le spectre tracé, 

Merci d'avance pour vos réponses


Sujet déplacé : ce topic rassemble toutes les questions sur MCNP.
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: Re: MCNP : the Godfather   Mer 5 Fév 2014 - 20:27

Bonsoir,
Il suffit de pondérer ton tally par E.
Par exemple E8 0 1024i 3.072 pour un 1024 canaux.
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
xln10
Spectateur
Spectateur



MessageSujet: Re: MCNP : the Godfather   Mer 2 Juil 2014 - 10:52

Bonjour a tous,

Enfin un Forum francais sur MCNP. Je vais tenter ma Chance!

Je simule actuellement un detecteur a scintillation NaI(Tl). Je souhaite obtenir le spectre du Cs137 pour commencer.
Ma Simulation est plutot complete: géometrie interne du detecteur, environnement...
Mes résultats sont jolis par leur forme. Mais le compton continuum est trop Faible par Rapport a mes essais tandis que le photopic est trop élevé. En jouant un peu Avec la densité de l'air ou du cristal, j'obtiens des résultats "satisfaisants". Mais voila, l'air a vu sa densité multipliée par 20.

Je suis tout ouie pour des pistes de reflections.

PS: je ne connais pas le vocabulaire technique... en francais!
Revenir en haut Aller en bas
GA_mcnp (x)
Ventriloque
Ventriloque



MessageSujet: Re: MCNP : the Godfather   Mer 2 Juil 2014 - 18:58

Bonjour,
Le mieux est de mettre en pj votre fichier d'entrée afin de regarder si il n'a pas des évidences à corriger.
Revenir en haut Aller en bas
http://mcnp.lanl.gov/
Gluonmou
Funambule
Funambule



MessageSujet: Re: MCNP : the Godfather   Ven 4 Juil 2014 - 16:12

Bonjour,

Etes vous sur des dimensions de votre cristal?
Un plateau Compton plus faible et un pic plus grand signe un détecteur plus efficace que prévu en absorption totale, c'est à dire globalement plus volumineux ou avec un rapport Volume/surface plus grand à volume identique

Cordialement

Gluonmou
Revenir en haut Aller en bas
xln10
Spectateur
Spectateur



MessageSujet: Re: MCNP : the Godfather   Mar 8 Juil 2014 - 10:40

Après verification, l'activité de la soruce a une erreur maximale de 3%. Apres Simulation, seul le photopic est trop elevé, le continuum est tres satisfaisant.

Mais treve de blabla, voici la bete (les annotations sont en allemand)

c BLOCK 1: CELLS ***********************************************
   1 0          101          imp:p=0 $ Au�er des Problems
   2 15 -.0012 -101 119 121 111 129 127 130 131 132 133 134 135 136 137 &
     138 139 140 141 147 150 151 imp:p=1 $ Au�er des Detektors
   3 15 -.0012 (104 109 110 -122 123)Sad-120 118 116 114:-115)Sad-113 110): &
     (-112 110)Sad-124 104)Sad-125 104) imp:p=1 $ Au�er des Housing    
  11 1 -3.667  -102          imp:p=1 $ Probe aus NaJ (Dichte = 3,667 g/cm�)
  12 2 -2.7    -103 144      imp:p=1 $ Alu extrabeschichtung
  13 2 -2.7    -104 103      imp:p=1 $ Kristall Housing (AlMg3)
  14 3 -2.23   -105 106      imp:p=1 $ Borosilicate glass
  15 4 -8.7    -107 105      imp:p=1 $ Mu Metall
  16 5 -.93    -108 107      imp:p=1 $ PE Housing
  17 2 -2.7    -109 108      imp:p=1 $ Housing (AlMg3)
  18 15 -.0012 -106 126 143 142 imp:p=1 $ Luft in Szintillator
  19 5 -0.93  (-110 146)Sad-110 -145)   imp:p=1 $ Plastik Sockel (PE)
  20 6 -8.47   -111 112 113  imp:p=1 $ Sockel (CuZn39Pl3)
  21 2 -2.7    -114 115      imp:p=1 $ Scheibe (Alu)
  22 7 -2.7    -116 117      imp:p=1 $ Tr�ger (Alu)
  23 2 -2.7    -118          imp:p=1 $ Scheibe (Alu)
  24 2 -2.7    -119 120      imp:p=1 $ Hintere Haube
  25 2 -2.7    -121 122      imp:p=1 $ Vornere Haube
  26 8 -1.38   -123 124 125  imp:p=1 $ Kappe (PET)
  27 9 -3.5    -126          imp:p=1 $ Dynode aus Stahl
  28 10 -1     -127          imp:p=1 $ Tisch aus Holz $ dichte 1 statt .75
  29 11 -2.25  -129:-151     imp:p=1 $ Wand und Boden aus Concrete MCNP
  30 2 -2.7 (-130:-131:-132) 140.1 141.1 imp:p=1 $ Unter Tr�ger
  31 2 -2.7   -133 140.1 141.1 121.1 imp:p=1 $ Tr�ger
  32 10 -1 (-134:-135) 136.1 imp:p=1 $ Drehtisch $ dichte 1 statt .75
  33 13 -8.07  -136          imp:p=1 $ Drehtisch Achse
  34 9 -5.5   -137:-138:-139 imp:p=1 $ Tisch Struktur
  35 9 -7.87   -140:-141     imp:p=1 $ Schrauben
  36 9 -7.87   -142          imp:p=1 $ Kollimator
  37 14 -4.5   -143          imp:p=1 $ 1. Dynode
  38 12 -3.58  -144 102      imp:p=1 $ Reflektor (MgO)
  39 9 -7      -146 145 -110 imp:p=1 $ Ring in Sockel
  40 9 -7      -147          imp:p=1 $ Buchse
  41 16 -1.8   -117          imp:p=1 $ Leiterblatt
  42 17 -1.873 -148          imp:p=1 $ Quelle
  43 10 -1     -149 148      imp:p=1 $ Papier um die Quelle
  44 5 -.93    -150 149      imp:p=1 $ Plastik um die Quelle
 
c BLOCK 2: SURFACE CARDS ***************************************
 101 RPP -100 100 -100 200 -100 100    $ Problem Zylinder
 102 RCC -2.54 0 0 5.08 0 0 2.54       $ Kristall Zylinder
 103 RCC -2.7737 0 0 5.3137 0 0 2.7737 $ Extra beschichtung (0,2337cm dick)
 104 RCC -2.8245 0 0 5.3645 0 0 2.8245 $ Housing Zylinder (0,0508cm dick)
 105 RCC 2.54 0 0 13.0251 0 0 2.54     $ Au�er Glass
 106 RCC 2.64 0 0 12.8251 0 0 2.44     $ Inner Glass (1mm dick)
 107 RCC 2.54 0 0 13.0251 0 0 2.5908   $ Mu Metall (0,0508cm)
 108 RCC 2.54 0 0 13.0251 0 0 2.8702   $ PE Housing
 109 RCC 2.54 0 0 13.0251 0 0 2.921    $ Alu Housing (0,0508cm)
 110 RCC 15.5651 0 0 4 0 0 2.8         $ Plastik Sockel
 111 RCC 18.5651 0 0 1.9 0 0 3.25      $ CuZn39Pb3 Sockel
 112 RCC 18.5651 0 0 1.9 0 0 1.7       $ Inner CuZn39Pb3 Sockel
 113 RCC 18.5651 0 0 1.5 0 0 2.85      $ Inner CuZn39Pb3 Sockel
 114 RCC 20.4651 0 0 0.4 0 0 2.95      $ Alu Scheibe
 115 RCC 20.4651 0 0 0.4 0 0 1.5       $ Alu Scheibe
 116 RPP 20.8651 29.0951 -2.8 2.8 -0.75 0.75 $ Alu Tr�ger
 117 RPP 20.8651 29.0951 -2.3 2.3 -0.75 0.75 $ Inner Alu Tr�ger
 118 RCC 29.0951 -0 0 0.4 0 0 2.95     $ Alu Scheibe
 119 RCC 20.4651 0 0 10.05 0 0 3.25    $ Hintere Haube
 120 RCC 20.4651 0 0 9.75 0 0 3        $ Inner hintere Haube
 121 RCC -3.2745 0 0 21.8396 0 0 3.25  $ Vordere Haube             
 122 RCC -3.1745 0 0 21.7396 0 0 3     $ Inner vordere Haube
 123 RCC -3.1745 0 0 1 0 0 3           $ PE Kappe  3,5mm
 124 RCC -3.1745 0 0 1 0 0 2.55        $ Inner PE Kappe
 125 RCC -2.8245 0 0 0.65 0 0 2.8245   $ Inner PE Kappe
 126 RPP 6.2 10.7 1 2 -1 1             $ Dynode
 127 RPP -100 100 -25.5 -22.5 -50 50   $ Tisch
c 128 RPP -140 -130 -100 205 -50 50     $ W�nde
 129 RPP -100 100 -100 -90 -100 100    $ Boden
 130 RPP 20 24 -6.5 -5 -13.25 13.25    $ Tr�ger Lineal
 131 RCC 22 -16.5 -11.25 0 10 0 2      $ Tr�ger Bein Links
 132 RCC 22 -16.5 11.25 0 10 0 2       $ Tr�ger Bein Recht
 133 RPP 20.5 23.5 -5 5 -7.25 7.25     $ Tr�ger
 134 RCC 0 -16.5 0 0 -3 0 40           $ �bere Drehtisch
 135 RCC 0 -19.5 0 0 -3 0 45           $ Untere Drehtisch
 136 RCC 0 -16.5 0 0 -6 0 3            $ Drehtisch Achse
 137 RPP -100 -98 -30.5 -25.5 -50 50   $ Tisch Struktur 1
 138 RPP 98 100 -30.5 -25.5 -50 50     $ Tisch Struktur 2
 139 RPP -98 98 -28.5 -25.5 -1 1       $ Tisch Struktur 3
 140 RCC 22 -8.5 -5 0 16 0 .75         $ Schrauben 1
 141 RCC 22 -8.5 5 0 16 0 .75          $ Schrauben 2
 142 RCC 6.04 0 0 .06 0 0 2            $ Kollimator
 143 RPP 6.2 7.2 -.75 .75 -1 1         $ 1. Dynode
 144 RCC -2.5585 0 0 5.0985 0 0 2.5585 $ Luft zwischen kristall und reflektor
 145 CX 2                              $ Ring in Sockel
 146 CX 2.5                            $ Ring in Sockel au�er
 147 RCC 30.5151 0 0 .86 0 0 .95       $ Buchse
 148 RCC -50.01 0 0 .02 0 0 3.51       $ Quelle
 149 RCC -50.025 0 0 .05 0 0 3.51      $ Papier um die Quelle
 150 RCC -50.0365 0 0 .075 0 0 4       $ Plastik um die Quelle
 151 RPP -100 100 195 200 -100 100     $ Dach
c BLOCK 3:DATA CARDS ******************************************
SDEF erg=d1 par=2 pos=-50 0 0 axs=1 0 0 ext=d3 rad=d2 WGT=4.31e8
SI1 L .66166
SP1 D 1.2e5
si2 0 3.51
sp2 -21 1
si3 0 0.05
sp3 -21 0
mode p
Phys:p 1 1 0 0 1 0
CUT:p j .01
nps  1e8
c
f18:p 11
e18 .06 861i .8
ft18 geb -.01 .07 -.08
c
m1   11000 0.45 $ Na
     53000 0.45 $ I
     81000 0.1  $ Tl (0,1% mole)
m2   13000 -0.97     $ Al (AlMg3)
     12000 -0.03     $ Mg
m3   5000 -0.040064  $ Glass Borosilicate
     8000 -0.539562
     11000 -0.02819
     13000 -0.011644
     14000 -0.377220
     19000 -0.00332
m4   28000 -0.77     $ Ni (Mu Metall - NiFe15Cu5Mo3)#(-111 112 113)#(-114 115)#(-116 117)
     26000 -0.15     $ Fe
     29000 -0.05     $ Cu
     42000 -0.03     $ Mo
m5   1000 -0.143716  $ PE ohne Bore
     6000 -0.856284
m6   29000 -0.58 $ CuZn39Pb3
     30000 -0.39
     82000 -0.03
m7   13000 -0.92 $ AlMgSi
     14000 -0.04
     12000 -0.04
m8   1000 -0.041960 $ PET (1,38)
     6000 -0.625016
     8000 -0.333024
m9   6000 -.005 $ Stahl
     26000 -.995
m10  1000 -.059642 $ Holz 0,75
     6000 -.497018
     7000 -.00497
     8000 -.427435
     12000 -.001988
     16000 -.004970
     19000 -.001988
     20000 -.001988
m11  1000 -.00453 $ Concrete MCNP
     8000 -.5126
     11000 -.01527
     13000 -.03555
     14000 -.36036
     18000 -.125
     20000 -.05791
     26000 -.01378
m12  8000  -.39696 $ Reflektor aus MgO
     12000 -.60304
m13  26000 -.000868 $ Messing
     29000 -.665381
     30000 -.325697
     50000 -.002672
     82000 -.005377
m14  29000 -1      $ Kupfer
m15  6000 -.000124 $ Luft
     7000 -.755268
     8000 -.231781
     18000 -.012827
m16  1000 -.25 $ Epoxy
     6000 -.625
     8000 -.125
m17  1000 -.037764
     6000 -.5625144
     8000 -.2997216
     55000 -.1 $ Cesium
Revenir en haut Aller en bas
xln10
Spectateur
Spectateur



MessageSujet: Re: MCNP : the Godfather   Mar 8 Juil 2014 - 15:27

J'ai trouvé!!!! Trop de Tallium dans mon cristal!
Merci tout de meme pour les pistes de reflections
Revenir en haut Aller en bas
 

MCNP : the Godfather

Voir le sujet précédent Voir le sujet suivant Revenir en haut 
Page 4 sur 4Aller à la page : Précédent  1, 2, 3, 4

 Sujets similaires

-
» MCNP : the Godfather
» Accès MCNP et MERCURAD
» Les protections neutroniques
» MCNPX, comment mieux le connaître ?
» Calcul de criticité 2

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-