Numerical solutions of three-dimensional, incompressible and unsteady Navier-Stokes equations for constant diameter swirling pipe flows are used to study vortex breakdown, including the detailed flow structures in the bubble domain and the "tail" behind the bubble during the vortex breakdown, and a comparison is made between the numerical solutions and the experimental results.