Blindajes Neutrones

Universidad Pedagógica y Tecnológica de Colombia
Escuela de Fı́sica
CALCULO MONTE CARLO DE BLINDAJES ESRUCTURALES
241
PARA EL TRANSPORTE DE FUENTES DE Am −9 Be
por: Alejandra Marı́a Garcia Rodriguez

para optar al grado de Fı́sico
Director: PhD Segundo Agustı́n Martı́nez Ovalle
2017
Ciencias Básicas
Índice general
1. Fı́sica de Neutrones 8
1.1. Generación de neutrones . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.1.1. Reacciones con partı́culas alfa (α, n) . . . . . . . . . . . . . . . . 9
1.1.2. Reacciones con partı́culas gamma (γ, n) . . . . . . . . . . . . . . 9
1.1.3. Fisión espontánea . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.4. Reacciones nucleares . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.2. Clasificación de los neutrones . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3. Interacción neutrón materia. . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3.1. Dispersión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3.2. Absorción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2. Blindajes de neutrones 15
2.1. Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.1. Agua . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.2. Polietileno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.3. Water Extended Polyester . . . . . . . . . . . . . . . . . . . . . . 16
2.1.4. Grafito . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.5. Acero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.2. Moderación de neutrones . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3. Atenuación de neutrones . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.4. Sección eficaz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.4.1. Sección eficaz geométrica . . . . . . . . . . . . . . . . . . . . . . 24
2.4.2. Secciones eficaces . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3. Método Monte Carlo 28

3.1. Ecuación de transporte . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.2. Código MCNPX . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.2.1. Estructura del fichero de entrada de MCNP . . . . . . . . . . . . 31
3.2.2. Especificación de la geometrı́a . . . . . . . . . . . . . . . . . . . . 31
3.2.3. Tarjeta de Datos . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2
Índice general 3
4. Reglamento Para El Transporte Seguro de Material Radiactivo 37

4.1. Conceptos Básicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.2. Limites de Actividad y Clasificación . . . . . . . . . . . . . . . . . . . . 38
4.2.1. Clasificación de los Bultos y Embalajes . . . . . . . . . . . . . . 39
4.3. Requisitos y controles para el transporte . . . . . . . . . . . . . . . . . . 40
4.4. Marcado, rotulado y etiquetado . . . . . . . . . . . . . . . . . . . . . . . 42
5. Resultados 45
6. Conclusiones 56
Índice de figuras
1.1. Esquema Dispersión Elástica . . . . . . . . . . . . . . . . . . . . . . . . 12

1.2. Esquema Dispersión Inelástica. . . . . . . . . . . . . . . . . . . . . . . . 12
2.1. Esquemas Sistema Laboratorio (SL) y Sistema Centro de masa (SCM). 18

2.2. Esquema velocidad del neutrón después de la colisión. . . . . . . . . . . 20
2.3. Sección eficaz principal para el Hidrógeno. . . . . . . . . . . . . . . . . . 25
2.4. Sección eficaz principal para el Oxigeno. . . . . . . . . . . . . . . . . . . 26
2.5. Sección eficaz principal para el Carbono. . . . . . . . . . . . . . . . . . . 26
2.6. Sección eficaz principal para el Nitrogeno. . . . . . . . . . . . . . . . . . 26
3.1. Visualización celdas de la geometrı́a del contenedor. . . . . . . . . . . . 32
4.1. Etiquetado Categorı́a I - BLANCA [5]. . . . . . . . . . . . . . . . . . . . 42

4.2. Etiquetado Categorı́a II - AMARILLA [5]. . . . . . . . . . . . . . . . . . 43
4.3. Etiquetado Categorı́a III - AMARILLA [5]. . . . . . . . . . . . . . . . . 43
5.1. Espectro de emisión del 241 Am −9 Be [28]. . . . . . . . . . . . . . . . . . 45

5.2. Contenedor real de transporte para una fuente de 241 Am −9 Be. . . . . 46
5.3. Geometrı́a del contenedor bajo simulación . . . . . . . . . . . . . . . . . 46
5.4. Esquema de la simulación . . . . . . . . . . . . . . . . . . . . . . . . . . 48
5.5. Fluencia de neutrones en los puntos 1 y 9 . . . . . . . . . . . . . . . . . 49
5.6. Fluencia de fotones en los puntos 1 y 9. . . . . . . . . . . . . . . . . . . 50
5.7. Fluencia de neutrones sobre los detectores 1, 2 y 3. . . . . . . . . . . . . 51
5.8. Fluencia de Fotones sobre los detectores 1, 2 y 3. . . . . . . . . . . . . . 51
5.9. Fluencia de neutrones en los puntos 2 y 5. . . . . . . . . . . . . . . . . . 52
5.10. Fluencia de fotones en los puntos 2 y 5. . . . . . . . . . . . . . . . . . . 53
5.11. Esquema de la simulación . . . . . . . . . . . . . . . . . . . . . . . . . . 55
4
Índice de cuadros
1.1. Fuentes de neutrones(α, n) [9]. . . . . . . . . . . . . . . . . . . . . . . . 9

1.2. Clasificación de neutrones [12]. . . . . . . . . . . . . . . . . . . . . . . . 11
2.1. Propiedades moderadoras de diferentes materiales [20]. . . . . . . . . . . 22

2.2. Secciones eficaces de captura. . . . . . . . . . . . . . . . . . . . . . . . . 27
3.1. Ejemplo input celdas y superficies. . . . . . . . . . . . . . . . . . . . . . 32

3.2. Especificación de la fuente. . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.3. Especificación de los tallys. . . . . . . . . . . . . . . . . . . . . . . . . . 34
4.1. Clasificación del numero de las Naciones Unidas para los bultos Tipo A 38
4.2. Categorı́as de los bultos y contenedores. . . . . . . . . . . . . . . . . . . 42
5.1. Especificaciones de los materiales . . . . . . . . . . . . . . . . . . . . . . 47

5.2. Valores de tasa de dosis equivalente para los detectores 1 - 10 . . . . . . 54
5
Introduccion
Las fuentes de neutrones de 241 Am −9 Be son ampliamente utilizadas en muchas

aplicaciones como, el análisis por activación [1], fuentes de calibración [2], la exploración
petrolı́fera, entre otras, porque son relativamente económicas, compactas, portables y
razonablemente constantes en el tiempo [3].
Sin embargo es necesario blindar este tipo de fuentes evitando ası́ los efectos adver-
sos de las radiaciones en el personal ocupacionalmente expuesto y el publico, ya que
algunos de los inconvenientes en el diseño del blindaje radican en la naturaleza de las
radiaciones (neutrones, fotones) y el espectro energético de los neutrones que abarcan
energı́as desde los pocos keV hasta mas de los 10 MeV [4], la selección de los materiales
forma parte fundamental en el diseño del blindaje, permitiendo ralentizar los neutrones
energéticos y absorber los neutrones térmicos mediante dispersiones, que pueden ser
elásticas o inelásticas, como resultado de dichas interacciones se producen partı́culas
secundarias que deben ser tenidas en cuenta en el blindaje.
Por lo anterior es aconsejable discutir las propiedades relevantes de los materiales

mas comúnmente utilizados en los contenedores de transporte para fuentes de neutro-
nes, la principal caracterı́stica de estos materiales son sus núcleos ligeros que permiten
las colisiones elásticas con los neutrones, la baja densidad de los materiales que reducen
el peso y el volumen del contenedor garantizando su capacidad de transporte.
El transporte de materiales radiactivos esta regulado por el Organismo Internacio-

nal de Energı́a Atómica (OIEA) a través del reglamento para el transporte sin riesgos
de materiales radiactivos desde su versión vigente de 1973 [5], y determina que para
este tipo de contenedor los valores de tasa de dosis en superficie, no deben superar los
2000 µSv/h, la elección de los materiales y el diseño del contenedor aseguran que el
blindaje utilizado permite obtener tasas de dosis dentro de los valores estipulados por
la normativa internacional para el transporte de este tipo de material radiactivo.
El diseño del blindaje se estudia mediante el transporte de los neutrones y su inter-

acción con la materia usando el método Monte Carlo y el código MCNPX, es importante
6
Índice de cuadros 7
resaltar que los materiales para el blindaje de neutrones siguen siendo objeto de estudio
y sigue abierta la oportunidad de producir materiales como el WEP establecida como
patente en este campo de investigación [6].
Capı́tulo 1
Fı́sica de Neutrones
Rutherford y sus colegas emplearon partı́culas α de gran energı́a, procedentes de

las sustancias radiactivas, como proyectiles para bombardear los núcleos atómicos con
el propósito de producir desintegraciones nucleares artificiales, este bombardeo según
descubrió rutherford provoca la emisión de partı́culas de protones de todos los elementos
con numero atómico inferior a 19. Al investigar de modo mas riguroso los efectos del
bombardeo del boro y del berilio con partı́culas α se llego a resultados nuevos, Bothe
y Becker (1930) descubrieron la producción de energı́as altamente penetrantes que se
creyeron podı́a ser γ. Curie y Joliot (1932) observaron la producción de protones de
alta energı́a. Chadwick (1931) pudo también demostrar que los rayos emitidos por
el berilio bombardeado provocaban el movimiento muy rápido de los átomos de las
sustancias sobre las que incidı́an (He, Li, Be, C, O y N) resultados que no podı́an ser
explicados si se trataran de partı́culas γ. Finalmente Chadwick (1932) demostró que
para explicar la energı́a de los rayos procedentes del berilio bombardeado estaban, en
realidad, formados por partı́culas de masa próxima a la del protón pero con carga nula.
Como la nueva partı́cula tenia carácter eléctricamente neutro y su masa era cercana al
protón se identifico como el neutrón [7].
1.1. Generación de neutrones
Se pueden generar haces de neutrones por medio de diferentes reacciones nucleares,

sin embargo la mayor producción de neutrones ocurren en las reacciones de fisión que se
realizan de forma regulada en reactores nucleares y que se observan también de forma
espontánea en algunos radionúclidos muy pesados (Z no menores de 90); sin embargo
se producen neutrones por medio de fuentes portátiles que se basan en reacciones (alfa
- neutrón) y (gamma - neutrón).
8
Capı́tulo 1. Fı́sica de Neutrones 9
1.1.1. Reacciones con partı́culas alfa (α, n)
La reacción responsable del descubrimiento del neutrón puede ser usada para pro-
ducir fuentes de neutrones portátiles usadas en laboratorios, en este caso el isótopo del
berilio estable 9 Be tiene un neutrón relativamente débil (energı́a de enlace 1.7 MeV).
Si una partı́cula α con decaimiento (5 - 6 MeV) choca contra un núcleo de 9 Be, un
neutrón sera liberado con un valor Q de 5.7 MeV de la forma,
He42 + Be94 → C612 + n10 (1.1)
Durante mucho tiempo los emisores alfa mas utilizados fueron (226 Ra y 210 P o).
Las fuentes de Ra-Be eran buenas por su amplia duración, pero peligrosas por la alta
emisión de radiación gamma procedentes de los decaimientos del 226 Ra y las fuentes
de Po-Be eran buenas por su baja emisión de radiación gamma pero su corta vida
provocaba que alrededor de un año la fuente quedaba reducida en una décima parte.
Ahora existen opciones mucho mejores como las fuentes de 241 Am −9 Be que po-
seen una alta duración y poca emisión gamma, de este modo las fuentes de neutrones
recientes reúnen las ventajas de las fuentes antiguas pero sin tener sus desventajas.
La limitación de este tipo de fuentes (α − n) consiste en la producción de neutrones
con amplios espectros de energı́a [8].
Algunos ejemplos de fuentes de neutrones bajo la reacción (α, n) se muestran en la
tabla 1.1,
Fuente Vida Media

