CRéation d'une source sous MCNP

Bonjour à tous,

Je débute avec MCNP et je cherche à créer une source avec plusieurs RN de ce type :

90% Cs-137 + 10 % Co-60, activité totale de 3e10 Bq (pour ce cas précis, je ne m'intéresse qu'aux photons)

j'utilise sdef de la manière suivante :

sdef x=d2 y=d3 z=d4 par=2 erg=D1 cel=1      $ j'utilise la méthode de rejection sur la cellule 1
si1 L E1....En                                               $ Energies des photons pour le Cs137 et le Co60
sp1 P1...Pn                                                  $ probabilités correspondant aux En
si2 0 414
sp2 0 1
si3 0 55
sp3 0 1
si4 0 244
sp4 0 1

J'envisage d'utiliser ensuite un tally F5 avec les facteurs de conversion issu de la cipr 74 pour obtenir des débits de dose (d'ailleurs, pourquoi la version japonaise est téléchargeable gratuitement sur le site de icrp alors que les autres versions sont payantes ? --> y a t'il une possibilité de l'obtenir en anglais ou français de manière gratuite ?)

Mes questions sont les suivantes :

1. sur la ligne sp1 dois-je pondérer les probabilités par les ratios correspondant aux spectres (c'est à dire 0.9*la proba pour le Cs137 et 0.1*les proba pour les 2 principales émissions liées au cobalt) ? (il me semble que oui).
   
2. Comment prendre en compte l'activité de la source dans le fichier input ?
   

Merci d'avance de vos réponses !

Coucou,

pour ta question 1 : oui !

Pour la question 2 : il faut rajouter la carte WGT
sdef x=d2 y=d3 z=d4 par=2 erg=D1 cel=1 WGT= (activité x somme des probabilités d'émission)
il me semble que c'est un truc comme ça.

pour la conversion fluence dose, j'utilise ça :
c Tally 
F5:p ......... (remplacer par les bonnes valeurs)
c c Coefficient CIPR 74 fluence <-> µSv/h
DE5   LIN 0.01       0.0125     0.015      0.0175     0.02       0.025      
          0.03       0.04       0.05       0.06       0.08       0.1        
          0.125      0.15       0.2        0.3        0.4        0.5        
          0.6        0.8        1          1.5        2          3          
          4          5          6          8          10         
c
DF5   LOG 2.4624E-04 1.6551E-03 3.0508E-03 3.7423E-03 3.8053E-03 3.5459E-03
          2.9584E-03 2.3494E-03 2.0853E-03 1.9889E-03 2.1100E-03 2.4253E-03
          2.9183E-03 3.4711E-03 4.5977E-03 6.8012E-03 8.8452E-03 1.0761E-02
          1.2535E-02 1.5808E-02 1.8779E-02 2.5094E-02 3.0497E-02 3.9770E-02
          4.8042E-02 5.5537E-02 6.3080E-02 7.7876E-02 9.2791E-02




bonne chance !

Merci beaucoup pour cette réponse rapide !

Bonjour,

Pour mettre un Facteur Multiplicatif, il faut utiliser la carte FM

Exemple :

DE15  0.010 0.015 0.020 0.030 0.040 0.050 0.060 0.080 0.1 0.15 0.2 &
     0.3 0.4 0.5 0.6 0.8 1 1.5 2 3 4 5 6 8 10
c facteur de conversion en Sv par heure H*(10) CIPR74
DF15  2.196E-10 2.988E-09 3.780E-09 2.916E-09 2.304E-09 1.98E-09 1.836E-09 &
      1.908E-09 2.196E-09 3.204E-09 4.320E-09 6.480E-09 8.568E-09 1.0548E-08 &
      1.2384E-08 1.5768E-08 1.872E-08 2.484E-08 3.096E-08 3.996E-08 4.824E-08 &
      5.580E-08 6.336E-08 7.776E-08 9.216E-08

FM15 2000

Donc là tes résultats du tally 15 sont multipliés par 2000

Merci Grumpy,

Du coup dois-je en déduire que le 2000 de la carte FM correspond à l'activité (c'est à dire 2000 Bq) puisque le résultat du tally est donné pour une particule entrante ?

Pour être sûr de bien comprendre, est-ce que la solution via la carte FM pour le tally est équivalente à la solution avec le WGT de la carte SDEF (cf message de Samael22) ?

Merci encore pour toute cette aide !