#### Abstract

The granular flow is one of the principal issues for the design of pebble bed reactors. Particularly, the clogging phenomenon raises an important issue for pebble bed reactors. In this paper, we conduct experiments and discrete particle simulation of two-dimensional discharge granular flow from a conical hopper, to study the effect of the particle bed height and hopper angle on the clogging phenomenon. In general, the clogging probability increases with height and starts to saturate when is larger than a critical value. The experimental result trends are supported by discrete simulations. To understand the underlying physical mechanism, we conduct discrete particle simulations for various values, focusing on the following parameters: the statistical averaging of the volume fraction, velocity, and contact pressure of particles near the aperture during the discharge. We found that, among all relevant variables, the contact pressure of particles is the main cause of the increasement of *J* when increases. An exponential law between the pebble bed and clogging probability *J* has been established based on these observations and Janssen model. As for hopper angle , *J* shows an almost constant behavior for any rise in followed by a sudden regression at . Surprisingly, the effect of is most obvious for intermediate values of , where we observe a sharp increasement of clogging probability. The same trend is observed in the two-dimensional discrete simulation results.

#### 1. Introduction

When a group of discrete particles move downwards inside the hopper, an arch may be formed if the size of the orifice is just a few times larger than the particles, leading to the development of the clogging phenomenon. Clogging occurs just above hopper’s aperture, and being a state of static equilibrium, it blocks the overall flow of particles making the system impractical. The phenomenon is not as simple as it appears to be, with its scope quite conceivable ranging from people escaping the room [1, 2] and blockage of industrial granular flow [3] to traffic flow [4]. Clogging also may happen in the core area of pebble bed high-temperature gas-cooled reactor (HTGR). HTGR contains a hopper-shaped core that is being continuously filled with graphite spherical pebbles [5]. This probable formation of arch, shaped in a convex-like structure, is capable enough to jam the system of pebble flow. Once clogging occurs, the refueling stops, initiating stern consequences.

Despite there exist many previous studies, the mechanisms of the clogging process remain unclear. Most previous studies focus on the influence of orifice [6–10]. From the pioneering work of To et al. in 2001 [6], we know that magnitude of clogging probability changes greatly with the outlet size. Furthermore, the existence of a critical outlet size above which the flow may never get jammed remains debatable [11]. Subsequently, it has been experimentally and numerically shown that some methods can reduce the probability of clogging by, for example, putting an obstacle above the exit [12, 13] and applying a helical texture to the pipes’ inner wall [14]. Other variables that have been studied are orifice geometry [15, 16], particle shape [17], particle polydispersity [18], driving force [19, 20], and width of silo [21].

To the best of our knowledge, the role of height of particle bed in the clogging process has not been studied in previous work. The flow pattern of particles before clogging has the resemblance of the flowing particles discharged from a hopper. For the case of flowing discharge flow from a hopper, it has been observed empirically that when the height of particle bed , where *L* is the hopper width, the flow rate of particles will become independent of the material height *h*. Such observation accounts for the fact that an hourglass is able to function properly [22]. Probably due to this reason, this parameter has not been studied in previous works. At the exit of the fuel cycle of a typical Chinese 10 MW high-temperature gas-cooled test reactor (HTR-10), the height of the pebble bed does not meet the condition *h* > 2.5 *L*, and the exit of the fuel cycle channel resembles a conical hopper with angle below [23, 24]. There exist few works that focus on the influence of hopper angle on granular clogging for high particle bed [6, 25]. The existing work reveals that clogging probability decreases by increasing the hopper angle. The authors noticed a sharp regression in clogging probability for a value of tilt angle around [6, 25]. The influence of low hopper angle value () has been ignored. It is reasonable to investigate the probability of the clogging phenomenon under small values of these two parameters. In the present work, we carried out an investigation of the effect of both, the particle bed height and angle of particle , on the clogging phenomenon. First, the experimental procedure and discrete simulation are described. The simulation is performed on LMGC90 software, which is based on the NSCD (nonsmooth contact dynamics method) approach. Second, the influence of particle bed height on the clogging problem has been studied; the experimental result trends are supported by discrete simulations; in order to find the cause for the increment of when increases, we conduct discrete particle simulations for various values, focusing on the statistical averaging of the volume fraction, velocity, and contact pressure of the particles near the aperture during the discharge. An exponential law between the pebble bed and clogging probability has been established based on these observations and Janssen model. Finally, the influence of angle of hopper on the clogging problem has been discussed.

#### 2. Materials and Methods

##### 2.1. Experimental Procedure

The experimental set-up consists of a two-dimensional conical silo with the height of , width of 90 mm, and thickness of 3.5 mm. Firm support has been given to the hopper by an aluminum base plate of 4.5 mm in thickness. Front wall comes up with a 4 mm thick sheet, made of perspex to allow visualization. Hopper side walls have been fabricated with aluminum rods of rectangular cross section firmly bolted to the base plate. The base of the hopper is made up of two aluminum rods whose thickness is 3.5 mm just to match the thickness of lateral walls and contains a sleeve inside them so that they can be moved with the help of screws along their length, providing us with the possibility of varying the angle at will. The rectangular outlet is at the center of the bottom of the hopper, with length *D* = 18 mm and thickness *W*, which is the same as the thickness of the hopper. Note that the thickness of rectangular outlet *D* is four times larger than the particle diameter *d*. The apparatus is grounded to prevent the buildup of static electricity.

After filling a monolayer of monodisperse stainless steel disks of thickness along with a diameter of , this weak polydispersity of particles is introduced to avoid the effect of crystallization. With no external drifting force on particles to ensure their discharge course, gravitational force would be the only deciding factor to traverse them out of the hopper. In this paper, the experiment is performed to validate the results obtained by simulation. Following To et al. [6], we have only focused on the clogging probability for experiments, and other information will be obtained by discrete particle simulations subsequently. The clogging probability is defined as the ratio of clogging events divided by the number of total trial . is 200. We do not refill the container with particles. The main parameters in this experiment are the particle bed height , which is the perpendicular distance from the aperture to the level displaced by the particles inside the hopper and the angle of the hopper bottom (see Figure 1). The dimensions and parameters of test facility are presented in Table 1.

##### 2.2. Simulated System

###### 2.2.1. Discrete Particle Simulation Model and Settings

We performed the simulations by using the open source software LMGC90, which adopts the nonsmooth contact dynamics method (NSCD); this method is based on the concept of rigid spheres and on the implementation of Coulomb’s law without regularization. The particles are treated inelastic. Frictional coefficients that have employed to enable contact dissipation are being and , with the particles and the walls, respectively. A weak polydispersity is set to be to avoid the crystallization effects. We set the time step to . More simulation details can be found in Zhou et al. [26] and Radjai et al. [27]. We have developed a two-dimensional conical hopper of angle , orifice size , and width (see Figure 2). The orifice is located at the bottom center. We impose to ensure that silo width does not influence the flow appreciably. Particles with diameter are arranged into granular column by random deposition to diminish the gravitational potential in a closed hopper. First, we have studied the clogging probability by varying the particle bed height and angle . The number of time steps is set long enough to ensure that all particles must get discharged out of the hopper when there would be no occurrence of clogging. 1000 trails have been performed for each case. Second, in order to study the underlying physical mechanism, the statistical averaging of the volume fraction, velocity, and contact pressure of the particles near the aperture during the discharge (before clogging occurs) is measured for various *h*. The computational domain is periodic in the vertical direction to keep constant the number of particles. The vertical boundary of the computational domain has been set to be above and below the hopper, where is the maximum diameter of the particle. The total number of time step is set to , to ensure the accuracy of the statistical estimation of the average value of volume fraction, velocity, and stress [26]. From the discrete numerical representation of LMGC90, these averages are then computed to examine the processes at the scale of the outlet that could rule the clogging phenomenon at the hopper’s scale. The geometry parameters of silo without and with periodic boundary condition are reported in Table 2. Simulation details are reported in Tables 3 and 4 . Specifically, the number of particles and the number of time step for the case without boundary condition are reported in Table 3; the number of particles and the number of time step to calculate the statistical averaging for the case with boundary condition are reported in Table 4.