210 P o − Be 138 dı́as
210 P o −B 138 dı́as
241 Am − Be 432,7 años
226 Ra − Be 1600 años
226 Ra −B 1600 años
239 P u − Be 24100 años
Cuadro 1.1: Fuentes de neutrones(α, n) [9].
1.1.2. Reacciones con partı́culas gamma (γ, n)
Este proceso es similar a las reacciones (α, n) ya expuestas en este trabajo, se

hace la mezcla de un emisor gamma con berilio o con oxido de deuterio, sin embargo
este tipo de fuentes presentan una ventaja sobre las alfa porque producen neutrones
prácticamente monoenergéticos siempre y cuando se presenten partı́culas gamma que
superen la energı́a de ligadura del neutrón que para el caso del berilio es (1.7 MeV)
y deuterio (2.23 MeV). Este tipo de fuentes (gamma - neutrón) son conocidas como
fotoneutrónicas.
Esta reacción se representa de la forma,
γ +9 Be →8 Be + n (1.2)
La fuente fotoneutrónica mas empleada es el antimonio - berilio (124 Sb −9 Be),

aunque su emisión gamma es muy alta, su rendimiento neutrónico es de los mas altos
(1.9 x 105 n/s Ci), casi todos los neutrones poseen una energı́a de 22 keV y una octava
parte son de 418 keV [10].
1.1.3. Fisión espontánea

Una fuente común de neutrones se produce mediante la fisión espontánea, entre
los núclidos que presentan este mecanismo se toman aquellos cuyos valores de Z sean
iguales o mayores a 90, este tipo de reacción consiste en la división de un átomo pesado
en dos átomos mas pequeños (fragmentos de fisión), de masas comparables y en dos o
mas neutrones. Uno de los isótopos que presentan fisión espontánea es el 252 Cf con un
periodo de semi desintegración de 2.65 años, los neutrones son directamente generados
en el proceso de fisión con una velocidad cercana a 4 por cada fisión, la fisión ocurre
en solo el 3 % de los decaimientos (el decaimiento α representa el resto), la velocidad
de producción de los neutrones es de 2.3 x 1012 neutrones/s por gramo de 252 Cf o 4.3
x 109 n/s por Ci de 252 Cf [11].
235 U →94 Sr +139 Xe +1 n +1 n
1.1.4. Reacciones nucleares

Hay por supuesto muchas reacciones nucleares que pueden producir neutrones, di-
chas reacciones requieren un acelerador que permita producir un haz de partı́culas para
iniciar la reacción y, por lo tanto, no son tan convenientes como el tipo de fuentes que
se dan mediante los decaimientos radiactivos anteriormente descritos. Sin embargo,
seleccionando cuidadosamente la energı́a incidente y el ángulo en el que observamos
el neutrón emitido, podemos obtener un haz razonablemente mono energético de casi
cualquier energı́a deseada [11].
Algunas reacciones que podrı́an utilizarse son:
3H + d →4 He + n Q = + 17.6 MeV
7 Li + p →7 Be + n Q = - 1.6 MeV
2H + d →3 He + n Q = + 3.3 MeV
1.2. Clasificación de los neutrones

La clasificación de los neutrones se realiza según la necesidad de detalle que se re-
quiera de su espectro energético como se observa en la tabla 1.2. Para las aplicaciones
más comunes, es habitual clasificarlos englobando los neutrones denominados epitérmi-
cos, lentos e intermedios, cabe destacar el grupo de neutrones térmicos, definiéndolos
como los neutrones que se encuentran en equilibrio térmico con el ambiente. Este grupo
neutrónico se caracteriza por poseer energı́as del orden de magnitud de la energı́a de
agitación térmica de los átomos que le rodean y su incapacidad para ceder más energı́a
al entorno[12].
Clasificación Energı́a (eV ) Velocidad (m/s)

Neutrones Frı́os 0.005 9, 66.102
Neutrones Térmicos 0.025 2, 2.103
Neutrones Epitérmicos 1 1, 4.104
Neutrones Intermedios 104 1, 4.106
Neutrones Rápidos 106 1, 4.107
Neutrones Ultrarrapidos 108 1, 4.108
Cuadro 1.2: Clasificación de neutrones [12].
1.3. Interacción neutrón materia.

Cuando un neutrón se aproxima al núcleo de un átomo diferentes fenómenos pue-
den suceder, el cambio en la trayectoria del núcleo sin modificarse su numero o masa
atómica (esta reacción es llamada dispersión), el cambio en las velocidades del neutrón
(dispersión) y el núcleo o la captura del neutrón por parte del núcleo (absorción).
1.3.1. Dispersión
Frecuentemente las partı́culas pueden ser dispersadas por impactos con los núcleos
atómicos de dos modos diferentes, uno es por dispersión elástica y el otro es por dis-
persión inelástica.
Dispersión elástica
En este caso tanto la energı́a cinética total como el momento total antes del choque,
son iguales, respectivamente, a la energı́a cinética total y al momento total después del
impacto.
Si un neutrón rápido choca contra una partı́cula ligera como un protón, o un deuterón
el núclido puede acoger una cantidad considerable de energı́a, si por el contrario el
neutrón choca con un núclido pesado, las leyes de la conservación permiten observar
que se puede otorgar muy poca energı́a cinética al núcleo pesado, y que el neutrón
tomara una nueva dirección, con casi toda su energı́a inicial.
Figura 1.1: Esquema Dispersión Elástica
De tal forma, si se desea frenar los neutrones rápidos hasta velocidades relativamen-
te pequeñas, mediante la dispersión elástica, los materiales deben ser preferiblemente
compuestos por átomos de bajo peso atómico como el caso del hidrogeno [13].
Hay dos tipos generales de colisiones elásticas entre neutrones y núcleos atómicos,
el primero implica la captura del neutrón para formar un núcleo compuesto, el segundo
no presenta formación de núcleos compuestos.
Dispersión Inelástica.
Durante este proceso el neutrón choca contra un núcleo que puede ser pesado, el
neutrón pierde una cantidad considerable de energı́a y el núcleo queda en un estado
excitado, el neutrón abandona el núcleo con una nueva dirección y con menor valor de
su energı́a inicial. Los núcleos excitados, al retomar su estado normal emite radiación γ.
Figura 1.2: Esquema Dispersión Inelástica.
Los choques inelásticos obedecen la conservación del momento dando lugar al re-
troceso de la partı́cula sin embargo la energı́a cinética no se conserva. En estos eventos
el neutrón impacta al núcleo y luego de la colisión el núcleo es llevado a un estado de
excitación y regresa a su estado base emitiendo un rayo γ, en este caso no se conserva

la energı́a cinética, aunque si se conserva la energı́a total [14].
XZA + n10 → (XZA ) ∗ +n10 (1.3)
(XZA )∗ → XZA + γ
1.3.2. Absorción
Son reacciones nucleares en las que el neutrón se integra en la estructura nuclear del
núcleo blanco, dando lugar a un isótopo distinto del mismo elemento. La naturaleza del
isótopo resultante, si es estable o inestable, puede dar lugar a cadenas de desintegración
posteriores, pudiendo resultar en elementos estables diferentes al inicial[11]. Entre las
reacciones de absorción se estudian los procesos de:
Captura radiactiva
Captura no radiactiva
Captura de activación
Fisión
Captura radiactiva
En las reacciones de captura radiactiva, o reacciones (n, γ), el núcleo compuesto

excitado emite su excedente de energı́a en forma de radiación gamma, denominada ra-
diación gamma de captura, dejando finalmente el núcleo compuesto en su estado de
mı́nima energı́a (estado fundamental). El proceso se representa simbólicamente median-
te la ecuación,
XZA + n10 → [XZA+1 ]∗ → XZA+1 + γ (1.4)
Donde XZA es el blanco, de número atómico Z y numero másico A. El producto de

la reacción de captura radiactiva XZA+1 es un isótopo del blanco X, puesto que tiene
número atómico igual y numero másico una unidad mayor, el núclido X A+1 puede o
no ser radioactivo según las circunstancias[12].
Este tipo de interacción de absorción predomina el estudio necesario para la elabo-
ración de este trabajo.
Captura no radiactiva
Se forma un isótopo del núcleo blanco sin emitirse radiación por quedar el producto
en estado fundamental. La reacción de este tipo viene dada por la ecuación (1.5). Para
que se pueda dar este tipo de captura es necesario que la energı́a del neutrón incidente
sea del orden de magnitud de la diferencia de energı́a entre los dos isótopos [12].
XZA + n10 → XZA+1 (1.5)
Captura de activación
El producto resultante de este proceso de absorción es inestable y, de acuerdo con

sus caracterı́sticas (vida media y mecanismo de desintegración) alcanzara un estado
fundamental como otro elemento diferente al de partida[12]. Según dichas caracterı́sticas
se puede clasificar como:
Captura con desintegración β
Captura con emisión de partı́culas α
Captura con emisión de protones
Fisión
El neutrón es absorbido por un núcleo pesado que después se divide en dos frag-
mentos o productos de fisión con energı́as cinéticas del orden de 100 MeV, liberándose
además un número variable de neutrones [12].
Capı́tulo 2
Blindajes de neutrones
Entre las medidas defensivas contra la irradiación se encuentran tres elementos

fundamentales: distancia, tiempo y blindaje. De estos solo uno es intrı́nseco al diseño de
la instalación, un diseño adecuado asegura que las tasas de dosis que los trabajadores
ocupacionalmente expuestos y el publico reciben se encuentren dentro de los limites
seguros.
El blindaje radiológico es un material que se ubica en medio de la fuente emisora de
radiación y un punto determinado, con fin de reducir las radiaciones emitidas desde la
fuente y que llegan al punto [15], en el caso de la radiación neutronica los mejores blin-
dajes son aquellos que tienen materiales de núcleos ligeros, por dos razones: primera,
al ser pequeña la masa de los núcleos, el choque elástico absorbe una buena propor-
ción de la energı́a del neutrón incidente, por lo tanto se necesitan menos choques, es
decir menor espesor de material para lograr termalizar los neutrones, como ya se vera
en la tabla 2.1. Segunda para el caso de los núcleos ligeros las secciones eficaces de
captura de neutrón son muy pequeñas, o bien el producto es estable, ası́ que no se pre-
sentara activación del blindaje.Entre los blancos ligeros que cumplen las propiedades
antes descritas se encuentran los núclidos de hidrogeno, carbono, oxigeno, nitrógeno y
silicio[8].
2.1. Materiales
La elección de los materiales para el contenedor de transporte de una fuente de

neutrones 241 Am −9 Be debe permitir la atenuación y moderación del espectro de neu-
trones emitidos por la fuente, estos materiales tambien deben ser poco densos logrando
ası́ reducir tanto volumen y peso como sea posible, garantizando las propiedades de
transporte del contenedor, dentro de los posible materiales utilizados dentro del blin-
daje de neutrones se tienen:
15
Capı́tulo 2. Blindajes de neutrones 16
2.1.1. Agua
Debido a su alta densidad en átomos de hidrogeno, el agua es usualmente utilizada
en el blindaje de neutrones, los átomos de hidrogeno disminuyen las energı́as de los
neutrones mediante dispersión elástica, suficiente para moderar los neutrones a sus
estados térmicos donde son rápidamente capturados. El agua tiene la ventaja de ser
económica y manejable, no es toxica y en la mayorı́a de estados es estable, sin embargo
a altas temperaturas se evaporiza. La emisión de partı́culas γ por la captura en los
blindajes con agua puede ser disminuida con la presencia de boro [15].
2.1.2. Polietileno
Este plástico es un hidrocarburo puro (CH2 )n con una densidad 0.92 g/cm3 , el po-
lietileno en realidad contiene un 18 % en volumen más hidrógeno que el agua, mientras
el agua contiene 6, 7x1022 átomos de hidrógeno/cm3 , el polietileno tiene aproximada-
mente 8x1022 átomos de hidrógeno/cm3 . El polietileno no se activa, pero se ablanda a
110 °C y soporta la combustión [16].
2.1.3. Water Extended Polyester

