Abstract: Gas SGS turbulent kinetic energy equations were derived based on a single-phase flow derivation method under the framework of a gas large-eddy simulation model and a particle second-order moment model (LES-SOM model), and a sub-grid scale (SGS) turbulent kinetic energy model was proposed for gas-solid two-phase flow. The particle second-order moment model was used considering anisotropic characteristics of particle velocity fluctuation. The turbulent behavior of gas-solid two-phase flow in a circulating fluidized bed was simulated using the LES-SOM model. The concentration distribution and particle velocity predicted agree better with experimental data than those of the simulations without considering SGS turbulent kinetic energy model. Different SGS turbulent models were compared and the turbulent characteristics of gas are better simulated by the SGS turbulent kinetic energy model than those of the Smagorinsky model. With the increase of gas velocity, the SGS turbulent kinetic energy and dissipation are increased and the non-uniform characteristics of radial distribution are enhanced.