The stream function for a plane potential incompressible doublet flow is supposed to be
psi = -(Q/2pi)*arctan((2arsin(theta)/(r^2-a^2))
Why?
Q = superficial flow rate (volume flow rate/unit length)
r, theta = polar coords
2a = length of the doublet.
Thanks, people.