Fully resolved simulations of particulate and aggregative fluidization systems are performed successfully with the so-called combined lattice Boltzmann method and time-driven hard-sphere model(LBM-TDHS).In this method,the discrete particle phase is described by time-driven hard-sphere model,and the governing equations of the continuous fluid phase are solved with lattice Boltzmann method.Particle-fluid coupling is implemented by immersed moving boundary method.Time averaged flow structure of the simulated results show the formation of core-annulus structure and sigmoid distribution of voidage in the axial direction,which are typical phenomena in fluidization systems.Combining the results of the simulation,the energy consumption N_(st)for suspending and transporting solids is calculated from the direct numerical simulation(DNS)of fluidization,and the stability criterion N_(st)/N_T=min proposed in EMMS/bubbling model is verified numerically.Furthermore the numerical results show that the value of N_(st)/N_T in particulate fluidization is much higher than that in aggregative fluidization,but N_(st)/N_T=min is effective for both particulate and aggregative fluidization. Fully resolved simulations of particulate and aggregative fluidization systems are performed successfully with the so-called combined lattice Boltzmann method and time-driven hard-sphere model (LBM-TDHS). In this method, the discrete particle phase is described by time-driven hard -sphere model, and the governing equations of the continuous fluid phase are solved with lattice Boltzmann method. Particle-fluid coupling is implemented by immersed moving boundary method. Time averaged flow structure of the simulated results show the formation of core-annulus structure and sigmoid distribution of voidage in the axial direction, which are typical phenomena in fluidization systems. Combining the results of the simulation, the energy consumption N_ (st) for suspending and transporting solids is from from direct numerical simulation (DNS) of fluidization, and the stability criterion N_ (st) / N_T = min proposed in EMMS / bubbling model is verified numerically.Furthermore the numerical results show that the value of N st / N_T in particulate fluidization is much higher than that in aggregative fluidization, but N st / N_T = min is effective for both particulate and aggregative fluidization.