Resumo. Neste trabalho será realizado um estudo de estabilidade da solução totalmente discretizada do problema de condução de calor em regime transiente analisando a influência do parâmetro de estabilização e do ∆t. O problema de condução de caloré caracterizado por equações parabólicas e a variável espacial será discretizada usando um método de elementos finitos híbrido estabilizado combinado com aproximações de diferenças finitas, método de Euler implícito, para variável temporal. O método híbrido estabilizado de elementos finitos consiste no acoplamento de problemas locais, onde a solução da variável primalé obtida, com um problema global para os multiplicadores de Lagrange, identificado como traço da temperatura, e tem a continuidade imposta de forma fraca. A metodologia de implementação utilizada para a resolução do problemaé a denominada Condensação Estática que tem como vantagem o fato de ser mais eficiente do ponto de vista computacional. Os resultados numéricos mostrados comprovam taxasótimas de convergência na norma L 2 (Ω) e não apresentam oscilações espúrias para tempos muito pequenos. Palavras-chave. Elementos Finitos, Métodos Híbridos, Análise de Estabilidade, Equação do Calor
IntroduçãoProblemas de condução de calor em regime transiente são comumente representados por equações diferenciais classificadas como parabólicas. As abordagens mais conhecidas são baseadas em formulações semidiscretas de elementos finitos para a aproximação espacial combinadas com esquemas de diferenças finitas para a aproximação temporal. O método de Galerkin clássico, queé usualmente definido de maneira que a aproximação espacial seja contínua entre os elementos da discretização,é bastante empregado para resolver numericamente essa classe de problemas. Contudo, quando esteé o escolhido e o passo do tempoé reduzido com um tamanho de malha fixo, oscilações espaciais espúrias aparecemà medida que o tempo aumenta, poluindo a solução nos tempos iniciais [5]. Com o intuito de 1 daianasb@lncc.br 2 jkfi@lncc.br 3