The purpose of this research is to develop the modeling and simulation of a frictional and compressible powder flow subjected to oscillatory shear stress in a device that consists of a rectangular shaped cavity with a free surface, two oscillatory moving walls, and three non-moving walls. The finite differences method under the explicit scheme is used to approximate the real solution to the resulting equations from the mathematical modeling of the studied system. The simulation results were compared qualitatively with experimental results, and the solution method was verified by comparing the simulation results of a viscous glycerin flow with the results of the software STAR-CCM+ under the same conditions of the powder flow. The results showed that the powder velocity profiles were coherent with the experimental results. Furthermore, the solution method was validated finding that the 83.5% of the compared data presented absolute errors less than 0.5 cm/s. The complete velocity profiles obtained from the simulation help to have a deep understanding of the granular material behavior subjected to oscillatory shear stress.