Resumo: O objetivo deste trabalho e reduzir o tempo de CPU (Central Processing Unit) necessario para resolver problemas difusivos bidimensionais, discretizados com malhas anisotropicas. Os modelos matematicos considerados referem-se a tres problemas bidimensionais lineares de conducao de calor, governados pelas equacoes de Laplace e Poisson, com condicoes de contorno de Dirichlet. O metodo de diferencas finitas e usado para discretizar as equacoes diferenciais com esquema de diferenca central (CDS) de segunda ordem. Os sistemas de  equacoes algebricas sao resolvidos usando-se os metodos Gauss-Seidel lexicografico e redblack, associados ao metodo multigrid geometrico com esquema de correcao (CS) e ciclo V. Foram resolvidos problemas com anisotropia geometrica, diversas malhas e razoes de aspecto. O numero de iteracoes internas (ƒË) foi verificado em um intervalo de 1 a 3.000. A analise do numero de niveis foi realizada utilizando-se o numero maximo de niveis ( ) maximo L e .1 maximo L , . 2, . 3 maximo maximo L L , .4 maximo L . Sao feitas comparacoes entre diversos algoritmos de engrossamento: engrossamento padrao (EP), semi-engrossamento (SE), semiengrossamento completo (SEC),  emi-engrossamento seguido de engrossamento padrao (SEEP) e engrossamento padrao seguido de semi-engrossamento (EP-SE). Tambem sao realizadas comparacoes entre alguns operadores de restricao: injecao, meia ponderacao e ponderacao completa. Sao propostos tres tipos de restricao para problemas anisotropicos: meia ponderacao geometrica, ponderacao geometrica completa e ponderacao parcial. O processo de prolongacao utilizado e a interpolacao bilinear. Tambem foi investigado o efeito sobre o tempo de CPU causado por: numero de pontos na malha (N); numero de iteracoes internas no solver (v); e numero de malhas (L). Verificou-se que: o algoritmo SE-EP e o mais rapido entre os cinco algoritmos testados; e confirmou-se que, para problemas isotropicos e anisotropicos, o solver Gauss-Seidel red-black com restricao por ponderacao parcial resulta em menor tempo de CPU em relacao ao Gauss-Seidel lexicografico

Oliveira, Fabiane de

Orientador : Carlos Henrique MarchiCo-Orientador: Prof. Dr. Marcio Augusto Villela PintoTese (doutorado) - Universidade Federal do Paraná, Setor de Tecnologia, Programa de Pós-Graduação em Engenharia Mecânica. Defesa: Curitiba,26/05/2010Bibliografia: fls. 149-153Resumo: O objetivo deste trabalho é reduzir o tempo de CPU (Central Processing Unit) necessário para resolver problemas difusivos bidimensionais, discretizados com malhas anisotrópicas. Os modelos matemáticos considerados referem-se a três problemas bidimensionais lineares de condução de calor, governados pelas equações de Laplace e Poisson, com condições de contorno de Dirichlet. O método de diferenças finitas é usado para discretizar as equações diferenciais com esquema de diferença central (CDS) de segunda ordem. Os sistemas de equações algébricas são resolvidos usando-se os métodos Gauss-Seidel lexicográfico e red-black, associados ao método multigrid geométrico com esquema de correção (CS) e ciclo V. Foram resolvidos problemas com anisotropia geométrica, diversas malhas e razões de aspecto. O número de iterações internas (?) foi verificado em um intervalo de 1 a 3.000. A análise do número de níveis foi realizada utilizando-se o número máximo de níveis (Lmáximo) Lmáximo - 1, Lmáximo - 2, Lmáximo - 3, Lmáximo - 4. São feitas comparações entre diversos algoritmos de engrossamento: engrossamento padrão (EP), semi-engrossamento (SE), semiengrossamento completo (SEC), semi-engrossamento seguido de engrossamento padrão (SE-EP) e engrossamento padrão seguido de semi-engrossamento (EP-SE). Também são realizadas comparações entre alguns operadores de restrição: injeção, meia ponderação e ponderação completa. São propostos três tipos de restrição para problemas anisotrópicos: meia ponderação geométrica, ponderação geométrica completa e ponderação parcial. O processo de prolongação utilizado é a interpolação bilinear. Também foi investigado o efeito sobre o tempo de CPU causado por: número de pontos na malha (N); número de iterações internas no solver (v); e número de malhas (L). Verificou-se que: o algoritmo SE-EP é o mais rápido entre os cinco algoritmos testados; e confirmou-se que, para problemas isotrópicos e anisotrópicos, o solver Gauss-Seidel red-black com restrição por ponderação parcial resulta em menor tempo de CPU em relação ao Gauss-Seidel lexicográfico.Abstract: The purpose of this work is to reduce the CPU (Central Processing Unit) time necessary to solve two-dimensional diffusive problems, discretized with anisotropic grids. The mathematical models considered are related to three two-dimensional linear heat diffusion problems, governed by Laplace and Poisson equations, with Dirichlet boundary conditions. The differential equations are discretized by the finite difference method, with second order approximations, given by central differences scheme (CDS) for derivatives. The systems of equations are solved with the lexicographical and red-black Gauss-Seidel methods, associated to geometric multigrid with correction scheme (CS) and V-cycle. For problems with geometric anisotropic, several grids and aspect ratios are considered. The inner iterations are verified for the interval between 1 and 3.000. The analyses for the number of levels are accomplished using the maximum level number (Lmáximo) and Lmáximo - 1, Lmáximo - 2, Lmáximo - 3, Lmáximo - 4. Comparisons among several semi-coarsening algorithms are made: standard coarsening (EP), semi-coarsening (SE), full semi-coarsening (SEC), semi-coarsening followed by standard coarsening (SE-EP) and standard coarsening followed by semi-coarsening (EP-SE). Comparisons are also made among some restrictions schemes: injection, half weighting and full weighting. Three restriction schemes for anisotropic problems are proposed: geometric half weighting, geometric full weighting and partial weighting.  The prolongation process used is the bilinear interpolation. The CPU time changes caused by the following effects are also studied: the number of nodes of a grid (N); the number of inner iterations of the solver (v); and the number of grids used (L). It was verified that: the algorithm SE-EP is the fastest of the five algorithms tested; for isotropic and anisotropic problems, red-black Gauss-Seidel solver with partial weighting restriction results in smaller CPU time compared to lexicographical Gauss-Seidel

