Efetuou-se um estudo numérico do escoamento entre cilindros concêntricos e excêntricos girantes com efeitos de dissipação viscosa presente. Assumiu-se que o cilindro interno estava girando com uma velocidade angular constante e que o externo estava fixo. Assumimos a condição de contorno de não deslizamento para velocidade em ambas as paredes. A temperatura do fluido nas paredes dos cilindros foi assumida constante e seu campo assumimos ser governado por efeitos de convecção e condução, com uma fonte de energia proveniente da dissipação viscosa do escoamento. A viscosidade do fluido foi assumida variável com a temperatura. A discretização espacial foi feita por elementos finitos através de funções de base lineares, mais uma função estabilizadora para velocidade no caso dos cilindros excêntricos. O método das características foi usado no problema excêntrico. Resultados computacionais ilustram a eficiência da técnica adotada.