É bem sabido que métodos clássicos de elementos finitos e diferenças finitas para o problema de Helmholtz apresentam efeito de poluição, que pode deteriorar seriamente a qualidade da solução aproximada. Controlar o efeito de poluição é especialmente difícil quando são utilizadas malhas não uniformes. Para malhas uniformes com elementos quadrados são conhecidos métodos (p. e. o QSFEM, proposto por Babuska et al) que minimizam a poluição. Neste trabalho apresentamos inicialmente dois métodos de elementos finitos de Petrov-Galerkin com formulação relativamente simples, o RPPG e o QSPG, ambos com razoável robustez para certos tipos de distorções dos elementos. O QSPG apresenta ainda poluição mínima para elementos quadrados. Em seguida é formulado o QOFD, um método de diferenças finitas aplicável a malhas não estruturadas. O QOFD apresenta grande robustez em relação a distorções, mas requer trabalho extra para tratar problemas não homogêneos ou condições de contorno não essenciais. Finalmente é apresentado um novo método de elementos finitos de Petrov-Galerkin, o QOPG, que é formulado aplicando a mesma técnica usada para obter a estabilização do QOFD, obtendo assim a mesma robustez em relação a distorções da malha, com a vantagem de ser um método variacionalmente consistente. Resultados numéricos são apresentados ilustrando o comportamento de todos os métodos desenvolvidos em comparação com os métodos de Galerkin, GLS e QSFEM. / It is well known that classical finite elements or finite difference methods for Helmholtz problem present pollution effects that can severely deteriorate the quality of the approximate solution. To control pollution effects is especially difficult on non uniform meshes. For uniform meshes of square elements pollution effects can be minimized with the Quasi Stabilized Finite Element Method (QSFEM) proposed by Babusv ska el al, for example. In the present work we initially present two relatively simple Petrov-Galerkin finite element methods, referred here as RPPG (Reduced Pollution Petrov-Galerkin) and QSPG (Quasi Stabilized Petrov-Galerkin), with reasonable robustness to some type of mesh distortion. The QSPG also shows minimal pollution, identical to QSFEM, for uniform meshes with square elements. Next we formulate the QOFD (Quasi Stabilized Finite Difference) method, a finite difference method for unstructured meshes. The QOFD shows great robustness relative to element distortion, but requires extra work to consider non-essential boundary conditions and source terms. Finally we present a Quasi Optimal Petrov-Galerkin (QOPG) finite element method. To formulate the QOPG we use the same approach introduced for the QOFD, leading to the same accuracy and robustness on distorted meshes, but constructed based on consistent variational formulation. Numerical results are presented illustrating the behavior of all methods developed compared to Galerkin, GLS and QSFEM.
| Identifer | oai:union.ndltd.org:IBICT/oai:agregador.ibict.br.BDTD_LNCC:oai:lncc.br:65 |
| Date | 02 March 2009 |
| Creators | Daniel Thomas Fernandes |
| Contributors | Fréderic Gerard Christian Valentin, Alexandre Loureiro Madureira, Eduardo Gomes Dutra do Carmo, Fernando Alves Rochinha, Saulo Pomponet Oliveira, Gustavo Benitez Alvarez, Abimael Fernando Dourado Loula |
| Publisher | Laboratório Nacional de Computação Científica |
| Source Sets | IBICT Brazilian ETDs |
| Language | Portuguese |
| Detected Language | Portuguese |
| Type | info:eu-repo/semantics/publishedVersion, info:eu-repo/semantics/doctoralThesis |
| Format | application/pdf |
| Source | reponame:Biblioteca Digital de Teses e Dissertações do LNCC, instname:Laboratório Nacional de Computação Científica, instacron:LNCC |
| Rights | info:eu-repo/semantics/openAccess |
Page generated in 0.0027 seconds