arxiv:cond-mat/ v1 [cond-mat.stat-mech] 12 Sep 2000

Size: px

Start display at page:

Download "arxiv:cond-mat/ v1 [cond-mat.stat-mech] 12 Sep 2000"

Mervyn Watts
5 years ago
Views:

1 arxiv:cond-mat/ v1 [cond-mat.stat-mech] 12 Sep 2000 An Experimental Study of a Granular Gas Fluidized by Vibrations Éric Falcon 1, Stéphan Fauve 1, and Claude Laroche 2 1 Laboratoire de Physique Statistique, École Normale Supérieure, 24, rue Lhomond, Paris Cedex 05, France 2 Laboratoire de Physique, École Normale Supérieure de Lyon, 46 allée d Italie, Lyon Cedex 07, France Abstract. We report experimental results on the behavior of an ensemble of inelastically colliding particles, excited by a vibrated piston in a vertical cylinder. When the particle number is increased, we observe a transition from a regime where the particles have erratic motions (granular gas ) to a collective behavior where all the particles bounce like a nearly solid body. In the gaslike regime, we measure the pressure at constant volume, and the bed expansion at constant external pressure, as a function of the number N of particles. We also measure the density of particles as a function of the altitude, and find that the atmosphere is exponential far enough from the piston. From these three independent measurements, we determine a state equation between pressure, volume, particle number and the vibration amplitude and frequency. 1 Introduction Once fluidized a vibrated granular medium looks like a gas of particles that can be described using kinetic theory of usual gases. However, granular gases basically differ from ordinary ones mostly due to the inelasticity of collisions, i.e., nonconservation of energy. While over the years many attempts based on kinetic theory [1] have been made to describe such dissipative granular gases, no agreement has been found so far both with experiments [2,3] and numerical simulations [3,4,5,6], for the dependence of the granular temperature, i.e., the mean kinetic energy per particle, on the parameters of vibration [7,8,9,10]. The aim of this study is to guess possible gas-like state equations for such a dissipative granular gas and to observe new kinetic behaviors which trace back to the inelasticity of collisions, e.g., the tendency of such media to form clusters. Although this feature has probably been known since the early observation of planetary rings [11] and although various cluster types in granular flows have been observed numerically [12], there exist only a few recent laboratory experiments. One experiment, with a horizontally shaken 2 D layer of particles, displayed a cluster formation, but the coherent friction force acting on all the particles was far from being negligible [13]. We performed a similar experiment by exciting a 3 D granular medium with a vertically vibrated piston. We did observe clustering, but we could not rule

2 2 Éric Falcon et al. out a lock-in mechanism involving the time scale connected with gravity and the period of vibration[14]. We thus repeated this experiment in a low-gravity environment, where inelastic collisions were the only interaction mechanism, and observed a motionless dense cluster that confirms that the inelasticity of collisions alone can generate clustering [15]. This paper is devoted to the study of the low-density situation, where clustering does not occur. It is organized as follows. Section 2 is devoted to the presentation of our experimental setup. The experimental results are presented in Secs We report in Sec. 3 (resp. Sec. 4) the measurements of pressure (resp. volume) of a gas of spherical particles excited by a vibrating piston and undergoing inelastic collisions. At constant external driving, we show that the pressure passes through a maximum for a critical number of particles before decreasing for large N. The density of particles as a function of the altitude is studied in Sec. 5, where we observe an exponential atmosphere far enough from the piston. From measurements of Secs. 3 5, we show in Sec. 6 that the dependence of the granulartemperature, T, on the piston velocity, V, is of the form T V θ, where θ is a decreasing function of N. Finally, we discuss our results in the light of previous works [2,3,4,5,6,7,8,9,10] and give our conclusions in Sec Experimental Setup The experiment consists of a transparent cylindrical tube, with an inner diameter of 60 mm, filled from 20 up to 2640 stainless steel spheres, 2 mm in diameter, roughly correspondingto 0 up to about 5 particle layersat rest. An electrical motor, with eccentric transformer from rotational to translational motion, drives the particles sinusoidally with a 25 or 40 mm amplitude, A, in the frequency range from 6 to 20 Hz. A lid in the upper part of the cylinder, is either fixed at a given height, h (constant volume experiment) or is mobile and stabilized at a given height h m due to the bead collisions (constant pressure experiment). The heights h and h m are defined from the lower piston at full stroke. When the vibration is strong enough and the number of particles is low enough, the particles display ballistic motion between successive collisions like molecules in a gas(see Fig. 1a).When the density ofthe medium is increased, the gaslike state is no longer stable but displays the formation of a dense cluster bouncing like a nearly solid body (see Fig. 1b). This paper is devoted to the study of the gaslike state. Its aim is to determine experimentally a state equation between pressure, volume, particle number and the vibration parameters (amplitude and frequency).

