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] Calcul d'isodose

Voir le sujet précédent Voir le sujet suivant Aller en bas 
AuteurMessage
Samael22
Trapéziste
Trapéziste



MessageSujet: [MCNP] Calcul d'isodose   Lun 28 Sep 2015 - 9:21

Me voici devant un nouveau problème. 

Énoncé: 
Une traversée (pour passer des câbles) va être réalisé dans une cellule. Un calcul MCNP doit prouver que la traversée n'engendre pas un débit de dose >25µSv/h.

J'ai donc placé des détecteurs tally F5 un peu de partout autour de la traversée. leur résultat est < 25 µsV/h, ce qui est déjà un bon point mais mon calcul ne satisfait pas le client.
Celui-ci voudrait un calcul représentant les isodoses. 

Je n'ai aucune idée de comment faire ça. Avez-vous des idées ?
Peut-être avec le logiciel Visu ?

Merci beaucoup !
Revenir en haut Aller en bas
xln10
Chauffeur de salle
Chauffeur de salle



MessageSujet: Re: [MCNP] Calcul d'isodose   Lun 28 Sep 2015 - 10:22

Bonjour,

tu as deja envisagé une représentation graphique sur MCNP grace au Meshtally?
-> représentation en Couleur (pour des présentations PPT clown )
-> possiblité de "cibler" un Point de ta Simulation et connaitre ses caracteristiques

Par contre c'est MCNP: tout en 2D (Tranche)
Revenir en haut Aller en bas
Samael22
Trapéziste
Trapéziste



MessageSujet: Re: [MCNP] Calcul d'isodose   Lun 28 Sep 2015 - 12:39

Je n'ai qu'une connaissance limité de MCNP. Je m'en servais jusqu'à lors pour calculer des dépôts d'énergie.
Je n'ai aucune idée du fonctionnement du Meshtally.
Pourrais-tu m'en dire un peu plus ?
(de mon coté je vais lire le "manuel")
Revenir en haut Aller en bas
xln10
Chauffeur de salle
Chauffeur de salle



MessageSujet: Re: [MCNP] Calcul d'isodose   Mar 29 Sep 2015 - 9:51

Il faut que tu places cette commande dans le paragraphe Tally.

La structure est assez simple, il s'agit de décomposer l'espace à étudier en Petits cubes (carrés ou rectangles). Le mieux est de le décomposer à l'aide de la fonction i, pour avoir des espaces de meme largeur... On peut aussi utiliser des coordonnées cylindriques ou spheriques (RMESH, CMESH, SMESH)
Pour un calcul plus rapide, focalise toi sur la cible, choisi une Resolution basse (Taille des carrés plus grande) en premier lieu, tu pourras afiner si tes resultats te satisfont.

Il y a 4 types de tally (tout bien expliqué au chapitre 5.6.23.- dans mon manuel)
Je ne connais que tres peu les tally autre que TOTAL (energy deposited)

Petit exemple:
TMESH $ Balise indiquant qu'on utilise le Meshtally
RMESH3 TOTAL $ TAlly type 3, Total d'energie déposée 
CORA3 -200 598i 100 $ coordonnées des intervales étudiés selon x [-200;100], 600 intervales
CORB3 -150 598i 150 $ coordonnées des intervales étudiés selon y [-150;150], 600 intervales
CORC3 -100 100 $ coordonnées des intervales étudiés selon z [-100;100], 1 intervale
ENDMD $ balise de fin du Meshtally

La ou ca se complique c'est pour "ploter" le Resultat:
D'abord lancer le calcul (mcnp irx Name=blabla.i)
une fois fini lancer le Dessin (mcnp z r=blabla.ir)
dans la commande, entrer "plot"
la plot window apparait.
cliquer sur xy donne le plan de coupe x=y=0
cliquer sur "tal" a droite pour activer le tally
puis cliquer sur "cel" a coté de Color jusqu'a voir apparaitre le nom de ta tally.
cliquer sur redraw pour afficher limage
(desactiver les Labels si besoin en cliquant sur "cel" en face de Label)

Petit tuto pas tres propre, mais jespere que ca te suffira, bon courage
Revenir en haut Aller en bas
Samael22
Trapéziste
Trapéziste



MessageSujet: Re: [MCNP] Calcul d'isodose   Mar 29 Sep 2015 - 10:39

En tout cas cheers  merci beaucoup. Je vais essayer ça demain ou après demain. Je te tiens au courant !
Revenir en haut Aller en bas
Samael22
Trapéziste
Trapéziste



MessageSujet: Re: [MCNP] Calcul d'isodose   Ven 16 Oct 2015 - 14:24

Me revoilà, finalement je n'ai pu essayer qu'aujourd'hui. 
Alors le  calcul a bugé avant la fin.

Le message d'erreur : 
bad trouble in mcrun ....
source particule n° ....
starting ....
Detected zero or negative mesh volume


après quoi j'ai quand voulu essayer de ploter.
ce message d'erreur apparaît : xlib: unable to connect to server display
puis ca revient à 
plot>


En tout cas merci beaucoup pour votre réponse qui m'a permis de lancer un début de quelque chose.
Revenir en haut Aller en bas
 
[MCNP] Calcul d'isodose
Voir le sujet précédent Voir le sujet suivant Revenir en haut 
Page 1 sur 1
 Sujets similaires
-
» Calcul proteines ?
» pour le calcul de la rétro , le .38 et .48 cumulatif?
» calcul des apss 2012-2013....c koi là...????
» aide au calcul de salaire
» CFE : chiffre d'affaires / calcul ajustement à l'année / congé de maternité

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: