# Robust experimental study of avalanche precursory events based on reproducible cycles of grain packing destabilizations

Luc Oger, Renaud Delannay, Yves Le Gonidec

► To cite this version:

Luc Oger, Renaud Delannay, Yves Le Gonidec. Robust experimental study of avalanche precursory events based on reproducible cycles of grain packing destabilizations. Powders and grains 2021, Jul 2021, Buenos Aires, Argentina. pp.03023, 10.1051/epjconf/202124903023. hal-02987121v2

HAL Id: hal-02987121

https://hal.science/hal-02987121v2

Submitted on 30 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Robust experimental study of avalanche precursory events based on repro-
ducible cycles of grain packing destabilizations

Luc Oger¹,*, Renaud Delannay¹,**, and Yves Le Gonidec²,***

¹Univ. Rennes, CNRS, IPR [(Institut de Physique de Rennes)]-UMR 6251, F-35000 Rennes, France

²Univ. Rennes, CNRS, Géosciences Rennes-UMR 6118, F-35000 Rennes, France

**Abstract.** Quasi-periodic collective displacements of grains at the free surface of a tilted grain packing consti-
tute precursors of granular avalanches. Laboratory experiments are commonly performed by slowly tilting the
packing from 0° to the maximal stability angle θ<sub>A</sub>. In these conditions, the number of precursors is too small
to assess reproducible and robust statistical analyses of the precursor activity. To go beyond this limitation,
we have developed a specific experimental protocol consisting in tilting the packing with successive oscillation
cycles. We use a high resolution optical camera and process the images of the packing free surface to identify
precursory events during many consecutive cycles of a single packing. We observe the same behavior for all
half-cycles, forth and back: appearance of the first precursors after the same variation of inclination, exponen-
tial evolution of the weak surface activity for the first precursors and linear growth of stronger surface activity
for the following ones. The experimental protocol provides both reproducible precursor measurements based
on large sample statistical inferences and a quasi-stationary state after one full-cycle. This approach is very
promising to highlight the effects of external parameters, including humidity and packing geometry.

# 1 Introduction

Understanding the triggering mechanisms of granular
flows is a subject of great interest in many domains.
For instance, in pharmaceutical and agricultural domains,
grains can be stored, pilled or displaced; in natural do-
mains, snow avalanches or giant landslides require pre-
diction to prevent the recurrence of such dramatic events.
At small scales, laboratory experiments are performed to
study gravity destabilization of a packing of glass beads by
gently inclining it, assuming a quasi-static behavior. The
quasi-static regime is satisfied if the inclination process al-
lows the non-displacement of the grains except when local
shear forces reach a critical value related to the Coulomb's
friction law. The destabilization is generally character-
ized by grain rearrangements at the free surface, based on
optical measurements [1–3] but many other experimental
techniques have also been successfully used (see [4] for a
comprehensive presentation).

Nerone et al. [5, 6] have shown that the size of the
grain rearrangements at the free surface increases with
the inclination angle θ for freshly prepared packings and
that quasi-periodical large collective events can be ob-
served few degrees before the maximum stability angle
associated to the avalanche θA. They called precursors
such events and observe that their size, i.e. the number
of grains involved, grows exponentially with respect to
θ. A better understanding of the precursor activity is re-

quired to predict the onset of the avalanche. The avalanche
triggering and the irruption of precursors are controlled
by several physical parameters, such as the packing di-
mensions (height [7], length [2] and width [8–11]), the
density or volume compaction of the packing [2, 7, 12–
14], the tilting regime [4, 15, 16] and the humidity in-
side the packing [17]. Most of the previous experi-
ments [1–3, 5, 7, 13, 17–19] dedicated to precursor iden-
tification have been performed by direct tilting experi-
ments of a grain packing from the horizontal position to
the avalanche.