3 An Experimental Study of a Granular Gas Fluidized by Vibrations 3 Fig.1. Transition from a dissipative granular gas to a dense cluster: (a) N = 480; (b) N = 1920, respectively corresponding to roughly 1 and 4 particle layers at rest. The parameters of vibration are f = 20 Hz and A = 40 mm. The driving piston is at the bottom (not visible), the inner diameter of the tube being 52 mm 3 Pressure Measurements Time averaged pressure measurements have been performed as follows. Initially, a counterwheight of mass 46 g balances the lid mass. The piston drives stainless steel spheres in erratic motions in all directions (see Fig. 1a). Particles are hitting the lid all the time, so that to keep it at a given height h we have to hold the lid down by a given force, Mg, where M is the mass of a weight we place on the lid and g the acceleration of gravity. At a fixed h, i.e. at a constant volume, Fig. 2 shows the time averaged pressure P = Mg/S exerted on the lid as a function of the number N of beads in the container, for different frequencies of vibration, S being the area of the tube crosssection. At constant external driving, i.e. at fixed f and A, the pressure passes through a maximum for a critical value of N roughly corresponding to 0.8 particle layers at rest. This critical number is independent of the vibration

4 4 Éric Falcon et al. 80 Pressure, P (Pa) N Fig.2. Mean pressure P as a function of N. From the lower ( marks) to the upper (hexagrams) curve, vibration frequency f varies from 10 to 20 Hz with a 1 Hz step. For all these experiments, h h 0 = 5 mm and A = 25 mm, h 0 being the bed height at rest. One single layer of particles at rest corresponds to N = 600. Lines join the data points frequency. A further increase of the number of particles leads to a decrease in the mean pressure since more and more energy is dissipated by inelastic collisions. Note that gravity has a small effect in these measurements that are performed for V 2 gh, where V = 2πfA is the maximum velocity of the piston. For N such that one has less than one particle layer at rest, most particles perform vertical ballistic motion between the piston and the lid. Thus, the mean pressure increases roughly proportionally to N. When N is increased such that one has more than one particle layer at rest, interparticle collisions become more frequent. The energy dissipation is increased and thus the pressure decreases. 4 Volume Measurements We now consider the bed expansion under the influence of collisions on a circularwire mesh lid placed on top ofthe beads leavinga clearanceofabout 0.5 mm between the edge of the lid and the tube. Due to the bead collisions, the lid is stabilized at a given height h m from the piston at full stroke. Although the lid mass is roughly 50 times smaller than the total mass of the beads, the lid proves to be quite stable and remains horizontal. The expansion, h m h 0, of the bed is displayed in Fig. 3 as a function of N for different vibration frequencies. h 0 is the bed height at rest. At fixed f, the expansion passes through a maximum for a critical value of N roughly corresponding to 0.6 particle layers at rest. This critical number is independent of the vibration frequency. When N is further increased, the expansion decreases showing,

