Streamlines for creeping flow 
past a solid sphere
Author's Data: Housam BINOUS
Department of Chemical Engineering
National Institute of Applied Sciences and Technology
Tunis, TUNISIA
Email: binoushousam@yahoo.com

U = 10; R = 1;

r = (Sqrt[x^2 + y^2])

θ = ArcSin[y/r]

ContourPlot[
      U R^2 Sin[θ]^2 (1/2 (r/R)^2 - 3/4 (r/R) + 
            1/4 (R/r)), {x, -5, 5}, {y, - 5,  5}, 
      Contours -> 30, 
      ColorFunction -> ((Hue[ArcTan[2.45 - 2.45 #]] &)), 
      DisplayFunction -> Identity]

Show[%, Graphics[Disk[{0, 0}, 1]], 
      DisplayFunction -> $DisplayFunction]