El water extended polyester es una resina que consiste en partes iguales de monóme-
ro de estireno y resina de poliéster combinada con hasta 65-70 % en peso de agua con
una densidad de masa promedio de 1.1 g/cm3 , el wep en comparación con el agua posee
una cantidad de hidrogeno significativa y logra superar a la lucite [17].
2.1.4. Grafito
Dado que está disponible en forma muy pura, exhibe pocos rayos gamma secundarios
a diferencia del agua, posee poco poder de atenuación pero es ampliamente utilizado
como moderador, también se encuentra en escudos térmicos y ocasionalmente como
material de blindaje primario en reactores [15].
2.1.5. Acero
El acero es utilizado en los blindajes de neutrones por sus propiedades estructurales
garantizando que los componentes del blindaje se encuentren contenidos de la mejor
manera posible, el acero se utiliza además por sus propiedades atenuadoras frente a los
rayos γ, en aceros inoxidables su estabilidad fisica no presenta compromisos soportando
temperaturas elevadas y ofreciendo una resistencia a la corrosión [15].
Los materiales que se utilizan en los blindajes de fuentes radiactivas deben tener
las siguiente caracterı́sticas adicionales:
Alta resistencia al fuego o temperaturas elevadas.
No debe ser toxico, ni permitir la emisión de gases tóxicos al elevarse la tempe-

ratura.
El material debe ser durable en el tiempo y permanecer lo mas inoloro posible.
No deben producirse deformaciones hasta alrededor de los 150°C.
El material debe ser fácilmente removible, reemplazado o reparado.
2.2. Moderación de neutrones
El proceso de moderación permite reducir las velocidades de los neutrones mediante

procesos de dispersión por choque elástico producida por los núcleos ligeros y disper-
sión por choque inelástico producida por los núcleos de masa media y elevada, durante
estos choques los neutrones pierden parte de su energı́a hasta que alcanzan el equilibrio
térmico con el medio y son capturados por los núcleos del material presente en el medio
[18].
Para la mayorı́a de los fines prácticos el proceso importante en la moderación se da

para neutrones con energı́a superior a 1 MeV.
Resolviendo las ecuaciones establecidas en el principio de la conservación de la

energı́a y la cantidad de movimiento del neutrón antes y después de la colisión puede
determinarse la cantidad de energı́a perdida en una colisión, para ello se utilizaran dos
sistemas de referencia, el primero es el sistema laboratorio en el que se considera al
núcleo en reposo antes de la colisión y el neutrón incidente se aproxima al mismo; el
segundo sistema de referencia que se va a utilizar es el denominado centro de masas o
centro de gravedad, donde se considera que el centro de masa del conjunto constituido
por el núcleo y el neutrón esta en reposo, y tanto el neutrón como el núcleo se aproximan
a el [18].
Figura 2.1: Esquemas Sistema Laboratorio (SL) y Sistema Centro de masa (SCM).
En la figura 2.1 se puede observar la relación entre los sistemas de referencia (labo-
ratorio y centro de masa) en el primero antes de la colisión el neutrón de masa m se
mueve hacia el núcleo de masa M con una velocidad V0 y una energı́a E0 . En el segundo
se observa Vcm como la velocidad del neutrón y Vc la velocidad del núcleo.
De la definición del centro de masa
m1 r1 + m2 r2
Rcm = (2.1)
m1 + m2
La velocidad del centro de masa respecto al sistema laboratorio,
m1 ṙ1 + m2 ṙ2
Vcm = Ṙcm = (2.2)
m1 + m2
Como el núcleo en el sistema laboratorio se encuentra en reposo V2 = 0,
m1 V0 + m2 (0)
Vcm = (2.3)
m1 + m2
m
Vcm = V0 (2.4)
m+M
Por adición Galileana,
M
Vcm = V0 − Vcm = V0 (2.5)
m+M
Como en el sistema centro de masas el momento total es nulo Pcm = 0
mVcm = −M Vc (2.6)
m m
Vc = Vcm = V0 (2.7)
M m+M
Se pueden presentar 3 casos durante las colisiones de los neutrones con los núcleos,
en el primero se evalúa una colisión rasante donde el ángulo de dispersión φ = 0
V = Vcm + Vc (2.8)
La cantidad de energı́a perdida por el neutrón es despreciable, por lo que se puede

considerar que :
E = E0
la velocidad del neutrón después de la colisión sera,
M m
V = V0 + V0 (2.9)
m+M m+M
Para el caso donde la colisión genera un ángulo φ = 180°, el neutrón retrocede en
sentido opuesto sobre su dirección original en este caso hay una perdida de energı́a
durante la colisión.
V = Vcm − Vc (2.10)

M m M −m
V = V0 − V0 = V0 (2.11)
m+M m+M m+M
La perdida de energı́a durante la colisión sera la razón entre la energı́a inicial y final
1
E mV 2
= 21 2
(2.12)
E0 2 mV0
2
V2

E M −m
= 2 = (2.13)
E0 V0 m+M
Para colisiones que generan un ángulo de dispersión φ entre 0 y 180°, la velocidad
del neutrón después de la colisión se calcula en función del ángulo aplicando la ley del
coseno en la solución de triángulos oblicuángulos.
V 2 = Vcm
2
+ Vc2 − 2Vcm Vc cos φ (2.14)
2 2
2 M m Mm
V = V0 + V0 −2 cos φ (2.15)
m+M m+M m+M
De la misma forma que en el proceso anterior se observa la perdida de energı́a
durante la colisión,
1
E mV 2
= 21 2
(2.16)
E0 2 mV0
E M 2 + m2 + 2M m cos φ
= (2.17)
E0 (m + M )2
M
si se tiene en cuenta la relación m = A entre la masa del núcleo y la masa del
neutrón, la ecuación anterior toma la forma,
E A2 + 1 + 2A cos φ
= (2.18)
E0 A2 + 1
Conviene expresar la relación entre las energı́as en función de la cantidad
(A − 1)2
r= (2.19)
(A + 1)2
Se puede expresar la máxima energı́a que pierde un neutrón en una sola colisión trans-
formando la ecuación (2.18)
E 1+r 1−r
= + cos φ (2.20)
E0 2 2
De este modo la máxima perdida de energı́a se da cuando φ = 180 asi cos φ = −1

y E = rE0 ; por el contrario, si φ = 0 el cos φ = 1 y E = E0 .
Relacionando los ángulos de dispersión φ del sistema centro de masa y el ángulo θ del
sistema laboratorio,
Figura 2.2: Esquema velocidad del neutrón después de la colisión.
Puede relacionarse el ángulo de dispersión φ, según el ángulo del sistema centro de

masa y con el ángulo de dispersión en el sistema laboratorio como puede observarse en
la figura 2.2; Con C y B definidos como:
C = Vcm cos φ + Vc
B = Vcm sen φ (2.21)
Remplazando los valores de Vcm y Vc en (2.21)

M m
C= V0 cos φ + V0
m+M m+M
M
B= V0 sen φ (2.22)
m+M
C
Se deduce que la cot θ = B por lo que,
cos φ + A1
cot θ = (2.23)
sen φ
A partir del valor de la cotangente puede obtenerse el coseno del ángulo mediante
la relación:
cot θ
cos θ = 1 (2.24)
(1 + cot2 θ) 2
Al dar a cot θ el valor definido por la ecuación (2.23), se puede deducir que
1 + A cos φ
cos θ = 1 (2.25)
(1 + A2 + 2A cos φ) 2
de esta forma se obtiene la relación entre los ángulos [7]

Se dice que la dispersión es el numero de la fracción total superficie que pasa a
través de un area asociada con un rango de ángulo dφ y φ, el area absoluta de un anillo
de radio sen φ y ancho dφ en una esfera con radio de una unidad es 2π sen φdφ. La
1
fracción de área total de la esfera es 4π esto es 2 sen φdφ = −d(cos φ)/2. El numero
diferencial de dispersión a través de esta área es:
−d(cos φ)
dn = n (2.26)
2
Ahora el promedio del coseno del ángulo generado en la dispersión es,

R
cos θdn
cos θ = (2.27)
n
Es necesario involucrar las ecuaciones (2.25) (2.26) y ası́ dejar la integral en términos
de una variable común cos φ, la integral de la ecuación (2.27) sobre la variable cos φ
para -1 a 1 es,
2
cos θ = (2.28)
3A
La ecuación (2.28) permite observar el ángulo promedio de la dispersión [19].
Cuando el valor de A es grande los núcleos dispersores son pesados esto quiere decir
que cos θ es pequeño y la dispersión, según el sistema laboratorio es prácticamente
isótropa; cuando A es pequeño como ocurre con los núcleos ligeros, se dispersan mas
neutrones en la dirección original de su movimiento que hacia la la opuesta.
Ahora puede calcularse la perdida media de energı́a por colisión que se representa
abreviadamente por ξ, es conveniente utilizar esta magnitud en el calculo de la mode-
ración de los neutrones dada su independencia frente a la energı́a de los mismos. como
la ecuación (2.20) es una función lineal del cos φ y todos los valores son igualmente
E
posibles, se deduce que lo mismo ocurre con los valores de E0 .
De esta forma la probabilidad P dE para un neutrón cuya energı́a inicial es E0 tenga
una energı́a comprendida entre E y E + dE después de una colisión esta dada por,
dE
P dE = (2.29)
E0 (1 − r)
Por definición,
E0
ξ = ln E0 − ln E = ln (2.30)
E
De esta forma,
Z Z
E0 E0 dE
ξ= ln P dE = ln (2.31)
E E E0 (1 − r)
Resolviendo la integral,
2 2
A−1 A−1
r A+1 ln A+1
ξ =1+ ln r = 1 + 2 (2.32)
1−r

A−1
1− A+1
ξ puede transformarse en,
(A − 1)2

A−1
ξ =1− ln (2.33)
2A A+1
El promedio de ξ sera,
2
ξ= 2 (2.34)
A+ 3
Se presenta en la tabla 2.1 las propiedades moderadoras de diferentes materiales
Elemento A θ ξ N° de choques
Hidrógeno 1 48° 1 18
Deuterio 2 70.5° 0.7261 25
Berilio 9 86° 0.2078 87
Carbono 12 87° 0.1589 115
Uranio 238 89° 0.0084 2160
Cuadro 2.1: Propiedades moderadoras de diferentes materiales [20].

Donde A es el numero atómico, θ corresponde al promedio del coseno del ángulo,

ξ es el valor logarı́tmico de la tasa de cambio de la energı́a y se presenta el numero
de colisiones necesarias para reducir la energı́a de un neutrón de 2 MeV a energı́as
térmicas, 0,025 eV [20].
2.3. Atenuación de neutrones

Se define la atenuación de neutrones como el medio por el cual el número de neutro-
nes inicial disminuye durante el proceso de interacción pero sus energı́as no cambian,
a diferencia de la moderación de neutrones donde las energı́as de estos varı́an durante
el proceso de interacción, en la moderación el número de neutrones inicial se mantiene
durante todo el proceso [18]. Dichas posibilidades (atenuación y moderación) son go-
bernadas principalmente por las probabilidades estadı́sticas donde un neutrón puede o
no interactuar por unidad de distancia recorrida, esta probabilidad denominada µ es
llamada el coeficiente de atenuación lineal que depende de la energı́a del neutrón y el
material perforado, este coeficiente tiene dimensiones inversas a la longitud[18].
EL coeficiente de atenuación esta relacionado con la sección eficaz σ como:
µ = σn (2.35)
Donde n, es el numero de átomos por unidad de volumen del medio en el cual se

mueve el neutrón, n se calcula como:
ρNA
n= (2.36)
M
Donde ρ es la densidad del medio, NA es el numero de avogadro, M es el peso