5 An Experimental Study of a Granular Gas Fluidized by Vibrations 5 10 Expansion, h m h 0 (cm) N Fig.3. Maximal bed expansion, h m h 0, as a function of N, for various frequencies f of vibration. From the lower ( marks) to the upper ( marks) curve, f varies from 7 to 20 Hz with a 1 Hz step and A = 25 mm. One single layer of particles at rest corresponds to N = 600. Lines join the data points as for pressure measurements, an increase in dissipated energy by inelastic collisions. Note that the height h m of the granular gas is much larger than for the pressure measurements of Fig. 2. Consequently, gravity is important here. As already found experimentally in 1 D [3] and 3 D [16,17] and numerically [3,18], we observe that, at N fixed, the granular medium exhibits (see Fig. 4a) a sudden expansion at a critical frequency corresponding to a bifurcation similar to that exhibited by a single ball bouncing on a vibrating plate [3,16]. Moreover, Fig. 4a shows that this critical frequency depends on the number of layers, n. When n increases above 0.4, a transition from the 1 D like behavior to a 3 D one is observed: the expansion at the critical frequency becomes less abrupt and tends to increase regularly with f (see Figs. 4a b). 5 Density Measurements Time averaged density measurements at a given height z are performed by means of two closely coupled coils, δz = 5 mm in height, and 64 mm in inner diameter, the cylindrical tube now being 52 (resp. 62) mm in inner (resp. outer) diameter. An 1.5 khz a.c. voltage is applied to the primary coil, the turns ratio of the transformer being roughly equal to 2. Steel spheres moving accross the permanent magnetic field of the primary coil, generate an inductive voltage variation accross the secondary one. This root mean square a.c. voltage U is a function of the mean inductance and mutual variations which are proportional to the mean number of particles in the

6 6 Éric Falcon et al. Expansion, h m h 0 (cm) (a) Expansion, h m h 0 (cm) (b) f (Hz) f (Hz) Fig.4. Maximal bed expansion, h m h 0, as a function of f, for various numbers of particles N. (a) From the lower ( marks) to the upper ( marks) curve, N = 20, 40, 80, 120, 240 and 360. (b) From the upper ( marks) to the lower ( marks) curve, N = 480, 600, 720, 960, 1200, 1440, 1680, 1920, 2160 and For all these experiments A = 25 mm. One single layer of particles at rest corresponds to N = 600 volume delimited by the sensor at altitude z from the piston surface at full stroke. We have calibrated the sensor with steel spheres at rest and we have checkedthat U N.Wehavealsofoundthattheeffectofspheresoutsideof the sensor volume decays exponentially with the distance to the sensor, with a 10±1mmdecaylengthindependentofN,forourrangeofN.Particlesdensity as a function of altitude is shown in Figs. 5a c for 3 different total numbers of particles and for various frequencies. For each N, the profile density in log linear axes displays a decay (at low f), a plateau (at intermediate f) or a dip (at high f) near the piston and an exponential decay in the tail at high altitude, whatever f. As for an isothermal gas, the atmosphere is found to be exponential far enough from the piston, but on very different length scales, i.e., few cm (resp. km) for our experiment (resp. for air). Such a dense upper region supported on a fluidized low-density region near the piston has been also reported numerically[18] and predicted theoretically[19]. Although the dip in the density profiles at the bottom was already observed in a 2 D granular gas experiment, non-negligible coherent friction force acting on all the particles did not allow determination of the granular temperature dependence on the piston velocity V from exponentional Boltzmann distributions fitted to tails of density profiles [2]. We can fit an exponential curve to the tail of the profile density. From the decay rate, ξ, in the fitted exponential and using kinetic theory [2], we can extract the dependence of the granular temperature T on the piston velocity. In Fig. 6, we plot 1/ξ (which is proportional to T) against V with log-log-axis together with power law fits of the form T V θ where θ is n dependent. Note that this power law, being observed only on a small range of velocities, one cannot rule out another functional behavior. In

