This work investigated the influence of thickness-direction boundary conditions on the flow characteristics of granular material in a quasi-two-dimensional silo using the discrete element method (DEM). Two types of boundary conditions were considered in the thickness direction: wall conditions and periodic boundary conditions. The simulation results indicate that under wall conditions, velocity waves propagate upward, manifested by the formation of bubble-like sub-flow zones in the velocity field, and the particle motion in the upper bed region exhibits a clear stick–slip feature. In contrast, under periodic boundary conditions, particle motion displays a resonant mode. Further statistical analysis reveals that, despite the distinct macroscopic motion mode under the two boundary conditions, the probability distributions of particle vertical fluctuating velocities share similar characteristics: both exhibit fat-tailed and asymmetric features and deviate from Gaussian distribution. Additionally, under wall conditions, the horizontal distributions of particle vertical velocity conform to the kinematic model throughout the bed, whereas under periodic boundary conditions, the horizontal distributions in the upper bed region display plug flow characteristics. In summary, the results of this work demonstrate that thickness-direction boundary conditions play a crucial role in determining the flow characteristics of granular assembly in silos.
Hu et al. (Thu,) studied this question.