Descripción e implementación de un precondicionador multinivel algebraico para esquemas de elementos finitos H1-conformes
DOI:
https://doi.org/10.15359/ru.34-2.4Palabras clave:
Métodos de elementos finitos, esquemas H1 conformes, aproximaciones de bajo orden, técnicas multiniveles, implementación computacional, MATLAB®Resumen
En este artículo se presenta, en forma detallada, aspectos sobre la implementación del Método de Elementos Finitos (FEM, por sus siglas en inglés), para resolver una ecuación de Poisson con condiciones de frontera homogéneas. El objetivo de este trabajo es clarificar los detalles de esta implementación, tales como la construcción de los algoritmos, creación de experimentos numéricos y los resultados acerca de estos. Por ello, se describe el problema continuo y se muestra un enfoque FEM clásico para resolverlo. Después, se establece una técnica multiniveles para la resolución eficiente del sistema lineal correspondiente, que describe e incluye algunos diagramas para explicar el proceso y presenta los códigos de la implementación en MATLAB®. Finalmente, se realiza una validación de los códigos con varios experimentos numéricos. Los resultados muestran un comportamiento adecuado del precondicionador debido a que el número de iteraciones del método PCG no se incrementa, incluso cuando el tamaño de la malla se reduce.
Referencias
Beck, R. (1999). Graph-based algebraic multigrid for Lagrange-type finite elements on simplicial meshes. Berlin, Germany: Konrad Zuse Zentrum für Informationstechnik. Recuperado de https://pdfs.semanticscholar.org/c719/beb475ff33cab3f235c361f7e8d76e9f820c.pdf
Beirão da Veiga, L.; Brezzi, F.; Cangiani, A.; Manzini, G.; Marini, L. D. & Russo, A. (2013). Basic principles of virtual element methods. Mathematical Models and Methods in Applied Sciences, 23(01), 199-214. doi: https://doi.org/10.1142/S0218202512500492
Burden, R. L.; Faires, J. D. & Burden, A. M. (2015). Numerical Analysis (10th ed.). United States: Cengage Learning.
Carstensen, C. & Klose, R. (2002). Elastoviscoplastic finite element analysis in 100 lines of Matlab. Journal of Numerical Mathematics, 10(3), 157-192. doi: https://doi.org/10.1515/JNMA.2002.157
Castillo, P. E. & Sequeira, F. A. (2013). Computational aspects of the Local Discontinuous Galerkin method on unstructured grids in three dimensions. Mathematical and Computer Modelling, 57(9-10), 2279-2288. doi: https://doi.org/10.1016/j.mcm.2011.07.032
Chen, K. (2005). Matrix preconditioning techniques and applications. United Kingdom: Cambridge University Press. doi: https://doi.org/10.1017/CBO9780511543258
Ciarlet, P. G. (2002). The finite element method for elliptic problems. United States: SIAM. doi: https://doi.org/10.1137/1.9780898719208
Cockburn, B. & Shu, C. W. (1998). The local discontinuous Galerkin method for time-dependent convection-diffusion systems. SIAM Journal on Numerical Analysis, 35(6), 2440-2463. doi: https://doi.org/10.1137/S0036142997316712
Johnson, C. (2009). Numerical solution of partial differential equations by the finite element method. United States: Courier Corporation.
Ruge, J. W. & Stüben, K. (1987). Algebraic multigrid. In: S.F., McCormick. (Ed.), Multigrid Methods (pp. 73-130). United States: SIAM. doi: https://doi.org/10.1137/1.9781611971057.ch4
Saad, Y. (1996). ILUM: a multi-elimination ILU preconditioner for general sparse matrices. SIAM Journal on Scientific Computing, 17(4), 830-847. doi: https://doi.org/10.1137/0917054
Saad, Y. (2003). Iterative methods for sparse linear systems (2th ed.). United States: SIAM. Recuperado de https://www-users.cs.umn.edu/~saad/IterMethBook_2ndEd.pdf
Shewchuk, J. R. (1996). Triangle: engineering a 2D quality mesh generator and Delaunay triangulator. In Workshop on Applied Computational Geometry (pp. 203-222). Heidelberg, Germany: Springer. doi: https://doi.org/10.1007/BFb0014497
Stüben, K. (2001). A review of algebraic multigrid. In: C., Brezinski, L., Wuytack. (Eds.), Numerical Analysis: historical Developments in the 20th Century (pp. 331-359). Elsevier. doi: https://doi.org/10.1016/B978-0-444-50617-7.50015-X
Trottenberg, U.; Oosterlee, C. W. and Schuller, A. (2000). Multigrid. United States: Academic Press.
Wilkinson, J. H. (1994). Rounding errors in algebraic processes. United States: Courier Corporation.
Descargas
Publicado
Cómo citar
Número
Sección
Licencia
Los autores que publican en esta revista están de acuerdo con los siguientes términos:
1. Los autores garantizan a la revista el derecho de ser la primera publicación del trabajo al igual que licenciado bajo una Creative Commons Attribution License que permite a otros compartir el trabajo con un reconocimiento de la autoría del trabajo y la publicación inicial en esta revista.
2. Los autores pueden establecer por separado acuerdos adicionales para la distribución no exclusiva de la versión de la obra publicada en la revista (por ejemplo, situarlo en un repositorio institucional o publicarlo en un libro), con un reconocimiento de su publicación inicial en esta revista.
3. Los autores han afirmado poseer todos los permisos para usar los recursos que utilizaron en el artículo (imágenes, tablas, entre otros) y asumen la responsabilidad total por daños a terceros.
4. Las opiniones expresadas en el artículo son responsabilidad de los autores y no necesariamente representan la opinión de los editores ni de la Universidad Nacional.
Revista Uniciencia y todas sus producciones se encuentran bajo una Licencia Creative Commons Atribución-NoComercial-SinDerivadas 4.0 Unported.
No existe costo por acceso, revisión de propuestas ni publicación para autores y lectores.