7 An Experimental Study of a Granular Gas Fluidized by Vibrations (a) 10 2 (b) Density (arb. unit) 10 1 Density (arb. unit) Height (cm) Height (cm) (c) Density (arb. unit) Height (cm) Fig.5. Mean density as a function of the height, for various frequencies f of vibration and 3 numbers of particles (a) N = 480: f = ( ) 5, (+) 6, ( ) 7, ( ) 8, ( ) 9, ( ) 10, ( ) 11, ( ) 12 and ( ) 13 Hz; (b) N = 720: f = ( ) 5, (+) 5.6, ( ) 7.1, ( ) 8.6, ( ) 10.1, ( ) 12.2, ( ) 14 and ( ) 15.2 Hz; (c) N = 1440: f = ( ) 7, (+) 8, ( ) 10, ( ) 11.4, ( ) 12.8, ( ) 15 and ( ) 17 Hz. For all these experiments, A = 40 mm. One single layer of particles at rest corresponds to N = 480. Lines join the data points particular, the faster increase of T at high velocity is not significant because of imprecision on the exponential decay of the density (see Fig. 5). 6 Towards a State Equation In order to use the above measurements to determine a state equation, we have to find the appropriate dependence T = T(V,n) of the granular temperature as a function of the vibrating velocity V and the number of particle layers n. It is known that for a fixed number n of granular layers at rest, one has in the low-density limit PΩ T [20], where P is the mean pressure and Ω the volume. Taking into account this law for small densities and our observation T V θ(n) from density measurements in an exponential atmosphere, we have plotted log(p) and log(h m h 0 ) as functions of log(v). On

8 8 Éric Falcon et al. T (arb. unit) V (m s 1 ) Fig. 6. Log-log plot of granular temperature versus V for various numbers of layers: ( ) 0.8, ( ) 1, ( ) 1.2, ( ) 2, ( ) 2.4 and ( ) 3. Experiments ( ), ( ) and ( ) (resp. ( ), ( ) and ( )) are performed for A = 25 mm (resp. A = 40 mm). Power law fits of the form V θ are dispayed in solid lines θ (n) n Fig. 7. Evolution of the exponent θ as a function of the number of layers, n, from pressure ( ), volume ( ) and density ( ) measurements. Fits are θ(n) = 1 tanh(n n c) with n c = 2 (straight line), 2.3 (dot-dashed line) and 3.5 (dashed line) the reported frequency range, these curves are straight lines, the slopes of which give θ(n). The behavior of θ(n) for the experiments at constant volume in Sec. 3 (resp. constant pressure in Sec. 4) is displayed in Fig. 7 with mark (resp. mark) together with the one in mark extracted from exponential density profiles of Sec. 5. The three curves, obtained with different experimental conditions and independent measurements have the same shape which could be simply fitted by θ = 1 tanh(n n c ) where n c = 3.5 (resp. 2.3) and 2.

9 An Experimental Study of a Granular Gas Fluidized by Vibrations 9 We can now use the observed law T V θ(n) to scale the pressure and bed expansion measurements of Figs. 2 and 3. The results are displayed in Figs. 8a b and show a rather good collapse of all the data on a single curve. We have thus shown that the law, PΩ T, together with T V θ(n), provide a correct empirical state equation for our dissipative granular gas in the kinetic regime. As shown earlier, this regime is limited at high density by the clustering instability [13,14,15] and on the other side, for a fixed too small number of particles, when the gas suddenly contracts on the piston below a critical frequency (see Fig. 4a and [3,16,18]) (a) (b) P / f θ (Pa s θ ) 10 1 (h m h 0 ) / f θ (cm s θ ) n n Fig.8. (a) Mean pressure P from Fig. 2 rescaled by f θ as a function of n. θ(n) = 1 tanh(n n c) with n c = 3.5. (b) Maximal bed expansion, h m h 0, from Fig. 3 rescaled by f θ as a function of n. θ(n) = 1 tanh(n n c) with n c = Conclusion The aim of this study was an experimental determination of the state equation of dissipative granular gases. It is known that for a fixed number of granular layers at rest, one has in the low-density limit PΩ T [20], where P is the mean pressure, Ω the volume and T the granular temperature. However, the dependence of T on the vibration amplitude, A, and frequency, f, of the piston and on the number of particles, n, is still a matter of debate [7,8,9,10]. Kinetic theory [2,9] or hydrodynamic models [7] show T V 2 n 1, whereasnumerical simulations [3,4,5,6] or experiments [2,3] givet V α n β, with 1.3 α 2 and 0.3 β 1, where V = 2πfA is the maximum velocity of the piston. These previous conflicting results may be explained by the V dependence of T that we have found from our measurements: T V θ(n), with θ continuously varying from θ = 2 when n 0, as expected from kinetic theory, to θ 0 for large n. We emphasize that we have only considered the dependence of the granular temperature on the vibration velocity V. For the

