Open-access Incipient motion of tire rubber grains bed

Movimento incipiente do leito de grãos de borracha de pneu

ABSTRACT

In the present work, it was verified that the use of rubber grains from tires, obtained by recycling, could be a viable alternative for use in mobile-bed river models. In this way, it was experimentally verified whether the rubber grains met the incipient motion criteria proposed by Shields. So, a mobile-bed of tire rubber grains was placed inside a rectangular channel, with 3 m long, 14 cm wide and with a slope of 0.001 m/m. Water flowed over the mobile-bed with different discharge increments until the incipient motion of the particles at the bottom was verified. Video image analysis was used to establish a criterion for the incipient motion. Two diameter classes of rubber particles were evaluated (3.00 mm to 2.00 mm and 2.00 mm to 1.70 mm). The experimental critical values of the Shields parameter of the rubber particles were satisfactory because the values were within the dispersion range of the experiments by Shields diagram.

Keywords:
Incipient motion; Rubber grains; Physics model

RESUMO

No presente trabalho verificou-se que a utilização de grãos de borracha de pneus, obtidos por reciclagem, pode ser uma alternativa viável para utilização em modelos de leito móvel de rios. Desta forma, foi verificado experimentalmente se os grãos de borracha atendiam aos critérios de movimento incipiente propostos por Shields. Então, foi colocado um leito móvel de grãos de borracha de pneu dentro de um canal retangular, com 3 m de comprimento, 14 cm de largura e inclinação de 0,001 m/m. A água escou sobre o leito móvel com diferentes incrementos de descarga até que se verificasse o movimento incipiente das partículas no fundo. A análise de imagens de vídeo foi utilizada para estabelecer um critério para o movimento incipiente. Foram avaliadas duas classes de diâmetro de partículas de borracha (3,00 mm a 2,00 mm e 2,00 mm a 1,70 mm). Os valores críticos experimentais do parâmetro Shields das partículas de borracha foram satisfatórios porque os valores estavam dentro da faixa de dispersão dos experimentos do diagrama de Shields.

Palavras-chave:
Movimento incipiente; Grãos de borracha; Modelos físicos

INTRODUCTION

According to the ASCE (ASCE Task Committee on Hydraulic Modeling, 2000), the physical modeling method for incipient bottom movement commonly derives from the Shields parameter (dimensionless shear stress), it also is termed particle mobility number (Yalin, 1989), flow intensity (Einstein, 1942), particle Froude number Fr* and densimetric Froude number (Sharp, 1981). After applying a dimensional analysis, there is a relation of the incipient movement of a particle that follows the critical value of the Shields parameter, τ*c, with others parameters as defined by Equation 1,

τ * c = τ o γ s γ f d s = ρ u * 2 γ s γ f d s = f t a n , u * d s ϑ f , l i f t d r a g , s h a p e (1)

Where τo is the shear stress at the bottom, ds is the average sediment diameter, γs is the specific weight of the sediment, γf is the specific weight of the fluid, ρ is the mass density of the fluid, u* is the shear rate, and νf is the kinematic viscosity of the fluid. Julien (2010) mentions that Duboys (1879) was first derived the Shields parameter, and defined its proportionality with tan (ϕ). Shields (1936) defined the relationship to the grain shear Reynolds number, Re*=u*dsϑf.

Later studies have allowed establishing a new relation between the Shields parameter and dimensionless grain diameter (Julien, 2010),

d * = d s G 1 g ν f 2 1 3 (2)

Where G=γsγf

Thus, the grain shear Reynolds number can be replaced by the dimensionless grain diameter, resulting in the in the modified Shields diagram, Julien (2010) proposed the equation, considering the angle of repose (ϕ):

τ * c 0.3 e d * 3 + 0.06 tan ϕ 1 e d * 20 (3)

Therefore, considering the angle of repose in the incipient movement is more important in very angular grains (crushed material) than in naturally rounded grains.

