Return to search

New enriched element methods for unsteady reaction-advection-diffusion models / Novos métodos de elementos finitos enriquecidos aplicados a modelos de reação-advecção-difusão transientes

Several problems in physics and engineering are modeled by reaction-advection-diffusion (RAD) equations. However, when the diffusive terms are small compared with the other ones, these problems can become difficult to solve numerically. Besides, formulating the unsteady version of these models in a semi-discrete fashion, it can be interpreted that the overall diffusivity gets smaller as the time step decreases. To overcome these drawbacks, this thesis considers the development of Galerkin (or Petrov-Galerkin) finite element methods based on approximation spaces enriched by residual-free bubbles (RFB) or multiscale functions. Beginning with the unsteady reaction-diffusion problem, new methods using multiscale functions are presented which improve the solutions in the reaction-dominated regime and/or when small time steps are adopted. They also give rise to a general concept of stabilizing unsteady problems differently along the time. In the following, it is shown that switching RFB by suitable multiscale functions in the elements connected to the outflow boundaries of the domain increases the accuracy of the solutions in this region for RAD problems with advection. Next, this methodology is further studied for systems of RAD equations. In a final contribution, an extension of the RFB method is introduced for the shallow waters equations. All these methods are tested through benchmark problems and compared with stabilized methods presenting stable and accurate results. / A modelagem de vários problemas físicos e de engenharia envolve a solução de problemas de transporte do tipo reação-advecção-difusão (RAD), porém, estes podem tornar-se singularmente perturbados quando os termos difusivos são pequenos comparados aos demais. Além disso, ao adotar formulações semi-discretas em problemas transientes, observa-se que diminuir o passo de tempo tem um efeito de redução da componente difusiva. Para superar estas dificuldades, esta tese considera o desenvolvimento de métodos de elementos finitos de Galerkin (ou Petrov-Galerkin) baseados em espaços de aproximação enriquecidos por funções bolhas livres do resíduo (RFB) ou funções multiescala. Começando pelo problema de reação-difusão transiente, novos métodos utilizando funções multiescala são apresentados, os quais melhoram as soluções no regime reativo-dominante e/ou quando pequenos passos de tempo são adotados. Com estes métodos, discute-se também o conceito de estabilização variável ao longo do tempo para problemas transientes. Na seqüência, verifica-se que utilizar funções multiescala nos elementos conectados às fronteiras de saída de fluxo do domínio e RFB nos demais elementos aumenta a precisão das soluções nesta região em problemas de RAD com advecção dominante. A seguir, esta metodologia é estudada para sistemas de RAD. Como contribuição final, estende-se o método RFB para o modelo de águas rasas. Todos estes métodos são submetidos a testes de robustez e comparados com métodos estabilizados, apresentando resultados estáveis e precisos.

Identiferoai:union.ndltd.org:IBICT/oai:agregador.ibict.br.BDTD_LNCC:oai:lncc.br:14
Date20 December 2005
CreatorsJairo Valões de Alencar Ramalho
ContributorsMarcio Arab Murad, Eduardo Gomes Dutra do Carmo, Alexandre Loureiro Madureira, Fréderic Gerard Christian Valentin, Leopoldo Luis Cabo Penna Franca
PublisherLaboratório Nacional de Computação Científica
Source SetsIBICT Brazilian ETDs
LanguagePortuguese
Detected LanguagePortuguese
Typeinfo:eu-repo/semantics/publishedVersion, info:eu-repo/semantics/doctoralThesis
Formatapplication/pdf, application/pdf, application/pdf, application/pdf
Sourcereponame:Biblioteca Digital de Teses e Dissertações do LNCC, instname:Laboratório Nacional de Computação Científica, instacron:LNCC
Rightsinfo:eu-repo/semantics/openAccess

Page generated in 0.0023 seconds