(OverVector[U])·∇(OverVector[U]) = -(OverVector[U])×∇×(OverVector[U])+∇((U^2)/2)
in cylindrical coordinates

Your tasks:
 ☺ After initializing the "kernel" by executing UDelU.nb, add a section here that finds (OverVector[U])·∇(OverVector[U]) in cylindrical coordinates.  (Easy?)