Desarrollo de un código neutrónico de difusión 2D y 3D estacionario por el Método de Volúmenes Finitos

Abstract

El objetivo del trabajo es el desarrollo de un código neutrónico modal de difusión en 2D y 3D estacionario utilizando el Método de Volúmenes Finitos, a partir de códigos libres y que se pueda aplicar a reactores de cualquier geometría. Actualmente, los métodos numéricos que más se utilizan en lo códigos de difusión proporcionan buenos resultados en malla estructurada, pero su aplicación a malla no estructurada no es fácil y puede presentar problemas de convergencia y estabilidad de la solución. Respecto a la malla no estructurada, su uso está justificado por su fácil adaptación a geometrías complejas y por el desarrollo de códigos acoplados termohidráuliconeutrónico, así como el desarrollo de códigos fluidodinámicos (CFD) que incentivan el desarrollo de un código neutrónico que tenga la misma malla que la de los códigos fluidodinámicos, que en general suele ser no estructurada. Por otra parte, el refinamiento de la malla y su adaptación a geometrías complejas es otro estímulo de cara a conocer con más detalle lo que ocurre en el núcleo del reactor. Finalmente, el código se ha validado con la simulación de un reactor homogéneo y otro heterogéneo para 2D y 3D