The mobile-bed river models are useful when sediment transport is import, for examples drop structures, local scours, erosion below spillways, locks and dams, reservoir sedimentation, etc. The mobile-bed river modeling process consists in establishing the similarity criteria to four parameters: Froude number; Manning-Strickler coefficient; critical values of Shields parameter; dimensionless particle diameter; referring to the prototype and the model; this means Frp=Frm, np=nm, *c p= *c m and d*p= d* m (the subscripts p and m designate the prototype and the model, respectively). The ds of the model may not follow the geometric scale established, since diameters smaller than 0.22 mm show interactions between flow and sediment, causing changes in the behavior of sediment transport on the model scale (Heller, 2011). To maintain the similarity of dimensionless sediment diameter, it is sometimes necessary to reduce grain density and increase model diameter. The use of lighter materials (low density) and larger diameters depends on the physical and chemical characteristics of the material. Materials present problems (Julien, 2002), for example: polystyrene floats after some time in water, nut shells have little durability, PVC presents hydrophobia, acrylic and river sand have dirt mixed in them, coal has dirt and non-homogeneity of specific density, Bakelite is porous and tends to rot, etc. In the present work, it was verified whether the use of rubber grains from tires, obtained by recycling, could be a viable alternative for use in mobile-bed river models (Martinez et al., 2023; Teixeira et al., 2020). In this way, it was experimentally verified whether the rubber grains met the incipient motion criteria proposed by Shields.

MATERIAL AND METHODS

The material used is crushed tire rubber. It was washed with detergents in tanks to remove impurities from the rubber, then it was left to rest, and part of the rubber that floated in the water was removed, leaving the rubber with more density than water. After drying in open air, the rubber was sieved to obtain particle size range of retained material between the diameters 1.18 mm to 3.00 mm, representing approximately 65% of the material's particle size distribution. In the other particle size ranges, there was not enough material accumulation to make the experiment of the incipient motion. The sand was collected from the reservoir of the Small Hydroelectric Power Plant of Salto de Paraopeba, from the banks, the adduction channel and in deposits formed close to the dam. The samples selected for analysis and particle size classification were those that represented the material from the adduction channel. Afterwards, the material went through the washing and drying process again, and separated into the particle size range between diameters of 0.125 mm to 0.50 mm, obtaining almost 80% of the material. In the other particle size ranges, there was not enough material accumulation, preventing its use in experiments at the incipient movement.

The angle of repose (ϕ) of sand and tire rubber was obtained experimentally using a funnel completely filled with the sample material. Then, the funnel was raised at a constant speed without contact with the deposited material and the funnel exit, thus forming a material deposition cone (Method C1444-00 - American Society of Testing and Materials, 2001; Chik & Vallejo, 2005). This process was performed three times for each particle size class of materials used in experiments. After the formation of the deposition cone, a photographic camera was used. It was adjusted to the horizontal plane of the stacked material. A graduated ruler was positioned orthogonally to the horizontal plane to record and measure results. The rest angle was determined by the software AUTOCAD (2016) as seen in Figure 1.

Figure 1
Determination of the angle of repose: a) design of the deposition cone; b) sand cone class S1.

The experiments were carried out for the rubber particle size classes R1 (1.70 mm to 2.00 mm) and R2 (2.00 mm to 3.00 mm), and for sand particle size classes S1 (0.125 mm to 0.250 mm) and S2 (0.250 mm to 0.500 mm). The particle size class of rubber within the range of 1.18 mm to 1.7 mm was not considered in the experiment, as it presented very low flow rates that were difficult to reproduce in the pumping system.

The tests to determine the incipient movement of the tire rubber particles were conducted in a rectangular section of an acrylic channel with 0.14 m of base, 0.20 m high, and 3 m long. The original slope of the bed in the reservoir prototype was 0.001 m/m. The rubber grains were used to model the movement of the natural sediment. There was a need to keep the rubber particles humidified in containers separated by particle size classes to fill the pores of the material with water, thus avoiding a change in density with the incorporation of air. It is remarkable that, without this technique, it would not be possible to determine the beginning of its movement. To avoid the effects of excavation of the bottom due to the entrance and exit of the flow in the flume, a fixed bottom of approximately 15 cm in length was glued to the entrance and another one was fixed at the bottom with 5 cm in length. Both bottoms had a roughness similar as that of the tested material, and a thickness of 1 cm. At the exit, an acrylic sill was placed at a height lower than the flow sheet to obtain the largest stretch of uniform flow in the canal. The analysis material was disposed along the channel and distributed evenly using a spatula and a trowel. The thickness of the material in the canal bed was 1 cm.

