Comparing non-local granular fluid continuum models for silo discharge: Toward clogging prediction

Abstract

Non-local constitutive theories have received increasing attention in continuum descriptions of granular flows. However, these models have not been systematically compared for silo discharge within a unified numerical framework. We address this gap with two-dimensional finite-volume method (FVM) simulations of silo discharge using the Basilisk platform. We first validate our FVM implementation of the dynamic non-local granular fluidity (NGF) model against the material point method results of Dunatunga & Kamrin (J. Fluid Mech., 2022, 940, A14), obtaining quantitative agreement. Second, we relate the discharge rate $Q$ to the outlet-to-particle size ratio $D/d$ and the non-local amplitude $A$. From the simulated $Q$, we then evaluate the clogging probability $J(D/d, A)$ within the probabilistic framework of Janda et al. (Europhys. Lett. 84 (4), 44002). The predicted $J$ decays exponentially with $D/d$, consistent with the experimental trend of Janda et al. (Europhys. Lett. 84 (4), 44002). Rather than directly predicting flow arrest, our approach captures the continuous probabilistic transition. Finally, within the same numerical framework and using identical values of $A$, we compare several non-local constitutive models, including several linearised variants that we derive. Almost all models predict a reduction in the discharge rate with increasing $A$, yet significant quantitative differences are observed among the models. The results are further classified into groups according to their predicted flow behaviour, revealing close correspondences among certain formulations. Notably, using the non-local amplitudes reported in the literature [Bouzid et al. (Phys. Rev. Lett. 111, 238301), Henann & Kamrin (Proc. Natl Acad. Sci. USA 110(17))] yields near-zero discharge rates. Ill-posed issues are discussed. The implementation of all models is open-sourced and computationally efficient.