Kiesgen de Richter et al. [20] have stated that the
destabilization activity of the grain packing depends
strongly on the history of the system. Information on that
past history, and in particular the packing fabrication, is
stored in the initial network of the physical contacts that
exist between the grains. To study the influence of the sys-
tem history on the precursor activity, the authors applied
a series of consecutive forth-and-back tilting cycles: the
box was slowly inclined from the horizontal to an "ageing"
angle smaller than the avalanche angle and then inclined
back to the horizontal. The authors observed a decrease of
the overall precursor activity during the first cycles. Af-
ter four cycles, the activity reaches a limit state that cor-
responds to a strongly reduced intensity. It is worth not-
ing that, using acoustical methods, the authors observed
‘mirror-type’ rearrangements in the bulk of the pile. The
ascending and descending parts of the cycles form a kind
of anti-symmetric hysteresis loop which remains during up
to 15 cycles. This behavior has motivated the present work
where we implement a new experimental protocol. The

<footnote>*e-mail: luc.oger@univ-rennes1.fr</footnote>

<footnote>**e-mail: renaud.delannay@univ-rennes1.fr</footnote>

<footnote>***e-mail: yves.legonidec@univ-rennes1.fr</footnote>
specific protocol consists in continuously tilting the packing with successive oscillation cycles between $\theta_C < \theta_A$ and $-\theta_C$. As shown in the following, this original approach allows providing a large number of precursors that improves the accuracy of the analyses.

In section 2, we present the experimental setup and the techniques used to monitor the surface activity. Section 3 describes the results and presents the advantages of the successive oscillation cycles, defined as 'full-cycle', to characterize the precursor activity.

# 2 Description of the reproducible cycle experiments

## 2.1 Grain packing fabrication

To achieve the so-called "same history" of a packing fabrication [21], we have followed the approach introduced by [20] which consists in a four-steps protocol.

In the first step, we position a 1 cm square mesh grid at the bottom of a rectangular aluminum container with length $L_b=440$ mm, width $W_b=200$ mm and height $H_b=185$ mm. Note that the dimensions of the container are adjustable ($L_b=240, 440, 640$ or $1000$ mm and $W_b=150, 200$ or $300$ mm). In the second step, we fill the container with monodisperse glass beads with a diameter $D = 2$ mm. The number of grain layers is about 90, which is high enough to avoid any influence on the maximal stability angle of the packing, the well-known wall effect in granular packings [22]. Note that this number is much larger than the limit of 10 layers observed by Aguirre et al. [13]. In the third step, we slowly pull out the grid, an operation which results in a nearly loose packing. In the fourth step, we flatten the free granular surface of the packing by pushing horizontally a rigid metallic bar along the $L_b$ direction.

The average initial packing fraction (ratio between the volume of sphere material with the volume of the whole packing) evaluated in [20] is 0.594 ± 0.001 but the last step of the packing fabrication may compact the superficial layer of grains by only few tenths but enough to result in a "scrapped surface" of the packing.

## 2.2 Experimental setup and successive tilting cycles principle

Preliminary developments of an automatic tilting experimental setup have been designed to introduce original results in characterizing precursors and avalanches of granular piles [4, 19, 23]. For the present work, this experimental setup has been significantly improved to reduce artifacts, including mechanical vibrations and electrostatic effects between the glass beads, to consider different sizes of the packing container, to synchronize several physical measurements not discussed here, and to perform large range of inclinations positives or negatives (Fig. 1).

The experimental setup, placed in a humidity and temperature controlled room, is constituted of an optical table managed by four pneumatic attenuators and mounted on damping pads in order to reduce external vibrations that

may modify the dynamics of the packing. Both the aluminum rectangular box that contains the grains, a camera used to scan the free surface and lights used to illuminate the free surface of the packing are fixed on an inclinable plate that can rotate to a maximum angle of 30° around an horizontal axis made of Teflon. The box is positioned so that the horizontal axis is directed according to $W_b$ and is located at the center of the $L_b$ length of the box. The plate inclination is controlled by the use of a linear actuator and the tilting speed can be set between 1.7 and 14 °/mn in a LabVIEW™ program that also records the inclination angle. In the present work, the experiments have been performed at the minimum speed.

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Photography of the experimental setup.</figcaption></figure>