For sand, five experiments were carried out for each of the particle size class (S1 and S2). For rubber, 19 experiments were carried out for the particle size class R1 and 21 experiments for the particle size class R2. In laboratory, the sand mass density of 2,850 kg/m3 and the tire rubber mass density of 1,140 kg/m3 were determined (Associação Brasileira de Normas Técnicas, 1986).

To evaluate the incipient movement of particles, preliminary experiments were carried out to observe flow and sediment behavior along the channel. The results showed that the movement process, for flows much lower than Shields' theoretical flow, already occurred close to the downstream boundary condition. After the increase in flow, the movement of particles occurred from downstream to upstream, quickly and with each increase in flow. However, particle movement was limited to approximately one meter from the downstream boundary condition, and upstream of this point the channel remained without movement. From this stage onwards, with the increase in flow, the movement in these two regions of the channel did not change significantly until a certain flow range, from which movement began in the rest of the channel.

Based on what was previously observed, three regions were established in the channel: Region A, characterizing the zone of influence of the downstream boundary condition; Region B, corresponding to the range for observing and defining the incipient movement of particles, acquiring images and measuring hydraulic depths; and, Region C, subject to the influence of the upstream boundary condition, Figure 2 illustrates the location of these regions for analysis of incipient movement. In addition, in Region B was defined an incipient motion zone, between 90 cm and 165 cm of the total length of the canal, where the flow had already fully developed without the influence of boundary conditions. The video camera used in the tests was positioned at 120 cm from the channel entrance, and it was adjusted so as to form an approximate angle of 30º to 45º with the channel bed, then the camera focused the bottom of the channel in the measure zone of the incipient motion.

Figure 2
Sketch of the channel and location of the regions of influence of the canal, image acquisition, and evaluation of the beginning of movement.

In the zone of incipient movement, a sub-zone called the movement measurement zone was defined, with an area of 280 cm2, in this sub-zone the video images were analyzed. As a criterion for incipient movement, it was established that the frequency of movement between particles must have occur in a time interval of less than 1s. If this criterion was not met, the flow rate was increased until particle movement was close to the established movement frequency criterion. The analysis of the images considering this criterion made it possible to observe the behavior of the incipient movement of the material for each flow rate, as well as predict whether, for flows close to the theoretical value, the entire section of the channel had already started moving, whether the bed still remained stable or whether found himself on the verge of moving.

RESULTS AND DISCUSSIONS

Table 1 shows the values of angle of repose. By using the graph developed by Simons (1957), which relates rest angles of non-cohesive materials, angle of repose of the sand were determined considering it as a very rounded sediment. In the case of rubber, an angular shape was considered.

Table 1
Average rest angles of rubber and sand.

The values obtained for the angle of repose of S1 and S2 were close to the values determined in the graph of Simons (1957). The formation of the pile cone for classes R1 and R2 showed an upper curvature (Figure 3), because the angle of repose values were higher than those determined by the Simons graph.

Figure 3
Cone formation for rubber, classes R1 and R2.

The tables 2, 3 and 4 show the experimental results of the incipient movement for the sand and tire rubber, respectively. Q is the rate flow, yh is the flow height, Rh is the hydraulic radius, τc exp is the critical shear stress at the bottom, τ*c exp is the experimental critical value of the Shield parameter, and slope of 0.001 m/m.

Table 2
Results of the incipient movement of sand S1 (0.125 mm to 0.250 mm) and S2 (0.250 mm to 0.500 mm).
Table 3
Results of the incipient movement of rubber R1 (2.00 mm to 1.70 mm).
Table 4
Results of the incipient movement of rubber R2 (2.00 mm to 1.70 mm).

Experiments demonstrated that determiner the critical values of the Shields parameter may vary significantly depending on what degree of bed mobility is set as the threshold of incipient motion (Shvidchenko & Pender, 2000).

