É proposta neste trabalho uma nova metodologia para resolução das equações governantes de fluidos viscoelásticos, baseada no método dos volumes finitos, usando o arranjo co-localizado para as variáveis e malhas estruturadas. São utilizadas aproximações de alta ordem para os fluxos lineares e não-lineares médios nas interfaces dos volumes, e para os termos não-lineares que surgem da discretização das equações constitutivas. Nesta metodologia, os valores médios das variáveis nos volumes são usados durante todo o procedimento de resolução, e os valores pontuais são obtidos ao final, através da deconvolução dos valores médios. O sistema de equações discretizadas é resolvido de forma simultânea, pelo método de Newton. A metodologia é exemplificada para um problema clássico em mecânica de fluidos computacional, o escoamento stick-slip, usando como equação constitutiva o modelo de Oldroyd-B. As soluções obtidas apresentaram boa precisão, sendo livres de oscilações mesmo em regiões de grandes gradientes das variáveis.
Fluidos viscoelásticos; simulação; método dos volumes finitos