molecular.
Las expresiones (2.35) y (2.36) permiten obtener los coeficientes de atenuación de
cada componente atómico del medio en el que se mueve un neutrón a diferentes energı́as.
El proceso de atenuación para los neutrones monoenergéticos es de tipo exponencial,
de tal forma que la ecuación que expresa el numero de neutrones que atraviesan un
espesor x de material es:
N (x) = N0 e−µx (2.37)
Donde N (x) representa el numero de neutrones que atraviesan dicho espesor sin
interactuar y N0 el numero de neutrones incidentes.
La distribución de probabilidad de la longitud del camino recorrido por el neutrón
sin interactuar durante el espesor del material viene dada,
P (x) = µe−µx (2.38)

La longitud que recorre un neutrón sin interactuar se denomina recorrido libre medio
y se denota λ, y viene dada por el reciproco de la sección eficaz macroscópica de forma
que:
1
λ= (2.39)
Σ
2.4. Sección eficaz

2.4.1. Sección eficaz geométrica
Un gran número de experimentos sobre las colisiones entre partı́culas atómicas
indican que los núcleos se pueden estimar como un conglomerado de esferas muy cer-
canas las unas a las otras, si el tamaño aproximado del radio medio del nucleón es
ro = 1, 2x10−13 cm.
De esta forma el radio del núcleo se puede calcular mediante la siguiente ecuación [21],
√
3
R = ro A
Para este caso la medida obtenida sirve como una aproximación al tamaño del blanco,
conocida como sección eficaz geométrica y esta dada por
σg = πR2
2.4.2. Secciones eficaces

Cuando los neutrones inciden sobre algún material, todos los núcleos del blanco no
se comportan siempre como si sus secciones eficaces fueran constantes. Por ejemplo en
la captura de un neutrón, previa a una desintegración radiactiva, puede parecer que el
tamaño nuclear es completamente diferente de la sección eficaz geométrica calculada.
Para los neutrones rápidos, una sección eficaz nuclear puede ser relativamente pequeña,
mientras que los neutrones lentos pueden tener secciones eficaces grandes.
Por lo tanto la sección eficaz no es el área del blanco en sentido literal, es mas un
concepto simbólico que busca expresar la probabilidad de interacción, ası́ la probabili-
dad de interacción entre un núcleo y un neutrón se llama sección eficaz microscópica σ
y se refiere al área de blanco efectiva.
Si un haz de neutrones incide sobre un cm3 del material del blanco, la sección eficaz
será nσ, donde n es el numero total de átomos en el cm3 , este producto es la sección
eficaz macroscópica y se designa como,
X
= nσ
De esta forma la sección eficaz total sera la suma de las diferentes secciones efica-
ces, cada una presente en los diversos procesos que alteran la energı́a o el momento del
neutrón.
σt = σs + σc
Donde σs es la sección eficaz producida por la dispersión y σc la sección eficaz para la

captura radiactiva [13].
Como se observo antes, los neutrones interaccionan con la materia mediante dife-
rentes procesos que se rigen por los valores de los coeficiente de atenuación de cada
material y que dependen de sus secciones eficaces, en las figuras 2.3,2.4 y 2.5 se mues-
tran las secciones eficaces de los diferentes elementos utilizados en los materiales que
sirven como blindaje.
Figura 2.3: Sección eficaz principal para el Hidrógeno.
En la simulación que realizamos en este trabajo se han utilizado secciones eficaces

tomadas de la referencia[], para cada uno de los diferentes elementos con energı́as de
hasta 10 MeV. La sección eficaz total es la suma de los diferentes procesos de interacción
para el caso del hidrogeno, como se observa en la figura 2.3, la sección de dispersión
elástica predomina a lo largo de las diferentes energı́as ası́ como la producción de fotones
gamma, y se observa la ausencia de la sección de absorción.
Figura 2.4: Sección eficaz principal para el Oxigeno.
Figura 2.5: Sección eficaz principal para el Carbono.
Figura 2.6: Sección eficaz principal para el Nitrogeno.
Dentro de los elementos que constituyen los medios de interés el oxigeno y el carbono
forman parte fundamental en la constitución de cada uno de los materiales utilizados
en el diseño del contenedor de transporte. Las figuras 2.4 ,2.5 y 2.6 muestran la sección
eficaz total de cada elemento, la suma de la sección eficaz de dispersión elástica cuyo
valor es superior a las secciones de absorción y producción gamma. Estas secciones
eficaces determinan que la probabilidad de interacción que predomina en cada elemento
es la dispersión elástica es clave durante la moderación y la atenuación de neutrones.
En la tabla 2.2 se presentan las secciones eficaces de captura de unos pocos elementos
[13].
Elemento Masa Atómica (g/mol) σc (barns)

Hidrogeno 1,01 0,33
Deuterio 2 0,00046
Carbono 12,01 0,0032
Berilio 9,02 0,010
Aluminio 27 0,23
Hierro 55,85 2,5
Cobre 63,54 3,6
Cuadro 2.2: Secciones eficaces de captura.

Capı́tulo 3
Método Monte Carlo
El método de Monte Carlo (MC) es un método numérico que permite resolver

problemas matemáticos mediante la simulación de variables aleatorias, las palabras
((variables aleatorias)) suelen emplearse para subrayar el hecho de no conocer el valor
que tomará esta variable en un caso concreto dado, pero sabemos qué valores puede
tomar y también las probabilidades de unos u otros valores, en base a estos datos
no podemos predecir con exactitud el resultado de una prueba relacionada con esta
variable aleatoria pero sı́ podemos prever con gran seguridad los resultados de un gran
número de pruebas, cuanto mayor sea el número de prueba, mayor exactitud tendrán
las predicciones; se considera como fecha de nacimiento del metódo el año 1949 cuando
apareció el artı́culo titulado (the Monte Carlo method) [22] la creación de este metódo
suele asociarse a los nombres de los matemáticos norteamericanos J. Von Neumann y
S. Ulam [23].
Los cálculos MC han sido ampliamente utilizados en todas las disciplinas involu-
cradas en análisis cuantitativos sin embargo es el transporte de radiación una de las
aplicaciones del método MC mas utilizadas, el transporte de radiación se refiere a la
difusión de partı́culas tales como electrones, fotones, neutrones, iones e incluso átomos
neutros dentro de un medio como el aire, el agua, una estrella, la atmósfera de la tierra
y hasta el cuerpo humano. El punto es que pequeñas partı́culas como los neutrones en
un reactor o los rayos X en el cuerpo interactúan con el medio, el problema es poder ob-
tener información sobre el campo de las partı́culas por ejemplo, el numero de partı́culas
en un determinado volumen o la energı́a transferida por las partı́culas al medio.
3.1. Ecuación de transporte

La ecuación de transporte de radiación puede ser explicada mediante la siguiente
relación,
(a) = − (b) − (c) + (d) + (e) (3.1)
28
Capı́tulo 3. Método Monte Carlo 29
Estos 5 términos tienen las siguientes interpretaciones (a). la tasa de aumento del
numero de partı́culas en un volumen V con unas energı́as y direcciones especificas, (b).
el numero neto de tales partı́culas que fluyen de el V a través de una superficie por
unidad de tiempo, (c). el numero total de partı́culas que sufren colisiones en el V por
unidad de tiempo (estas partı́culas cambian sus energı́as y direcciones o simplemente
son absorbidos), (d). el numero de partı́culas secundarias que son producidas por las
interacciones en el V, (e). numero de partı́culas que son producidas en el V por unidad
de tiempo por fuentes independientes al campo de radiación.
La relación de balance puede ser expresada como:
Ω.∇φ(r, E, Ω) + µ(r, E)φ(r, E, Ω) =

Z ∞ Z
dE 0 dΩ0 µ(r, E 0 → E, Ω0 → Ω)φ(r, E 0 , Ω0 ) + S(r, E, Ω). (3.2)
0 4π
Esta ecuación integro diferencial es conocida como la ecuación de transporte de

Boltzmann y lleva el nombre de Ludwig Boltzmann quien fue la primera persona en
derivar la ecuación hace mas de un siglo, esta ecuación sirve como una descripción
precisa del campo de radiación en todos los problemas de transporte para fotones y
neutrones, sin embargo es suficientemente compleja para ser resuelta analı́ticamente
excepto para casos simples, el método monte carlo permite obtener la solución para
esta ecuación en los problemas reales del transporte de partı́culas.
Esta ecuación o una variante de la misma es la base de muchos de los grandes
códigos de transporte que se utilizan rutinariamente hoy en dı́a, esta ecuación es bas-
tante general y es aplicable a cualquier geometrı́a, para fotones y neutrones y para la
interacción de las partı́culas con el medio, de la ecuación 3.2 tenemos,
El coeficiente de interacción µ(r, E) y µ(r, E 0 → E, Ω0 → Ω) que depende del

material que compone el medio y las interacciones de interés.
dE 0 que es el diferencial de la energı́a
El flujo de densidad angular φ(r, E, Ω)
Ω.∇φ que expresa las coordenadas utilizadas en el sistema
S(r, E, Ω) que es la distribución de la fuente de radiación primaria, normalmente

incluye la presencia de radionucleidos.
En una simulación de transporte de partı́culas neutras de Monte Carlo, primero se

especifica la geometrı́a del sistema, normalmente por geometrı́a combinatoria, aunque
se pueden usar otros enfoques. Luego, mediante el muestreo de muchas distribuciones,
se puede obtener una simulación completa de una partı́cula a medida que migra a
través del espacio de fase, el muestreo de la fuente es necesario para escoger el punto
de partida inicial, la dirección y la energı́a de las partı́culas.
La partı́cula es rastreada hasta que es absorbida o deja el lı́mite del problema. A
medida que una partı́cula se mueve a lo largo de su trayectoria, se actualizan varias
cuentas para que, después de muchas historias, se pueda estimar alguna propiedad
deseada del campo de radiación.
La mayorı́a de las simulaciones de Monte Carlo emplean una o más técnicas de
reducción de la varianza.
3.2. Código MCNPX
Generalmente uno de los códigos de transporte de radiación mas utilizado es MCNPX

un acrónimo para Monte Carlo N-Partı́cula y el X es la versión del código, esta ver-
sión a diferencia de otras es capaz de simular el transporte de 35 partı́culas entre las
que se puede resaltar los protones, particulas alfa, neutrones entre otras y es capaz del
transporte enlazado entre muchas de estas partı́culas (neutrón, fotón) donde los fotones
son producido por interacciones de neutrones, (neutrones, fotones, electrones) (fotones,
electrones) o (electrones, fotones).
MCNP es soportado por la mayorı́a de sistemas operativos entre Microsoft Win-
dows, Mac OS X, Linux para 32 y 64 bits, adicional a esto MCNP puede ser compilado
mediante multi procesadores, los principales autores para MCNPX son Jhon S. Hen-
dricks, Gregg W. Mckinney, Michael L. Fensin, Franz X. Gallmeie, Michael R. James,
Denise B. Pelowitz, Laurie S. Waters and M. William Johnson.[Exploring monte carlo
methods]
Por la naturaleza general de este código, puede ser utilizado en muchos campos de
estudio donde es inherente la necesidad de conocer las interacciones de la radiación con
la materia, las aplicaciones comunes incluyen pero no se limitan a la fı́sica médica, al
estudio de reactores, a los ciclos de combustible nuclear, la radiografı́a, la producción de
isótopos, el desarrollo de detectores y el análisis ambiental. MCNPX tiene aplicaciones
adicionales en los aceleradores de partı́culas y la fı́sica de altas energı́as.
MCNP permite transportar neutrones con energı́as desde 0 a 150 MeV para ciertos
isótopos, para fotones permite el transporte con energı́as desde 1 keV a 100 GeV y para
electrones las energı́as puedes ser de 1 keV hasta 1 GeV. MCNPX extiende el rango de
energı́a de las partı́culas anteriores utilizando modelos fı́sicos de sección transversal.
y transporte de neutrones, fotones y electrones con la materia, la versión utilizada