The Figure 4 presents the experimental results of Delft Hydraulics (Rijn, 1989) . natural non-cohesive sediment: curve A. occasional particle movement at some locations; curve B. frequent particle movement at some locations; curve C. frequent particle movement at many locations; curve D. frequent particle movement at nearly all locations; curve E. frequent particle movement at all locations; curve F. permanent particle movement at all locations; curve G. general transport (initiation of ripples). and Shields curve (Shields, 1936). In fact. Shields did not derive a model or an equation (Miedema. 2010). and published his findings as a graph preserving the representation of the data by a wide band. but Rouse (1939, 1949) subsequently fitting a line through the data (Buffington, 1999).

Figure 4
The experimental result of Delft Hydraulics [after Rijn, 1989] and incipient movement results for sand and tire rubber sediments.

The results of class sand. S1 (0.125 mm to 0.250 mm) and S2 (0.250 mm to 0.500 mm). were between curves A and D. these results are expected. as our criterion of incipient movement is almost similar to those adopted for curves A and D.

On the other hand. the results of the rubber R1 (2.00 mm to 1.70 mm) and R2 (2.00 mm to 1.70 mm) were between curves E and G. These results were due to the lower density of the rubber sediment. and experimentally a strong sensitivity of the increase in the movement of the rubber sediments with the increase in shear stress (increase in flow rate) was observed. quickly changing from an occasional movement (similar to curve A) to a general movement (similar to curve G).

Another criterion of incipient movement was adopted by Graf & Pazis (1977) for natural non-cohesive sediment. where in an area the number of moving particles is counted. 1. 10. 100 and 1000. Figure 5. The results for the sand class. S1 (0.125 mm to 0.250 mm) and S2 (0.250 mm to 0.500 mm). were outside curves 1 and 10. One explanation is that the experiments by Graf and Pazt were carried out for coarser sediment. The results for rubber. R1 (2.00 mm to 1.70 mm) and R2 (2.00 mm to 1.70 mm). were between curves 100 and 1000. but it must be considered that the experiment by Graf and Pazt was made with natural sediment and the rubber experiment used a lower density material. so a lower density sediment offers less resistance to the shear stress of flow

Figure 5
The experimental result of Graf & Pazis (1977) and incipient movement results for sand and tire rubber sediments.

Comparing the incipient movement of the rubber sediment with other sediments. Graf (2010) presents materials with different specific gravity and authors. Figure 6: Amber (1.05. Shields.); Lignite (1.27. Shields); Glass beads (2.49. Vanoni); Barite (4.25. Shields); Granite (2.70. Shields); Steel shot (7.90. White); Sand in air (2.10. White); Sand (2.51. White) and de others sand (2.65. Casey. Vanoni. Gilbert. Kramer. and USWES).

Figure 6
Shields’ diagram. Dimensionless shear stress vs Shear Reynolds number (Graf, 2010).

The results plotted in the Shields diagram (Graf, 2010) were compared with the results of sand and rubber sediments. It can be seen that the results of S1 and S2 were below the Shields curve. in the case of the results of rubber sediments (R1 and R2) values are around the Shields curve.

One of the objections to the use of the Shields diagram that the dependent variables appear in both ordinate and abscissa parameters. Depending on the nature of the problem. the dependent variable can be critical shear stress or grain size (Yang, 1996).

A second parameter is used (Vanoni, 1975). that only depends on the grain diameter. d*=dϑ0.1G1gd1/2. Afterwards. Julien (2010) proposed. d*=dsG1gνf213.

Figure 7 is the Shields diagram modified by Julien. The critical values of the Shields parameter for sand (S1 and S2) were below the Shields curve. the critical values of the Shields parameter for rubber (R1 and R2) were slightly above the Shields curve. but between the critical values of Shields and other authors.

Figure 7
Modified Shields diagram (Julien, 2010).

Using Equation 3 (Julien, 2010). the theoretical critical values of the Shields parameter for sand and rubber were calculated (Table 5). The critical value of the Shields parameter for sand presents values within the data cloud (Figure 7). but the critical values of the Shields parameter for rubber presented values (0.036 and 0.0409) inside the range of τ*c.

