Granular dampers in microgravity: sharp transition between modes of operation

When a container filled with granular material is subjected to sinusoidal vibration in microgravity, dependent on the amplitude of the oscillation, the granulate may exhibit one of two distinct dynamical modes: at low amplitude, a gas-like state is observed, where the particles are relatively homogeneously distributed within the container, almost independent of the phase of the oscillation. In contrast, for large amplitude, collective motion of the particles is favoured, termed collect-and-collide regime. Both regimes are characterized by very different dissipation characteristics. A recent model predicts that the regimes are separated by a sharp transition due to a critical amplitude of the vibration. Here we confirm this prediction of a sharp transition and also the numerical value of the critical amplitude by means of experiments performed under conditions of weightlessness.


Granular dampers
When granular particles collide with one another or with container walls, they dissipate mechanical energy. This fundamental property which, in fact defines the term granular particle originates from the fact that grains are macroscopic or thermodynamic bodies with many internal degrees of freedom. In contrast, molecular particles with only a few internal degrees of freedom cannot dissipate energy but quickly reach thermodynamic equilibrium where on average, the energy is evenly distributed among all degrees of freedom. The dissipative properties of granular matter do not only give rise to a rich physical phenomenology but are also of significant interest for vibration damping in technical applications, ranging from medicine [1,2], through recreational sports [3] to construction [4][5][6][7].
In their most basic form, granular vibration dampers comprise an assembly of particles housed in a sealed container or cavity. When exposed to oscillations of an adequately high vibrational energy ( Ē ∝ 2 A 2 , where and A are, respectively, the angular frequency and amplitude of the vibration), the particles will overcome the sedimenting force of gravity, collide with one another and the walls of the housing container and, consequently, dissipate energy. This dissipation attenuates the vibration and, thus, damps the vibrating system [8][9][10][11][12][13][14]. In order to isolate and better understand these phenomena, it is often valuable to study the dynamics of these systems in a zero-gravity environment, where complexities due to gravitational acceleration, and the resultant spatial heterogeneities 1 and strong particle-wall frictional interactions induced thereby are limited. It should be noted, however, that the simplified, gravity-free systems studied here are not without direct real-world application -see, for example, the Russi an EVA hamme r, used for work outside the MIR space station.

Vibrated granular systems in microgravity
The dynamics of vibrated granular media in microgravity have been the subject of academic study for more than two decades. Early work by Falcon et al. [15,16] showed that a dilute granular medium in microgravity will maintain a relatively homogeneous density distribution, whereas a denser system will typically collapse into a single dense cluster. Low-gravity environments have also proven a valuable means of studying fundamental phenomena and properties of granular media, such as velocity-dependent restitution coefficients [17,18], wave propagation [19], and the coupling of rotational and translational energy [20][21][22][23], an important problem in the field of granular physics, where equipartition cannot be assumed [24,25]. Indeed, the simplicity of microgravity systems makes them ideal for investigating various non-classical behaviours exhibited by granular gases [26,27], in particular their well-known non-Gaussian velocity distributions [28,29]. More recently, significant interest has developed in the study of the phase transitions of vibrated dry granular media in microgravity [30][31][32][33][34], and it is one such transition which forms the focus of the present work.