para este trabajo es MCNPX [24].
3.2.1. Estructura del fichero de entrada de MCNP
El código permite crear un archivo con diferentes especificaciones, este archivo con-
tiene información sobre el problema en áreas como:
La especificación de geometrı́a.
La descripción de los materiales y la selección de las evaluaciones transversales.
La ubicación y caracterı́sticas de la fuente de neutrones, fotones o electrones.
El tipo de respuestas o cifras deseadas.
Cualquier técnica de reducción de la varianza utilizada para mejorar la eficiencia.
3.2.2. Especificación de la geometrı́a
Las geometrı́as en MCNP trabajan bajo una configuración tridimensional arbitraria

de materiales definidos por el uso de celdas geométricas delimitadas por superficies de
primer y segundo grado y toros elı́pticos de cuarto grado. Las celdas están definidas
por las intersecciones, uniones y complementos de las regiones limitadas por las superfi-
cies. Las superficies se definen proporcionando coeficientes a las ecuaciones superficiales
analı́ticas o, para ciertos tipos de superficies, puntos conocidos en las superficies. MCNP
tiene una geometrı́a más general que la que está disponible en la mayorı́a de los códigos
de geometrı́a combinatoria [25].
La especificación de la geometrı́a que se realizo para esta simulación utiliza las
superficies generadas por macrocuerpos que van desde conos troncados, esferas, cilindros
entre otros.
Se realizaron un total de 25 superficies y 25 celdas que generaron la geometrı́a del
contenedor figura 3.1 y un volumen mundo que es una esfera con radio de 200 cm, sin
embargo la simulación se llevo a cabo bajo 36 celdas y superficies, dentro de estas se
ubicaron 10 detectores a lo largo de toda la superficie del contenedor.
1 3 -7.96 -1 4 -2 -17
2 3 -7.96 -2 1 3 4 7 -17
18 2 -1.1 -18 8 9 10 11 22 CELDAS
999 0 998
1 RCC 0 0 -1 0 0 0.06 1.5
2 RCC 0 0 -1 0 0 19.9 1.9 SUPERFICIES
18 TRC 0 0 18.88 0 0 10.6 23.5 9.7
998 SO 100
Cuadro 3.1: Ejemplo input celdas y superficies.
Figura 3.1: Visualización celdas de la geometrı́a del contenedor.
A continuación se presenta un pequeño ejemplo del esquema input o archivo de

entrada con las especificaciones de la geometrı́a de la simulación para este trabajo y
como se pudo ver anteriormente se limita a la producción de superficies que están
delimitadas por un cierto numero de celdas.
La estructura de las celdas están definidas por valores que corresponden al numero
que identifica la celda, el numero que indica el material que llena la celda, si el numero es
0 entonces la celda estará vacı́a, la densidad del material y las superficies que delimitan
la celda, mientras que las superficies se definen utilizando signos positivos y negativos
para indicar su orientación respecto al origen del sistema cartesiano.
Parámetros Comando Ejemplo

Posición POS = X Y Z 000
Celda CEL = N° de Celda inicial 1
Energı́a ERG = Energı́a fuente 14 MeV
Peso WGT = Peso inicial 1
Tiempo TIME = Tiempo inicial 0
Partı́cula PAR = Partı́culas emitidas 1 para N, 2 para P y 3 para E
Cuadro 3.2: Especificación de la fuente.
3.2.3. Tarjeta de Datos

Dentro de la tarjeta de los datos en el archivo de entrada de MCNP, se encuentra
la especificación de la fuente, los materiales, la fı́sica involucrada en la simulación y los
tallys necesarios para el desarrollo del problema.
Especificación de la fuente
La fuente en el código MCNP se especifica mediante el comando SDEF, y permite

proporcionar ciertos parámetros para la fuente a utilizar. Dentro de los parámetros se
encuentran:
A continuación se presenta un ejemplo del archivo de entrada para la especificación

de la fuente
SDEF POS = 0 0 0 CEL = 1 ERG = 14 WGT = 1 TME = 0 PAR = 1
Especificación de los materiales
Los datos en esta sección especifican la composición isotópica de los materiales

en las celdas y las secciones eficaces que se utilizarán en cada librerı́a, que depende
principalmente de las posibles interacciones de la radiación con la materia, el comando
en el archivo de entrada es (m) y puede tener restricciones 1 ≤ m ≤ 99999 numero que
acompañara al comando y que dependerá de la celda utilizada.
los materiales ingresan en el archivo de entrada como un bloque de datos de la forma
ZZZAAAID donde ZZZ representa el numero atómico, AAA representa el numero de
masa atómica mientras que ID hace referencia a la librerı́a de secciones eficaces que se
van a utilizar.
Especificación del modo
MCNP puede transportar la radiación bajo diferentes modos,

Comando Descripción del tally Unidades *Unidades

F1:N,P,E Corriente en una superficie Partı́culas MeV
F2:N,P,E Flujo promedio en una superficie Partı́culas/cm2 MeV/cm2
F4:N,P,E Flujo promedio en una celda Partı́culas/cm2 MeV/cm2
F5a:N,P Flujo en un detector Partı́culas/cm2 MeV/cm2
F6:N,NP,P Energı́a depositada en una celda MeV/g jerks/g
F7:N Energı́a depositada por fisión MeV/g jerks/g
F8:P,E,PE Distribución de energı́a de un pulso pulso MeV
+F8:E Carga depositada carga N/A
Cuadro 3.3: Especificación de los tallys.
MODE
N : Solo transporte de neutrones.
N P : Transporte de neutrones y transporte de fotones inducidos por neutrones.
P : Solo transporte de fotones.
E : Solo transporte de electrones.
P E : Transporte de fotones y transporte de electrones inducidos por fotones.
N P E : Transporte de neutrones, fotones inducidos por neutrones y el transporte

de electrones.
Hemos utilizado durante la simulación el modo N P que como se vio ya transporta

neutrones y fotones que son inducidos por la interacción de los neutrones con la materia.
Especificación de los tallys
Los tallys son utilizados dentro del fichero de entrada para el input en MCNP
ofreciéndole al usuario el tipo de información que desea obtener del calculo MC, es
decir, el flujo en un determinado punto, la energı́a depositada, entre otros, el comando
que se utiliza para el registro de los tallys es (Fn), y n varia dependiendo la información
que se desea conocer.
Existen 7 tipos de tallys disponibles para neutrones, 6 tipos de tallys para fotones
y cuatro tipos de tallys para electrones disponibles para el código MCNP. Dentro de
los tallys básicos para los cálculos MC en MCNP se muestran a continuación.
Para el desarrollo de este trabajo se utilizaron 3 tipos de tallys: F4, F5 y F6, que
nos permitieron obtener los valores de flujo de neutrones y fotones en cada uno de
los puntos de interés, que son detectores ubicados a lo largo de toda la superficie del
contenedor, ası́ como los datos de energı́a depositada en dichos puntos junto con los
valores de tasa de dosis que serán tratados en el capitulo de resultados.
Transporte de neutrones en MCNP
Cuando una partı́cula (que representa cualquier número de neutrones, dependiendo

del peso de la partı́cula) choca con un núcleo, tiene lugar la siguiente secuencia:
1. Se identifica el nucleı́do de colisión:

Si hay n nucleı́dos diferentes que forman el material en el cual ocurrió la colisión, y
si ξ es un número aleatorio en el intervalo unitario (0,1), de esta forma el k-esimo
nucleı́do es elegido a colisionar si
k−1
X n
X k
X
Σti < ξ Σti ≤ Σti (3.3)
i=1 i=1 i=1
Donde Σti es la sección eficaz macroscópica total del nucleı́do i-esimo.
2. Generación opcional de fotones:

Los fotones se generan mediante la interacción de neutrones / fotones y si el nucleı́do
de colisión tiene una sección eficaz de producción de fotones distinta de cero. El
número de fotones producidos es una función del peso de los neutrones, el peso de
la fuente de neutrones, los lı́mites de peso de fotones, la sección eficaz de producción
de fotones, la sección eficaz total de neutrones y la importancia de las celdas, la
creación de fotones bajo la importancia de su peso es,
Wn σγ
Wp = (3.4)
σT
Donde Wp es el peso de los fotones, Wn el peso de los neutrones, σγ es la sección

eficaz en la producción de fotones y σT es la sección eficaz total.
3. Absorción:
La absorción se trata de una de dos maneras: analógica o implı́cita, los términos
absorción y captura se usan indistintamente para los nucleı́dos no fisionables,
Absorción analógica: En este proceso las partı́culas son capturadas por la pro-
babilidad σa /σT , donde σa es la sección eficaz de absorción y σT la sección eficaz
total, todas las partı́culas que han sufrido una absorción analógica, depositan toda
su energı́a y su peso en la celda donde ocurre la colisión.
Absorción implı́cita: Durante esta absorción los neutrones que sobreviven cam-
bian su peso Wn por un nuevo peso Wn , de la forma

σa
Wn = 1− ∗ Wn (3.5)
σT
4. Dispersión elástica e inelástica:

Si las condiciones no dan para un tratamiento térmico, las partı́culas experimentaran
dispersiones elásticas e inelásticas, la probabilidad de que ocurra una dispersión
elástica esta dada por,
σel σel
= (3.6)
σin + σel σT − σa
Donde σel es la sección eficaz de dispersión elástica, σin es la sección eficaz de

dispersión inelástica, σa corresponde a la sección eficaz de absorción y σT la sección
eficaz total, σT = σel + σin + σa .
Si la sección eficaz elastica y total se ajustan al tratamiento termico del gas libre
con energias termicas entonces la selección de una colisión inelástica se hace con la
probabilidad restante
σin
Pin = (3.7)
σT − σa
5. Tratamiento térmico:
El tratamiento térmico denotado como S(α, β) esta completamente representado
por la dispersión inelástica térmica de moléculas y cristales solidos, por medio de
dos procesos, el primero mediante secciones eficaces de dispersión inelastica y un
acople entre energia y angulo en representación de la ley de dispersión térmica. El
segundo por medio de dispersión elástica sin cambios en la energı́a de neutrones de
salida para sólidos con sección eficaz σel Y un tratamiento angular derivado de los
parámetros de la red. El tratamiento de dispersión elástica se elige con probabilidad,
σel
Pel = (3.8)
σel + σin
Este tipo de tratamiento por dispersión térmica permite considerar las dispersiones
para moléculas multiatómicas[25].
Capı́tulo 4
Reglamento Para El Transporte

Seguro de Material Radiactivo
La Organización Internacional de Energı́a Atómica está autorizada por sus estatutos

a “establecer o adoptar, normas de seguridad para proteger la salud y reducir al mı́nimo
el peligro para la vida y la propiedad”; Mediante sus disposiciones de reglamentación
de la seguridad nuclear y radiológica. La seguridad no es un fin en sı́ misma, sino un
requisito indispensable en la actualidad y en el futuro. Los riesgos relacionados con
la radiación ionizante deben evaluarse y controlarse sin restringir indebidamente la
contribución de la energı́a nuclear al desarrollo equitativo y sostenible.
4.1. Conceptos Básicos

Embalaje: Comprende de uno o mas recipientes que cumplen la función de con-
tención y otras funciones de seguridad.
Bulto: Por bulto se entenderá el producto completo de la operación de embalaje,

que consiste en el embalaje y su contenido preparados para el transporte.
Contenedor: Por contenedor se entenderá un elemento de equipo de transporte

de carácter permanente y por tanto, suficientemente resistente para permitir su
empleo repetido; especialmente destinado a facilitar el transporte de mercancı́as,
por una u otra modalidad de transporte.
A1 y A2 : Por A1 se entenderá el valor de la actividad de los materiales radiactivos