10 10 Éric Falcon et al. dependence on the number of particles, another term of the form n β should exist, that we cannot determine from these measurements. References 1. J. T. Jenkins, S. B. Savage, A theory for the rapid flow of identical, smooth, nearly elastic, spherical particles. Fluid Mech. 130, (1983); C. S. Campbell, Rapid granular flows. Ann. Rev. Fluid Mech. 22, (1990) 2. S.Warr,J.M. Huntley,G.T.H.Jacques,Fluidization ofa2 Dgranularsystem: Experimental study and scaling behavior. Phys. Rev. E 52, (1995) 3. S. Luding, E. Clément, A. Blumen, J. Rajchenbach, J. Duran, Studies of columns of beads under external vibrations. Phys. Rev. E 49, (1994) 4. S. Luding, H. J. Herrmann, A. Blumen, Simulations of 2 D arrays of beads under external vibrations: Scaling behavior. Phys. Rev. E 50, (1994) 5. S. Luding, Granular materials under vibration: Simulations of rotating spheres. Phys. Rev. E 52, (1995) 6. H. J. Herrmann, S. Luding, Modeling granular media on the computer. Continuum Mech. Thermodyn. 10, (1998) 7. J. Lee, Scaling behavior of granular particles in a vibrating box. Physica A 219, (1995) 8. S. McNamara, S. Luding, Energy flows in vibrated granular media. Phys. Rev. E 58, (1998) 9. V. Kumaran, Temperature of a granular material fluidized by external vibrations. Phys. Rev. E 57, (1998) 10. J. M. Huntley, Scaling laws for a 2 D vibro-fluidized granular material. Phys. Rev. E 58, (1998) 11. P. Goldreich, S. Tremaine, The dynamics of planetary rings. Ann. Rev. Astron. Astrophys. 20, (1982) 12. M. A. Hopkins, M. Y. Louge, Inelastic microstructure in rapid granular flows of smooth disks. Phys. Fluids A 3, (1991); S. McNamara, W. R. Young, Inelastic collapse and clumping in a 1 D granular medium. Phys. Fluids A 4, (1992); I. Goldhirsch, G. Zanetti, Clustering instability in dissipative gases. Phys. Rev. Lett. 70, (1993) 13. A. Kudrolli, M. Wolpert, J. P. Gollub, Cluster formation due to collisions in granular material. Phys. Rev. Lett. 78, (1997) 14. É. Falcon, S. Fauve, C. Laroche, Cluster formation, pressure and density measurements in a granular medium fluidized by vibrations. Eur. Phys. J. B 9, (1999) 15. É. Falcon, R. Wunenburger, P. Évesque, S. Fauve, C. Chabot, Y. Garrabos, D. Beysens, Cluster formation in a granular medium fluidized by vibrations in low gravity. Phys. Rev. Lett. 83, (1999) 16. C. E. Brennen, S. Ghosh, C. R. Wassgren, Vertical oscillation of a bed of granular material. J. Appl. Mech. 63, (1996) 17. M. L. Hunt, S. Hsiau, K. T. Hong, J. Fluid Eng. 116, (1994) 18. Y. Lan, A. D. Rosato, Macroscopic behavior of vibrating beds of smooth inelastic spheres. Phys. Fluids 7, (1995) 19. D. A. Kurtze, D. C. Hong, Effect of dissipation on density profile of 1 D gas. Physica A 256, (1998) 20. S. McNamara (1998) Private Communication

Velocity statistics in excited granular media

CHAOS VOLUME 9, NUMBER 3 SEPTEMBER 1999 Velocity statistics in excited granular media W. Losert, D. G. W. Cooper, and J. Delour Department of Physics, Haverford College, Haverford, Pennsylvania 19041 A.