###### 2.2.2. Postprocessing Method

In this paper, we present an extensive numerical study of the micromechanical properties displayed by granular media flowing through an orifice at the bottom of the hopper. We describe the volume fraction, velocity, and stress of granular media over a confined region near hopper’s aperture. Following Zhou et al. [26], we measured the profile of the particle volume fraction , the horizontal velocity , and the vertical velocity at the orifice. The measured region is near the orifice, which is the dashed red box in Figure 3, with a height *d*. It is the average of the statistics performed over the set of snapshots of a computation during a period of discharge without clogging. For each snapshot, one considers the statistics to be homogeneous over a rectangular domain centered on the measurement location *x* (see blue box in Figure 3). The width of this rectangle is , which is small with regard to the typical variation of the profiles. The particle volume fraction, , is the average over the integration volume (in space and time) of the indicator function, whose value is over the spatial extent of particles and otherwise. According to the ergodicity theorem [28], we can derive the relative accuracy of the statistical estimation of the average value as follows [29]:where is the average outlet velocity. The velocity profile at the outlet of a particle population is obtained as the ensemble average of the individual particle velocity weighted by their volume intersection with the integration domain. It is therefore the space and time average of the product of the indicator function by the individual velocity divided by . Since standard deviation of the velocity distribution is small, the convergence of statistics toward the average value is more rapid than for . The detail of the postprocessing technique can be found in the study by Zhou et al. [26]. By using an averaging process popularized by Goldhirsch et al. [30, 31], we compute the stress tensor from the interparticle forces in the granular medium. The measured region is the yellow box in Figure 2. Let us consider at time *t* a set of *N* particles indexed by , having the mass , whose centers of mass are and which move with a speed . The macroscopic stress field can be decomposed into a contact stress field related with particle-particle interactions and a kinetic stress field which is induced by the fluctuating motion of grains, which carry their momentum from one place to another. By choosing a common step function, the definition of the stress tensor is equivalent towhere is the fluctuating velocity of particle . The second sum is performed over the contacts *c* inside the volume . Here, is the interaction force modulus between the two particles in contact in *c*. The vector is called the branch vector and is equal to the vector connecting the centers of two particles and forming the contact if and belong to the volume but is equal to the portion of this vector included in *V* if one particle is outside the volume *V*. The trace of the stress tensor is defined as the contact pressure. Figure 4 shows the contact pressure of a typical case with , which is compatible with the contact pressure field obtained in Rubio-Largo’s work [32].

#### 3. Results and Discussion

##### 3.1. Influence of Particle Bed Height

We studied the effects of height on clogging probability . Figure 5 plots with respect to under various hopper angles , where Figure 5(a) and Figure 5(b) show the experiment results and simulation results, respectively.

**(a)**

**(b)**

From Figure 5(a), we can see that the clogging probability increases with height and starts to saturate when is larger than . With , decreases sharply for any particle bed height . For smaller and greater heights , apart from that clogging decreases sharply at , no appreciable difference in clogging can be seen by transforming the inclination of hopper walls. The influence of hopper angle will be further discussed in the next section. The same trend has also been observed in simulation results (see Figure 5(b)); the clogging probability increases with height and starts to saturate when is larger than .

In order to explore the underlying physical mechanism of the influence of particle bed , we have performed the discrete element simulation with periodic condition. The parameters for the simulation are also found in Tables 1 and 3. First, following Zhou et al. [26], we investigated the volume fraction and velocity profiles near the outlet. Figure 6(a) shows the measured particle volume fraction at the outlet, as a function of the horizontal position *x* normalized by the diameter of the outlet for various heights . We can observe that as the height increases, the data points tend to overlap with each other. Average values of volume fraction are plotted versus the particle bed height normalized by the hopper width in Figure 6(b). With different values, the average values of volume fraction remain constant. The error bars represent the relative accuracy given by equation (1) which is around . Regardless of any effect induced by pebble bed height , we observed a relatively constant behavior in volume fraction with a minute surge at elevated values. The difference between the maximum and minimum values is percent of the minimum value. In the case of flowing particles without clogging, the variation of volume fraction at the center of the outlet depends on the outlet size and on the particle diameter, exhibiting a dilatancy for small aperture in order to maintain the flow of the material [26, 33]. From the results displayed in Figure 6, it can be concluded that the change of height of particles bed does not influence the dilatancy at the outlet. Namely, the volume fraction does not contribute to the variation of clogging probability when varies.

**(a)**

**(b)**

In the same way, the profiles of the horizontal component of velocity, , and the profiles of the vertical component of the velocity, , for various column heights are plotted, respectively, in Figures 7(a) and 7(b). The horizontal velocity depends linearly on , and it increases.

**(a)**

**(b)**

**(c)**

As we move away on either side from the central axis of the hopper at the outlet, these profiles show that the particle flow has a horizontal component towards the center. These values of horizontal velocity are independent of the pebble bed height . Figure 7(b) reveals that the vertical velocity profiles are quasiparabolic and depend on the particle bed height. The evolution of the average value of vertical velocity versus the particle bed height normalized by hopper width is plotted in Figure 7(c). We can observe that, as the pebble bed heights increase, the vertical velocities would also rise. It has been observed experimentally that when the driving force increases, the clogging probability decreases. Arevalo et al. [19] found that the physical magnitude behind clogging reduction by increasing gravity or particle density is the averaged kinetic energy per particle. Assuming that the kinetic energy of the particles is related to the total kinetic energy, the higher the kinetic energy, the longer it takes to the system to stabilize, and during this time is when the particles can keep flowing, it seems reasonable that increasing the kinetic energy decreases the clogging probability. In our case, when increases, the kinetic energy increases but the clogging probability increases; this suggests that the physical magnitude behind this clogging increment by increasing the particle bed height is not the kinetic energy.

Next, we investigated the results of contact pressure at a region close to the outlet (see the yellow box in Figure 2); at the vertical central line of the region mentioned above is plotted against vertical position normalized by the diameter of the outlet for various heights (see Figure 8(a)). We can observe that the contact pressure increases with the vertical position and saturates at a certain position. The contact pressure at the position is plotted against in Figure 8(b); one can see that it depends strongly on the particle bed height and starts to saturate when is larger than . We can conclude that, among all relevant variables, the contact pressure of particles is the main cause of the increasement of when increases.

**(a)**

**(b)**

In 1895, Janssen proposed a model in which they measured the pressure at the bottom of a silo filled with corn [34]. They noticed that as long as the corn keeps being filled, pressure saturates, which is different from the case of liquids. In light of Janssen’s model, we will establish a naive relationship between the clogging probability and particle bed height . First of all, we will study the contact pressure on the upper edge of the square shown as the yellow box in Figure 2 on the position . According to Janssen’s model, the vertical stress at is given bywhere and are constants which depend on the geometry and material of particles and hopper. According to Janssen’s assumption, the horizontal normal stress is proportional to the vertical normal stress, and the contact pressure can be written as follows:where , , and are fitting parameters. We adjusted the simulation data of the contact pressure near the orifice by equation (4) by using the least-squares method (see black dotted line in Figure 8(b), where the contact pressure at location is plotted versus the particle bed height). The agreement is found to be fairly good, which validates equation (4).