The experimental setup allows performing multi-cycle experiments composed of an initial quarter-cycle $C_0$ from 0° to $\theta_C < \theta_A$, corresponding to a classical tilting experiment with stop before reaching avalanche, followed by a series of successive half-cycles indexed by $i$ from $+\theta_C$ to $-\theta_C$ ($C_i^+$) and from $-\theta_C$ to $+\theta_C$ ($C_i^-$). After $N$ full-cycles, the inclination of the last half-cycle $C_i^-$ increases up to the maximal stability angle $\theta_A$.

# 3 Measurement of the precursor activity

## 3.1 Optical detection of the precursory events

During the experiment, the optical camera takes one image of the free granular surface per second with a resolution of 2448 × 2050 pixels on 8 bits gray level. This total surface is noted $S_0$. Positioned 1 m above the free granular surface and set on an aperture speed of 1/15000 s, the camera allows imaging a surface of 230 × 230 mm² located at the center of the granular surface and the grain 2D projection corresponds to about 300 pixels. The image acquisition is performed by the LabVIEW™ program and is thus synchronized with the inclination angle measurements as a function of the recording time.
Four stripes of LEDs illuminate the packing free surface to produce grain shadows. Changes in the shadows surface between two consecutive images are used to detect grain displacements: an ImageJ [24] script, initially developed by Kiesgen de Richter et al. [2], consists in subtracting two consecutive images, thresholding the result to remove the noise and then quantifying the surface S of the displaced grain shadows as a function of the inclination angle. Based on the ratio S/S₀ hereafter called the surface activity, the detection of a precursory event is similar to the approach used in previous studies [2, 5, 17, 19] and corresponds to a surface activity close to unity.

## 3.2 Precursor activity in the classical range 0° to $\theta_C$

The initial quarter-cycle $C_0$ of our multi-cycle protocol is similar to classical avalanche experiments but the inclination is stopped before the avalanche occurs, at about 28°: just after its fabrication, the grain packing is slowly tilted from the horizontal $\theta = 0°$ to $\theta_C = 24°$. The image processing allows detecting several main peaks of activity which appear quasi-periodically (Fig. 2, blue curve). They can thus be associated to precursory events. The red dashed line is drawn in a semi-logarithm scale to ensure the Nerone et al.'s assumption of an exponential increase [1]. Obviously, this assumption is only valid for the very first and weak amplitude precursory events. At larger inclination angles ($\theta > 18°$), the precursor activity is getting higher but increases more slowly. It is difficult to know, without a detailed analysis, if the first small precursors, already observed by Nerone et al. [1], are of the same nature as the large ones that follow them.

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Surface activity measured for a 2 mm beads packing and $\theta_C = 24°$ during the initial quarter-cycle $C_0$. The dashed line is an exponential approximation of the activity behavior.</figcaption></figure>

Weak events can also be observed between the precursor events (Fig. 2) and are discarded in the present study where we chose to identify precursory events when at least half of the free surface of the grain packing has moved. We thus characterize the dynamic of the precursor activity after the angle of the first precursor defined by $S/S_0 > 0.5$. In the present $C_0$ experiment which starts at $\theta = 0°$, we measure a variation $\theta_0 \approx 16.5°$ of the inclination angle before the first large precursor and detect $M_0=5$ precursors. Variations in $\theta_0$ and $M_0$ may exist from one experiment to another not only because of the use of different grain packings but also because of the noise in

the image processing. When only considering $C_0$ experiments to characterize the precursor activity with the angle $\Delta\theta$ that separates two consecutive precursors, a large number of realizations is required to increase the total number of precursory events and get more representative statistical analyses of the precursor activity. Such an approach using different initial packings has strong limitations, in particular when we aim at quantifying the effects of external conditions (humidity, temperature, inclination speed, angle $\theta_C$, ...) on the precursor activity.

