ACELERACIÓN DE ALGORITMOS DE RESOLUCIÓN DE SISTEMAS LINEALES DE MATRICES DISPERSAS PARA CÓDIGOS DE COMBUSTIBLES NUCLEARES

MATIAS LOZA PERALTA; MARTIN LEMES; ALEJANDRO SOBA

Congreso; ENIEF; 2019

La resolución de ecuaciones diferenciales por métodos discretos como el de elementosfinitos requiere la resolución de grandes sistemas de ecuaciones lineales en general dispersas. Tal es elcaso de las obtenidas en el desarrollo de códigos de combustibles nucleares, especialmente cuando setrabaja en tres dimensiones. Para acelerar el cálculo minimizando el tiempo de resolución de dichossistemas es necesario programar solvers de ecuaciones lineales eficientes, basados en métodositerativos entre los que destacan aquellos que utilizan espacios de Krylov, ya que satisfacen lascondiciones de velocidad y eficiencia requeridas. Así mismo y para optimizar aún más estosalgoritmos, es posible incluir herramientas de paralelización, sea mediante métodos que utilizanmemoria distribuida (MPI) o memoria compartida (Openmp). En este trabajo se presentan losresultados obtenidos durante el desarrollo de varios solvers iterativos, el CG (gradientes conjugados),el BiCG (gradientes biconjugados) y el BiCGSTAB (gradientes biconjugados estabilizado). Los tressolvers se utilizaron para resolver problemas originados dentro del Código de combustibles nuclearesDIONISIO, usando distintos tipos de precondicionadores y aplicados a sistemas lineales simétricos yasimétricos, comparándose tiempo y eficiencia de los mismos.