Modes of operation of granular dampers
Previous experiments [35][36][37] revealed two very distinct dynamical modes of granulate in a granular damper when subjected to sinusoidal motion, A(t) = A cos( t + x ) , in weightlessness: at low driving amplitudes, the granulate is found to exist in a gas-like state, where the particles near-homogeneously fill the available space. In contrast, at large amplitudes the particles move coherently in synchrony with the external vibration which suggests the modeling of the granulate's motion as a single completely inelastic quasi-particle oscillating between the walls of the container. The current mode of operation is independent of the frequency but depends only on the amplitude of the vibration. The threshold amplitude, A 0 , discriminating the modes depends only on the physical free space inside the container, i.e. the space not occupied by particles.
A similar transition between a gas-like and a solid-like state has also been observed under gravity [38]. However, unlike in the microgravity experiments described in [35][36][37] which suggest a single phase boundary between the two states, these experiments showed strong hysteresis [38]. Other experiments under Earth's gravitational conditions also demonstrated more complex dynamics, such as a liquid-like 'sloshing' motion [39] that is not observed in microgravity experiments.
The previous microgravity experiments discussed above focused on the highly dissipative collective regime, while the region near the threshold was studied insufficiently. In the current paper we address this issue, exploring in detail the damping properties in the gas regime and near the critical amplitude. A few measurements for the collective regime have been repeated for calibration.
In order to understand the key rôle of the threshold separating the two regimes, we first need to understand the main differences between said regimes [35,36]: at low driving amplitudes, the particles fill the whole container nearly homogeneously with the exception of a gradient which forms where the walls of the container sweep the space. Particles in the zones which are swept by walls once per period get accelerated in the process and successively lose their kinetic energy in the subsequent cascade of dissipative particle-particle and particle-wall collisions. Consequently, at low amplitudes in each period only a small fraction of all particles comes into contact with the wall and, therefore, the gas state is characterized by a relatively low energy dissipation rate.
In contrast, at large driving amplitude, particles are observed to undergo collective motion: during the inward stroke of a wall, the particles are gathered on the inwardaccelerating wall due to a granular collapse [40,41]. At the point of inflection of the motion, the bulk of particles leaves the wall as an agglomerate with only a small variation in velocity between its constituent particles. Although these residual differences in the particles' velocities cause some small dispersion of the collective while streaming through the container [42], the dispersed collective will be collected again by a granular collapse, when it reaches the opposite container wall at a phase when the wall is accelerating inwards, that is, towards the approaching particles. For this regime the term collect-and-collide was coined [35]. A corresponding model description [36] explained the behavior of the granular damper including its two modes of operation and the corresponding dissipative properties in quantitative agreement with experiments. This model also explains why the amplitude of a spring with an attached granular damper decays nearly linearly in time [43], which was a long standing question. Under certain conditions, a granulate consisting of as few as 4 particles shows the collective behavior described above [37].
Under weak assumptions, the collect-and-collide and gas regimes are separated by a threshold amplitude of the oscillation: (1) where L g is the free space inside the container which the cluster has to travel through between opposite walls of the container. Interestingly, the theory predicts that the value of A 0 is independent of the frequency of the oscillation and also material properties, in agreement with experimental results and numerical particle simulations [35][36][37]. For details of the derivation of Eq. (1) see [36]. Since the energy dissipation rate of both regimes is very different, the value of A 0 can be determined experimentally by measuring the power which is dissipated by the granular damper as a function of the vibration amplitude.

Experiment
The experimental setup, sketched in Fig. 1, follows the same design as that employed in [36]: a cuboidal box constructed from 4 mm-thick polycarbonate plates is filled with granular material and mounted on a linear bearing driven by a linear actuator. A force sensor is mounted between the wagon of the linear rail and the box, allowing measurement of the force F(t) needed to impose sinusoidal motion on the box, where A and are, respectively, the amplitude and the angular frequency of the motion. A total of four such boxes are housed in the experimental rack. Each pair of boxes is mounted on a single linear rail and driven in anti-parallel orientation, such that the total momentum of each pair adds up to zero in order to minimise the vibrations exerted upon the experimental rack. Further, a Hall-effect based position encoder is mounted on each linear rail to measure the instantaneous position x(t) of the box. Each box is monitored by a camera ( 640 × 480 pixels at 240 fps) from the side (same perspective as the sketch in Fig. 1) and each pair of boxes is monitored by a high-speed camera ( 1280 × 1024 pixels at 500 fps) from the top.
To exclude the sedimenting effect of gravity, this experiment has been performed in an aircraft performing parabolic (2) x(t) = A cos t + x flights, where in each parabola approximately 22 s of microgravity is achieved. Prior to each parabola, the apparatus is set to oscillate at desired amplitudes and a frequencies. At the same time, the recording of the force and the positional data is initiated at sample rate 10 kHz , and the cameras are triggered to record throughout the interval of microgravity. After the completion of each parabola, the oscillation is stopped and, if necessary, the boxes are swapped. This way we were able to study a variety of granular systems with different particle numbers and/or physical or geometrical properties. In order to obtain maximal information regarding the transition from the collect-and-collide to the gas regime, we focused on the range of amplitudes near the expected transition. The measurements presented below, were obtained with three distinct systems as described in Table 1, all containing the same filling volume, indicated by the empty space length, L g .
Due to the box size and the amount of filling described by the empty space length, L g , we expected the transition between the gas to the collective regime at amplitude A 0 = L g ∕ ≈ 29.9 mm . Therefore, the range of amplitude in the experiment was chosen 20 mm ≤ A ≤ 40 mm , to ensure the observation of the transition between the regimes.

