A robust MPS model to simulate granular flows is developed.

A constitutive law to describe granular flows is validated by MPS model.

Glass bead and sand flows are extensively modeled, and surface profiles and front wave movements are simulated well.

The surface flowing and quasi-static internal structure in granular flows are clearly produced

The simulation results show that initial column ratio plays a significant role in granular spreading.