Abstract: The study of rheological properties of two-phase flow, especially of turbulent flow, is one of the key issues of gas-solid two-phase flow of slender particles. The study of rheological properties is more difficult due to lack of efficient coupling relation between slender particles and turbulent flow. On the basis of the three-dimensional force and motion model of slender particles in the light of Discrete Element Method（DEM） and rigid dynamics respectively, a three-dimensional two-way coupling model between slender particles and turbulent flow was founded based on the coupling correlation between Lagrangian time scales and ?-ε model, and numerical analysis on a cold-state two-phase flow of slender particles was implemented using the model. It is found that the axes of most slender particle are nearly parallel to the direction of main stream of flow during the fluidization; the volume fraction distribution of flow field is every uneven due to interaction between slender particles; at height direction, on the one hand the residence of slender particles results in reduction of section surface of turbulent flow thereby results in rising of turbulent kinetic energy, and on the other hand slender particles residing in local turbulent flow will consume some turbulent kinetic energy thus lead to decline of local turbulent kinetic energy; at horizontal direction local turbulent kinetic energy will decline where some slender particles reside; the residence of slender particles result in the increment of pressure gradient of the turbulent flow field.