FEMFFUSION: Un código de elementos finitos para la simulación de sistemas nucleares

FEMFFUSION es un código neutrónico escrito en C++ que resuelve la ecuación multigrupo del transporte neutrónico usando la aproximación de la difusión y la aproximación de esféricos armónicos simplificados, también conocida como las ecuaciones SPN. El código utiliza para la discretización espacial de las ecuaciones elementos finitos continuos de Galerkin siendo capaz de resolver todo tipo de geometrías y dimensiones del problema. El código está construido sobre la librería de código abierto Deal.ii que le suministra, además del soporte, avances en el método de elementos finitos. Esta librería incluye muchas técnicas de desarrollo propio para la resolución de ecuaciones en derivadas parciales, problemas de álgebra lineal y estrategias de ciencia computacional, pero además es compatible y posee interfaz con otras librerías de código abierto punteras en estos campos. En nuestro caso, además del Deal.ii, nos basamos de las estructuras de PETSc que es una herramienta con técnicas muy sofisticadas para la resolución en paralelo de problemas de álgebra lineal. Por otra parte, utilizamos técnicas desarrolladas por la librería SLEPc junto con otras metodologías propias para la resolución de problemas de autovalores. Por último, destacar que el código permite el uso de técnicas matrix-free que no necesitan montar ni guardar las matrices involucradas en el sistema, para mantener demandas de la memoria computacional en valores razonables y aumentar la velocidad de cálculo.

Más información: