Hui Zee
Then
,
Teruyo
Sekiguchi
and
Ko
Okumura
*
Physics Department and Soft Matter Center, Ochanomizu University, Tokoyo 112-8610, Japan. E-mail: okumura@phys.ocha.ac.jp
First published on 6th August 2020
The rising motion of an obstacle in a vibrated granular medium is a classic problem of granular segregation, and called the Brazil nut (BN) effect. Identification of the controlling vibration parameters of the effect is a long-standing problem. The simple possibility that the BN effect can be characterized solely by vibration velocity has recently been pointed out. The issue has become controversial over the long history of research, with only a few systems providing evidence for this simple possibility. Here, we investigate the rising motion of an obstacle in a vertically positioned one-layer granular bed under continuous vibrations. We find the rising motion is composed of two distinct regimes, and the first and second regimes are both governed, in terms of vibration parameters, solely by the vibration velocity. We further demonstrate simple scaling laws that well describe the two regimes. Our results support the emergent simple possibility for the controlling parameters of the BN effect and suggest that this feature could be universal. We propose two possible mechanisms of convection and arch effect for the two distinct regimes and demonstrate that these mechanisms explain the scaling laws followed by our experimental data.
However, as for the controlling vibration parameters of the BN effects, our understanding has become controversial. It was believed that characterization of the convection-driven rising motion requires at least two vibration parameters: the acceleration and frequency.8,9 However, several years ago, it was clearly shown that the convection-driven rising motion is characterized, in terms of vibration parameters, solely by the vibration velocity, in a wide range of experimental parameters.10 Subsequently, a detailed study on the granular convection induced by vibration revealed that the convection velocity is well characterized by the vibration velocity,11 supporting the recent study.10 However, despite a long history of research, the vibration-velocity governed BN effect has been reported only in a few cases,10,12 and the problem of the controlling vibration parameters has become controversial (e.g., well-known MRI studies8,13 did not demonstrate that their data can be characterized solely by vibration velocity). Therefore, providing experimental data that can be characterized solely by vibration velocity in different systems is an emergent important issue for settling this long-standing issue.
In order to elucidate the controlling vibration parameters in the BN effect, we investigate the rising motion in a one-layer granular bed under continuous vibration. As a result, we find that the dynamics are divided into first and second regimes, and both are well characterized solely by the vibration velocity through simple scaling laws. In addition, we propose possible mechanisms of convection and arch effect for the two regimes, and demonstrate that these mechanisms explain the scaling laws followed by the experimental data.
We used a cell with side walls slightly down-facing with an angle of about 5° (specified below) throughout this study. A slight side-wall angle was necessary to observe stable rising motion of the disk. When the side walls were vertical, the rising motion was less smooth and less reproducible. Dependence on the angle may originate from a weak friction on the side walls (their surfaces are smooth) and a relatively weak polydispersity of our beads (the bead diameters were approximately in the range 0.9–1.15 mm), which tends to create lattice structures. With a slight down-facing angle, lattice structures tend to be broken and the movement of beads near the side walls are restrained as in the case when friction on the side wall is strong. Note that the importance of friction on the vertical side wall for convection has been pointed out, and it has also been demonstrated that the introduction of an angle in the opposite up-facing direction can reverse the direction of convection.5 However, in the present case, the angle was smaller and introduced to conduct a reproducible study of the BN effect. In this sense, the slight side-wall angle contributes to enhance the effect of polydispersity and friction near the side walls.
A cell of thickness t = 1.2 mm was made from two transparent acrylic plates of thickness 3 mm, separated with spacers of thickness t. The cell height H was fixed to 140 mm. The top and bottom widths of the cell, W1 and W2 (see Fig. 1a), were slightly different with a finite angle α defined by tanα = (W2 − W1)/(2H). The width of the cell W was defined as W = (W1 + W2)/2 for convenience. The cell width W was either 55, 80, 85, or 105 mm. In the case of W = 85 mm, W1 and W2 were set to W1 = 80 and W2 = 90 mm, while W1 and W2 for different W were determined such that the angle α was the same with that for W = 85 mm, i.e., α was fixed to a value α = arctan[(90–80)/(2 × 140)] (∼5°).
The cell contained a metal disk (diameter D = 12, 15, or 18 mm, with thickness 1 mm), which plays the role of an obstacle or an intruder, and one layer of small balls of aluminum oxide with average diameter d = 1 mm (AL-9, AS ONE Corp.), which play the role of grains or a granular medium. The balls filled the cell with the obstacle to a depth H0 = 100 mm (or 80 nm in a few cases). The density of disks (7.7–7.9 g cm−3) are larger than that of the small alumina balls (3.95 g cm−3).
The cell containing the obstacle and grains was mounted vertically on a vibration generator system (m060/MA1-CE, IMV Corp.), which was controlled by a multi-function generator (WF1948, NF Corp.) and monitored by an acceleration meter (VM-1970, IMV Corp.). Digital images were obtained with a CCD camera (STC-MB33USB, SENTECH Co., Ltd), which was synchronized with the vibration system through the function generator.
The vibration generation system can produce continuous sinusoidal waves with control of the angular frequency ω and amplitude A. The sinusoidal wave can be characterized by the vibration velocity v = ωA, acceleration a = ω2A, and frequency f = ω/(2π).
![]() | (1) |
![]() | (2) |
![]() | ||
Fig. 2 (a and b) The rising velocity in the first regime U1 and that in the second regime U2 as a function of the vibration acceleration a for different vibration frequencies, 12.5 Hz (squares), 20 Hz (triangles), 30 Hz (circles) and 40 Hz (diamonds). Data for (D,W) = (12,55), (12,85), (15,85), (18,85) and (18,105) are represented by black, red, yellow, blue, and magenta symbols, respectively. (c) All the data in (a) replotted with renormalized axes. The line with slope 2 represents eqn (1) with k1 = 0.88 ± 0.022. (d) All the data in (b) replotted with renormalized axes. The line with slope 1 represents eqn (2) with k2 = 0.85 ± 0.016. |
The following information could be gleaned from the values of k1 and k2. (I) The fact that k1 and k2 are both close to one confirms that the orders of magnitude predicted by eqn (1) and (2) are consistent with our experimental data. (II) The fact that the standard deviations for k1 and k2 thus obtained are only a few per cent (they are 2.51 and 1.85%, respectively) suggests that our experimental data and analyses are of high quality.
Details for obtaining the values and error bars for Fig. 2c and d are explained as follows. We prepared n sets of data under the same experimental conditions, i.e., for a set of (D,W,f,v). For each set of data, we first determined the second region exploiting its clear linearity with a slope (corresponding to U2), and the data in the remaining region was fit by a straight line with another slope (corresponding to U1). This second fit was conducted by selecting a linear region as wide as possible with the criterion that the coefficient of determination became below 0.98 (this coefficient was determined by where zP and zA are the predicted and average values, respectively, for m data of zi). This well-defined method to determine U1 was used based on the following observations. (1) Rising motion was unstable near the starting point with less reproducibility. (2) However, before the second region, which is clearly linear, another linear region tended to appear and its slope was reproducible although its width was not. These features are visible in Fig. 1b if we compare the slopes and the experimental data. Note that k1 and k2 used in Fig. 1b are the same values obtained from Fig. 2c and d. In this way, we obtained n different sets of (U1, U2) for a single set of (D,W,f,v). If n was larger than 10 we employed the standard deviation as the error bar. Otherwise, the highest and lowest positions of the error bar were determined by the maximum and minimum values, respectively. Some data in Fig. 2 and the data in Fig. 1b were obtained at H0 = 80 mm but the U1 and U2 obtained at H0 = 100 and 80 mm were indistinguishable, as demonstrated by a good agreement between the data and the solid and dashed lines in Fig. 1b.
![]() | ||
Fig. 3 (a) Illustration of two arches repeatedly formed at the bottom of a layer during vibration. The width of the arch is determined by the distance between the side walls of the cell. This illustration superimposes two cases, a wider cell and a narrower cell at the same vibration velocity when the arch becomes most significant (such a moment, corresponding to the second rightmost snapshot in Fig. 3b below, repeatedly appears with a given frequency). The height at the center is the same for the two cells because it is set by v2/g as explained in the text. (b) Sequential snapshots, separated by 10 ms, of a cell of W = 85 mm under vibration at f = 12.5 Hz and a = 30 mm s−2 with an intruder of D = 15 mm. The arch-shaped void is created at the bottom of the granular layer, which is most visible in the photograph second from the right. (c) Void filling mechanism suggested in previous studies: black particles continuously fill up the void created beneath the intruder during shaking, which results in the rise of the intruder. Figure (c) is reproduced with permission from ref. 7, Copyright, American Physical Society, 1992. |
For simplicity, we consider the case a ≫ g to gain physical insight into this sequential motion (in fact, the ascent can be observed only when a > g, as also reported in ref. 10 and 14). When a sinusoidal wave is applied to the system, the normal force acting on the system (total mass M) of one grain layer and the intruder from the cell bottom is given by N = −Mg − Mα, as long as the system holds contact with the cell bottom, if we neglect the friction effect near the side walls. Here, α = −acos
ωt is the acceleration corresponding to the movement of the position of the base plate of the cell described by z = A
cos
ωt. This means, if a ≫ g, at the moment when the cell moves upwards from the central position of vibration, z = 0, the normal force N becomes zero and the system starts a parabolic motion under gravity with the initial speed comparable to v = ωA. The maximum height of the parabolic motion scales as h with v2 ≃ gh (this relation is exact in Newtonian mechanics, in the absence of side walls with the constraint that all the particles are on the vertical plane; thus the relation tends to hold in the central region if the wall distance is large and friction with the front and back walls is small, as in the present case). Furthermore, if the maximum height h is considerably larger than the vibration amplitude A (this implies ω is large in the light of the relation a = ω2A), the maximum height of the layer bottom relative to the cell bottom is comparable to h, which is the same under a fixed v. In summary, the maximum height tends to scale as h ≃ v2/g.
However, small beads close to the side walls cannot move (relative to the cell) because of the friction effect near the side walls, while small beads around the center of the cell repeatedly go up with the intruder under vibration till the maximum height comparable to h ≃ v2/g (relative to the cell bottom) as estimated above. This leads to the formation of a velocity gradient in the granular layer in the direction of cell width, which results in an arch-shaped void at the bottom of the layer as shown by the snapshots in Fig. 3b. During the formation of an arch-shaped void, some beads close to the bottom of the layer enter into the void due to the velocity gradient, which results in a shear force. As a result, the height of the intruder relative to the cell bottom when it is next pushed to the cell bottom becomes higher than before. This is the mechanism of ascent by void filling.
This mechanism of the ascent by void filling is consistent with the dependence of the ascent velocity U1 on the cell width W and on the vibration velocity v, predicted in eqn (1). As long as a ≫ g, for a fixed v, the maximum height of the arch tends to be the same (comparable to h given above), meaning that the velocity gradient increases as W decreases. Accordingly, we expect that U1 increases as W decreases, which qualitatively explains why U1 scales with 1/W in eqn (1). For a fixed W, since h tends to scale with v2, if the void were completely filled, U1 would scale with ωv2 by noting that the period of vibration scales as 1/ω. However, the void filling is much less effective, as observed in Fig. 3b, because the filling is induced by a shear force and thus takes a finite time. This implies the effectiveness of void filling increases with the period ≃1/ω. Thus, expect the ascent per vibration may scale not with ωv2 but this factor multiplied with the effectiveness factor 1/ω, which results in U1 scaling with v2, in accordance with eqn (1). Note that the number of beads falling off from the arch into the void will increase as the life time of the transient arch increases, i.e., the period of vibration increases (as shear tends to increase the number). In this sense, the effectiveness of void filling increases with 1/ω (as W decreases the effectiveness reducing shear).
According to eqn (1), U1 is independent of the disk diameter D. This is in contrast with the ascent ascribed to the void filling mechanism in previous studies.2,3,7 The void considered in previous studies is not the arch-shaped void but a void formed beneath the disk as illustrated in Fig. 3b. This implies the ascent velocity increases with D. However, in the present case, the size of the arch-shaped void is independent of the size of the disk (this is because the width is set by W and the height is set by v2/g as explained above), and, thus, it is natural that U1 is independent of D.
For later discussion, we here confirm experimentally that the width of the downward flow developed near a side wall is comparable to the size of the grain, which is consistent with previous reports.5 This is shown in Fig. 4b quantifying the width of the downward flow near the wall. To obtain this plot, we took snapshots of a vibrated cell with a single tracer particle and superposed only a selected set of snapshots for a given duration, where we selected only the snapshots in which the particle moved down compared with its previous snapshot. An example of such superimposed snapshots is given in Fig. 4c. We estimated the width from the maximum width (as a function of height) as indicated in the figure. The width thus obtained is given as a function of time for two cells with different width as in Fig. 4b. As shown in the plot, the width saturates with time and the saturated value, which is several times the grain size, is independent of the cell width.
The independence of the width of the downward flow near the side wall from the cell width, confirmed in Fig. 4b is physically natural because such a downward flow is caused by the friction effect near the side walls, and this effect should be the same for cells of different width as long as the side walls are well separated compared with the width of the downward flow. In addition, this implies that the flux of the downward flow is independent from the cell width W for a given v, and that the shear force which governs the flow is characterized not by v/W but by v/d.
Considering that the flow flux is conserved in the rolling motion, we expect that upwards velocity in the absence of the disk (and away from the disk) scales with a dimensionless factor d/W as illustrated in Fig. 5a. However, alongside the disk, it scales with d/(W − D) as suggested in Fig. 5b. In other words, there are two factors d/W and d/(W − D), and both govern the rising velocity of the disk. Considering that the former is a decreasing function of W and the latter introduces a factor that increases with D, and expecting that the expression is independent of the smallest scale, the bead diameter d, a simplest possibility incorporating the two factors is a dimensionless factor D/W. This explains why U2 scales with D/W in eqn (2).
![]() | ||
Fig. 5 (a) The effect of cell width on upward flow velocity. (b) The effect of an intruder's diameter on upward flow velocity. |
This D dependence of U2 is in contrast with results reported in previous studies:5,10 they reported that the velocity of ascent due to convection is independent of the size of the intruder. This may because the previous cases were more insensitive to the effect; they both considered not a two-dimensional case as ours, but a three-dimensional case, and they considered cases in which the ratio of the intruder diameter to grain diameter is larger (the ratio is ≃3.4–306 in ref. 5 and ≃8–13 in ref. 10 while it is ≃5–7 in our case).
This journal is © The Royal Society of Chemistry 2020 |