Implementing discrete element method for large-scale simulation of particles on multiple GPUs

Numerical study of granular flow in a slit funnel with a novel structure to avoid particle clogging

To solve the problem of particle clogging in slit funnels and to obtain a stable discharge flow rate, we proposed a new funnel structure, namely the slit baffle funnel. We conducted a systematic investigation using the discrete element method (DEM) to study the effects of funnel half-angle θ, outlet width W, and baffle height H on flow rate and flow pattern. We found that the proposed structure could effectively avoid particle clogging and guarantee a continuous and stable flow rate with small outlet width. Under the condition of H >3 d, a bigger flow rate was obtained at a smaller funnel half-angle. This new funnel structure could be applied to solve clogging problems associated with granular matter in the slit geometry in mining, agriculture, food, and pharmaceuticals.

Flow-induced surface crystallization of granular particles in cylindrical confinement

An interesting phenomenon that a layer of crystallized shell formed at the container wall during an orifice flow in a cylinder is observed experimentally and is investigated in DEM simulation. Different from shear or vibration driven granular crystallization, our simulation shows during the flow the shell layer is formed spontaneously from stagnant zone at the base and grows at a constant rate to the top with no external drive. Roughness of the shell surface is defined as a standard deviation of the surface height and its development is found to disobey existed growth models. The growth rate of the shell is found linearly proportional to the flow rate. This shell is static and served as a rough wall in an orifice flow with frictionless sidewall, which changes the flow profiles and its stress properties, and in turn guarantees a constant flow rate.

Inclined granular flow in a narrow chute

In this paper we presents a detailed description of granular flow down a flat, narrow chute using discrete element method simulations, with emphasis on the influence of sidewalls on the flow. The overall phase diagram is provided and it is found that there are four flow regimes (no flow, bulk flow, surface flow, and gas flow). The H̃_stop curve is very complicated and quite different from that in the case without sidewalls. The effective friction coefficient [Formula: see text] increases with pile height H̃ and a surface flow occurs when the inclination angle [Formula: see text] exceeds a critical value. The profile of [Formula: see text] shows that the [Formula: see text] rheology is valid in boundary layers. Furthermore, [Formula: see text] increases with the velocity of particles and there is a saturation to a nonzero value in static heap. For small H̃, the static heap vanishes and there is a bulk flow. A similarity between basal particles and sidewall particles indicates a universal role of the boundaries. In this bulk flow, there is a transition of the velocity profile with wall friction [Formula: see text]. When [Formula: see text] is large, the velocity is linear and decreases with increasing height. With small [Formula: see text], the velocity is nonlinear and the flow rate is roughly proportional to H̃^3/2.

Preliminary research on flow rate and free surface of the accelerator driven subcritical system gravity-driven dense granular-flow target

A spallation target is one of the three core parts of the accelerator driven subcritical system (ADS), which has already been investigated for decades. Recently, a gravity-driven Dense Granular-flow Target (DGT) is proposed, which consists of a cylindrical hopper and an internal coaxial cylindrical beam pipe. The research on the flow rate and free surface are important for the design of the target whether in Heavy Liquid Metal (HLM) targets or the DGT. In this paper, the relations of flow rate and the geometry of the DGT are investigated. Simulations based on the discrete element method (DEM) implementing on Graphics Processing Units (GPUs) and experiments are both performed. It is found that the existence of an internal pipe doesn’t influence the flow rate when the distance from the bottom of the pipe to orifice is large enough even in a larger system. Meanwhile, snapshots of the free surface formed just below the beam pipe are given. It is observed that the free surface is stable over time. The entire research is meaningful for the design of DGT.