Repositório Digital Institucional da UFPR

Efeito de malhas anisotrópicas bidimensionais sobre o desempenho do método Multigrid Geométrico

Orientador : Carlos Henrique MarchiCoorientador: Prof. Dr. Marcio Augusto Villela PintoTese (doutorado) - Universidade Federal do Paraná, Setor de Tecnologia, Programa de Pós-Graduação em Engenharia Mecânica. Defesa: Curitiba,26/05/2010Bibliografia: fls. 149-153Resumo: O objetivo deste trabalho é reduzir o tempo de CPU (Central Processing Unit) necessário para resolver problemas difusivos bidimensionais, discretizados com malhas anisotrópicas. Os modelos matemáticos considerados referem-se a três problemas bidimensionais lineares de condução de calor, governados pelas equações de Laplace e Poisson, com condições de contorno de Dirichlet. O método de diferenças finitas é usado para discretizar as equações diferenciais com esquema de diferença central (CDS) de segunda ordem. Os sistemas de equações algébricas são resolvidos usando-se os métodos Gauss-Seidel lexicográfico e red-black, associados ao método multigrid geométrico com esquema de correção (CS) e ciclo V. Foram resolvidos problemas com anisotropia geométrica, diversas malhas e razões de aspecto. O número de iterações internas (?) foi verificado em um intervalo de 1 a 3.000. A análise do número de níveis foi realizada utilizando-se o número máximo de níveis (Lmáximo) Lmáximo - 1, Lmáximo - 2, Lmáximo - 3, Lmáximo - 4. São feitas comparações entre diversos algoritmos de engrossamento: engrossamento padrão (EP), semi-engrossamento (SE), semiengrossamento completo (SEC), semi-engrossamento seguido de engrossamento padrão (SE-EP) e engrossamento padrão seguido de semi-engrossamento (EP-SE). Também são realizadas comparações entre alguns operadores de restrição: injeção, meia ponderação e ponderação completa. São propostos três tipos de restrição para problemas anisotrópicos: meia ponderação geométrica, ponderação geométrica completa e ponderação parcial. O processo de prolongação utilizado é a interpolação bilinear. Também foi investigado o efeito sobre o tempo de CPU causado por: número de pontos na malha (N); número de iterações internas no solver (v); e número de malhas (L). Verificou-se que: o algoritmo SE-EP é o mais rápido entre os cinco algoritmos testados; e confirmou-se que, para problemas isotrópicos e anisotrópicos, o solver Gauss-Seidel red-black com restrição por ponderação parcial resulta em menor tempo de CPU em relação ao Gauss-Seidel lexicográfico.Abstract: The purpose of this work is to reduce the CPU (Central Processing Unit) time necessary to solve two-dimensional diffusive problems, discretized with anisotropic grids. The mathematical models considered are related to three two-dimensional linear heat diffusion problems, governed by Laplace and Poisson equations, with Dirichlet boundary conditions. The differential equations are discretized by the finite difference method, with second order approximations, given by central differences scheme (CDS) for derivatives. The systems of equations are solved with the lexicographical and red-black Gauss-Seidel methods, associated to geometric multigrid with correction scheme (CS) and V-cycle. For problems with geometric anisotropic, several grids and aspect ratios are considered. The inner iterations are verified for the interval between 1 and 3.000. The analyses for the number of levels are accomplished using the maximum level number (Lmáximo) and Lmáximo - 1, Lmáximo - 2, Lmáximo - 3, Lmáximo - 4. Comparisons among several semi-coarsening algorithms are made: standard coarsening (EP), semi-coarsening (SE), full semi-coarsening (SEC), semi-coarsening followed by standard coarsening (SE-EP) and standard coarsening followed by semi-coarsening (EP-SE). Comparisons are also made among some restrictions schemes: injection, half weighting and full weighting. Three restriction schemes for anisotropic problems are proposed: geometric half weighting, geometric full weighting and partial weighting.  The prolongation process used is the bilinear interpolation. The CPU time changes caused by the following effects are also studied: the number of nodes of a grid (N); the number of inner iterations of the solver (v); and the number of grids used (L). It was verified that: the algorithm SE-EP is the fastest of the five algorithms tested; for isotropic and anisotropic problems, red-black Gauss-Seidel solver with partial weighting restriction results in smaller CPU time compared to lexicographical Gauss-Seidel

Universidade Federal do Paraná

https://acervodigital.ufpr.br/bitstream/1884/24274/1/Tese_V5_17_agosto.pdf

Efeito de malhas anisotrópicas bidimensionais sobre o desempenho do método Multigrid Geométrico

Abstract

Similar works

Full text

Available Versions

Repositório Digital Institucional da UFPR

Universidade Federal do Paraná