Summary: | Este trabalho apresenta a aplicação de técnicas de Decomposição de Domínio e Processamento Paralelo na solução de grandes sistemas de equações algébricas lineares provenientes da modelagem de fenômenos eletromagnéticos pelo Método de Elementos Finitos. Foram implementadas as técnicas dos tipos Complemento de Schur e o Método Aditivo de Schwarz, adaptadas para a resolução desses sistemas em cluster de computadores do tipo Beowulf e com troca de mensagens através da Biblioteca MPI. A divisão e balanceamento de carga entre os processadores são feitos pelo pacote METIS. Essa metodologia foi testada acoplada a métodos, seja iterativo (ICCG), seja direto (LU) na etapa de resolução dos sistemas referentes aos nós internos de cada partição. Para a resolução do sistema envolvendo os nós de fronteira, no caso do Complemento de Schur, utilizou-se uma implementação paralisada do Método de Gradientes Conjugados (PCG). S~ao discutidos aspectos relacionados ao desempenho dessas técnicas quando aplicadas em sistemas de grande porte. As técnicas foram testadas na solução de problemas de aplicação do Método de Elementos Finitos na Engenharia Elétrica (Magnetostática, Eletrocinética e Magnetodinâmica), sejam eles de natureza bidimensional com malhas não estruturadas, seja tridimensional, com malhas estruturadas. === This work presents the study of Domain Decomposition and Parallel Processing Techniques applied to the solution of systems of algebraic equations issued from the Finite Element Analysis of Electromagnetic Phenomena. Both Schur Complement and Schwarz Additive techniques were implemented. They were adapted to solve the linear systems in Beowulf clusters with the use of MPI library for message exchange. The load balance among processors is made with the aid of METIS package. The methodology was tested in association to either iterative (ICCG) or direct (LU) methods in order to solve the system related to the inner nodes of each partition. In the case of Schur Complement, the solution of the system related to the boundary nodes was performed with a parallelized Conjugated Gradient Method (PCG). Some aspects of the peformance of these techniques when applied to large scale problems have also been discussed. The techniques has been tested in the simulation of a collection of problems of Electrical Engineering, modelled by the Finite Element Method, both in two dimensions with unstructured meshes (Magnetostatics) and three dimensions with structured meshes (Electrokinetics).
|