The present work is a 2-D numerical simulation of pulsed fluidized bed with immersed tubes using DEM-LES coupling method. The pulsed inflow of gas phase is modeled as U-0(1 + sin(2 pi ft)), in which four pulsating frequencies of f = 5, 10, 15 and 20 of velocity inflow are used. The discrete element method (DEM) simulation for particle motion coupled with the large eddy simulation (LES) for gas phase is used. The fluidized bed with five immersed tubes of staggered arrangement and six immersed tubes of in-line arrangement is simulated, respectively. It is found that the pressure drop, the mean drag force and the mean pressure gradient force experienced by particles are forced oscillated. The different effects of pulsed fluidization on the circumferential distribution of particle-tube collision on the outer surface of tubes at different pulse frequencies and modes of arrangement of immersed tubes are numerically analyzed. Finally, it is found that the pulsed motion of fluid with high frequency leads to suppression of particle fluctuating motion. (C) 2009 Elsevier Ltd. All rights reserved.