Resolución la ecuación de difusión en tres dimensiones mediante técnicas de multigrilla

Abstract

Se ha desarrollado un programa para la solución de la ecuación de difusión en tres dimensiones. Se usaron métodos iterativos, entre ellos SIPSOL, y se aceleró la convergencia con métodos de multigrilla. Se programó una serie de subrutinas que permiten implementar el algoritmo básico de dos redes, los algoritmos ”V” y ”W” entre múltiples niveles de redes, y el algoritmo de multigrilla completa. El programa trabaja con fuentes y difusividades dependientes de la posición. Se realizaron las pruebas básicas para validar el programa.

A program for the solution of the three dimensional diffusion equation has been developed. Iterative methods are used, among them SIPSOL. Convergence is accelerated with multigrid methods. A number of subrutines were developed, they implement the basic two meshes algorithm, ”V” and ”W” algorithm between various levels of meshes, and the full multigrid algorithm. The program works with localized sources and position dependent, diffusivities. The standard tests to validate the program were satisfactorily performed.