Oh man, you are lucky I'm in the middle of a numerical run right now.
From Viscous Fluid Flow(3rd ed.), Frank White, Chapter 4, Section 10.6.
If a round jet emerges from a circular hole with sufficient momentum, it remains narrow and grows slowly, the radial changes \partial /\partial r being much larger than axial changes \partial / \partial x
Continuity:
\frac{\partial u}{\partial x} + \frac{1}{r}\frac{\partial}{\partial r}(rv) = 0
x momentum
u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = \frac{v}{r}\frac{\partial}{\partial r}\left(r\frac{\partial u}{\partial r}\right)
Schlichting reasoned that the jet thickness grew linearly, so the similarity variable is r/x. he defined the stream function
\psi(r,x) = \nu x F(\eta)\,\, \eta = \frac{r}{x}
From which the axisymmetrical velocity components are:
u = \frac{1}{r}\frac{\partial \psi}{\partial r} = \frac{\nu F'}{r}
v = -\frac{1}{r}\frac{\partial \psi}{\partial x} = \frac{\nu}{r}(\eta F' - F)
Substitution into the x-momentum equation gives the following third-order non-linear differential equation.
\frac{d}{d\eta}\left( F'' - \frac{F'}{\eta}\right) = \frac{1}{\eta^2} (FF'' - \eta F'^2 - \eta FF'')
The boundary conditions are F(0) = F'(0) = F'(infinity) = 0. The exact solution is:
F = \frac{(C\eta)^2}{1 + (C\eta /2)^2}
Where C is a constant determined from the momentum of the jet
J = \rho \int^\infty_0 u^2 2\pi r\,dr = \frac{16\pi}{3}\rho C^2 \nu^2
C = \left( \frac{3J}{16\pi \rho \nu^2}\right)^{1/2}
The axial jet velocity is then:
u = \frac{3J}{8\pi \mu x}\left( 1 + \frac{C^2 \eta^2}{4}\right)^{-2}
The term in parenthesis is the shape of the jet profile. The jet centerline velocity drops off as x^{-1}. The mass flow rate across any axial section of the jet is:
\dot{m} = \rho \int^{\infty}_0 u2\pi r\,dr = 8\pi \mu x
IF the jet is laminar, and we assume a simple plane laminar jet, we find that the velocity distribution is:
u(x,y) = u_{max} \sech^2 a\eta Or:
u(x,y) = u_{max}\sech^2 \left[ 0.2752 \left(\frac{J\rho}{\mu^2 x^2}\right)^{1/3} y \right]
Where J is the momentum flux. At this point, if we define width of the jet as twice the distance y where u = 0.001u_{max} then:
\mbox{Width} = 21.8 \left(\frac{x^2 \mu^2}{J \rho}\right)^{1/3}
That may be slightly more helpful on the width of the jet question. On the effect of the surface...no idea.