Abstract: Abstract: The clogging of labyrinth channel emitters is of great concern to many researchers. It is great significance to carry out the flow analysis in the labyrinth channel, to master a more comprehensive understanding of two-phase flow mechanism, and to achieve the visualization of the fluid flow and particle movement inside the labyrinth channel for the clogging of the emitters. In this paper, a high quality structured mesh is used on the basis of reference based on the ICEM CFD software. The applicability of a series of k-ε turbulence models is systematically investigated and the particle movement law is further analyzed based on the Lagrangian method for discrete phase model. Compared with the experimental data cited form the literature Liming Yu et al. (2009), this paper, the various forces of single-particle were comprehensively considered, and the change law of the trajectory of single-particle with different densities and diameters were predicted, the variation of the migration velocity between single-particle and its surrounding fluid under different diameters and the variation law of the drag force of single-particle were analyzed, et al. The adaptability of 3 turbulence models (standard k-ε, RNG k-ε and Realizable k-ε) in the labyrinth channel emitter were analyzed, the average errors predicted by RNG k-ε and Realizable k-ε model are 8.96% and 8.11% respectively, while the average error predicted by standard k-ε model is only 2.32%, the standard k-ε model is more suitable for simulating of the pure continuous-phase complex turbulent flows inside a labyrinth channel emitter. The computed trajectory characteristics of single-particle showed that the stochastic trajectory model, which considered gravity, buoyancy, virtual mass force and drag force, can be employed to capture precisely the movement of the particles in labyrinth channel. The paths of magnesium particle, sand particle and aluminum particles passing through the labyrinth channel were 43.97 mm, 44.25 mm and 50.79 mm, respectively, while the paths of 65, 100 and 150 μm sand particles passing through the labyrinth channel were 34.04, 44.25 and 49.06 mm respectively, so the change of particle density had less effect on the particle trajectory than that of particle diameter. With the increase of particle density, the fluid flow field was greatly disturbed by the particle, the mainstream region and low velocity region begin to mix and the separation of the particle motion and the liquid flow was obvious. With the increase of the particle diameter, the velocity deviation between singe-particle and the surrounding fluid increased continuously, the disorderly phenomenon of the particle-phase and liquid-phase velocities peak emerges, the changes of the velocity amplitude and the mean velocity was decreased, the drag force of single-particle was increased, and the trajectory of single-particle becomes chaotic. Compared the numerical results of literature Yu Liming et al (2009), the average error of the particle movement velocity was only 2.34%, and the single-particle trajectories were more in accordance with the objective law. With the increase of particle diameter or density, the particle velocity decreased totally, the change of particle velocity amplitude from large value to small value were in order of minimum velocity, average velocity and maximum velocity. As the particle diameter or density increases, the particle velocity decreased greatly when the particle flowed through the vortex zone, and the following behaviors of the particles became worse, and the particles possibly deposit on the inside of the corner of the vortex zone under the action of the inertial force, which probably leads to the clogging in the labyrinth channel of the emitter. The analyses can provide a reference for two-phase flow mechanism and anti-clogging design for the emitter.