This monograph presents a methodology for simulating non-isothermal viscoelastic incompressible fluid flows where the viscosity and the relaxation time of the fluid are temperature-dependent. The viscoelasticity of the fluid is modeled by the Oldroyd-B and PTT (Phan-Thien-Tanner) models, where the temperature-dependent parameters are modeled by the WLF (Williams-Landel-Ferry) formulation. The numerical methodology used to solve the non-isothermal model is based on the MAC method for viscoelastic fluid flows via projection method. In this methodology, the Navier-Stokes equations and the Oldroyd-B and PTT constitutive equations are discretized by the finite difference method on a staggered grid. The numerical method was verified by simulation two-dimensional non-isothermal Poiseuille flow. Finally, the numerical methodoly was apllied for solving the 4 : 1 contraction problem in order to analyze the influence of parameters on the vortex dynamic.Esta dissertação apresenta uma metodologia para a simulação de escoamentos incompressíveis viscoelásticos não-isotérmicos, onde a viscosidade e o tempo de relaxação do fluido são dependentes da temperatura. A viscoelasticidade do fluido é modelada pelas equações constitutivas Oldroyd-B e PTT (Phan-Thien-Tanner), onde os parâmetros dependentes da temperatura são modelados pela relação WLF (Willians-Landel-Ferry). A metodologia numérica empregada para resolver o modelo não-isotérmico é baseada no método MAC para escoamentos viscoelásticos via método de projeção. Nesta metodologia, as equações de Navier-Stokes e as equações constitutivas Oldroyd-B e PTT são discretizadas pelo método de diferenças finitas em uma malha deslocada. A metodologia foi verificada na simulação do escoamento não-isotérmico bidimensional entre placas paralelas Poiseuille Flow. Finalmente, a metodologia numérica foi aplicada para resolver o problema da contração 4:1, onde são analisados os efeitos das variação de parâmetros na dinâmica dos vórtices.