en forma especial y que se utiliza para determinar los lı́mites de actividad para los
requisitos del Reglamento de la OIEA, por A2 se entenderá el valor de la actividad
de los materiales radiactivos, que no sean materiales radiactivos en forma especial
37
Capı́tulo 4. Reglamento Para El Transporte Seguro de Material Radiactivo 38
y que se utiliza para determinar los lı́mites de actividad para los requisitos del
Reglamento de la OIEA.
Uso exclusivo: Por uso exclusivo se entenderá el empleo exclusivo, por un solo
remitente, de un medio de transporte o de un gran contenedor, respecto del cual
todas las operaciones iniciales, intermedias, finales de carga, descarga y expedi-
ción sean efectuadas de conformidad con las instrucciones del remitente o del
destinatario.
Índice de transporte: Es un número asignado a un bulto, sobreenvase o contenedor

sin embalar, que se utiliza para controlar la exposición a las radiaciones.
Sistema de contención: Es el conjunto de componentes del embalaje, que están

destinados a contener los materiales radiactivos durante el transporte.
Arreglos especiales: Por arreglos especiales se entenderá aquellas disposiciones,

aprobadas por la autoridad competente, en virtud de las cuales podrá ser trans-
portada una remesa que no satisfaga todos los requisitos aplicables del Reglamen-
to de la OIEA.
4.2. Limites de Actividad y Clasificación

Cada material radiactivo deberá tener un numero especificado por las Naciones
Unidas y que dependerá del tipo de bulto al que corresponda, en este trabajo el tipo
de bulto utilizado es el Tipo A, en la tabla 4.1 se muestran los números otorgados y la
descripción correspondiente.
Numero N.U Descripción

UN 2915 No en forma especial, no fisionables o fisionables exceptuados(b)
UN 3327 No en forma especial, fisionables
UN 3322 En forma especial, no fisionables o fisionables exceptuados(b)
UN 3333 En forma especial, fisionables
Cuadro 4.1: Clasificación del numero de las Naciones Unidas para los bultos Tipo A
Para cada clasificacion y limitacion que el reglamentación de la OIEA expide se

utilizan los valores de A1 y A2 que dependen del radionucleido que contenga el material
radiactivo que en este caso es Americio - 241.
Radionucleido A1 (TBq) A2 (TBq)
Am-241 1∗ 101 1 ∗ 10−3
4.2.1. Clasificación de los Bultos y Embalajes

La actual versión del reglamento de la OIEA establece 5 tipos de embalajes o bultos,
Exceptuado
Industrial
Tipo A
Tipo B
Tipo C
Esta clasificación está relacionada con la actividad y la forma fı́sica del material ra-
diactivo contenido en el embalaje. La OIEA también establece estándares de desempeño
como requisitos de diseño y procedimientos de prueba, para cada tipo de embalaje; per-
mitiendo establecer posibles riesgos potenciales derivados del transporte del material
radiactivo.
Embalajes o Bultos exceptuados
Los embalajes exceptuados son aquellos con niveles de radiación suficientemente

bajos que no simbolizan peligros potenciales y por lo tanto no demandan pruebas en
la integridad hermética. Muchos de estos son bultos vacı́os que contuvieron materiales
radiactivos durante determinados periodos de tiempo y que pueden hacer uso de ins-
trumentos o artı́culos manufacturados con uranio natural, uranio empobrecido o torio
natural.
Embalajes o Bultos Industrial
Los embalajes industriales son usados para transportar dos tipos de materiales:
Material de Baja Actividad Especı́fica (BAE), un ejemplo de embalaje BAE son

los desechos hospitalarios de productos radiofarmacéuticos.
Objetos Contaminados Superficialmente (OCS), las maquinarias del Ciclo de

Combustible o las partes de los reactores nucleares, cuyas superficies han sido
contaminadas por fluidos refrigerantes o agua procesada, son consideradas OCS.
Estos tipos de materiales son intrı́nsecamente seguros debido a la poca dispersión del
material radiactivo almacenado y su baja actividad.
Los Embalajes Industriales (IP) se subdividen en tres categorı́as designadas co-
mo IP-1, IP-2 y IP-3, las cuales difieren de acuerdo al grado de rutina y condiciones
normales del transporte que ellas puedan resistir.
Embalajes o Bultos Tipo A
Los embalajes Tipo A son usados para el transporte de pequeñas, pero significativas,
cantidades de material radiactivo. Ya que se supone que este tipo de embalaje puede
ser, en teorı́a, dañado en un grave accidente y que una porción de su contenido se
puede salir, la cantidad de radionucleidos que ellos pueden contener está limitada por el
Reglamento de la OIEA. En tal caso, estos lı́mites aseguran que los riesgos de radiación
o contaminación externa sean muy bajos.
Embalajes o Bultos Tipo B
Los embalajes Tipo B se usan para el transporte de material altamente radiactivo.

Estos embalajes deben resistir las mismas condiciones normales de transporte como los
embalajes de Tipo A, pero como su contenido excede los lı́mites del Tipo A, es necesario
especificar la resistencia adicional a la radiación o al material radiactivo en caso de
daños accidentales. Este tipo de embalaje debe ser capaz de resistir las condiciones
de accidente esperadas, sin que se agriete el contenedor o haya un incremento en la
radiación a niveles que pongan en peligro al público en general y a aquellos involucrados
en las labores de rescate o de limpieza.
Embalajes o Bultos Tipo C
La Edición de 1996 del Reglamento de Transporte de la OIEA, incorporó un requi-

sito para un embalaje de diseño más robusto, el Embalaje Tipo C, para el transporte
aéreo de material de mayor radiactividad por vı́a aérea. Los embalajes Tipo C deben
satisfacer todos los requisitos adicionales de los embalajes Tipo A y la mayorı́a de los
requisitos adicionales de los embalajes Tipo B [26].
4.3. Requisitos y controles para el transporte

La OIEA establece estándares de desempeño como requisitos de diseño y procedi-
mientos de prueba, para cada tipo de bulto, este enfoque clasificatorio del bulto y por
consiguiente el embalaje está relacionado con el peligro potencial que pudiera surtir de
las operación del transporte de este tipo de material radiactivo.
De acuerdo al reglamento para el transporte seguro de materiales radiactivos nues-
tro contenedor esta sujeto a las disposiciones de un bulto o embalaje tipo A.
Aunque se da por supuesto que los bultos de este tipo no permanecerán intactos
en caso de accidente, han de ser capaces de resistir, sin ninguna pérdida o dispersión
del contenido radiactivo y sin ningún aumento del nivel máximo de radiación en la
superficie externa del embalaje, lo que el reglamento denomina condiciones normales

de transporte[27].
Determinación del ı́ndice de transporte
El IT de un bulto, sobreenvase o contenedor será la cifra obtenida aplicando el

siguiente procedimiento:
Se determinará el nivel de radiación máximo en unidades de milisievert por hora

(mSv/h) a una distancia de 1 m de las superficies externas del bulto, sobreenvase
o contenedor.
El valor determinado se multiplicará por 100 y la cifra obtenida será el IT.
Siguiendo el procedimiento anterior para la fuente de 241 Am −9 Be utilizada en este

trabajo,
El nivel de radiación de la fuente sobre la superficie del contenedor es de 2 mSv/h,
utilizando la ecuación 4.1 para hallar el nivel de radiación a un metro de distancia
X1 D1 = X2 D2 (4.1)
X1 D1
D2 =
X2
Remplazando por los valores iniciales
(0,01m)(2mSv/h)
D2 =
1m
D2 = 2mSv/h
Multiplicando el valor del nivel de radiación a un metro por 100 obtenemos un IT

de 2
Categorı́as de los bultos o empaques
Los bultos se clasificarán en las categorı́as I-BLANCA, II-AMARILLA o III-AMARILLA,

de conformidad con las condiciones especificadas en la tabla 4.2.
Para la fuente de 241 Am −9 Be que utilizamos en el trabajo el valor del ı́ndice de
transporte es 2 y el nivel de radiación en superficie es 2 mSv/h podemos clasificar
nuestro contenedor mediante la tabla 4.2 dentro de la categorı́a III-AMARILLA.
Índice Trans. Nivel de Radiación máximo en Categorı́a

superficie
0 Hasta 0,005 mSv/h I-BLANCA
Mayor que 0 pero no Mayor que 0,005 mSv/h pero no mayor II-
mayor que 1 que 0,5 mSv/h AMARILLA
Mayor que 1 pero no Mayor que 0,5 mSv/h pero no mayor III-
mayor que 10 que 2 mSv/h AMARILLA
Mayor que 10 Mayor que 2 mSv/h pero no mayor que III-
10 mSv/h AMARILLA(b)
Cuadro 4.2: Categorı́as de los bultos y contenedores.
4.4. Marcado, rotulado y etiquetado

Marcado
Todo bulto deberá llevar la marca de las Naciones Unidas, el tipo de bulto al que
pertenezcan (un bulto del Tipo A llevará marcada la inscripción “TIPO A”) y para los
bultos cuya masa bruta exceda de 50 kg llevará marcada su masa bruta permitida de
manera legible y duradera en el exterior del embalaje.
Etiquetado
Todo bulto, sobreenvase y contenedor deberá llevar las etiquetas que correspon-
dan a los modelos aplicables figuras 4.1, 4.2 y 4.3 , ası́ mismo todas las etiquetas no
relacionadas con el contenido deberán retirarse o cubrirse.
Figura 4.1: Etiquetado Categorı́a I - BLANCA [5].
La figura 4.1 muestra la etiqueta para la categorı́a I-BLANCA. El color de fondo

de la etiqueta será blanco, el trébol y los caracteres y lı́neas impresos serán negros y la
barra que indica la categorı́a será roja.
Figura 4.2: Etiquetado Categorı́a II - AMARILLA [5].
La figura 4.2 muestra la etiqueta para la categorı́a II-AMARILLA. El color de

fondo de la mitad superior de la etiqueta será amarillo y el de la mitad inferior blanco,
el trébol y los caracteres y lı́neas impresos serán negros y las barras que indican la
categorı́a serán rojas.
Figura 4.3: Etiquetado Categorı́a III - AMARILLA [5].
La figura 4.3 muestra la etiqueta para la categorı́a III-AMARILLA. El color de

fondo de la mitad superior de la etiqueta será amarillo y el de la mitad inferior blanco,
el trébol y los caracteres y lı́neas impresos serán negros y las barras que indican la
categorı́a serán rojas [5].
Capı́tulo 5
Resultados
Para realizar este trabajo la fuente de 241 Am −9 Be se simulo mediante su espectro

de emisión caracterı́stico reportado por IAEA 318 [28], en la Figura 5.1 se muestra el
espectro usado en el cálculo,
Figura 5.1: Espectro de emisión del 241 Am −9 Be [28].
Una vez caracterizada la fuente, se procedió a elaborar la simulación recreando la

geometrı́a construida bajo los parámetros de un contenedor real como se muestra en la
fı́gura 5.2, la simulación se llevo acabo empleando el código MCNPX.
45
Capı́tulo 5. Resultados 46
Figura 5.2: Contenedor real de transporte para una fuente de 241 Am −9 Be.
El contenedor tiene tres secciones diferentes, la parte central consta de un cilindro

que contiene dentro una pequeña cámara de aire donde se alberga la fuente de 241 Am−9
Be, en las extremidades del contenedor se ubican dos conos truncados que se unen
mediante un cilindro. Compactando el contenedor se ubica una capa de acero inoxidable
sobre toda la superficie como se puede observar en la fı́gura 5.3
Figura 5.3: Geometrı́a del contenedor bajo simulación

El volumen es de unos 60 litros, el contenedor debe ser lo suficientemente pequeño

