A coupled approach of discrete element method and pore network model (DEM-PNM) is developed for simulating particle–fluid flow. By this approach, the particle/solid flow is described at a particle scale by DEM and the fluid flow is described at a pore scale by PNM. The coupling scheme between DEM and PNM and the boundary conditions for realizing particle–fluid flows under different conditions are explicitly introduced. The capability is examined in three systems: a dynamic cell, a gas-fluidized bed, and a packed bed with lateral gas injection. The simulated results by the DEM-PNM approach agree with those results by either DEM-LBM (lattice Boltzmann method) approach or experiments. In particular, the simulated reasonable flow patterns in gas fluidization and lateral gas injection demonstrate the applicability of the proposed approach to complex particle–fluid flow. However, the current model presents less accuracy in dilute flow region and does not consider some important factors, which should be further improved in the future.