Curso MCNP Aula2011a

Fazer download em ppt, pdf ou txt
Fazer download em ppt, pdf ou txt
Você está na página 1de 30

Introduo ao Cdigo de Transporte de Radiao MCNP

[email protected]

O Mtodo de Monte Carlo (MMC)


MMC um mtodo matemtico estocstico, isto , baseado em funes de probabilidade; MMC no sinnimo de MCNP MCNP um cdigo computacional que utiliza MMC.

[email protected]

O Mtodo de Monte Carlo (MMC)


Provavelmente a primeira utilizao do MMC que se tem notcia foi realizada em 1772 por Buffon; Laplace sugeriu, em 1786, que poderia ser calculado por MMC Supe-se que Kelvin tenha usado MMC para calcular as integrais temporais da energia cintica que aparecem na teoria cintica dos gases.
[email protected]

O Mtodo de Monte Carlo (MMC)


A primeira aplicao do MMC no transporte de radiao foi realizada na dcada de 1940, em Los Alamos-EUA, para traar movimentos de nutrons em materiais fissionveis (Desenvolvimento da Bomba Atmica Projeto Manhattan). O MMC pode ser aplicado a qualquer tipo de problema, desde que seja estabelecida uma funo de probabilidade coerente.

[email protected]

O Mtodo de Monte Carlo (MMC)


Exemplo: Clculo de rea

S=?

SQ=LxL

S/SQ =P =NS/NT

P Probabilidade de uma pedra cair na rea S, NS numero de pedras que caem na rea S NT o nmero total de pedras lanadas
[email protected]

O cdigo Computacional MCNP


MCNP: Monte Carlo N-Particle; Cdigo de propsito geral (Aplicao em metrologia, medicina, militar, etc); MCNP no faz milagre: para extrair uma resposta coerente com a realidade deve-se primeiro conhecer profundamente o problema fsico antes de se realizar uma simulao; importante, quando possvel, validar os resultados da simulao com medidas experimentais;

[email protected]

Arquivo de entrada do MCNP (Estrutura)


O arquivo de entrada contm a descrio da geometria, da fonte, dos materiais e das respostas desejadas para o problema simulado. Ele deve ser organizado conforme a estrutura abaixo:
Ttulo do problema, somente uma linha (Opcional, mas fortemente recomendvel) Descrio das clulas . . Linha em branco Especificao das superfcies descritas no carto de clulas . . Linha em branco Especificao da fonte, dos materiais utilizados no problema, das respostas desejadas (tally) e das instrues de execuo.

Bloco de Clulas* Bloco de Superfcie Bloco de Dados

*Clula a regio que representa o volume formado por uma ou varias superfcies.
[email protected]

Arquivo de entrada do MCNP (Caractersticas)


Cada linha de comando deve possuir at 80 colunas, caso ultrapasse, utilize & e continue na linha seguinte ou escreva diretamente a partir da quinta coluna; Uma linha de comentrio iniciada com a letra C e no est limitada 80 colunas; O Comentrio de fim de linha feito com o caractere $; No h limite de nmero linhas O arquivo pode ser feito no wordpad ou bloco de notas, entretanto, deve-se retirar a extenso .txt; O Texto no deve conter acentos ou ; O nome do arquivo de entrada no pode exceder a 8 caracteres (letras {sem acentos} e nmeros)

[email protected]

Exemplo de arquivo de entrada


Exemplo introducao ao MCNP C bloco de celulas (associacao de superficies) 1 1 -1 (-100) imp:N=1 $ esfera de agua com raio 5cm 2 0 (100 -200) imp:N=1 $ universo externo vacuo 3 0 (200) imp:N=0 $ fim do transporte C bloco de superficie 100 SO 5 $ esfera centrada na origem 200 SO 10 $ esfera concentrica a 100 e de raio 10cm C bloco de dados SDEF POS=0 0 0 ERG=1 $ Fonte pontual monoenergetica (1 MeV) C Resposta F1:N 200 $ Neutrons que atravessam a superficie 200 M1 1001 2 8016 1 $ Descrio do material (H-2 O-1) NPS 10000 $ numero de eventos limitador de execuo

[email protected]

Representao geomtrica de superfcies no MCNP


Superfcies bsicas so designadas atravs de cdigos (mnemnicos); Planos, cilindros e cones possuem extenso infinita; Dimenses das superfcies: centmetros; As Superfcies so delimitadoras de regio/volume/clula;

[email protected]

Tipo

Representao de algumas superfcies no MCNP


representao no MCNP PX (ou PY,PZ) SO SX (ou SY,SZ) CX (ou CY, CZ) C/X (ou C/Y,C/Z) Equao X-D=0 X2+Y2+Z2=R2 (X-X)2+Y2+Z2=R2 Y2+Z2=R2 (Y-Y)2+(Z-Z)2=R2

Parmetros de entrada

Plano Normal ao eixo X Esfera Centrada na origem Esfera Centrada sobre o eixo X Cilindro Sobre o eixo X Cilindro Paralelo ao eixo X

