Solução Numérica do Modelo Giesekus para Escoamentos Através de uma Contração Planar Bidimensional 4:1

Abstract

Este trabalho apresenta um método numérico, que utiliza diferenças finitas, para simular escoamentos viscoelásticos bidimensionais governados pela equação constitutiva Giesekus [2]. A velocidade e a pressão são calculadas por um método implícito, enquanto que a equação constitutiva Giesekus é resolvida por um método Runge-Kutta de segunda ordem. A convergência do método é verificada por meio de refinamento de malha. O código resultante é aplicado para simular o escoamento em uma contração planar 4:1. Resultados numéricos utilizando vários números de Reynolds e Weissenberg são apresentados.