Este trabalho apresenta o desenvolvimento de um método numérico para simular escoamentos viscoelásticos com superfícies livres de um fluido governado pelo modelo de Oldroyd-B. As equações governantes para um fluido Oldroyd-B são consideradas. A derivada temporal é aproximada por um método de segunda ordem. Uma formulação para o cálculo do tensor de tensão extra nos contornos rígidos é apresentada. As equações governantes são resolvidas pelo método de diferenças finitas numa malha deslocada utilizando variáveis primitivas. O método numérico descrito neste trabalho foi implementado no ambiente de simulação Freeflow-2D. Resultados numéricos demonstrando que o método numérico empregado neste trabalho aplicado a vários escoamentos bidimensionais de um fluido Oldroyd-B são apresentados.