Table 5
Critical Valore of Shields Parameter by Julien´s equation.

When critical values of Shields parameters and dimensionless particle diameters in the model and in the prototype match (Sharp, 1981). it is possible to establish a scale relation between particle size and particle apparent weight.

ρ s m ρ w a t e r ρ s p ρ w a t e r = d s p d s m 3 (4)

Where dsm = sediment diameter in the model; dsp = sediment diameter in the prototype; ρwater = water density; ρsm = sediment density in the model. in this case rubber – 1.140 kg/m3; and ρsp= density of the natural sediment in the prototype. whose values can vary between 2.000 kg/m3 and 2.900 kg/m3 depending on the origin of the material under study. Considering the model diameter as of class R1 (average diameter of 1.85 mm). the sediment range of the prototype is 0.77 mm to 0.93 mm. For class R2 (average diameter of 2.5 mm). the sediment range of the prototype is 1.048 mm to 1.25 mm.

CONCLUSIONS

The experiments to determine the incipient movement of sand particles (S1 and S2) presented critical values of Shields parameters (dimensionless shear stress) slightly below the Shields and Julien curves. this is due to the incipient movement criterion adopted in this study. which tends to underestimate the critical value of dimensionless shear stress. when considering that incipient movement occurs when there is an interval of one second between the movement of one particle and another.

The experimental critical values of Shields parameters for rubber sediments R1 (3.00 mm a 2.00 mm) and R2 (2.00 mm to 1.70 mm to mm) were within the dispersion range of the experiments by Shields diagram. Then. it was experimentally verified that the rubber grains meet the incipient movement criteria proposed by Shields. Finally. it is possible to confirm the potential use of rubber sediments in the study of sediment transport in mobile-bed river models.

ACKNOWLEDGEMENTS

To the Minas Gerais Research Support Foundation (FAPEMIG) and Companhia Energética de Minas Gerais (CEMIG) for funding the research.