Data processing
We obtain the dissipation rate of the granular dampers by averaging the time dependent power needed to set the damper in vibration, P(t), over time. Note that the momentary value of P(t) may be positive or negative, depending on the phase of the oscillation. We obtain  where Eq. (2) was used. Expanding the force in harmonics with coefficients F i and i , we write w h e r e w e e x p l o i t t h e i d e n t i t y 2 cos(x) sin(y) = sin(x + y) − sin(x − y) . Averaging over one period of oscillation, we obtain the energy supplied to the system by the driving mechanism. Since we are in stationary state, this is identical to the energy dissipation rate, where t * is an arbitrary time and we took into account that the integral vanishes for all terms in Eq. (4), except the (n − 1)-term for n = 1 , corresponding to the frequency of the excitation.
Equation (5) thus allows us to determine the energy dissipation rate from the lowest order Fourier expansion of the measured driving force.

Results
To compute the energy dissipation rate of the granular damper, we measure the driving force as a function of time during the interval of weightlessness, compute the Fourier transform, and select the respective amplitude, F 1 , and phase, 1 , corresponding to the frequency of the driving. Then together with the amplitude and phase of the excitation, A and x , we determine P diss via Eq. (5). Figure 2 shows the dissipated power of samples No. 1-3, described in Table 1 as functions of the vibration amplitude.
For amplitude A ≲ 30 mm the system is in the gas state where the energy dissipation rate is low, ( P diss ≲ 7 mW ) while for A ≳ 30 mm the system is in the collect-andcollide regime, characterized by larger dissipation rate, P diss ≳ 30 mW . The transition amplitude agrees well with the prediction made by the model developed in [36], A 0 ≈ 29.9 mm , see Eq. (1). The transition between the regimes can also be seen from the images of the cameras. For sample 2, Fig. 3 illustrates the abrupt transition between the regimes through series of snapshots from one period of the vibration taken for amplitudes, A = 29.5 mm (left column) and A = 31.5 mm , close to the transition point, A 0 .
For a more quantitative analysis of the transition between the dynamical regimes and the corresponding energy dissipation rates, we compare the results from the current experimental campaign presented in Fig. 3 with earlier data obtained during the 54th ESA parabolic flight where the energy dissipation rate of an equally sized granular damper was studied, but filled with varying masses of 4 mm steel beads (for details see [36]). For the comparison, we scaled the data in the same way as in the presentation of [37], namely the amplitude of the vibration was scaled by  Table 1 as functions of the vibration amplitude L g , such that the value A∕A 0 discriminates the amplitudes corresponding to the gas state from the amplitudes corresponding to the collect-and-collide state. Second, the energy dissipation rate was scaled by the maximum energy dissipation rate which the system can possibly dissipate, P max = 2 mA 2 3 similar to [36]. Figure 4 shows the current results superimposed by the results from the 54th ESA parabolic flight, scaled in the described way. We notice a very good quantitative agreement between the data from [36] and the current data. Deviations such as the point from sample No. 4 at A∕A 0 ≈ 0.67 may be attributed to the coarse grained representation due to the sparse data in the range A ≈ A 0 available in the earlier data set. While both sets of data obtained in independent experiments support one another, the new set clearly supports our conjecture made in [36], regarding the numerical value of the transition amplitude, A 0 , and the fact that the transition is not smooth but occurs abruptly when A approaches A 0 .

Summary
When a granular damper is vibrated under conditions of microgravity, the energy dissipation rate depends strongly on the amplitude of the vibration, A, but shows no significant dependence on the frequency. When slowly increasing the amplitude, at the transition amplitude A 0 , the granulate leaves the gas state characteristic for A ≲ A 0 and assumes a different dynamical state called collect-and-collide regime characterized by collective motion of the particles. The transition between the dynamical states corresponds to an abrupt change of the energy dissipation rate. These statements were the main results of an experiment performed during the 54th ESA parabolic flight in 2011 and published in [36]. Since the theoretical model for A 0 was not available at this time, the amplitude was swept over a wide range of values such that the interesting range in the neighborhood of A 0 was covered only by a few data points, for details see [36]. Therefore, one could still have doubts about the statements above, thus, to refine the data in the range A ≈ A 0 with the aim to characterize the transition between the regimes in more detail was the aim and motivation of the current study. The current experiments obtained during the 31st DLR parabolic flight campaign resolve the region around the transition amplitude, A ≈ A 0 in much greater detail than previous measurements. As a result, the data support our conjecture made in [36], regarding the numerical value of the transition amplitude, A 0 , and the fact that the transition is not smooth but occurs abruptly when A approaches A 0 . Fig. 4 Normalised energy dissipation rate as a function of the normalized amplitude of the vibration: Symbols connected by lines show current data while data represented only by symbols correspond to systems studied during the 54th ESA parabolic flight campaign in 2011 [36]