## 3.3 Multi-cycle precursor activity from $\theta_C$ to $-\theta_C$

This section presents the first results of the multi-cycle experiment composed of $N=4$ full-cycles with $\theta_C = 24°$ (Fig. 3a).

<figure><img src="image_4.png" /><figcaption>**Figure 3.** (a) Surface activity $S/S_0$ measured as a function of the recording time during a multi-cycle experiment composed of 4 full-cycles (blue line). During the first quarter-cycle and each following half-cycle, the inclination angle $\theta$ changes linearly with the time (red lines). (b) Activity of the $C_1^+$ half-cycle following $C_0$ as a function of the inclination angle.</figcaption></figure>

The experiment starts with the initial quarter-cycle $C_0$ discussed in the previous section, represented here as a function of the acquisition time. It is followed by the first half-cycle $C_1^+$ (Fig. 3b) during which the first "large" precursor appears at $\theta \approx 8°$. The variation of inclination angle necessary to observe this first large precursor is thus $\theta_1^+ = \theta_C - 8° \approx 16°$, similar to $\theta_0$ measured in the previous case $C_0$. Then, we observe a long series of precursors whose activity increases up to the opposite extreme angle $-\theta_C$. Interestingly, the activity increases rapidly for the weak activity precursors preceding $8°$ and then linearly at a much slower rate for the precursors of stronger surface activity, after $8°$. This is symmetric to what we observed in figure 2 and confirms the mirror-type effect observed in [20]. In this first half-cycle $C_1^+$, with $\theta_1^+ \approx 16°$, the number of precursors $M_1^+ = 16$ is already much larger than $M_0$. In the second half-cycle $C_1^-$, we observe a symmetric behavior with $\theta_1^- \approx 19°$, and a number of precursors $M_1^- = 13$. So, in practice, we can average inter-precursor angles $\Delta\theta$ on each group of precursors during $C_1^+$ half-cycle or the summation of the inter-precursors angles during each half-cycle up to the given full-cycle $N$.

As a result, the total number $M_N = \sum_{i=1}^N (M_i^+ + M_i^-)$ of precursory events detected for a single grain packing can be very high. The results show a good reproducibility of the precursor activity after the initial quarter-cycle $C_0$ and no more than two full-cycles are required to get a representative grain packing behavior. Thanks to the robust statistical approach allowed by the new multi-cycle protocol,
only few full-cycles are enough to obtain a good accuracy of the results. In practice, four full-cycles following the $C_0$ quarter-cycle will be used in our future work.