REFERENCES

  • American Society of Testing and Materials – ASTM. (2001). Standard method for measuring the angle of repose of free-flowing mold powders. Method C1444-00. In American Society of Testing and Materials – ASTM. Annual Book of ASTM Standards (pp. 694-695). West Conshohocken: ASTM International.
  • ASCE Task Committee on Hydraulic Modeling - ASCE. (2000). Hydraulic Modeling: Concepts and Practice. ASCE Manual and Reports on Engineering Practice (N° 97). Reston: ASCE.
  • Associação Brasileira de Normas Técnicas – ABNT. (2016). NBR 6457 Amostras de solo - Preparação para ensaios de compactação e ensaios de caracterização Rio de Janeiro: ABNT.
  • AUTOCAD. (2016). For Students and Educators. USA: Autodesk. Free software.
  • Buffington, J. M. (1999). The legend of A. F. Shields. Journal of Hydraulic Engineering (New York, N.Y.), 125(4), 376-387. http://doi.org/10.1061/(ASCE)0733-9429(1999)125:4(376)
    » http://doi.org/10.1061/(ASCE)0733-9429(1999)125:4(376)
  • Chik, Z., & Vallejo, L. E. (2005). Characterization of the angle of repose of binary granular materials. Canadian Geotechnical Journal, 42(2), 683-692. http://doi.org/10.1139/t04-118
    » http://doi.org/10.1139/t04-118
  • Duboys, M. P. (1879). Etudes du régime du Rhône et de l’ation exercée par les eaux sur un lit à fond de graviers indéfiniment affouillable. Annales des Ponts et Chaussées, 5(18), 141-195.
  • Einstein, H. A. (1942). Formula for Transportation. Trans. ASCE, 107, 561-577.
  • Graf, W. H. (2010). Hydraulics of sediment transport (2nd ed.). New York: McGraw-Hill.
  • Graf, W. H., & Pazis, G. C. (1977). Les Phenomenes de Deposition et d’Erosion dans un Canal Alluvionnaire. Journal of Hydraulic Research, 15, 151-165. http://doi.org/10.1080/00221687709499653
    » http://doi.org/10.1080/00221687709499653
  • Heller, V. (2011). Scale effects in physical hydraulic engineering models. Journal of Hydraulic Research, 49(3), 293-306. http://doi.org/10.1080/00221686.2011.578914
    » http://doi.org/10.1080/00221686.2011.578914
  • Julien, P. Y. (2002). River mechanics (Vol. 56). Cambridge: Cambridge University Press. http://doi.org/10.1017/CBO9781139164016
    » http://doi.org/10.1017/CBO9781139164016
  • Julien, P. Y. (2010). Erosion and sedimentation Cambridge: Cambridge University Press. http://doi.org/10.1017/CBO9780511806049
    » http://doi.org/10.1017/CBO9780511806049
  • Martinez, C. B., Saliba, A. P. M., Viana, E. M. F., Tarqui, J. L. Z., Teixeira, E. D., Costa, M. E. F., & Bastos, A. S. (2023). Proposal to use new methods of design and construction in mobile riverbed physical models: the Salto do Paraopeba dam case. RBRH: Revista Brasileira de Recursos Hídricos, 28, e13.
  • Miedema S. A. (2010). Constructing the Shields curve. a new theoretical approach and its applications. See also other papers on the Shields curve. In WODCON XIX Beijing China.
  • Rijn, L. C. (1989). Sediment transport by currents and waves (handbook). Rep. H 461 Delft, Netherlands: Delft Hydraulics.
  • Rouse, H. (1939). Discussion of Laboratory investigation of flume traction and transportation. Transactions of the American Society of Civil Engineers, 104, 1303-1308.
  • Rouse, H. (1949). Fundamental principles of movement. Engineering Hydraulics. In Proceedings of the 4th Hydr. Conf. (pp. 786-804). New York: Wiley.
  • Sharp, J. J. (1981). Hidraulic modeling London: Butterworth & Co.
  • Shields, A. (1936). Application of similarity principles and turbulence research to bed-load movement. Mitt. der Preuss. Versuchsanstalt für Wasserbau und Schiffbau, 26(5-24), 47.
  • Shvidchenko, A. B., & Pender, G. (2000). Flume study of the effect of relative depth on the incipient motion os course uniform sediments. Water Resources Research, 36(2), 619-628. http://doi.org/10.1029/1999WR900312
    » http://doi.org/10.1029/1999WR900312
  • Simons, D. B. (1957). Theory and design of stable channels in alluvial material Ph. D. dissertation. Colorado State University, Colorado.
  • Teixeira, E. K. C., Rinco, A. V., Coelho, M. M. L. P., Saliba, A. P. M., Pinto, E. J. A., & Furtado, L. M. (2020) Methodology for physical modeling of reservoir sedimentation. RBRH: Revista Brasileira de Recursos Hídricos, 25, e27.
  • Vanoni, V. A. (1975). Sedimentation engineering Reston: ASCE. ASCE task Committee for the Preparation of the Manual on Sedimentation of the Sedimentation Committee of the Hydraulics Division.
  • Yalin, M. S. (1989). Fundamentals of hydraulics physical modeling. In R. Martins (Ed.), Recent advances in hydraulic physical modeling (NATO ASI Series Klower). London: Academic Publishers. http://doi.org/10.1007/978-94-009-2344-7_1
    » http://doi.org/10.1007/978-94-009-2344-7_1
  • Yang, C. T. (1996). Sediment transport. Theory and practice New York: McGraw-Hill.

Edited by

  • Editor in-Chief:
    Adilson Pinheiro
  • Associated Editor:
    Edson Cezar Wendland

Publication Dates

  • Publication in this collection
    25 Apr 2025
  • Date of issue
    2025

History

  • Received
    17 Dec 2024
  • Reviewed
    12 Feb 2025
  • Accepted
    20 Feb 2025
location_on
Associação Brasileira de Recursos Hídricos Av. Bento Gonçalves, 9500, CEP: 91501-970, Tel: (51) 3493 2233, Fax: (51) 3308 6652 - Porto Alegre - RS - Brazil
E-mail: rbrh@abrh.org.br
rss_feed Stay informed of issues for this journal through your RSS reader
Report error