Maybe something like this?
ParametricPlot3D[{r*Cos[phi], r*Sin[phi], Re[u1[r, phi, 3]]}, {r, 0, 3}, {phi, 0, 2 Pi}, MeshFunctions -> {#3 &}]
