Integración de la ecuación de la difusión neutrónica en geometrías multidimensionales. Aplicación a reactores nucleares. Cálculo de los Modos Lambda.

Abstract

AgradecimientosEn primer lugar, quiero dar las gracias a los directores del presente trabajo, Gumersindo Verd u y Jos e L u i s M ũ noz-Cobo, por su entusiasmo y la gran ayuda que me han prestado, sin los cuales, no hubiese sido posible la realizaci on del mismo.Tambi en quiero mostrar mi agradecimiento a Vicente Vidal por su colaboraci on y a yuda en la resoluci on de los problemas de c alculo num erico que han ido surgiendo.Debo asimismo agradecer a Ximo, Loles y a Ana su ayuda con las gr a cas y dibujos de la memoria, extendiendo el agradecimiento a todos los amigos que han estado conmigo en el Departamento de Ingenier a Qu mica y Nuclear durante la realizaci on del trabajo, por la ayuda que me han prestado y el excelente ambiente de trabajo del que he disfrutado durante este tiempo.Asimismo, quiero dar las gracias a Antonio Soria por sus comentarios y sugerencias, que han contribuido a mejorar en lo posible la versi on nal de la memoria.Por ultimo, quiero recordar aqu especialmente a mi familia, por su ayuda y apoyo constante.

Estudio de librerías paralelas de libre distribución y algoritmos paralelos iterativos multipaso para la resolución de sistemas de ecuaciones lineales dispersos. Aplicación a la ecuación de difusión neutrónica.

“…Para simplificar (2.11), se utiliza una aproximación multigrupo [47] [123], que consiste en obtener una ecuación para los neutrones cuya energía esté comprendida en cada intervalo [E g , E g+1 ], g = 1, ..., G − 1, donde G es el número de grupos de energía considerados. Esto es posible porque, en general, las secciones eficaces están en función de la energía de los neutrones [47].…”

“…Esto es posible porque, en general, las secciones eficaces están en función de la energía de los neutrones [47].…”

“…Se dice que un reactor está en estado crítico cuando la proporción a la que se producen los neutrones en su interior es igual a la proporción a la que se pierden [47] [98]. En estas condiciones el reactor se encuentra en estado estacionario.…”

“…Así pues, para estudiar el estado estacionario de un reactor dado, se puede forzar la criticidad del mismo de un modo artificial, multiplicando las secciones eficaces, relacionados con la producción de los neutrones por procesos de fisión, por un número λ (o bien dividiéndolas por un número k, λ=1/k) [47]. De esta forma, se espera que exista un número λ que satisfaga las ecuaciones (2.17) y (2.18), mismas que pueden reescribirse de la siguiente manera: …”

“…La ecuación constituye un sistema de ecuaciones en derivadas parciales, por lo que su resolución numérica precisa de métodos de discretización, tanto para la parte espacial como para la parte temporal. Mucha de la información mostrada en este capítulo puede encontrarse en forma más detallada en [81] [47] [117] [116].…”