como para estar dentro del rango portátil y su peso en promedio esta entre los 60 -
65 kilos, estableciéndose como una de las restricciones fı́sicas a tener en cuenta en la
selección de los materiales del blindaje.
La elección de los materiales utilizados durante la simulación se dio bajo ciertas
consideraciones, la primera de ellas involucro un análisis de los elementos con altos coe-
ficientes de absorción para rayos gamma µ(E, Z) producidos por la fuente y aquellos
fotones inducidos por la captura neutrónica, al igual que la sección eficaz de neutro-
nes σ(E, A). Los materiales con mayor conveniencia en estos aspectos deberı́an ser la
concentración de diferentes isótopos de elementos con altos coeficientes de atenuación
para neutrones térmicos y epitérmicos, ya que la fuente emite neutrones dentro de un
rango de 4.4 MeV en un 37.5 %.
Basados en las consideraciones fı́sicas antes mencionadas como el decremento en la
energı́a y el numero de colisiones necesarias para la moderación y posterior atenuación
de los neutrones, los mejores elementos son el hidrógeno, el agua pesada, el boro y
el grafito. Estos elementos tienen una sección eficaz de dispersión relativamente alta
y se pueden considerar como componentes de un blindaje primario. Por otro lado los
materiales con altos valores en su coeficiente de atenuación lineal para la atenuación
gamma son aquellos con un numero alto en su Z.
En la literatura es recurrente encontrar el uso de tres polı́meros (polietileno, poli-
estireno, y metacrilato de polimetilo) para blindar las fuentes de neutrones, todos ellos
basados principalmente en átomos de C y H. Por ello el polietileno es empleado dentro
de los materiales de la simulación y se utiliza como tapón para sellar la puerta de acceso
a la fuente de 241 Am −9 Be, la resina de WEP cuya composición se puede observar en
la tabla 5.1 cumple con las consideraciones antes hechas y su densidad de 1.1 g/cm3
permitió cumplir con las limitaciones del peso del contenedor.
Finalmente los materiales seleccionados como blindaje para la fuente de 241 Am −9
Be son (Aire, Water Extended Polyester (WEP), Acero inoxidable y Polietileno) su
distribución a lo largo del contenedor se puede observar en la figura 5.3. La composición
de cada uno estos materiales y sus densidades se muestran en la tabla 5.1 .
Material Densidad Composición

WEP 1.1 g/cm3 H(9.7 %) + C(25,3 %) + O(65 %)
Acero Inoxidable 7.96 g/cm3 Cr(18 %) + Fe(70 %) + Ni(10 %) + Mo(2 %)
Aire 0.00129 g/cm3 O(24.4 %) + C(2.5 %) + N(75 %) + Ar(0.5 %)
Polietileno 0.94 g/cm3 H(14.37 %) + C(85.63 %)
Cuadro 5.1: Especificaciones de los materiales

Durante la simulación se ubicaron 10 detectores localizados en los puntos de interés,

es decir 8 detectores Monte Carlo (MC) sobre la superficie del contenedor y los 2
detectores MC restantes a un metro de la superficie del blindaje, la ubicación de los
detectores se puede observar en la Figura 5.4.
Figura 5.4: Esquema de la simulación
En la figura 5.4 los números indican las posiciones de los detectores MC, los puntos
1-8 se encuentran sobre la superficie del contenedor permitiéndonos obtener los valores
de tasa de dosis equivalente que son limitadas por la normativa internacional para el
transporte de este material radiactivo (< 2000µSv/h), mientras que los puntos 9-10 se
ubican a un metro de la superficie externa del contenedor, donde se hallaran valores
de tasa de dosis permisibles para el publico, ya que se encuentran a distancias medias
donde podrı́a presentarse publico cuando ocurre el transporte.
La simulación se llevo acabo con un total de 450 millones de historias y un tiempo
de calculo aproximado de 45 dı́as, los datos de salida del calculo corresponden a una
fuente de 241 Am −9 Be de 18 Ci, el programa permitió conocer la fluencia de fotones
y neutrones en cada uno de los detectores, también obtuvimos el valor de tasa de
dosis equivalente H ∗ permitiéndonos evaluar el rendimiento del blindaje radiobiologico
debido a fotones y neutrones como se puede observar en la tabla 5.2, todos los resultados
obtenidos del calculo son producidos por partı́cula emitida desde la fuente y presentan
incertidumbres menores al 1 %.
Como resultado de dichos calculos se muestran las fluencia de neutrones y fotones

en el punto 1 y 9 como se puede observar en las Figuras 5.5 y 5.6, el primer punto se
ubica sobre la puerta del contenedor y el segundo punto a un metro de distancia del
primero.
Figura 5.5: Fluencia de neutrones en los puntos 1 y 9
En la Figura 5.5 el espectro superior esta relacionado con la posición 1 del detec-
tor MC, mientras que la curva inferior es la obtenida para la posición del detector 9.
Comparando las fluencias de neutrones de la gráfica anterior se observa que hay una dis-
minución de la fluencia en un factor de 2 sin embargo para bajas energı́as la atenuación
es mayor, para interpretar lo que se observa en la gráfica anterior se debe recordar que
el punto 9 se encuentra a un metro de distancia y el detector esta ubicado en medio de
aire, por lo tanto usando las gráficas 2.4 y 2.6 reportadas por el centro estadounidense
de datos nucleares (LANL), se observa la sección eficaz para el oxigeno y el nitrógeno,
en el oxigeno el factor dominante es la sección eficaz total de captura neutrónica que
alrededor de energı́as cercanas a 1 MeV es de 1 bar, por otro lado para el caso del
nitrógeno, podemos observar que la captura de neutrones en la región térmica tiene
valores ligeramente mayores a la captura total a la energı́a de 1 MeV, recordemos que
en la composición del aire el nitrógeno es 4 veces mas abundante que el oxigeno, es por
esto que se observan las diferencias entre los espectros de neutrones para el punto 1 y
9. La figura también permite apreciar neutrones dentro de todo el rango de energı́as.
Figura 5.6: Fluencia de fotones en los puntos 1 y 9.
Para la Figura 5.6 se observa que la forma espectral no presenta grandes diferencias
entre si, interpretamos este resultado debido a que tanto el material blindante como
el aire tienen una eficacia de absorción similar y que la diferencia entre el área de los
espectros esta relacionada con la ley del inverso al cuadrado de la distancia.
Las Figuras 5.7 y 5.8 muestran las fluencias de neutrones y fotones, en los dos casos
las graficas reflejan la interaccion de los neutrones con el material blindante, se observa
la representación de los detectores ubicados en los puntos 1, 2 y 3.
Las diferencias entre los espectros son consistentes, el espectro de neutrones permite
observar que se encuentran dos regiones que muestran una diferencia, la primera se
ubica aproximadamente en 1x10− 7 MeV y el otro esta sobre los 2.5 MeV donde el flujo
de neutrones es mayor sobre el punto 1, que es la puerta de acceso al contenedor, esto
quiere decir que los puntos 2 y 3 presentan mayor atenuación neutrónica en la zona de
altas energı́as, mientras que para el punto 1 se evidencias fugaz de neutrones rápidos.
El espectro de fotones tiene dos regiones que presentan una notable diferencia al
igual que el espectro de neutrones, para este caso el primer pico se ubica sobre los 0.2
MeV y los 2.5 MeV, el detector 2 evidencia lo mencionado en el espectro de neutrones
donde esta zona ofrece mayor atenuación de neutrones, al encontrarse mayor espesor
Figura 5.7: Fluencia de neutrones sobre los detectores 1, 2 y 3.
Figura 5.8: Fluencia de Fotones sobre los detectores 1, 2 y 3.

Figura 5.9: Fluencia de neutrones en los puntos 2 y 5.
de material para este punto y por consiguiente genera fotones inducidos de energı́as de
2,5 MeV de la captura neutrónica.
En las Figuras 5.9 y 5.10 se presentan los espectro de neutrones y fotones de los
puntos 2 y 5 que representan los valores maximos y minimos obtenidos durante el calculo
de las fluencias de los detectores MC ubicados sobre la superficie del contenedor. La
Figura 5.9 permite observar un espectro superior que esta relacionado con el punto de
ubicación 5, representa el máximo espectro obtenido durante el calculo Monte Carlo
sobre los detectores de la superficie, el espectro inferior se obtiene como resultado
del detector ubicado en la posición 2 que representa el valor con menor fluencia de
neutrones, los espectros de los demás detectores se ubican en la zona entre los dos
espectros con ligeras diferencias entre si.
La diferencia entre el área de los espectros se relaciona principalmente por la distancia
en la que se encuentran ubicados los detectores respecto a la posición de la fuente,
el punto 2 esta localizado a 23.8 cm, mientras que el punto 5 esta sobre los 18.8 cm.
Como es de esperarse a medida que el material blindante posea un mayor espesor la
atenuación frente a los neutrones sera mayor.
Al igual que en la Figura 5.9 se observa un espectro superior que esta relacionado con
el punto de ubicación 2, representa el máximo espectro de fotones obtenido durante el
calculo , el espectro inferior se obtiene como resultado del detector ubicado en la posición
Figura 5.10: Fluencia de fotones en los puntos 2 y 5.
5 que representa el valor con menor fluencia, los espectros no presentan cambios en su
forma y su diferencia de área se ve relacionada al espesor de material utilizado. Los
espectros de los demás detectores se ubican en la zona entre los dos espectros con ligeras
diferencias entre si.
Para evaluar la eficacia del contenedor como blindaje radiobiologico de transporte,

se calculo la tasa de dosis al rededor de la superficie del contenedor, recordando que los
valores máximos permisibles son (< 2000µSv/h), estos valores se obtuvieron a través
de los detectores MC que nos permitieron obtener del código el tally F5 (Fluencia en
un detector) valores que fueron operados por unos factores de calidad.
Los valores obtenidos para la tasa de dosis de cada uno de los detectores ubicados
en la simulación se presentan en el cuadro 5.2, en la primera columna de la tabla se
encuentran los detectores según su ubicación, en la segunda columna los valores de tasa
de dosis equivalente debido a neutrones, en la tercera columna el valor de tasa de dosis
debido a fotones por captura, en la cuarta columna el valor de tasa de dosis total que
se obtiene de la suma de las tasa de dosis de neutrones y fotones, en la ultima columna
se da la distancia entre los detectores y la fuente.
De- H ∗ 10 (Neutrones H ∗ 10 (Fotones H ∗ 10 (Total Distancia

tector µSv/h) µSv/h) µSv/h) (cm)
1 1739 81 1820 24.3
2 681 97 778 23.9
3 1062 119 1181 21.7
4 1177 120 1297 19.8
5 1833 150 1983 18.8
6 1178 120 1298 19.8
7 1062 119 1181 21.7
8 681 97 778 23.9
9 43 3 46 124.3
10 32 4 36 121.7
Cuadro 5.2: Valores de tasa de dosis equivalente para los detectores 1 - 10
Del Cuadro 5.2 se puede observar varias aspectos, el primero de estos es que el de-
tector 5 nos muestra que son requeridos al menos 19 cm de material WEP para obtener
tasas de dosis dentro de los limites estipulados por la normativa para el transporte.
La segunda que el blindaje utilizado en el contenedor de transporte cumple con los
requerimientos de la normativa internacional para el transporte de este material ra-
diactivo, los valores de tasa de dosis de cada uno de los detectores se encuentran dentro
del limite (< 2000µSv/h), demostrando que la elección de los materiales utilizados en
la simulación tabla 5.1, cumplen de manera adecuada con cada uno de los paráme-
tros establecidos para el blindaje radiobiologico, sin embargo estos valores deberı́an ser
optimizados a tasas de dosis aproximadas a 1000µSv/h, esto se puede lograr como el
cuadro nos permite observar agregando 3 cm al espesor del material blindante WEP es
decir con 22 cm, como lo refleja los cálculos de los detectores MC 3 y 7.
En la puerta de acceso a la fuente de 241 Am −9 Be se observa un valor de tasa
de dosis total de 1820 µSv/h y hay un factor de 20 entre la tasa de dosis equivalente
por neutrones y captura gamma, por lo tanto se puede decir que la composición de
los materiales seleccionados es altamente eficiente para este tipo de fuentes de emisión
binaria, sin embargo cabe resaltar que los valores presentados reflejan una fuga de
radiación que debe ser reevaluada con el fin de optimizar los valores de tasa de dosis
que aunque se encuentran dentro de los limites establecidos deben garantizar tanta
protección radiológica como sea posible.
En particular vale la pena mencionar que para el detector 9 se observa que la relación
entre la tasa de dosis de captura gamma y la debida a neutrones se redujo a un factor
de 15, esta diferencia entre los detectores 1 y 9 se presenta debido a que el detector 9
se encuentra en aire y el análisis es igual al realizado en la Figura 5.5.
Los detectores MC 9 y 10 muestran valores en tasa de dosis de 46 y 36 µSv/h

