Abstract: In this study, an investigation on particle hydrodynamic behaviors of gas-particle slug flow is modeled and numerically simulated. A novelty particle kinetic energy model at subgrid-scale level is developed to consider the effect of gas turbulence on particle movement and a four-way coupling approach is combined. Anisotropic particle dispersions and the multiphase interactions is revealed by using a second-order moment turbulence model and large eddy simulation is performed to solve the Euler-Euler two-fluid transport equations, Predictions are well agreed with the experimental data and reported results by discrete particle model built-in the multiphase flow model with interphase exchanges code. Bubble motions have great effects on the particle dynamics due to bubble entrainment. Flow patterns of particle circulation are generated between higher and lower concentrations. The largest frequency of power spectrum density at center of top surface is approximately 4.5 times larger than that of wall region. Bubblelike granular temperatures are much larger than smaller-scale particle granular temperature. At the top surface, maximum value of vertical bubblelike Reynolds stress is found at x/W = 0.28, it is 4.0 times larger than that of the lateral x-direction. In addition, the mean values, and the standard deviation of SGS particle kinetic energy at center are both approximately 3.5 times larger than those of z = 10.0 cm. For particle shear stress at middle section, maximum ratios of mean value and standard deviation at central regions to those of walls reach up to 8.0 and 11.0.