A method of solution to solve the compressible 3D Navier-Stokes Equations in cylindrical co-ordinates coupled to the continuity equation in cylindrical coordinates is presented in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow. A dimensionless parameter is introduced whereby in the large limit case a method of solution is sought for in the boundary layer of the tube. It is concluded that the total divergence of the flow can be expressed as the integral with respect to time of the line integral of the dot product of inertial and azimuthal velocity. The line integral is evaluated on a contour that is annular and traces the boundary layer as time increases in the flow. Finally it's explicit dependence on the gradient of a density function and Laplacian of pressure is shown with complete dependence on the vorticity and it's rate of change for a flow associated with the Navier Stokes equations.