Let $n(x,y,z)$ be the unit outward normal of the surface $S$ of the cylinder $x^2 + y^2 \leq 1 , 0 \leq z \leq 3$.
We have to compute $\int_{S} v.n dS = ?$ where $v = xz \hat{i} + 2yz\hat{j} + 3xy\hat{k}$.
the unit normal vector is given by the gradient of $S$, But how to find the gradient I think I would take gradient of $x^2 + y^2 -4$ that is $(2x,2y,0)$.
Now how do we evaluate $\int_{S} v.n dS = ?$ that is $\int_{x}\int_{y} (2x^2z + 4y^2z dx dy$ where $x,y$ vary such that $x^2 + y^2 \leq 4$ ?
The integral breaks up into a top and bottom caps and the side of the cylinder. For the side the normal vector is clearly $\hat \rho$, while it is $\pm \hat z$ for the caps. You might want to convert $\hat \rho$ to Cartesian.
Note that the points on the side of the cylinder are of the form $2(\cos\phi,\sin\phi,z_0)$ so a unit tangent vector is $(-\sin\phi,\cos\phi,0)$ and a normal would be $(\cos\phi,\sin\phi,0)$ which you may convert to $x,y$.
However, it is simplest here to compute the dot product in Cartesian and the convert to cylindrical since your surface element (for the side) is then simply $2d\phi dz$.