<figure><img src="image_5.png" /><figcaption>**Figure 4.** (a) Evolution of the inter-precursor angle $\Delta\theta$ averaged (<>on each half-cycle $C_i^{\pm}$ (disk with its error bar) and accumulated average from 0 up to $n$ (open diamond) versus the number $n$ up to 23 full-cycles. (b) Evolution of the accumulated data error bar versus $M_n = \sum_{i=1}^{n}(M_i^{+} + M_i^{-})$, the dashed line corresponds to the fit for the latest values.</figcaption></figure>

Figure 4 demonstrates that, after the first full-cycles, we can observe a quasi-stationary of the precursor measurements for the average inter-precursor angle $\Delta\theta$ between two consecutive precursory events. We can notice also the accuracy increase with the increase of $M_n$ after $n$ full-cycles.

## 4 Conclusion

The original multi-cycle experimental protocol designed to study precursor activity of grain packings allows us to provide precursor measurements based on large sample statistical inferences. In complement - not shown here - this quasi-stationary state already obtained after two full-cycles has also permitted us to demonstrate that this setup protocol don't seem to depend much on the initial packing fabrication. With only two full-cycles, we have highlighted this quasi-stationary state and, for accuracy and reproducibility improvement, four full-cycles will be used as standard protocol in our future experimental works. This new approach allows generating accurate and reproducible results required to measure and analyse in more details weak variations induced by changes in physical parameters such as humidity, packing geometry and inclination speed that motivate future works on grain packing precursor events.

**Acknowledgments** We would like to thank P. Chasle for all the developments of our complex LabVIEW™ code and Y. Robert and G. Pecheul for the mechanical realizations.

## References

[1] N. Nerone, M. Aguirre, A. Calvo, D. Bideau, I. Ippolito, Phys. Rev. E **67**, 011302 (2003)

[2] S. Kiesgen de Richter, G. Le Caër, R. Delannay, Journal of Statistical Mechanics: Theory and Experiment **2012**, p04013 (2012)

[3] M. Duranteau, Ph.D. thesis, Université Rennes 1 (2013)

[4] R. Delannay, M. Duranteau, V. Tournat, Comptes Rendus Physique **16**, 45 (2015)

[5] N. Nerone, M. Aguirre, A. Calvo, I. Ippolito, D. Bideau, Physica A **283**, 218 (2000)

[6] N. Nerone, S. Gabbanelli, Granular Matter **3**, 117 (2001)

[7] M. Aguirre, N. Nerone, A. Calvo, I. Ippolito, D. Bideau, Granular Matter **3**, 75 (2001)

[8] P. Boltenhagen, The European Physical Journal B **12**, 75 (1999)

[9] S. Courrech du Pont, P. Gondret, B. Perrin, M. Rabaud, Europhysics Letters (EPL) **61**, 492 (2003)

[10] W. Bi, R. Delannay, P. Richard, N. Taberlet, A. Valance, Journal of Physics: Condensed Matter **17**, S2457 (2005)

[11] J. Métayer, Ph.D. thesis, Université de Rennes I (2008)

[12] P. Evesque, D. Fargeix, P. Habib, M.P. Luong, P. Poirion, J. Phys. I France **2**, 1271 (1992)

[13] M.A. Aguirre, N. Nerone, A. Calvo, I. Ippolito, D. Bideau, Physical Review E **62**, 738 (2000)

[14] N. Gravish, D.I. Goldman, Physical Review E **90**, 032202 (2014)

[15] A. Kabla, G. Debrégeas, J.M. di Meglio, T.J. Senden, Europhysics Letters (EPL) **71**, 932 (2005)

[16] A. Jarray, V. Magnanimo, S. Luding, Powder Technology **341**, 126 (2019)

[17] L. Oger, C. El Tannoury, R. Delannay, Y. Le Gohidec, I. Ippolito, Y.L. Roht, I. Gómez-Arriaran, Physical Review E **101**, 022902 (2020)

[18] V. Gibiat, E. Plaza, P. De Guibert, The journal of the Acoustical Society of America **123**, 3142 (2008)

[19] M. Duranteau, V. Tournat, V. Zaitsev, R. Delannay, P. Richard, *Identification of avalanche precursors by acoustic probing in the bulk of tilted granular layers*, in *Powders and Grains 2013* (AIP, 2013), Vol. 1542, pp. 650–653

[20] S. Kiesgen de Richter, V.Y. Zaitsev, P. Richard, R. Delannay, G. Le Caër, V. Tournat, Journal of Statistical Mechanics: Theory and Experiment **2010**, P11023 (2010)

[21] D. Bideau, E. Guyon, L. Oger, in *Disorder and Fracture*, edited by J.C. Charmet, S. Roux, E. Guyon (Springer US, 1990), Vol. 204, pp. 255–268, NATO ASI series (series b: physics) edn., ISBN 13:978-0-306-43576-8

[22] L. Oger, J.P. Troadec, D. Bideau, J.A. Dodds, M.J. Powell, Powder Technol. **46**, 121 (1986)

[23] S. Kiesgen de Richter, Ph.D. thesis, Université de Rennes 1 (2009)

[24] W. Rasband, Imagej, U. S. National Institutes of Health, Bethesda, Maryland, USA, http://imagej.nih.gov, 1997-2018