The effect of particle aspect ratio and surface geometry on granular flows is assessed by performing numerical simulations of rod-like particles in simple shear flows using the discrete element method (DEM). The effect of particle surface geometry is explored by adopting two types of particles: glued-spheres particles and true cylindrical particles. The particle aspect ratio varies from one to six. Compared to frictionless spherical particles, smaller stresses are obtained for the glued-spheres and cylindrical particle systems in dilute and moderately dense flows due to the loss of translational energy, which is partially converted to rotational energy, for the non-spherical particles. For dilute granular flows of non-spherical particles, stresses are primarily affected by the particle aspect ratio rather than the surface geometry. As the particle aspect ratio increases, the effective particle projected area in the plane perpendicular to the flow direction increases, so that the probability of the occurrence of the particle collisions increases, leading to a reduction in particle velocity fluctuation and therefore a decrease in the stresses. Hence, a simple modification is made to the kinetic theory for granular flows to describe the stress tensors for dilute flows of non-spherical particles by incorporating a normalized effective particle projected area to account for the effect of particle collision probability. For dense granular flows, the stresses depend on both the particle aspect ratio and the surface geometry. Sharp stress increases at high solid volume fractions are observed for the glued-spheres particles with large aspect ratios due to the bumpy surfaces, which impede the flow. However, smaller stresses are obtained for the true cylindrical particles with large aspect ratios at high solid volume fractions. This trend is attributed to the combined effects of the smooth particle surfaces and the particle alignments such that the major/long axes of particles are aligned in the flow direction. In addition, the apparent friction coefficient, defined as the ratio of shear to normal stresses, is found to decrease as the particle aspect ratio increases and/or the particle surface becomes smoother at high solid volume fractions.

This review article focuses on the modeling of complex granular flows employing the discrete element method (DEM) approach. The specific topic discussed is the application of DEM models for the study of the flow behavior of nonspherical, flexible, or cohesive particles, including particle breakage. The major sources of particle cohesion—liquid induced, electrostatics, van der Waals forces—and their implementation into DEM simulations are covered. These aspects of particle flow are of great importance in practical applications and hence are the significant foci of research at the forefront of current DEM modeling efforts. For example, DEM simulations of nonspherical grains can provide particle stress information needed to develop constitutive models for continuum-based simulations of large-scale industrial processes.

Granular shear flows of flat disks and elongated rods are simulated using the Discrete Element Method. The effects of particle shape, interparticle friction, coefficient of restitution, and Young's modulus on the flow behavior and solid phase stresses have been investigated. Without friction, the stresses decrease as the particles become flatter or more elongated due to the effect of particle shape on the motion and interaction of particles. In dense flows, the particles tend to have their largest dimension aligned in the flow direction and their smallest dimension aligned in the velocity gradient direction, such that the contacts between the particles are reduced. The particle alignment is more significant for flatter disks and more elongated rods. The interparticle friction has a crucial impact on the flow pattern, particle alignment, and stress. Unlike in the smooth layer flows with frictionless particles, frictional particles are entangled into large masses which rotate like solid bodies under shear. In dense flows with friction, a sharp stress increase is observed with a small increase in the solid volume fraction, and a space-spanning network of force chains is rapidly formed with the increase in stress. The stress surge can occur at a lower solid volume fraction for the flatter and more elongated particles. The particle Young's modulus has a negligible effect on dilute and moderately dense flows. However, in dense flows where the space-spanning network of force chains is formed, the stress depends strongly on the particle Young's modulus. In shear flows of non-spherical particles, the stress tensor is found to be symmetric, but anisotropic with the normal component in the flow direction greater than the other two normal components. The granular temperature for the non-spherical particle systems consists of translational and rotational temperatures. The translational temperature is not equally partitioned in the three directions with the component in the flow direction greater than the other two. The rotational temperature is less than the translational temperature at low solid volume fractions, but may become greater than the translational temperature at high solid volume fractions.

Bulk (Ga1-x Zn x )(N1-x O x ) as a photocatalyst has received increasing attention as a potential solution for the energy shortage challenge; however, its catalytic performance is highly limited by its bulk form. To improve the photochemical potential, the nanoscale form of this multiple-metal oxynitrides is desirable. In this work, a new type of (Ga1-x Zn x )(N1-x O x ) nanostructure is obtained. Its composition can tuned to the full range (0.18 < x < 0.95). The (Ga1-x Zn x )(N1-x O x ) nanostructure exhibits excellent photocatalytic activity for overall water splitting, and the highest quantum efficiency of (Ga1-x Zn x )(N1-x O x ) is as high as 17.3% under visible light irradiation. Using this new type of (Ga1-x Zn x )(N1-x O x ) nanostructure, the narrowing of the bandgap for (Ga1-x Zn x )(N1-x O x ) is not only due to an increase in the valence band maximum, but it is also related to a decrease in the conduction band minimum.

Abstract