The usefulness of discrete element method simulation for studying fluidization of nanoparticle agglomerates is explored. Nanoparticle agglomerates were simulated by using solid particles of equivalent sizes and densities. Validity of the present simulation was assessed through comparisons of simulation results and experimental observations of bed expansion, characteristic fluidization behaviour, and dense-bed settling. The simulation was then used to investigate initial bed expansion and bed uniformity under particulate fluidization conditions. The role of inter-agglomerate interparticle force in fluidization of nanoparticle agglomerates was examined. A stability analysis originally developed for addressing the transition from particulate to bubbling fluidization for conventional particles was used for predicting the start of bubbling in fluidized beds of nanoparticle agglomerates.