A numerical method for solving steady 2D and axisymmetrical viscoelastic flow problems with an application to inertia effects in contraction flows