This work studied the influence of the particle shape and flow regime on the lump coal breakage in pneumatic conveying using CFD-DEM simulation. A variety of agglomerates with different sphericities were modelled by the parallel bond method to analyse the breakage characteristics of lump coal. The numerical parameters, simulation conditions and CFD-DEM simulation results were separately validated by experimentation. To demonstrate the lump coal breakage process, the mechanism energy variation in the coal agglomerate was analysed. The fragmentation degree of the coal agglomerate was positively correlated with the energy difference in collision. The integrality ratio of the coal agglomerate increased with the particle sphericity and swirling number in pneumatic conveying. A quasi-periodical downgrade in the coal agglomerate integrality ratio was shown in the swirling flow. A remarkable improvement in the coal agglomerate integrality ratio from approximately 0.85 to 0.975 was observed in swirling flow pneumatic conveying. The regression relationship of the particle sphericity and swirling number fit the Exp3p2 exponential function well. (C) 2017 Elsevier B.V. All rights reserved.