D (Posio do plano no eixo) R (Raio da esfera) X(pos. cent. esf. em X R (Raio da esfera) R (Raio do cilindro) Y(pos. eix. Cil. Em Y) Z e R

[email protected]

Planos sobre os eixos X,Y,Z (PX,PY,PZ)


PX 3 PX plano que intercepta o eixo X 3 Posio (cm) em X PY 5 PY plano que intercepta o eixo Y 3 Posio (cm) em Y PZ 6 PZ plano que intercepta o eixo Z 6 Posio (cm) em Z

[email protected]

Esferas na origem (SO)

Esfera de 1 cm de raio, centrada na origem: SO 1

[email protected]

Esferas sobre os eixos X,Y,Z (SX,SY,SZ)


Esfera centrada sobre o eixo X: SX 5 1 Esfera centrada sobre o eixo Y: SY 3 1 Esfera centrada sobre o eixo Z: SZ 4 1

[email protected]

Definio geral de esfera (S)

S XYZ R (Centro) S 3 4 1 1

[email protected]

Cilindros sobre os eixos X,Y,Z (CX,CY,CZ)

Cilindro centrado em X: CX R Cilindro centrado em Y: CY R Cilindro centrado em Z: CZ R

Ateno: cilindros tm extenso infinita


[email protected]

Cilindros paralelos aos eixos X,Y,Z (C/X,C/Y,C/Z)

Cilindro paralelo ao eixo X: C/X Y Z R Cilindro paralelo ao eixo Y: C/Y X Z R Cilindro paralelo ao eixo Z: C/Z X Y R

[email protected]

Associao de superfcies construo de clulas


Clulas so construdas com uma ou mais superfcies. Somente aps definir as superfcies possvel construir as clulas.

[email protected]

Regio definida por uma superfcie


- Uma regio estabelecida atribuindo-se sinais as superfcies (no bloco de clulas)

Plano
Sentido do vetor normal igual do sentido do eixo (sinal + da super. no carto de clulas) Sentido do vetor normal diferente do sentido do eixo (sinal - da super. no carto de clulas)

Cilindro ou Esfera (super. Fechada)


A regio de interesse a parte interna -> sinal - da super. no carto de clulas; A regio de interesse a parte externa -> sinal + da super. no carto de clulas

Planos
n x n x

Demais superfcies

Sentidos iguais: +

Sentidos opostos:

Regio interna: Regio externa: +


[email protected]

Clulas no MCNP
No MCNP, uma clula a unidade bsica de volume; Todo o espao tem que ser definido, at o universo externo ao transporte de radiao; As clulas so definidas considerando as superfcies definidas no bloco de superfcie; As superfcies que formam uma clula so associadas atravs de operadores: unio (:), interseo (espao), complemento (#);

[email protected]

Operadores de associao de superfcies


Clulas so construdas com uma ou vrias superfcies associadas

Interseo Unio Complemento

Espao entre os nmeros das Superfcies Dois pontos (:) entre os nmeros das Superfcies

Caractere # seguido do n da clula complementada

[email protected]

[email protected]

Sintaxe do carto de Clulas


N da clula N do material densidade do material descrio da geometria

0= Vcuo

(-) g/cm (+) densidade atmica ( at/cmx10-24)

80 colunas

[email protected]

Ordem das operaes no MCNP


Intersees so realizadas primeiro que unies -1:-3 -2 equivalente a -1:(-3 -2)

1 2 3 3 2 3

1 2 3

-1:-3 -2

(-1:-3) -2

[email protected]

Exemplo 1
1 numere todas as superfcies; 2 numere todas as clulas; 3 defina todas as superfcies; 4 defina todas as clulas incluindo a regio externa ao transporte de radiao;

[email protected]

Exemplo 2
Cilindro ao longo do eixo z

[email protected]

Exemplo 3
Cubo simples

1 4 6 3 2

[email protected]

exerccios
Modele duas esferas concntricas em (0,0,0) e de raios 0,1m e 0,5m. Modele uma geometria cilndrica com dimenses internas de R= 2cm e L=4cm, e dimenses externas de R=2,5cm e L=5cm. Modele uma esfera de 0,3m de dimetro inscrita em um cilindro de 20cm de raio e 40 cm de altura.
[email protected]

Exerccio 1: Tanque cilndrico 3m

Espessura parede: 5cm

1 3 5 Z Y 1 2

2 3 5 4 3

1m

[email protected]

Exerccio 1: Tanque cilndrico 3m

Espessura parede: 5cm

1 3 5 Z Y 1 2

2 3 5 4 3 Titulo: execicio 1 1 0 (-1 2 -3) 2 0 (3 -5 -1 4):(-2 4 -5) 3 0 -4:5:1 1 2 3 4 5 PZ PZ CZ PZ CZ 300 5 45 0 50

1m
Titulo: execicio 1 1 0 (-1 2 -3) 2 0 ((3:-2) -1)(-1 -5 4) 3 0 -4:5:1 1 2 3 4 5 PZ PZ CZ PZ CZ 300 5 45 0 50
[email protected]

Você também pode gostar