Estudo de cavidades acusticas usando o metodo de elementos finitos via Galerkin/minimos quadrados
| Autor: | Ahmida, Khaled Mohamed |
|---|---|
| Přispěvatelé: | Pavanello, Renato, 1959, Devloo, Philippe Remy Bernard, Iguti, Fernando, Universidade Estadual de Campinas. Faculdade de Engenharia Mecânica, Programa de Pós-Graduação em Engenharia Mecânica, UNIVERSIDADE ESTADUAL DE CAMPINAS |
| Rok vydání: | 2021 |
| Předmět: | |
| Zdroj: | Biblioteca Digital de Teses e Dissertações da Universidade Estadual de Campinas (UNICAMP) Universidade Estadual de Campinas (UNICAMP) instacron:UNICAMP |
| DOI: | 10.47749/t/unicamp.1996.118678 |
| Popis: | Orientador: Renato Pavanello Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica Resumo: A solução numérica da equação de Helmholtz, via o Método de Elementos Finitos ou Diferenças Finitas, possui uma característica dispersiva ao contrárío do que ocorre com a solução exata. Soluções discretas deste tipo são funções do número de onda discreto, que é dependente da freqüência. O Método de Elementos Finitos (MEF) via a formulação clássica de Galerkin pode ser aplicado na solução da equação de Helmholtz sem que exista um limite teórico para o número de onda a ser analisado. Todavia, quando utiliza-se o MEF via Galerkin para elevados números de onda, faz-se necessário o uso de malhas extremamente refinadas para se obter soluções com precisão satisfatória e dispersão mínima o que conduz ao custo computacional muitas vezes proibitivo. Como regra geral, procura-se resguardar uma resolução de malha da ordem de dez elementos por período o que conduz a um número de equações muito elevado a ser resolvido, na medida em que se deseja resolver problemas da equação de onda no domínio das médias e altas freqüências, como é o caso de interesse em acústica, tema deste trabalho. O método de Galerkin Mínimos Quadrados (GMQ), derivado a partir de uma modificação na forma integral do problema, pode eliminar a poluição ou o erro na solução por elementos finitos da equação de Helrnholtz em domínios unidimensionais. Em domínios bidimensionais, esta poluição numérica pode ser reduzida mas não eliminada. A aplicação do método GMQ na solução da equação de Helrnholtz, para elevados números de onda, permite o uso de malhas com resolução da ordem de somente quatro elementos por período, o que aumenta consideravelmente o alcance das soluções numéricas. Neste trabalho, o método GMQ é aplicado em problemas bidimensionais da equação de Helmholtz, baseando-se em resultados numérícos obtidos para exemplos unidimensionais. Exemplos de decaimento exponencial e propagação de onda para condições de contorno de Dirichlet, serão apresentados. Elementos lineares, triangulares e quadrilaterais, em malhas regulares e irregulares, serão utilizados para a obtenção dos resultados numéricos Abstract: The numerical solution of the Helrnholtz equation, as finite element or finite difference solution, does not preserve the nondispersive character of the exact solution. I?iscrete solutions of this type are functions of a discrete wave number that depends on the frequency. The Galerkin finite element method is capable of modeling increasingly higher wave numbers by refining the mesh. However, this may become prohibitively expensive, as long as, an acceptable resolution of ten elements per wave, according to mIe of thumb, is to be considered. This would result in a large number of equations to be solved. The method of Galerkin Least Square (GLS), here utilized, for the numerical solution of the Helrnholtz equation can possibly eliminate the dispersive numerical effects by modification of the variational model in one-dimensional problems. However, in twodimensional problems, it is not possible to eliminate the pollution in the finite element error but can still be reduced. A resolution of four elements per wave is required for good results. In this work, the GLS method is applied to two-dimensional problems described by Helrnholtz equation, using results obtained in one-dimensional problems. Exponential decay problems and wave propagation problems with Dirichlet boundary conditions are considered. Numerical results are obtained using linear triangular and quadrilateral finite elements Mestrado Mecânica dos Sólidos e Projeto Mecânico Mestre em Engenharia Mecânica |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|