respectivamente, estas diferencias demuestran mayores virtudes como blindaje para
neutrones al material WEP frente al Polietileno, ademas estos valores tomados a una
distancia promedio de 1.2 metros permiten obtener valores en tasa de dosis para publico
menores a 5 mSv/h cuando el contenedor se encuentre dentro del vehiculo de transporte.
Figura 5.11: Esquema de la simulación
Los resultados del Cuadro 5.2 pueden ser plasmados mediante un mapa de dosis,
como se observa en la Figura 5.11, donde los valores de tasa de dosis mas altos se
encuentran en las tonalidades rojas mas intensas y se establecen cerca de la fuente de
241 Am −9 Be alcanzando distancias de 5 cm sobre las zonas de material WEP mien-
tras que en la zona del polietileno alcanza distancias de 7 cm, los valores intermedios
aproximados a 2000µSv/h se observan en tonalidades verdes alcanzando distancias de
19 cm para material WEP y 23 cm para material Polietileno, los valores con tasas de
dosis optimizadas aproximadamente 1000µSv/h tienen las tonalidades de azul oscuro
con distancias de 22 cm en WEP mientras que el Polietileno no logra estos valores de
tasa de dosis y por eso en su region no se observa tonalidades azul oscuro, se manifiesta
una deformidad en el campo de radiacion que puede apreciarse sobre la zona del tapon
de polietileno es decir la puerta de acceso al contenedor.
Capı́tulo 6
Conclusiones
De los resultados obtenidos por el detector MC 5 podemos decir que son requeridos
al menos 19 cm de material WEP para obtener tasas de dosis dentro de los limites
estipulados por la normativa internacional para el transporte seguro de materiales
radiactivos.
De los valores del detector MC 1 se observa que para el material Polietileno son
necesarios 24 cm de espesor de material para obtener tasas de dosis dentro de los
limites estipulados por la normativa internacional para el transporte.
El blindaje utilizado en el contenedor de transporte cumple con los requerimientos

de la normativa internacional para el transporte de este material radiactivo, todos
los valores de tasa de dosis de los detectores se encuentran dentro del limite
( < 2000µSv/h), demostrando que la elección de los materiales utilizados en la
simulación tabla 5.1, cumplen de manera adecuada con cada uno de los parámetros
establecidos para el blindaje radiobiologico.
Aunque los valores de tasa de dosis se encuentren dentro de los limites estos valores
deberı́an ser optimizados a tasas de dosis aproximadas a 1000µSv/h garantizando
tanta radioprotección como sea posible, esto se puede lograr como el cuadro 5.2
nos permite observar agregando 3 cm al espesor del material blindante WEP es
decir con 22 cm de espesor de material.
En la puerta de acceso a la fuente de 241 Am −9 Be detector MC 1 se observa un

valor de tasa de dosis total de 1820 µSv/h y hay un factor de 20 entre la tasa
de dosis equivalente por neutrones y captura gamma, por lo tanto se puede decir
que la composición de los materiales seleccionados es altamente eficiente para
este tipo de fuentes de emisión binaria, sin embargo cabe resaltar que los valores
presentados reflejan una fuga de radiación que debe ser reevaluada con el fin de
optimizar los valores de tasa de dosis.
56
Bibliografı́a
[1] Shahriari M. Sohrabpour M. 2000. Borehole parametric study for neutron induced
capture γ-ray spectrometry using the MCNP code, Applied Radiation and Isotopes
52(1):127–135.
[2] Mowlavi A. Fayegh R. 2004. Determination of 4.438 MeV γ-ray to neutron emis-
sion ratio from a 241 Am–9 Be neutron source, Applied Radiation and Isotopes
60(6):959–962.
[3] Croft S. 1989. The use of neutron intensity calibrated 9Be (α, n) sources as 4438
keV gamma-ray reference standards, Nuclear Instruments and Methods in Phy-
sics Research Section A: Accelerators, Spectrometers, Detectors and Associated
Equipment 281(1):103–116.
[4] Nasrabadi M. Baghban G. 2013. Neutron shielding design for 241 Am–Be neutron
source considering different sites to achieve maximum thermal and fast neutron
flux using MCNPX code, Annals of Nuclear Energy 59:47–52.
[5] International Atomic Energy Agency. 2010. Regulations for the safe transport of
radioactive material, International Atomic Energy Agency.
[6] Korf C. 1972 water extended polyester resin. August 29. US Patent 3,687,883.
[7] Kaplan I. 1961. F’ısica nuclear, Madrid: Aguilar:172,508–514.
[8] Guzman M. 1989. Basic nucleonics, pages 170,216.
[9] Turner J. 2008. Atoms, radiation, and radiation protection, John Wiley y
Sons.:209–230.
[10] Atomic Energy of canada. Neutron sources and their characteristics, Otawa:AECL,
Technical Bulletin NS-1.
[11] Krane K. Halliday D. 1988. Introductory nuclear physics, 465:446–447.
[12] Koerting V. Queral C. 2012. Energia nuclear y ciclo de combustible.
57
Bibliografı́a 58
[13] White HE. WHITE H. 1969. Introducción a la fı́sica atómica y nuclear, 1:225–231.
[14] Espinosa G. Rodrı́guez R. 2016. fisica de reactores nucleares un enfoque conceptual.
[15] Marino R. 2004. Radiological shielding of low power compact reactor: calculation
and design.
[16] Schaeffer N. 1973. Reactor shielding for nuclear engineers, page 453, Radiation
Research Associates, Inc., Fort Worth, Tex.(USA).
[17] Vega Carrillo H. Manzanares E. Hernandez V. Gallego E. Lorentec A. 2006. Neu-

tron shielding performance of water-extended polyester.
[18] S. Ovalle. Simulación monte carlo de la interacción de neutrones con la materia.

Tesis de maestria., 2007.
[19] Murray L. Ellis R. 1957. Nuclear reactor physics, Physics Today 10:36.
[20] Glasstone S. Sesonske A. 1990. Ingenierı́a de reactores nucleares.
[21] Semat H. 1971. Fı́sica atomica y nuclear:(curso de introducción), page 475.
[22] Metropolis N. Ulam S. 1983. The monte carlo method.
[23] Sobol I Vega C. 1983. Método de Montecarlo.
[24] Pelowitz D. 2005. MCNPX user’s manual version 2.5. 0., Los Alamos National
Laboratory, Los Alamos.
[25] Brown F. et al. 2003. MCNP–A General Monte Carlo N-Particle Transport Code,
Version 5, Los Alamos National Laboratory, Oak Ridge, TN.
[26] World Nuclear Transport Institute. Tipos de Embalajes usados para el Transporte
de Materiales Radiactivos, Ficha tecnica N°2.
[27] White M. El transporte de materiales radiactivos en bultos del tipo A, exentos y

de baja actividad especı́fica.
[28] International atomic energy agency. 2001. Compendium of neutron spectra and
detector responses for radiation protection purposes, IAEA 318.

Blindajes Neutrones

Cargado por

Copyright:

Formatos disponibles

Blindajes Neutrones

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Blindajes Neutrones

Cargado por

Copyright:

Formatos disponibles

Universidad Pedagógica y Tecnológica de Colombia

CALCULO MONTE CARLO DE BLINDAJES ESRUCTURALES

por: Alejandra Marı́a Garcia Rodriguez

Director: PhD Segundo Agustı́n Martı́nez Ovalle

3. Método Monte Carlo 28

4. Reglamento Para El Transporte Seguro de Material Radiactivo 37

1.1. Esquema Dispersión Elástica . . . . . . . . . . . . . . . . . . . . . . . . 12

2.1. Esquemas Sistema Laboratorio (SL) y Sistema Centro de masa (SCM). 18

3.1. Visualización celdas de la geometrı́a del contenedor. . . . . . . . . . . . 32

4.1. Etiquetado Categorı́a I - BLANCA [5]. . . . . . . . . . . . . . . . . . . . 42

5.1. Espectro de emisión del 241 Am −9 Be [28]. . . . . . . . . . . . . . . . . . 45

1.1. Fuentes de neutrones(α, n) [9]. . . . . . . . . . . . . . . . . . . . . . . . 9

2.1. Propiedades moderadoras de diferentes materiales [20]. . . . . . . . . . . 22

3.1. Ejemplo input celdas y superficies. . . . . . . . . . . . . . . . . . . . . . 32

5.1. Especificaciones de los materiales . . . . . . . . . . . . . . . . . . . . . . 47

Las fuentes de neutrones de 241 Am −9 Be son ampliamente utilizadas en muchas

Por lo anterior es aconsejable discutir las propiedades relevantes de los materiales

El transporte de materiales radiactivos esta regulado por el Organismo Internacio-

El diseño del blindaje se estudia mediante el transporte de los neutrones y su inter-

Rutherford y sus colegas emplearon partı́culas α de gran energı́a, procedentes de

1.1. Generación de neutrones

Se pueden generar haces de neutrones por medio de diferentes reacciones nucleares,

1.1.1. Reacciones con partı́culas alfa (α, n)

He42 + Be94 → C612 + n10 (1.1)

Fuente Vida Media

Cuadro 1.1: Fuentes de neutrones(α, n) [9].

1.1.2. Reacciones con partı́culas gamma (γ, n)

Este proceso es similar a las reacciones (α, n) ya expuestas en este trabajo, se

La fuente fotoneutrónica mas empleada es el antimonio - berilio (124 Sb −9 Be),

1.1.3. Fisión espontánea

1.1.4. Reacciones nucleares

1.2. Clasificación de los neutrones

Clasificación Energı́a (eV ) Velocidad (m/s)

Cuadro 1.2: Clasificación de neutrones [12].

1.3. Interacción neutrón materia.

Figura 1.1: Esquema Dispersión Elástica

Figura 1.2: Esquema Dispersión Inelástica.

excitación y regresa a su estado base emitiendo un rayo γ, en este caso no se conserva

XZA + n10 → (XZA ) ∗ +n10 (1.3)

En las reacciones de captura radiactiva, o reacciones (n, γ), el núcleo compuesto

Donde XZA es el blanco, de número atómico Z y numero másico A. El producto de

XZA + n10 → XZA+1 (1.5)

El producto resultante de este proceso de absorción es inestable y, de acuerdo con

Captura con desintegración β

Captura con emisión de partı́culas α

Captura con emisión de protones

Entre las medidas defensivas contra la irradiación se encuentran tres elementos

La elección de los materiales para el contenedor de transporte de una fuente de

2.1.3. Water Extended Polyester

Alta resistencia al fuego o temperaturas elevadas.

No debe ser toxico, ni permitir la emisión de gases tóxicos al elevarse la tempe-

El material debe ser durable en el tiempo y permanecer lo mas inoloro posible.

No deben producirse deformaciones hasta alrededor de los 150°C.

El material debe ser fácilmente removible, reemplazado o reparado.

2.2. Moderación de neutrones

El proceso de moderación permite reducir las velocidades de los neutrones mediante

Para la mayorı́a de los fines prácticos el proceso importante en la moderación se da

Resolviendo las ecuaciones establecidas en el principio de la conservación de la