Next, we hypothesize that the clogging probability depends linearly with the contact pressure of the particles when the particles are static before opening the orifice. We can deduce that the relationship between clogging probability and pebble bed height *h* is written as follows:

represents the clogging probability when tends towards infinity. This value depends on other parameters (like orifice size) of the system, and there exist some empirical models of this value [11]. and are the fitting parameters. We plotted versus with the flat-bottom hopper () for experiment and simulation in Figure 9(a) and 9(b), respectively. We adjusted the experimental and simulation data by equation (5) by using the least-squares method (see black dashed line); the fitting parameters are and for experiment and and for simulation. We can see that in two cases, is an order of magnitude larger than . When tends to infinity, we obtain clogging probability in the experiments and in simulation. We can see that the particles are more likely to jam in simulation; this is perhaps due to the reason that the experiment is quasi-2D which differs from a 2D case.

**(a)**

**(b)**

##### 3.2. Influence of Hopper Angle

As we can observe from Figure 5 that clogging probability depends on hopper angle, here we replot Figure 5 to better illustrate the influence of hopper angle. Figure 10 plots with respect to under various particle beds , where Figures 10(a) and 10(b) show the experiment results and simulation results, respectively. From Figure 10(a), we can see that shows an almost constant behavior for any rise in followed by a sudden regression for . However, for intermediate values of , we have noticed a sharp increasement in clogging till . This trend has been perceived mostly at . At a particle bed height , has proved to be affected by at its maximum. A clogging maximum has been attained at that has not been affected by particle heights.

**(a)**

**(b)**

The same trend has been observed for simulation results in Figure 10(b). increases with for intermediate values of when a maximum value has been attained at and starts to decrease at . It has been revealed that a convincing correlation exists between the hopper angle and the nature of flow of particles. For lower angles, the nature of flow is of funnel type, at which particles along the central axis of the hopper move faster than the ones closer to the walls since frictional effects dominate so that only a portion of particles can move while others form a stagnant zone [35]. Rose and Tanaka [36] found that due to the presence of a stagnant zone near the opening, the flow rate in a hopper does not depend on if is less than a critical value related to the angle of stagnant zone. Before clogging occurs, the flow of particles in our clogging experiment has the resemblance of the discharge of flow from hopper. It is possible that in our experiment, is between and so that the clogging probability is almost constant for . For hopper angle , the clogging probability decreases sharply; this is reasonable since at the extreme value when , clogging should not occur when . However, the surge of clogging probability for intermediate values *h* is somewhat unexpected. It is probably due to the augmentation of contact pressure. Currently simulations are being carried out for further study of this effect.

#### 4. Conclusions

In this paper, discrete numerical simulations and experiments are conducted to study the influence of container’s geometry on the clogging phenomenon. The height of the pebble bed and the angle of the conical hopper are our varied parameters. The volume fraction, the velocity, and the contact pressure profiles through the aperture were measured to study the underlying physical mechanism. According to the results, some conclusions were drawn as follows:(1)The clogging probability increases with height and starts to saturate when is larger than a certain value, suggesting that flow becomes more susceptible to clog if we increase . This behavior is consistent for every hopper angle.(2)Based on the results of volume fraction near the orifice, we observe a weak variation of volume fraction when varying particle bed height , suggesting that the change of height of particles bed does not influence the dilatancy at the outlet. We can conclude that the volume fraction at the outlet does not contribute to the augmentation of clogging probability when increases.(3)The results suggest that the physical magnitude behind this clogging increment by increasing the particle bed height is not the kinetic energy.(4)The contact pressure near the orifice depends strongly on the particle bed height and starts to saturate when is larger than a certain value. It thus has an imperative effect on the probability of clogging. An exponential law between the pebble bed and clogging probability has been established based on these observations and Janssen model.

#### Data Availability

The data used to support the findings of this study are available from the corresponding authors upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This project was supported by the National Natural Science Foundation of China (Grant no. 11802094).