Use of boundary-fitted curvilinear coordinates for the numerical simulation of complex viscoelastic flows