Solid-liquid two-phase flow widely existed in daily life. It was widely used in petrochemical, medical, environmental protection and other fields. The common particles in solid-liquid two-phase flow were spherical, cylindrical and rectangular. According to the elongated shape characteristics and particle orientation of cylindrical catalyst particles relative to spherical particles, five cylindrical particles with diameters of 2 mm and different lengths were selected to study. The computational fluid dynamics coupled discrete element method (CFD-DEM) was used to simulate the motion behavior of cylindrical particles during the settlement of tubular containers. Cylindrical particles were constructed by the combination of spherical element method and multi sphere method. A cylindrical particle sedimentation test-bed was established, and the sedimentation process of cylindrical particles in tubular container was verified by high-speed camera. The results showed that the particles near the wall were affected by the wall effect. When cylindrical particles with the same diameter and length were released at different positions, the particles drifted to the center in the process of sedimentation. The settling velocity of particles near the wall was slower and the settling time was longer. When releasing the columnar particles at the same position of the circular tube by only changing the angle between the cylindrical particle's principal axis and the horizontal plane, the particles would eventually rotate to a state where the principal axis was parallel to the horizontal plane. In the process of sedimentation, the greater the angle between the particle principal axis and the horizontal plane, the greater the resistance on the windward surface of the particle, and the particle rotation time would also increase. The Stokes equation for the settlement of cylindrical particles was derived, and the resistance coefficient in the Stokes equation was modified through the experimental data. The resistance coefficient corrected by the experiment was introduced into user defined function (UDF), and the end velocity of particle settlement was simulated and calculated. When using the resistance coefficient of spherical particles, the error was 40.6%. After correction, the relative error was reduced from 50% to less than 17%.