The role of rolling resistance in the rheology of wizarding quidditch ball suspensions

To elucidate the effect of particle shape on the rheology of a dense, viscous suspension of frictional, non-Brownian particles, experimental measurements are presented for suspensions of polystyrene particles with different shapes in the same solvent. The ﬁrst suspension is made of spheres whereas the particles which compose the second suspension are globular but with ﬂattened faces. We present results from steady shear and shear-reversal rheological experiments for the two suspensions over a wide range of stresses in the viscous regime. Notably, we show that the rheology of the two suspensions is characterised by a shear-thinning behaviour, which is stronger in the case of the suspension of globular particles. Since the shear-reversal experiments indicate an absence of adhesive particle interactions, we attribute the shear thinning to a sliding friction coefﬁcient which varies with stress as has been observed previously for systems similar to the ﬁrst suspension. We observe that the viscosity of the two suspensions is similar at high shear stress where small sliding friction facilitates particle relative motion due to sliding. At lower shear stress, however, the sliding friction is expected to increase and the particle relative motion would be associated with rolling. The globular particles attain a higher viscosity at low shear stress than the spherical particles. We attribute this difference to a shape-induced resistance to particle rolling that is enhanced by the ﬂattened faces. Image analysis is employed to identify features of the particle geometry that contribute to the resistance to rolling. It is shown that the apparent rolling friction coefﬁcients inferred from the rheology are intermediate between the apparent dynamic and static rolling friction coefﬁcients predicted on the basis of the image analysis. All three rolling resistance estimates are larger for the globular particles with ﬂat faces than for the spherical particles and we argue that this difference yields the stronger shear thinning of the globular particle suspension.


Introduction
Non-Brownian suspensions made of relatively rigid particles are ubiquitous in industry (fresh concrete, civil engineering, rocket fuel, etc.) and in natural flows (mud, lava flows, submarine avalanches, etc.).This widespread occurrence has encouraged active research in the past years that has revealed great complexity in the behaviour of these systems, which are usually composed of particles with irregular shape.Notably, it has been shown that even the simplest suspension, a non-Brownian suspension made of relatively rigid, single-sized rough spheres (of radius a) with negligible colloidal forces (no adhesion), suspended in a density-matched (no effect of gravity) Newtonian fluid (of viscosity η 0 ) and sheared in a viscous creeping flow (no inertial effect), can exhibit a rich variety of rheological behaviours.The best known feature is the divergence of shear viscosity, η, when the solid volume fraction, φ, tends to a maximum value known as the jamming volume fraction, φ m .However, the range of complex rheological behaviours can also include the occurrence of a yield stress (Dagois-Bohy et al. 2015;Ovarlez et al. 2015), shear-thinning (Vázquez-Quesada, Tanner & Ellero 2016; Lobry et al. 2019) or shear-thickening behaviours (Barnes 1989;Mari et al. 2014;Guy, Hermes & Poon 2015;Comtet et al. 2017;Madraki et al. 2017;Madraki, Ovarlez & Hormozi 2018;Madraki et al. 2020), normal stress differences, irreversibility under oscillating shear (Pine et al. 2005;Blanc, Peters & Lemaire 2011a), shear-induced microstructure (Gadala-Maria & Acrivos 1980;Blanc et al. 2011aBlanc et al. , 2013) ) and particle migration (Phillips et al. 1992;Snook, Butler & Guazzelli 2016;Sarabian et al. 2019;Rashedi, Ovarlez & Hormozi 2020).
Owing to the complexity already present in the 'simplest system', suspensions made of spheres have been studied extensively for decades.In contrast, the role played by the particle shape has only started to be investigated recently and still suffers from a dearth of experimental data.Yet, many suspensions found in industry and in nature are composed of globular particles, which have an irregular compact form with a global aspect ratio close to 1 (see figure 1).These particles are predominantly convex due to erosion.The present paper describes an experimental work that aims at reducing this deficit by studying the rheology of a viscous non-Brownian frictional suspension made of globular particles (2a ∼ 40 μm) and comparing it with a suspension of spheres made of the same solid material and suspended in the same solvent.For this purpose, some polystyrene (PS) beads have been crushed, while others have not, in order to create two similar suspensions (described in § 2): one made of beads (see the first sketch from the left in figure 1) and the other made of particles with irregular globular shapes (see the third sketch from the left in figure 1).Since the recent works of Le et al. (2023) have shown that the rheology of a suspension depends strongly both on the type of particles and the solvent, it is important to note that both types of PS particles studied in the present paper are separately dispersed in the same suspending liquid (silicone oil).Therefore, the only difference between the two types of suspension studied in the present paper is the solid particle shape and we investigate the role of shape disentangled from other factors.
In the last decade, the central role played by direct solid contact in the flow properties of non-Brownian frictional suspensions has been revealed by Boyer, Guazzelli & Pouliquen (2011), who succeeded in applying a granular paradigm to describe the rheological behaviour of non-Brownian and non-colloidal spheres suspended in a Newtonian fluid in the dense regime, showing the key role played by solid contact interactions between particles, existing thanks to their asperities.Later, using a discrete-element method (DEM)-like approach Gallier et al. (2014) have extensively studied the influence of asperity height, h r , and sliding friction coefficient, μ s , between spheres on the rheology of suspensions.They have notably shown that μ s is a key parameter that governs the flow The rheology of wizarding quidditch ball suspensions (a) ( b) ( c) Figure 1.Schematic of different types of 2D-projected particle shapes (from the left to the right): a simple sphere (disc), a regular polyhedron (polygon) and a globular/crushed particle that possesses flat faces and spherical arcs.
properties of frictional suspensions of spheres in the concentrated regime (φ > 0.40).
Several numerical studies (Gallier et al. 2014;Mari et al. 2014;Wyart & Cates 2014;Peters et al. 2016;Singh et al. 2018) have then shown that μ s changes the value of the jamming volume fraction, φ m .For instance, Seto et al. (2013) and Mari et al. (2014) have shown that the proliferation of frictional contacts is known to be the cause of the discontinuous shear-thickening (DST) observed in highly concentrated suspensions of spheres when the shear stress is high enough to overcome repulsive interactions between particles and push them into contact.As a consequence, the authors have measured, in the case of spherical particles, a decay of φ m from 0.66 to 0.58 when μ s increases from 0 (frictionless case) to 1 (frictional), in qualitative agreement with the experimental values from the literature for frictional suspensions of spheres: φ m ∈ [0.54; 0.62] (Zarraga, Hill & Leighton 2000;Ovarlez, Bertrand & Rodts 2006;Boyer et al. 2011;Blanc, Peters & Lemaire 2011b;Blanc et al. 2018).Later, Peters et al. (2016) numerically found that φ m decreases from 0.7 to 0.56 for the same variation of μ s (0 μ s 1), in quite good agreement with these previous works.Moreover, recent experimental studies have directly measured the values of μ s by atomic force microscopy (AFM) measurements between pairs of PS beads suspended in silicone oil (Arshad et al. 2021;Le et al. 2023).They found that 0.1 μ s 4, which confirms the considered range of the values of μ s in the numerical studies.Shear-thinning is common in viscous non-Brownian suspensions (Gadala-Maria & Acrivos 1980;Zarraga et al. 2000;Dbouk, Lobry & Lemaire 2013;Vázquez-Quesada et al. 2016, 2017;Blanc et al. 2018;Gilbert, Valette & Lemaire 2022) and can have different physical origin, depending both on the physical properties of the suspension and the range of applied shear stress, Σ 12 (the indices 1, 2 and 3 referring to the flow, gradient and vorticity directions, respectively).By studying a non-Brownian suspension made of polyvinyl chloride (PVC) particles suspended in a 1,2-cyclohexane dicarboxylic acid diisononyl ester (DINCH, Newtonian oil), Chatté et al. (2018) have notably proposed the possible existence of two successive regimes of shear-thinning behaviour separated by a shear-thickening regime related to the frictionless-frictional transition.The first shear-thinning regime occurs at small stress, when the suspension remains frictionless since repulsion prevents direct solid particle contacts.This system can be actually seen as a suspension of 'soft' particles, composed of a 'hard core' (of diameter d = 2a) to which a frictionless jacket of thickness, ξ , is added.The gap 2ξ between neighbouring particles is determined by balancing the normal force F N induced by the applied stress with the colloidal repulsive force, f N .When Σ 12 (and therefore the normal force F N between particles) increases, ξ decreases, and so the apparent size of the particles decreases, a app = a + ξ( f N ), inducing a decay of the apparent volume fraction of the suspension and, in fine, a decay of η (Krieger 1972;Maranzano & Wagner 2001a).When the particle pressure increases more and overcomes the repulsive forces (F N f C N ), the particles enter increasingly frequently into direct solid contact thanks to their asperities and the suspension passes from a frictionless state to a frictional one.
Interestingly, Mari et al. (2014) have shown that the onset of this frictionless-frictional transition ('fft') occurs for a critical shear stress (and not a shear rate, γ ): σ fft in ≈ 0.3 × f C N /(6πa 2 ) for spheres, whose value is independent of φ as already observed in many experiments (Bender & Wagner 1996;Frith et al. 1996;Maranzano & Wagner 2001a,b;Lootens et al. 2005;Fall et al. 2010;Larsen et al. 2010;Brown & Jaeger 2012, 2014).The authors have also shown that the stress range over which thickening occurs remains constant.This has motivated us to control the applied shear stress in the present study, instead of the shear rate.Once the load F N is large enough (F N f C N ), the direct solid contacts between particles saturate since all the particles in the suspension have contacts with their neighbours: the system is in the frictional state.Mari et al. (2014) have measured the occurrence of this second regime at σ fft out ∼ f C N /a 2 .In the frictional state, if Σ 12 increases further, then a potential second shear-thinning regime can be observed.We want to emphasise that it is precisely this second shear-thinning regime (when the suspension is frictional) that will be explored in the present paper.The physical origin of this complex behaviour remains an open question.For instance, Acrivos, Fan & Mauri (1994) suggested that the apparent shear-thinning behaviour observed in Couette flow can be due to a difference of density, Δρ, between the solid particles and the suspending fluid.Indeed, solid particles heavier than the suspending fluid settle because of gravity and form a more concentrated layer.Then, shear-induced viscous resuspension (Gadala-Maria 1979;Acrivos, Mauri & Fan 1993;Zarraga et al. 2000;Saint-Michel et al. 2019;d'Ambrosio, Blanc & Lemaire 2021) tends to homogenise the suspension when Σ 12 increases, which induces an apparent decay of the viscosity.However, while this mechanism may arise in some experiments with Couette rheometers, it cannot explain the shear-thinning behaviour observed in other types of flow.For instance, in the case of a parallel plates geometry, the shear-induced viscous resuspension would tend to increase the viscosity.In addition, we show that Σ 12 in the present study is large enough so that gravity would not cause significant deviation from uniform volume fraction, so the effect of any shear rate dependence related to gravity is absent.
Lastly, numerical simulations (Lobry et al. 2019) and experimental studies (Chatté et al. 2018;Arshad et al. 2021;Le et al. 2023) have shown that the shear-thinning behaviour observed for concentrated viscous non-Brownian frictional suspensions (i.e.beyond the DST) could be related to a sliding friction between solid particles that varies with the normal force F N .Following the model from Brizmer, Kligerman & Etsion (2007), Lobry et al. (2019) have considered that the contact between particles is elastic and occurs only through a few hemisphere-like asperities.In these conditions and according to the Hertz theory, the elastic contact area A contact is proportional to where F T denotes the tangential force.This model is in good agreement with experimental works (Chatté et al. 2018;Arshad et al. 2021;Le et al. 2023) which have directly determined the decay of μ s with the normal force F N by conducting AFM measurements between pairs of particles.Arshad et al. (2021) and Le et al. (2023) have conducted AFM measurements to measure the pairwise friction between pairs of PS beads (d ≈ 40 μm) immersed in an aqueous liquid and silicone oil, respectively.Note that the system of suspension studied by Le et al. (2023) is the same as that studied in the present paper.
The different studies (Lobry et al. 2019;Arshad et al. 2021;Le et al. 2023) performed on suspensions of spherical particles have all converged to the following equation based on the works from Brizmer et al. (2007): where L c corresponds to the critical normal force which scales the saturation of μ s .In other words, the sliding friction coefficient becomes constant and equal to μ ∞ s when F N L c , because of an elastic to plastic transition of asperities deformation (Lobry et al. 2019).In the case of a contact between a perfectly smooth half sphere and a flat surface, Brizmer et al. (2007)  2).On the other hand, Le et al. (2023) measured for PS beads in silicone oil: μ ∞ s = 0.15 (m = 0.4).Note that, since the particles of the suspensions studied in the present paper are of the same chemical composition found in these studies from the literature (and even the same solvent for Le et al. 2023), we reuse (1.2) coupled with the latter constants to characterise the shear-thinning behaviour of the studied suspensions.Lobry et al. (2019) have numerically determined the relationship between the normal force applied on spherical particles and the shear stress: F N = 6πa 2 Σ 12 /1.69.Equivalently, a critical shear stress, Σ c , can be defined as L c = 6πa 2 Σ c /1.69, which allows one to obtain the following updated equation for the variable sliding friction coefficient: It is known in granular media that the two possible motions for a particle are sliding (characterised by μ s ) and rolling.The one offering the least resistance will be favoured but both can obviously occur at the same time in a sheared suspension (Estrada, Taboada & Radjai 2008).One can easily understand that the particle shape might have a significant effect on one or even both of these motions, depending on the contact between particles.A decade ago, the numerical simulations of Estrada et al. (2011) in granular media have shown that the way a non-spherical shape provides resistance to rolling can be essentially modelled by approximating the non-spherical particle (like a globular one) by a sphere 'equipped' with an apparent rolling resistance torque, Γ F N (see figure 2).This shape-induced rolling resistance would be therefore characterised by a rolling friction coefficient, μ r , defined from a Coulomb-type law: . Schema of rolling on a plane from the left to the right for a regular octagon (on the left) and a disc (on the right), thanks to a tangential force, F T , applied at the centre of mass G of the particles.The two particles have the same perimeter.The two systems can be considered equivalent if the disc is 'equipped' with a rolling resistance torque, Γ F N , directly related to a 'rolling resistance force', F t N defined from the normal force, F N , applied on G as F t N μ r F N .Here, μ r is the rolling friction coefficient (Estrada et al. 2008(Estrada et al. , 2011)).
of the two studied suspensions, in order to characterise the rheological behaviour of non-Brownian viscous suspensions made of frictional particles with irregular shapes and compare it with the rheology of a basic suspension made of spheres of the same material.In the second part, we then determine by an image analysis process (see § 4) the rolling friction coefficient, μ r , of the studied globular particles in order to compare the numerical predictions of φ m from the literature with our own experimental data.

Suspensions
In this paper, the rheological behaviour of two different non-Brownian viscous suspensions are investigated.The two suspensions are very similar: they are both made of the same PS particles (TS40, Microbeads) with a density measured as ρ p = 1.06 g cm −3 and sieved between 36 and 45 μm in order to reduce the initially large size distribution, dispersed separately in the same solvent, a Newtonian silicone oil (Sigma-Aldrich) of density ρ f = 0.97 g cm −3 and viscosity η 0 = 0.98Pa s measured at T = 23 • C. To prepare a given suspension, a known mass of solid particles is carefully mixed with a known mass of liquid.The air bubbles are then removed by putting the sample in an ultrasound bath.The suspension is finally gently stirred in order to resuspend the particles that would have settled during the degassing procedure.
The only difference between the two suspensions remains in the shape of the PS particles.For the first suspension, labelled S PS40 , the solid particles are spheres and to make the second suspension labelled C PS40 , the PS particles have been crushed by a process described in Appendix A. Figure 3 shows examples of these particles captured with a basic microscope: some spherical particles are presented in figure 3(a) whereas a sample of crushed particles is shown in figure 3(b).One can already note that the population of crushed particles is slightly heteroclyte, being composed of different shapes classified from simple spheres to more facetted particles and particles having both spherical and flat surfaces (see the rightmost schematic in figure 1).It is this appearance, combining spherical arcs and flat surfaces similarly to a quidditch ball (the so-called quaffle), which motivated us to choose the title for the present paper.Figure 4 displays an enlarged image of a sample of crushed PS particles, which allows one to better appreciate this heteromorphism.A quantitative study by image analysis has been conducted over a few hundred images captured with a microscope such as those presented in figure 3 in order to characterise the size distribution of the two types of particles, displayed in figure 5.One can observe that the spherical and crushed PS particles have roughly the same size, and both populations can be considered monodisperse with mean and standard deviation of the diameter of d S PS40 ≈ (42 ± 1) μm and d C PS40 ≈ (43 ± 4) μm.If the crushed particles appear slightly larger than the spherical particles, it is likely due to the fact that the diameter d is calculated from the projected area of the particle.is controlled by the static lower plate and is set at T = 23 • C for all the experiments.The gap is imposed at 1 mm h 2 mm, which allows one to have enough particles (20 h/d 50) to minimise phenomena of layering and sliding.The preference of working in a parallel rotating disc is led by the near absence of shear-induced particle migration in such a geometry (Chow et al. 1994;Merhi et al. 2005), which helps in keeping a homogeneous suspension across the gap.However, the drawback of this geometry is that the shear rate is not constant.Indeed, γ increases from 0 at the centre to γR = ΩR/h at r = R, with Ω the angular velocity of the upper rotating plate.In the case of a non-Newtonian behaviour, this variation can be problematic since the viscosity of the suspension, η, depends on the shear rate, γ .In order to take into account this experimental bias and deduce the correct values of η, we use the well-known Mooney-Rabinowitsch correction:

Rheometry experiments
where η app is the apparent viscosity deduced by the rheometer from the measurements of shear rate at the rim of parallel plates, γR , and applied torque, Γ , We studied the rheological behaviour of each suspension over a wide range of shear stress, Σ 12 ∈ [5, 100] Pa, and solid volume fraction, φ ∈ [0.43, 0.51].Note that we work with a volume-imposed geometry and being sure of the volume fraction φ present in the gap is critical for our experiments.It is very difficult to prepare a proper sample with a known volume fraction when φ is close to the jamming volume fraction, φ m .This could be due to the presence of air bubbles hard to remove, instantaneous shear-induced migration when the sample is poured into the gap or even a yield stress which may prevent the suspension from flowing into the gap by gravity.For these reasons, the maximum value for φ in the experiments was kept at 0.51.For each Σ 12 and each φ (in total, 50 combinations of (φ, Σ 12 )), a shear reversal experiment was performed.We encourage the readers to consult Blanc et al. (2018) for details on the protocol.Briefly, the suspension is simply sheared at a given constant Σ 12 .Once the steady state has been reached (η is constant), the flow direction is reversed while the value of Σ 12 is kept constant.Then, the suspension is sheared in this new direction until the steady value of η is retrieved.For each φ on both types of suspension, a series of shear-reversal experiments was performed on two independent samples.The results in the following correspond to the average of these two independent measurements.Within these conditions, the values of the Péclet and Reynolds numbers characterise the suspension as non-Brownian and its flow as viscous (inertial effects are negligible), respectively: with k B the Boltzmann constant.At the same time, note that the Stokes number is kept small throughout all the experiments: St = ( 1 18 )ρ p d 2 Σ 12 /η 2 < 10 −5 .It is thus expected that only viscous and contact forces govern the suspension behaviour.
The maximum shear stress (Σ 12 = 100 Pa) is set by the occurrence of edge fracture which is expected for a first normal stress of the order of the capillary pressure (Keentok & Xue 1999): N 1 ≈ 2γ air-oil /h, with the surface tension of silicone oil, γ air-oil ≈ 30 mN m −1 .Since the literature shows N 1 0.5Σ 12 , we obtain the following criterion to avoid edge fracture: Σ 12 120 Pa, a value close to experimental observations.On the other hand, the minimum stress (Σ 12 = 5 Pa) is chosen in such a way that the Shield number, denoted Sh, is large enough (Sh 1) to ensure that the particles do not settle due to the slight difference of density, ρ, between the solid and liquid phases, and that a vertical homogeneous suspension is maintained throughout the entire experimental procedure: Σ 12 ρgd 10 2 with ρ = ρ p − ρ f = (0.09 ± 0.02) g cm −3 . (2.4) We want to underline that AFM measurements found in the literature (Le et al. 2023) do not observe any repulsive forces before contact for PS particles (also from Microbeads) in a silicone oil (from Merck, ρ f = 0.95 g cm −3 , η 0 ≈ 20 mPa s at 25 • C), meaning that we are already in the frictional regime for the range of Σ 12 studied in the present paper, and that particles are in contact even when γ → 0. This will be confirmed later by the measured values of φ m and the comparison with the literature (Gallier et al. 2014;Mari et al. 2014;Peters et al. 2016).
To conclude this section on the rheometry, we want to emphasise that the plate surfaces are smooth and we made sure that there was no wall slip phenomenon by measuring the viscosity of the suspensions at the largest volume fraction for different gap size.A viscosity found to be independent of the height of the upper plate indicates that there is no detectable wall slip (Yoshimura & Prud'homme 1988).

Rheological measurements
In this section, we aim to characterise the rheological behaviour of the suspension made of crushed PS particles (C PS40 ) and compare it with our measurements of the rheology of the suspension made of spherical PS particles (S PS40 ), which is more common in the literature.The values of viscosity measured on S PS40 within the explored range of Σ 12 are in quite good agreement with other previous works present in the literature (Blanc et al. 2018;Lobry et al. 2019;Le et al. 2023) and conducted on an identical system (i.e.PS spheres of size close to 40 μm dispersed in silicone oil).It appears in figure 6 that C PS40 exhibits a rheological behaviour which is broadly similar to that which characterises S PS40 .In particular, we observe for both suspensions that: (i) as expected, η r increases with φ for a given Σ 12 ; (ii) η r decreases with Σ 12 for a given φ, qualifying the non-Newtonian behaviour in the range of applied shear stress (Σ 12 ∈ [5-100] Pa) for both suspensions as shear-thinning; (iii) as expected, the decay of η r with Σ 12 is steeper (meaning the shear-thinning behaviour is more pronounced) at large φ.
On the other hand, the primary distinction between the suspensions is that the shear-thinning behaviour is stronger for C PS40 compared with S PS40 for a given φ. Figure 7 displays the normalised difference of relative viscosity between the two suspensions, (η , as function of the applied shear stress Σ 12 .One can then easily observe that the suspension made of crushed particles is more viscous than the suspension made of spherical particles at low shear stress, whereas the viscosities of the two suspension are nearly the same at high shear stress. According to the literature (Coussot & Piau 1994;Schatzmann, Fischer & Bezzola 2003;Sosio & Crosta 2009;Mueller, Llewellin & Mader 2010;Vance, Sant & Neithalath 2015), we can quantify the non-Newtonian behaviour of such suspensions by fitting the , as function of the applied shear stress Σ 12 .Each colour labels the solid volume fraction φ: 0.43 (blue), 0.45 (orange), 0.47 (green), 0.49 (red) and 0.51 (purple).experimental measurements by a power law (coloured straight lines in figure 6): where K and n are the consistency factor and the shear-thinning index, respectively.Their values resulting from the fits of the experimental data in figure 6 are displayed as functions of φ in figure 8.We observe in figure 8(a) that K increases with φ as expected.This reflects the increase of the viscosity with volume fraction.On the other hand, we observe in figure 8(b) that n decreases with φ, which accounts for the more pronounced shear-thinning behaviour at large φ.One can also note that n is systematically smaller in the case of C PS40 at a given φ, which reflects the more pronounced shear-thinning behaviour for the suspension made of crushed particles.More precisely, we observe that the relative variation of n over the range of studied φ is roughly twice as large for C PS40 than for S PS40 ( n/ n ∼ 0.2 for crushed particles whereas n/ n ∼ 0.1 for spheres).Regarding the consistency factor, K, it is interesting to see that apparently K C PS40 ≈ K S PS40 at a given φ.However, any further interpretation of this comparison in K can be difficult since its units are not exactly the same between the two suspensions because From figures 6 and 8, it can be seen that the more pronounced shear-thinning behaviour which characterises the suspension C PS40 compared with the same suspension made of spheres (S PS40 ) results from the observations that η at small Σ 12 .To conclude this section, we discuss why we have not considered the existence of a yield stress for either suspension.It is true that it is more relevant to characterise the rheological behaviour for some non-Brownian suspensions by using the Herschel-Bulkley (H-B) law:  the presence of weak adhesive forces between solid particles which would lead to particle aggregation.Thus, the value of τ c may be understood as the minimum stress required to break these aggregates.Furthermore, it is expected that the crushed particles, which have some flat faces, favour Van der Waals interactions since they offer a much larger contacting surface between particles compared with spheres, leading to a higher yield stress.In view of this, we have also fitted our experimental measurements in figure 6 by (3.2).The results have shown that the impact of the third fitting parameter τ c on K and n is negligible, since we found τ c < 1 Pa for both suspensions and all explored φ.For S PS40 , one can note that this is in good agreement with the works of Le et al. (2023) who measured τ c = 0.3 Pa for a very dense suspension made of PS beads having a size of 40 μm and concentration φ = 0.55 in a silicone oil (the same system as studied in the present paper).
The largest volume fraction studied in the present work being φ = 0.51, one can expect that the values of τ c for S PS40 are even smaller than this value within the range of studied φ.Thus, we can advance with enough confidence that the minimum applied shear stress in our study (Σ 12 = 5 Pa) is at least 10 times larger than τ c for S PS40 and C PS40 .We confirm by some measurements from the shear-reversal experiments that adhesive forces do not play a significant role in the rheological behaviour of the studied suspensions within the applied range of shear stress Σ 12 .

A stress-dependent jamming volume fraction
We want to recall that the shear-thinning regime observed for a frictional non-Brownian suspension is common and has already been observed extensively in the literature for suspensions made of spheres (Gadala-Maria & Acrivos 1980;Zarraga et al. 2000;Dbouk et al. 2013;Vázquez-Quesada et al. 2016, 2017) or even facetted (sugar) particles (Blanc et al. 2018).As explained in the introduction of the present paper, the physical origin of this complex behaviour remains an open question.Some recent works, including an experimental study from Chatté et al. (2018) and numerical simulations from Lobry et al. (2019), have demonstrated that the shear-thinning behaviour for frictional spheres could come from a decay of the sliding friction coefficient, μ s , when the shear stress, Σ 12 , increases, which induces an increase of the jamming volume fraction, φ m (Wildemuth & Williams 1984;Zhou, Uhlherr & Luo 1995;Blanc et al. 2018;Lobry et al. 2019;Gilbert et al. 2022).The introduction of a stress-dependent jamming fraction φ m (Σ 12 ) is thus very useful to describe accurately the complex rheological behaviour of a suspension.points correspond to the experimental data and, for each applied Σ 12 , the variation of the reduced viscosity, η r , with the volume fraction, φ, is fitted by a Maron-Pierce-type law: (3.3) Note that the parameter α 0 in (3.3) is used in order to get an accurate fit of our experimental data.If one were to apply (3.3) over the full range of particle volume fractions, α 0 would need to be 1 in order that η r = 1 when φ → 0. However, this fit only works in the dense regime, typically for φ 0.3 in the case of frictional spherical particles, and hence α 0 can have a value different from 1 in order to describe the variation of η r with φ accurately within this regime (Lobry et al. 2019).In our case, a very good fit for each applied shear stress (plotted as coloured lines in figure 9) is obtained for α 0 = 0.85 for both suspensions, a value not too far from the one used in the original equation of Maron & Pierce (1956): The value chosen here is also in good agreement with the numerical simulations of Lobry et al. (2019) who have found that 0.65 α 0 1 when 0 μ s 2, which are the typical values of μ s for common materials such as PS (Arshad et al. 2021;Le et al. 2023), polymethyl methacrylate, glass and rubber.In figure 9(b), we can observe that it satisfactorily fits the experimental data for crushed particles.
The good fit obtained with this given value of α 0 for both suspensions is not so surprising as it is known that the values of η r when φ → φ m are controlled primarily by the value of φ m (Blanc et al. 2018;Lobry et al. 2019).Figure 10 displays the variation of φ m with Σ 12 , determined by (3.3) with α 0 = 0.85 for the suspensions S PS40 made of PS beads (blue circle) and C PS40 made of crushed particles (orange squares).As expected, we observe that φ m increases with Σ 12 for both suspensions.This increase is larger for C PS40 when compared with S PS40 , which illustrates the more pronounced shear-thinning behaviour for the suspension made of crushed particles.To be precise, we observe that φ within the smaller end of the Σ 12 range whereas φ the largest Σ 12 values.This mirrors the previous observation made from the viscosity comparison between the two suspensions.
It is easy to understand that the values of φ m determined by (3.3) might depend slightly on the value of α 0 , and hence we have also fitted our experimental data of η r (φ) by (3.3) with the two known extreme values of α 0 (Lobry et al. 2019): α 0 = 0.65 and α 0 = 1 (the corresponding curves are not plotted in figure 9 in order to keep the graph clear).The confidence areas plotted in figure 10 for each suspension represent this influence of α 0 on φ m .Thus, one can observe that the values of φ m for S PS40 are between 0.560 ± 0.005 and 0.595 ± 0.009 in very good agreement with Lobry et al. (2019), whereas they are between 0.550 ± 0.004 and 0.595 ± 0.009 for C PS40 .It is quite satisfying that these values of φ m for both types of suspension are globally in very good agreement with the literature when non-Brownian frictional (μ s / = 0) suspensions are considered (Zarraga et al. 2000;Ovarlez et al. 2006;Boyer et al. 2011;Mari et al. 2014;Peters et al. 2016;Singh et al. 2018;Lobry et al. 2019;Singh et al. 2020).Furthermore, it is noteworthy that the primary observations that the jamming volume fraction φ m is smaller in the case of non-spherical particles at low shear stress whereas the rheological behaviours of non-spherical and spherical particles are characterised by the same φ m at high shear stress is not altered by the value of α 0 within its known range.
The rest of the paper focuses on finding a physical mechanism to explain the observed rheological difference between the suspension made of crushed particles and the suspension made of spheres.

Physical origin of the stronger shear-thinning regime for crushed particles
In this section, we want to understand the physical origin of the higher viscosity in the suspension of crushed particles (C PS40 ) for small shear stress, as well as the reason that the viscosity of the two types of suspension are similar when Σ 12 is increased.Since the only difference between the suspensions is the shape of particles present in them, it is obvious that this difference in viscosity is related to it.Two different possible physical origins will be thus investigated.First, we show in § 3.2.1 that it is unlikely that the small The rheology of wizarding quidditch ball suspensions remaining adhesion between particles (which is expected to be stronger for the crushed particles at a given shear stress) explains this observation.Second, we discuss in § § 3.2.2 and 3.2.3whether changes in viscosity can be explained by a variable sliding friction between particles coupled with a rolling resistance of particles related to the particle shape itself.To end this section, we study in § 3.2.4 the rheological behaviour of the suspensions in a frictionless case in order to confirm some assumptions of the considered model.

Shear reversal experiments and absence of adhesion
As mentioned previously in this paper, shear-thinning behaviour of a suspension is common and can have different possible physical origins depending on the studied system (Gadala-Maria & Acrivos 1980;Zarraga et al. 2000;Dbouk et al. 2013;Vázquez-Quesada et al. 2016, 2017;Blanc et al. 2018;Chatté et al. 2018;Lobry et al. 2019;Gilbert et al. 2022).One of them is adhesion.Weak adhesive forces exist between solid particles that would lead to particle aggregation.In this scenario, two main features would appear.First, a suspension would exhibit a yield stress τ c (Brown et al. 2010), which may be understood as the minimum stress needed to break these aggregates.Second, they would exhibit shear-thinning behaviour, related to the fact that increasing Σ 12 would break more and more aggregates, which would produce as a result a decrease in the viscosity of the suspension.Furthermore, this explanation would be suitable to explain the highest viscosity at low Σ 12 for crushed particles while the viscosity of the two types of suspensions (S PS40 and C PS40 ) would tend to be similar at large Σ 12 .As already mentioned previously in § 3.1.1,flat surfaces of crushed particles favour particle adhesion.Potential aggregates are then less likely to be destroyed in C PS40 than in S PS40 when both suspensions are sheared at a given small enough Σ 12 .This would make C PS40 more viscous than S PS40 when Σ 12 is small.In contrast, when Σ 12 τ c , all the aggregates are destroyed by shear, even in the case of C PS40 which then flows similarly to S PS40 .
The first flaw in this explanation has already been presented in § 3.1.1.Indeed, we have seen that the smallest value of Σ 12 that we apply to shear the suspension is at least 10 times larger than τ c .At Σ 12 = 10 Pa (second lowest value of applied shear stress), we have Σ 12 20 × τ c .Yet, a significant difference of viscosity between the two suspensions still remains at large φ, which raises doubt that adhesion could be the main physical origin of the stronger shear-thinning for C PS40 .For instance, at φ = 0.51 and Σ 12 = 10 Pa, η r = (110 ± 5) Pa s for C PS40 whereas η r = (80 ± 5) Pa s for S PS40 (see figure 9), which gives a difference of the order of 30 %.Nevertheless, we understand that this argument about Σ 12 and τ c alone is insufficient to support the statement that the adhesion is not mainly responsible for the more pronounced shear-thinning for C PS40 .It is indeed very difficult to estimate precisely when the adhesive forces can be neglected only from τ c .To go further, we have conducted a series of shear-reversal experiments on both types of suspension by following the procedure of Blanc et al. (2018).
A shear-reversal experiment may turn out to be very interesting.It is a very basic experiment (the suspension is simply sheared at a given Σ 12 in a given direction before the flow direction is reversed whereas Σ 12 is kept constant), characterised by a very specific transient response of η which has been observed in all shear reversal experiments (Gadala-Maria & Acrivos 1980; Blanc et al. 2011a) and in simulations (Ness & Sun 2016;Peters et al. 2016).Figure 11 displays an example of the transient response of η for S PS40 (in blue) and C PS40 (in orange) at Σ 12 = 10 Pa and φ = 0.51.As can be observed, a step-like drop of η occurs just after the shear reversal and the viscosity of the suspension reaches a minimum value, η min , at a strain γ = γ min (γ = 0 corresponds to the moment of reversal).This drop is then followed by a rebound of the viscosity which reaches the steady value, Roughly, η H ∼ η min /η 0 and η C ∼ (η s − η min )/η 0 .We refer readers to the numerical work of Peters et al. (2016) to better understand the physical origin of this result.In brief, the particles in contact tend to separate when the shear is reversed.The microstructure of the suspension is thus broken, which induces the drop of the viscosity.Progressively, the microstructure of the suspension is then rebuilt (mirroring the microstructure before the shear reversal since the flow direction has been reversed), which induces the rebound of η to its steady value.
In the present study, the transient viscosity induced by a shear reversal can be very interesting because, if a stress-dependent particle aggregation occurs, then it should also affect the values of η min and the characteristic strains, γ .Notably, Gilbert (2021) has studied the rheology of a non-Brownian frictional suspension composed of homemade soft PDMS particles (Young modulus, E PDMS = 1.8 MPa E PS ∼ 3 GPa) suspended in Span 80 (Newtonian liquid).By using the JKR theory (Johnson, Kendall & Roberts 1971), the author has observed for this suspension that adhesion plays a role if Σ 12 τ a ≈ 10 Pa.By doing shear reversal experiments, he has then shown (see figures 86-2 of Gilbert 2021) that η Σ 12 <τ a min η Σ 12 >τ a min , and that the characteristic deformation of the transient response for a shear reversal, γ s , was much larger than 10 (γ s ∼ 50 for φ = 0.4 and Σ 12 = τ a in the case of his suspension).
Figure 12  adhesive beads.First, η min at a given Σ 12 and a given φ.Second, η min is weakly dependent on Σ 12 for a given φ in the studied range of applied shear stress.Thus, it is likely that Σ 12 τ a for both suspensions in the present study and that adhesion can then be neglected.We also want to underline that, based on (3.4a,b), we can determine η H ≈ 5-6 for S PS40 concentrated at 45 % from the experimental data (roughly independent of Σ 12 ), which is in very good agreement with numerical simulations from Gallier et al. (2014) which shows η H ≈ η ∞ ≈ 5-6 for a non-Brownian viscous suspensions of (frictionless or frictional) spheres at φ = 0.45.Here η ∞ is the high-frequency dynamic viscosity (Van der Werff & De Kruif 1989).
Note that, in figure 12, the experimental data for η min of each suspension for Σ 12 45 Pa have been then fitted by a power law, based on (3.1) where now K ≡ K min and n ≡ n min , to quantify these observations of η min .The fitting parameters K min and n min resulting from this fit are presented in figure 13.The upper limit for the shear stress considered here for the fit (Σ max 12 = 45 Pa) is imposed due to the poor resolution of the measurement of η min when Σ 12 > Σ max 12 .Thus, the apparent plateau of η min observed at large Σ 12 has no physical meaning.It is an experimental artifact.Therefore, one can clearly see in figure 13 that η min for the suspensions S PS40 and C PS40 are characterised by the same rheology, as we observe that both K min and n min are independent of the considered suspension.In addition, S PS40 and C PS40 are both characterised by a Newtonian behaviour when η = η min : 0.94 n min 1 for both suspensions, when 0.43 φ 0.51 and 5 Σ 12 100 Pa.Moreover, the η min results indicate that hydrodynamic interactions are not significantly affected by particle shape.
Figure 14 displays the experimental values of characteristic strains, γ min (open symbols) and γ 0.5 (closed symbols) for the suspension made of PS spheres (blue) and that made of crushed particles (orange), as a function of φ when Σ 12 = 10 Pa.Whereas γ min corresponds to the accumulated strain from the moment of shear reversal to when η = η min , γ 0.5 is defined as the accumulated strain from the minimum state (η = η min )  to the moment when the viscosity has recovered 50 % of its reversal-induced deficit: The uncertainties of the experimental measurements for the characteristic strains are estimated at ±5 × 10 −2 .The experimental data are also compared with numerical (Pine et al. 2005;Peters et al. 2016) and experimental (Pine et al. 2005) results from the literature.We recall that Pine et al. (2005) have shown that a particle in a non-Brownian suspension subjected to oscillatory shear flow returns to its initial position at each oscillatory cycle consistent with Stokes flow reversibility as long as the strain amplitude does not exceed a critical value, denoted γ c .As explained by Peters et al. (2016), γ 0.5 corresponds to the strain necessary for spherical particles to form a significant amount of solid contacts which would then lead to displacements that violate Stokes flow reversibility.
One can observe that γ min decreases with φ for both suspensions and that γ Analogous to γ min , we observe that γ 0.5 decreases when φ increases, which is in good agreement with the literature.In addition, the experimental values for S PS40 at Σ 12 = 10 Pa (filled blue discs) are well-captured by the numerical simulations from Peters et al. (2016) (μ s = 0.5, h r /d = 10 −2 ).One can also note that γ for a given φ even though both are still of the same order and follow the same trend.This slight difference is interesting for two different reasons.First, as we have seen from the work of Gilbert (2021), having γ (and γ s ∼ 10 as can be observed in figure 11) is consistent with the inference that adhesion forces do not play a predominant role in the rheology of C PS40 compared with S PS40 , within the applied range of Σ 12 .Second, Peters et al. (2016) have explained that the force network is reestablished over a typical strain equal to γ 0.5 during a shear reversal experiment.According to this assertion, it would be a little harder for the particles in C PS40 to rearrange during the transient in order to rebuild the microstructure leading to contact forces (see the works from Peters et al. (2016) for details on the physical mechanism).We think this is related to the shape-induced rolling resistance and it could be interesting to study it using numerical simulations, because shear reversal gives access to the separate hydrodynamic and contact contributions to the stress.More generally, Peters et al. (2016) have studied the influence of μ s and h r /a on the values of characteristic strains and we think it could be interesting to also quantify the role played by μ r , if any, in the transient of a shear-reversal experiment.

Variable sliding friction coefficient
In the previous section, we have seen that adhesion cannot account for the shear-thinning behaviour of the two suspensions and that there is no evidence from shear reversal experiments of stronger adhesion in C PS40 than S PS40 .In this section, we show that, unlike adhesion, a variable sliding friction model allows us to explain the shear-thinning behaviour of the suspensions.
From the numerical works of Mari et al. (2014) and Gallier et al. (2014), it is well known that the jamming volume fraction, φ m , is strongly dependant on the sliding friction coefficient, μ s .In addition, as presented in the introduction of the present paper, the recent literature (Chatté et al. 2018;Lobry et al. 2019;Arshad et al. 2021;Le et al. 2023) relates the shear-thinning behaviour of a non-Brownian frictional suspension to a decay of μ s when the normal force F N between particles (directly proportional to Σ 12 ) increases: We recall that Σ c is a critical value which characterises the elastoplastic transition of asperities deformation (Lobry et al. 2019)   the model (Lobry et al. 2019) considers that the contact between two particles occurs at only one or two asperities (mono-asperity contact) and that the particle asperities are supposed to be close to hemispheres (for which m ∼ 1/3, Brizmer et al. 2007).Recent AFM measurements performed on PS beads (d ≈ 40 μm) suspended in an aqueous liquid (Arshad et al. 2021) or in silicone oil (Le et al. 2023) have given μ ∞ ≈ 0.2, Σ c ≈ 10 Pa and m ≈ 0.5.Figure 15(a) displays the variation of μ s with Σ 12 based on these values (−−).
Our main assumption is that (3.6) can describe the variation of μ s with Σ 12 in both suspensions.In addition to the form of the function, we assume that the values of μ ∞ s , Σ c and m are also identical for both types of particles: spheres and crushed.We understand that this statement is critical but several arguments tend to support it.We recall that sliding friction should depend on the local interaction of two surfaces.In the present study, the same PS particles (in size and material) constitute the two studied suspensions and, even The rheology of wizarding quidditch ball suspensions though the crushing process does change the radii of curvature of particles in some places, we assume that it does not significantly affect the topology of the asperities.Thus, as considered from Peters et al. (2016) for spheres, we assume that the solid contact between particles occurs through only a few asperities even to the crushed particles in C PS40 .In this scenario, the value of m determined by Arshad et al. (2021) and Le et al. (2023) for PS spherical particles (m ≈ 0.5) can be applied for the crushed ones.Moreover, the values of μ ∞ s and Σ c determined for PS beads (μ ∞ ≈ 0.2 and Σ c ≈ 10 Pa) by Arshad et al. (2021) and Le et al. (2023), depending on the properties of the solid particle material (Young's modulus E, Poisson's ratio ν, yield strength Y 0 ) and asperity height h r , can also be kept the same for the crushed PS particles.We want to underline that the assumption that μ ) is also consistent with the results displayed in figure 14 for the characteristic strain, γ min .Indeed, Peters et al. (2016) have shown the role played by μ s on γ min , and the experimental data from the present study tend to show that the values of μ s are between 0 and 1 for the studied suspensions and are very similar between the two.Lobry et al. (2019) have proposed the following phenomenological function φ m (μ s ) relating the jamming volume fraction to the sliding friction coefficient: where φ ∞ m and φ 0 m are specific values of φ m when the particles cannot slide (μ s and when the suspension is frictionless (μ s → 0), respectively.The expression contains a fitting parameter X p . Figure 15(b) redisplays the variation of the jamming volume fraction, φ m , with shear stress, Σ 12 (already shown in figure 10).In this new figure, the experimental data (represented as blue discs for S PS40 and orange squares for C PS40 ) are fitted by the model described by (3.6) and (3.7).Additionally, φ ∞ m and X p are left as free parameters while φ 0 m is set equal to 0.65, in good agreement with the literature when the frictionless (μ s = 0) regime is considered (Gallier et al. 2014;Mari et al. 2014;Gallier, Peters & Lobry 2018;Singh et al. 2018;Le et al. 2023).We show later in the present paper (in § 3.2.4) that it is also in very good agreement with the rheology of S PS40 and C PS40 sheared in the frictionless regime.
Then, one can observe in figure 15(b) that the variation of φ m with Σ 12 deduced from the fit exhibits two plateaus (−−), each located at extreme values of shear stress: the first when Σ 12 10 −1 Pa and the second when Σ 12 10 3 Pa.According to (3.6) (Lobry et al. 2019), the plateau when Σ 12 → +∞ is due to the saturation of μ s (plastic regime) when Σ 12 /Σ c 1 (see figure 15a).The other plateau predicted by the fit when Σ 12 → 0 is explained by the weak influence of μ s on the values of φ m when μ s is larger than 1 or 2, as demonstrated by the numerical works of Mari et al. (2014) Note that the estimated values of φ ∞ m and φ 0.2 m for S PS40 are in very good agreement with the literature (Fernandez et al. 2013;Gallier et al. 2014;Mari et al. 2014;Peters et al. 2016;Lobry et al. 2019;Le et al. 2023).Mari et al. (2014) , as expected.To sum up, we have observed by fitting the experimental data φ m (Σ 12 ) by (3.6) and (3.7) that the shear thinning behaviour of the two studied suspensions (S PS40 and C PS40 ) is induced by the same variable friction law, μ s (Σ 12 ).The main difference between the two is in φ ∞ m , whose value is smaller in the case of globular/crushed PS particles compared with the PS spheres.One can note that X p | C PS40 ∼ X p | S PS40 ∼ 2, which supports the statement about the sliding friction being the same for the two types of particles.Moreover, X p | S PS40 ≈ 2.3 is a value which is in good agreement with the literature (Lobry et al. 2019;Arshad et al. 2021;Le et al. 2023).

Geometry-related rolling resistance
A decade ago, Estrada et al. (2008Estrada et al. ( , 2011) ) have simulated rolling regular polygons and shown that the stress was the same as discs (with the same φ) equipped with a rolling friction coefficient, μ r (see the schema in figure 2).This would then mean that the geometric effect is a rolling resistance which, in the case of equivalent discs, can be obtained with a μ r .
More recently, in the frame of a study characterising the shear-thickening behaviour of suspensions made of hard spheres (for which μ s is kept constant), Singh et al. (2020) have numerically studied the role of torque-activated (or stress-activated) rolling resistance, which can be simply induced by the 'rough' particle shape of particles in real-life suspensions.Note that adhesive surfaces can also induce a resistance to rolling motion but we eliminated this physical origin in § 3.2.1.To this aim, the authors have simulated spherical particles with a rolling resistance characterised by a rolling friction coefficient, μ r .Singh et al. (2020) have then studied the role played by different combinations of μ r and μ s in determining the value of the jamming volume fraction, φ m .As shown in figure 16, which displays their result, Singh et al. (2020) demonstrated interestingly on the one hand that φ m depends weakly on μ r when μ s is small enough (typically, μ s 0.35).For instance, their results show that φ m decreases from 0.62 to 0.60 when μ r increases from 10 −3 (vanishing rolling resistance) to 10 (extremely strong rolling resistance), and μ s = 0.2 (see the blue curve in figure 16).Note that we determined in the present work: φ m = 0.61 ± 0.01 when μ s = 0.2 (see figure 15), which is in very good agreement with this observation.On the other hand, Singh et al. (2020) have predicted that φ m is strongly dependent of μ r when μ s 0.5.For instance, within the same range of rolling friction The rheology of wizarding quidditch ball suspensions coefficient (μ r ∈ [10 −3 − 10]), the authors showed that φ m decreases from 0.57 to 0.36 when μ s = 10 (see the purple curve in figure 16).Typically, this corresponds to the case where sliding is prevented and only rotation can occur (μ r μ s → ∞).In the frame of the present study, this latter result from the literature (Singh et al. 2020) is very interesting since, based on the assumption that μ r | C PS40 > μ r | S PS40 , it can explain the main observation obtained in the previous section: To go further, we plotted in figure 17 the variation of φ m with μ s , determined by our experimental data in the case of the suspensions S PS40 (blue solid line) made of spherical particles and C PS40 (orange dashed line) made of crushed particles.We compared these two different variations with the numerical simulations from Singh et al. (2020) for suspensions of spheres with two values of the rolling friction coefficient μ r = 0.03, and μ r = 0.10 (i.e. the two values of μ r that we determined in figure 16 with the corresponding value of φ ∞ m for the two suspensions).One can observe very good agreement between our experimental results and the numerical simulations from Singh et al. (2020), which tends to confirm the assessment already formulated by Estrada et al. (2008Estrada et al. ( , 2011) ) that the non-spherical globular particles can be approximated as spheres as long as the effect of their shape is reflected by a rolling resistance, characterised by μ r .In the case of the rheological measurements, it is then captured by the value of φ ∞ m .To sum up, the rheology of the suspensions (S PS40 and/or C PS40 ) is solely determined by μ r (induced by the non-spherical particle shape) when whereas it is nearly independent of shape when Estrada et al. (2008) have indeed demonstrated in the frame of a numerical study on granular material that the dominant mode of relative motion at the contacts (sliding or rolling) is that which minimises the coefficient of internal friction.This simply means that the particles prefer rolling if μ r μ s or sliding if μ r μ s .The case where μ r ∼ μ s is obviously more complex since it involves rolling and sliding motion at the same time.Thus, by considering the most extreme case where μ s = 10 in figure 16

Frictionless suspensions made with the same particles
We have briefly studied the rheology of the frictionless case (μ s = 0) of the two suspensions studied in the present paper, by dispersing the same PS particles present in S PS40 and C PS40 in an aqueous solution, labelled AQ0, and shearing the suspensions in a vane tool geometry.The aqueous solution is a mixture of deionised water with a small amount (less than 3wt%) of Triton-X-100 (surfactant, Sigma Aldrich) and sodium iodide.We encourage the reader to see the supplementary material of Madraki et al. (2020) for more details about this experimental procedure.Furthermore, the critical normal load f C N (occurrence of the frictionless-frictional transition) has been measured by AFM measurements by Madraki et al. (2020)   Figure 18 displays the experimental measurements of η r (coloured symbol) for spherical PS particles (blue discs) and crushed PS particles (orange squares) in aqueous solution AQ0, when Σ 12 ≈ 10 −2 Pa (frictionless case: Σ 12 σ fft in ↔ μ s = 0).As expected, the variation of the reduced viscosity η r with the volume fraction φ follows a Maron-Pierce law (coloured solid straight lines in figure 18, see (3.3)) (Peters et al. 2016;Lobry et al. 2019).Analogous to § 3.1.2,the fit using (3.3) has been done for α = 0.85 and a confidence area is displayed according to the fits of the experimental data when α = 0.65 and α = 1.The result of the fit of η s as function of φ gives φ m ≈ 0.66 ± 0.01 for the suspension made of spherical particles (blue solid line) whereas φ m ≈ 0.66 ± 0.02 for the suspension made of the non-spherical particles (orange dashed line).Several observations can be underlined from this result.First, in the case of frictionless spherical particles, the value of the jamming fraction is in good agreement with the literature (Mari et al. 2014;Gallier et al. 2018;Singh et al. 2018) and this confirms that the suspension is frictionless.Second, it confirms our previous choice to have assumed φ 0 m ≈ 0.65 in order to fit the experimental data for φ m (Σ 12 ) of the suspensions C PS40 and S PS40 by (3.6) and (3.7).Finally, φ 0 m ≈ φ m | crushed PS in AQ0 ≈ φ m | spheres PS in AQ0 is consistent with the numerical results of Singh et al. (2020) who found that φ m is independent of μ r when μ s → 0.
We mention that some literature (Donev et al. 2004;Baule & Makse 2014;Kallus 2016) shows that particles having a shape which deviates slightly from spheres are characterised by a random close packing concentration, φ RCP m , which is slightly larger than the known value for spheres, φ RCP m | spheres = 0.64.To the best of the authors' knowledge, there is no strong evidence in the literature showing that the jamming volume fraction φ 0 m (frictionless case: μ s = 0) and φ RCP m have to be equal, even though it is known that both have the same value in the case of spherical particles.In the present paper, we find that φ 0 m | crushed PS ≈ φ 0 m | spheres PS , but it is possible that a slight difference is hidden by the uncertainties.Nevertheless, we want to emphasise that, although the value of φ 0 m in (3.7) plays a significant role in the high-shear-stress regime (Σ 12 10 2 Pa), the shape-induced rolling resistance, related most strongly to φ ∞ m rather than φ 0 m , dominates the low-shear-stress regime (Σ 12 10 1 Pa).
In the second part of the present paper, we describe how we can determine a value of μ r for each type of particles (spheres and crushed), based on image analysis.The goal is to compare these new values with the ones predicted by the combination of the numerical works of Singh et al. (2020) based on shear rheology measurements coupled with the experimental data φ m (Σ 12 ) (see in § § 3.2.2 and 3.2.3)that we recall here: μ

Image analysis study
In this section, we focus on the direct determination of the value of the rolling friction coefficient, μ r , to be compared with the value, μ φ m r , inferred from our rheological measurements and the simulations of Singh et al. (2020) (see figure 16).Nevertheless, we note that the treatment of a non-spherical particle by a single rolling friction coefficient on a sphere is an approximation.It would not be exact for two reasons.One is that the resistance to rolling of the non-spherical particle would be different at different parts of the surface.The other is that the static rolling resistance one needs to overcome to initiate rolling could be larger than the time-averaged dynamic rolling resistance one needs to balance to maintain rolling.This difference was minimised by Estrada et al. (2011) by considering a uniform polygon.We want to determine how well either of these rolling friction coefficients helps to describe a more irregular but still compact particle rolling resistance.

Characterising quantities of particle shape
To the best of the authors' knowledge, a precise measurement of μ r between a pair of particles is much more difficult than the measurement of μ s , which can be done by AFM measurements (Chatté et al. 2018;Hsu et al. 2018;Arshad et al. 2021;Le et al. 2023).It is even more difficult for crushed particles with irregular shapes which require even more statistics.It is common in granular media to determine μ r by letting a particle roll over a slope (Agarwal et al. 2021).But the determination of μ r by this method can be very complicated or nearly impossible for small particles or particles with a large deviation from spherical shape.Because of these experimental limits, we have chosen here to use a novel The rheology of wizarding quidditch ball suspensions Figure 19.Sketch of a crushed particle (blue area), considered as an irregular convex polygon with a centre of mass G.The sides of the polygon/crushed particle are outlined in a darker blue.The vertices connecting each side are displayed as dark points.Here, L ij is the length of a side connecting two successive vertices, V i and V j .The radius a j is the length of the segment connecting G to the vertex V j .We use h ij to denote the height of the centroid G above the segment V i V j .The black point E ij is on the bisector of the segment V i V j and at the same height h ij from the segment V i V j as G.The eccentricity of the particle/polygon related to its side V i V j is defined as the ratio of the length of the vector e ij , denoted |e ij |, to the projected diameter, d. method introduced by Agarwal et al. (2021) and Tripathi et al. (2021), based on image analysis of static grains to calculate the rolling friction coefficient without considering any material properties of the particle.The basic principle of this novel method is to approximate the projected image of a given particle as a polygon that we can characterise by measuring: (i) the aspect ratio, a ratio , defined as the ratio of the longest 'height' (i.e. the length between the centre of mass, G, and a side V i V j ) of the polygon over the smallest one: a ratio = h max ij /h min ij ; (ii) the number of sides, n s ; (iii) the internal angle of each vertex, α j ; (iv) the length of each side, L ij ; (v) the eccentricity associated with each side, |e ij |/d.
A qualitative schema of an irregular polygon is displayed in figure 19 to help visualise the different characterising quantities that we aim to measure.Regarding the vector e ij , we want to underline that e ij = GE ij .As we show later, the parameter e ij , the horizontal component of e ij can be negative or positive depending on the relative positions of E ij , G and V j .The eccentricity is then defined as the magnitude of e ij .

Approximation of particles projected area as an irregular convex polygon
Figure 20 shows four examples of 2D approximations as irregular convex polygons for the images of particles composing the suspensions S PS40 (a,c) and C PS40 (b,d).The basic images are taken with a microscope (examples of basic photos shown in figure 3) with an approximate scale of 80 pixels per particle (projected) diameter.We recall that the projected diameter, d (see figure 5), for a crushed particle corresponds to the diameter of a sphere having the same projected area as the non-spherical particle.Note that, from the start, images with well-separated particles are captured, but if two or more particles are not distinct enough (see figure 3), they are simply not taken into account to compute μ r .Moreover, we want to emphasise that spherical particles, such as those in figure 20(a,c are also present in C PS40 .In the end, the resulting characterising quantities of particles presented above are determined for approximately 600 particles for each type of particle.
The image analysis process is described in Appendix B. The data for the physical particles are compared with results for 10 000 'reference' numerical spheres with similar diameters as the real particles, i.e. in the range 70 2a 90 px.This comparison allows us to examine the effect of the image resolution on the properties of the particles, which in all cases are approximated as polygons.

Image analysis results on characteristic quantities of particles shape
Once the coordinates (V i,x , V i,y ) of each vertex V i for a given polygon/particle are known, all the characteristic physical quantities for a captured polygon/particle (see schema in figure 19) can be determined.In particular, the area A p and the location (x G , y G ) of the The rheology of wizarding quidditch ball suspensions centre of mass G of each particle/polygon are determined as follows: In addition to the size distribution already shown in figure 5 where we observed that crushed and spherical particles have roughly the same size (d ∼ 40 μ m , with a slight larger degree of polydispersity for the crushed particles), figure 21 displays the distribution of the values of characteristic physical quantities determined for the 'reference' perfectly smooth spheres (in green), the real spherical particles in S PS40 (in blue) and the real crushed particles in C PS40 (in orange).One can observe that the particles from S PS40 (in blue) are mainly spheres since the differences from the reference data (in green) are small as indicated by: (i) the aspect ratio of the particles in S PS40 is close to 1 (a S PS40 ratio 1.2 with 90 % of a (iv) the eccentricity for the beads of S PS40 is very small (( The comparison of the 'reference' spheres and the spherical particles of the suspension S PS40 on α j , L ij and |e ij |/d allows us to characterise the slight deviation from perfect spheres, which is much less than the deviation of crushed particles from spherical shapes. At first glance, one can observe that the global shape of crushed particles does not deviate much from a sphere.In particular, Figure 21(a) shows that a ratio < 1.5 for crushed particles, with two-thirds of a C PS40 ratio 1.2.In addition, the crushed particles from C PS40 and the spheres from S PS40 are both globally approximated as polygons notably having: Figure 21.Distribution of the values of (a) polygonal aspect ratio a ratio = h max ij /h min ij , (b) number of segments per particle n s , (c) internal angle α i at each vertex V i , (d) relative length L ij /d and (e) eccentricity |e ij |/d measured across all the particles for the suspensions S PS40 (blue) and C PS40 (orange).For a ratio and n s , the statistics include roughly N ∼ 600 particles for each type of suspension, whereas for α j , L ij /d and |e ij |/d, the computation is done for all the vertices V i of all the polygons: n s × N ∼ 2 × 10 4 data for each suspension.Inset in (c,d,e): Logarithmic y-scale is used to highlight the smallest values of α j and the largest values of L ij /d and |e ij |/d, respectively.The data from real suspensions are compared with the results (in green) obtained by doing the same image analysis process on a numerical image of perfectly smooth spheres that we use as a reference.The rheology of wizarding quidditch ball suspensions However, significant differences between the two types of particles are brought out at the same time by figures 21(a), 21(c) and 21(e).On these three specific graphs, we observe as expected that the crushed particles from C PS40 are indeed characterised by: (i) an aspect ratio a ratio 1.2 for one-third of the particles; (ii) a larger portion of 'small' angles than in the case of particles from S PS40 (∼15 % of α To sum up, all these observations show, in fact, that the shapes of crushed particles in C PS40 do not deviate globally from a sphere.However, a small but non-negligible number of their sides are very different from spherical arcs, likely at least enough to induce the rheological differences between C PS40 and S PS40 observed in § 3.1.More precisely, these different measurements conducted to characterise the shape of particles tend to show that the rheological differences between C PS40 and S PS40 , if related to the particles shape, are mainly due to the three following quantities: a ratio , α i and (|e ij |/d).We show in the next section how these are all connected to each other and to μ r .| C PS40 ≈ 0.10 ± 0.01, which confirms the empirical proposition that the eccentricity can be used to estimate μ r .However, a drawback of this method to calculate μ r is that it is limited to particles whose shape does not deviate strongly from a sphere.For instance, it cannot be applied to regular polygonal particles (Estrada et al. 2011) for which we can expect obviously a higher rolling resistance than spheres despite the fact that their eccentricities are zero.Thus, we will follow and build upon the more fundamental approach of Estrada et al. (2011), in which μ r is derived based on the torque required for rolling which in turn is related to the particle shape parameters.

Determination of the rolling friction coefficient
Figure 22 displays a simple sketch of a crushed particle approximated here as an irregular convex polygon with centre of mass G and number of sides n s , rolling from the left to the right around one of its vertices (that named V j on the schematic in figure 22) as a result of a tangential force F T applied at the centroid G.We consider the conditions such that the irregular convex polygon/particle can only roll (μ s μ r , Estrada et al. 2008).As shown in figure 22, a normal force F N applied at G offers a resistance to the particle's rolling, and rotation occurs if and only if where Γ F T is the torque which tends to make the particle roll, and is defined as Here Γ F N is the rolling resistance torque and is defined as . Schema based on the first sketch drawn in figure 19 of a partial crushed particle (blue area), considered as an irregular convex polygon with a centre of mass G. Here, the vertices V i and V j are the two vertices considered initially in direct contact with another particle.The forces F T and F N are the tangential and normal forces applied at G, respectively.The force f T is the projection of F T on the tangent passing by G of the circle of centre V j and radius a j .In other words, f T is the component of F T that makes the polygon roll from left to right around V j by applying a torque Γ F T .Similarly, f N is the component of F N that offers a resistance for the polygon to roll around V j by applying a torque Γ F N .In this scenario, it is important to understand that the particle motion is from left to right and the particle can only roll (no sliding).In fact, a normal contact force opposing F N and a sliding friction force opposing F T are acting at the contact point V j to prevent it from sliding or moving vertically, but they are not represented here for simplicity.The angle θ ij corresponds to the angle between the vectors GV j and F N .
The forces f T and f N are the parts of the applied forces F T and F N , respectively, which contribute to the corresponding torques, and are defined as (see figure 22) The angle θ ij corresponds to the angle formed by the vectors GV j and F N when the polygon/particle rolls around its vertex V j from left to right.By coupling (4.4) and (4.5a,b), we obtain the following condition for the particle to roll around V j from left to right: Obviously, the value of θ ij evolves during the rotation of the particle and, as a result, so does the force required to make the particle roll.Figure 23 displays a qualitative sketch of the horizontal force F T that must be applied at the centre of mass G as function of the rotation angle ϕ, in order to make an irregular polygon/crushed particle (composed of five sides) roll over its entire perimeter.One can then observe that the resistance for the particle to roll around one of its vertices is locally maximum at the start of the rotation around the given vertex.
According to the literature (Estrada et al. 2008(Estrada et al. , 2011;;Singh et al. 2020), the rolling friction law between two grains of radii a 1 and a 2 defines the maximum torque transmitted by the contact from the rolling friction coefficient μ r as Γ max roll = μ r lF N , with l = a 1 + a 2 .By assuming that a given particle rolls around its vertex V j on a mirror particle in the studied suspensions S PS40 and C PS40 (consistent with suspensions roughly monodisperse and a ratio ∼ 1), we have l = 2a j which then leads to Γ max roll /a j = 2μ r F N .Thus, the applied tangential force F T to roll a sphere equivalent to a crushed particle would have to be greater or equal to Γ max roll /a j (see figure 2).The rheology of wizarding quidditch ball suspensions . Qualitative sketch of the variation of the ratio F T /F N based on (4.6) for an irregular convex polygon with five sides which rolls over its entire perimeter.Each coloured peak labels the rotation of the given pentagon around one of its vertices.Each peak can have a different height because the pentagon is irregular.The height of one peak corresponds to the configuration where the contact with the second particle is 'flat' (for instance, when V i and V j are both in contact with a second particle in figure 22), for which θ ij = θ C ij .Note that this corresponds to the position where the rolling resistance around the corresponding vertex is maximum.Then, after each coloured peak, F T /F N = 0 over an angle, ϕ, with . This is because F N no longer induces a rolling resistance once G has been 'vertically' aligned with the vertex/centre of rotation (θ ij = 0).Thus, it is no longer necessary to apply a force F T to continue rolling around the vertex in the same direction until the next vertex becomes the new contact point/centre of rotation of the particle.
rolling friction coefficient associated with the vertex V j when the particle rolls in a given direction (here from left to right), denoted μ ij r,s , can be then described as (see figure 22) where θ C ij is defined as the maximum possible value reached by θ ij when the particle rolls around a vertex V j in a given direction.As shown in figures 22 and 23, this occurs when the side V i V j of the polygon is in contact with the mirror particle.Moreover, we want to point out that, through the parameter θ C ij , the static rolling friction coefficient related to it, μ ij r,s , depends in fact on the length L ij of the segment V i V j , the parameter e ij and the height h ij (see figure 22) because It is important to note that the value of e ij is directly related to the vector e ij = GE ij , and can be positive or negative depending on the relative x position of E ij , G and V j (see figure 22) and the rolling direction: (i) if GE ij points in the direction opposite to rolling, then e ij < 0; (ii) if GE ij is the rolling direction, then e ij > 0.
Another example which shows the importance of the relative x position of these three points (E ij , G and V j ) is that if G was located to the right of V j in figure 22 (with the particle rolling from left to right), then μ ij r,s = 0 (as qualitatively shown in figure 23) because the force F N applied on G no longer induces a resistance torque.We observe that the rolling resistance is larger when: (i) the length of the side L ij is large; and/or (ii) e ij increases the value of (L ij /2) + e ij ; and/or (iii) the height h ij is small.These observations from (4.8) are actually quite intuitive.For instance, in the simple case of a regular (e ij = 0) 2D polygon with four sides, we can easily imagine that it is harder to roll a cube (tan θ C ij = 1) compared with a long rectangle laying on its small side (on the width, tan θ C ij → 0).On the other hand, the long rectangle on its long side will be much harder to roll (μ ij r 1 when θ C ij → 90 • ).As observed in figure 23, the irregular shape of 2D polygons/crushed particles implies that the static rolling friction coefficient μ r,s associated with a given irregular convex polygon/crushed particle is inhomogeneous in angular space.We have therefore chosen to define μ r,s of a given particle as equal to the maximum value of (4.9) with: The idea behind this choice follows the argument made by Estrada et al. (2008).We consider an irregular polygon laid on its side V i V j on a plane inclined with an angle θ ij .In order to make the polygon roll down the inclined plane (i.e. to change its side in contact with the inclined plane), the angle of the slope must be larger than a critical value: θ ij θ C ij .By rolling (without inertia), if the new critical angle θ C jk (associated with the new side V j V k in contact with the plane) is lower than the previous one (i.e.θ C ij ), then rolling continues.However, if a subsequent segment of the polygon has a higher value of critical angle, the polygon stops rolling.In determining μ r,s , we also consider the maximum resistance between the two possible directions of rotation.
Thus, the static rolling resistance is related to the torque required to initiate rolling (assuming the particle stopped rolling at its most resistant angle).However, it is important to understand that another rolling resistance can be related to the work required to maintain rolling at a constant angular velocity.Both should be important in different parts of a sheared suspension (and at different times at the same location).Analogous to Estrada et al. (2011), this second rolling resistance can be determined by calculating the total work required to roll a non-spherical particle over its entire perimeter, P p , and then balancing it with the total work of an equivalent sphere (of the same perimeter P p as the first one) with a resistance for rolling motion (i.e. a work balance instead of a torque balance).The resistance to rolling motion induced by the particle shape would be then characterised by a dynamic friction coefficient, μ r,d , instead of the static one, μ r,s .In figure 23, it would be then determined from the total (coloured) area under the curve F T /F N , instead of the maximum peak, and one can expect that μ r,d μ r,s .We describe the method of calculating μ r,d analogous to Estrada et al. (2011)  Figure 24.Distribution of (a) the static rolling friction coefficient, μ r and (b) the dynamic rolling friction coefficient, μ r,d , both determined by image analysis for the reference perfectly smooth spheres (in green), the spherical particles composing S PS40 (in blue) and the crushed particles from C PS40 (in orange).

(i) δy
In fact, one can note that μ r,d is related directly to the averaged particle dilatancy, which is in agreement with the work of Estrada et al. (2011).

Image analysis results for the rolling friction coefficient
Figure 24 displays the measured distribution of (a) the maximum static rolling friction coefficient μ r,s and (b) the dynamic friction coefficient μ r,d , based on (4.9) and (4.10), respectively.Before discussing these graphs, we emphasise that our geometrical approach to determine the coefficients μ r,s and μ r,d is correct and independent of the method of measurements or resolution.We point out that the measurements on the 'reference' perfectly smooth spheres (in green in figure 24) give a maximum resolution (i.e.lower limit for the values) on the order of 10 −1 for μ r,s and 10 −2 for μ r,d .These limits for μ r,s and μ r,d are induced by the discretisation process which divides a particle contour into a finite number of segments and could be reduced by increasing the resolution of the images.Nevertheless, the resolution is estimated to be sufficient here because the magnitudes of the experimental data appear to be higher than these limits.Moreover, it is expected that the value of μ r,s should not be affected at all by the resolution limit because it takes into account the maximum rolling resistance, unlikely induced by a spherical part of the particle.It is true that the resolution limit might play a role in the determination of the value of the dynamic rolling friction coefficient: the values of μ r,d are likely overestimated because the rolling resistance related to the spherical part of a particle would tend to decrease the value of μ r,d .This is the reason why the relative uncertainties on μ r,d are quite large.Moreover, as we show later, the value of μ r,d is already lower than the value of the rolling friction coefficient determined previously from φ ∞ m .Figure 24(a) displays the measured distribution of the static rolling friction coefficient, μ r,s .We observe that the values of μ r,s associated with the crushed particles (in orange) are globally larger and more broadly distributed than those associated with the spheres from S PS40 (in blue).This result is quite intuitive because, for a given particle, we consider only the maximum value of μ ij r,s to determine μ r,s .As it is unlikely that the spherical part of a crushed particle is taken into account following this, the difference from spheres is then emphasised.One can even note that the values of μ S PS40 r,s slightly differ from those determined for perfectly smooth spheres (in green in figure 24), which could come from real deformations of spheres composing S PS40 .The averaged value of μ r,s for each suspension is measured to be μ r,s | S PS40 ≈ 0.13 and μ r,s | C PS40 ≈ 0.2.Thus, on the one hand, we have μ Finally, the results of this study show that the globular/crushed PS particle geometry itself is enough to induce the rheological differences observed between C PS40 and S PS40 .To go further, it is quite satisfactory that μ r,d μ φ ∞ m r μ r,s for the two studied suspensions.We think that the experimental method described in the present paper to characterise the resistance to rolling motion induced by particle shape can be considered as another step to estimate the rolling friction coefficient for usage in DEM simulations, because it gives a framework for the value of μ r for real suspensions made of non-spherical particles.We recall that, in agreement with the works of Agarwal et al. (2021) and Tripathi et al. (2021) in dry granular media, a more accurate estimation can be obtained by considering the particle eccentricity defined by e/d.However, two limits of this ratio exist: it cannot be considered for particles with a regular polygon shape (e = 0) or a shape that deviates too much from a sphere.Therefore, the novelty of the present work is then to give a way of estimating a framework of μ r , particularly its upper (static) and lower (dynamic) bounds.

Concluding remarks
In this paper, two different but similar monodisperse suspensions have been sheared in a parallel plate rheometer in order to study their rheological behaviours and characterise their differences.More precisely, the main goal of the present paper was to study the influence of particle shape on the rheology of non-Brownian viscous frictional suspensions.Indeed, the rheology of suspensions composed of spherical particles has been studied extensively in the literature.However, understanding of the rheological behaviour of more complex suspensions composed of particles with irregular shape, which are more common in nature, remains more elusive.We have made two different suspensions composed of the same solid PS particles, separately dispersed in the same suspending Newtonian liquid.The only difference between the two lies in the shape of the particles present in each suspension: spheres in the first and crushed PS particles in the second.
We have then characterised the rheological behaviour of these different types of suspension by studying the variation of the jamming volume fraction with shear stress.Our main result shows that the suspension made of crushed particle is more viscous The rheology of wizarding quidditch ball suspensions than the suspension made of spheres at small shear stress whereas the viscosity of the two suspensions becomes equivalent at large shear stress.This results in a stronger shear-thinning behaviour for the suspension made of crushed particles.This observation is notably reflected by a jamming volume fraction smaller at low shear stress for crushed particles whereas it is of the same order of magnitude as that for spheres when shear stress increases.
To go further, we have tried to understand the physical mechanism behind this observation, obviously induced by the different particle shapes.The literature pointed out the influence of rolling resistance but two different origins could be related to it and have been proposed: changes of adhesive force strength with particles' local curvature and rolling resistance induced by locally normal contact forces acting at a non-spherical particle surface.
We have proposed two arguments which tend to demonstrate that adhesion is not important for the present rheological measurements.The first has been to show that the applied shear stress in the present study is much larger than the yield stress of suspensions (spheres and crushed).The second (and main) argument has been based on conducting shear reversal experiments and measuring the minimum value of viscosity, η min , and characteristic strains.The measurements of the same η min in the two types of suspensions and characteristic strains typical for non-adhesive particles was in contradiction with what could be expected if adhesive forces had played a significant role.
The second explanation relates the shear-thinning behaviour of both suspensions to a variable sliding friction, μ s , whereas the larger viscosity at low shear stress for the non-spherical particles is assumed to be related to the particle shape.This physical origin for the shear-thinning behaviour is notably supported in the case of the PS spherical particles by the literature with the works of Lobry et al. (2019) coupled with the AFM measurements on the same type of particles done by Arshad et al. (2021) and Le et al. (2023).Regarding the crushed particles, we have assumed a variation of μ s with Σ 12 identical to the spherical particles, based on the fact that the crushed particles and the spherical particles are from the same material and of the same size.Moreover, the global aspect of the crushed particles does not deviate much from spheres.AFM measurements performed on the crushed particles could validate this assumption in the future.
The recent numerical work of Singh et al. (2020) has shown that the rolling resistance of solid particles plays a predominant role determining the jamming volume fraction, φ m , when the sliding friction coefficient is large (μ s 0.5), but has almost no effect when μ s is small.A quick comparison on the suspensions made of PS spheres and crushed particles in the frictionless case (PS particles in an aqueous solution) has shown no rheological differences, which is consistent with the absence of impact of μ r when μ s → 0. We have shown that it is possible to fit the variation of the jamming volume fraction with shear stress for both types of suspensions by the same variable sliding friction model (Lobry et al. 2019), simply by predicting a smaller value of the jamming volume fraction for crushed particles when the shear stress tends to zero (i.e.sliding friction coefficient grows 'infinite').The obtained value of φ ∞ m (for which the predominant relative motion is rolling) from the fit (Lobry et al. 2019)  small non-spherical particles, we decided to use an image analysis process consisting of approximating particles as irregular convex 2D polygons to measure the characterising shape parameters such as aspect ratio, internal angle and eccentricity, and finally calculate the values of the static rolling friction coefficient associated with each side of each polygon, μ ij r .The static rolling friction coefficient of each particle, μ r,s , has then been defined as the maximum value of all μ ij r,s characterising each particle.On the other hand, we have also determined the value of the dynamic friction coefficient, μ r,d , from the work needed to roll the particle over a distance equal to its own perimeter, analogous to Estrada et al. (2011).Therefore, μ r,d can be then seen as an averaged value to characterise the shape-induced resistance of a particle to rolling motion based on the whole particle shape.In addition to the fact that the particle geometry of the two studied suspensions is enough to explain the rheological differences between the two, we have notably shown that the calculation of these two coefficient values (μ r,s and μ r,d ) gives a framework to estimate the value of μ r for usage in numerical simulations.
Interestingly, a very good agreement with the recent works of Agarwal et al. (2021) and Tripathi et al. (2021) has been found and we confirmed that the eccentricity, defined as the ratio e/d, gives a very good estimation of the value of rolling resistance for usage in DEM simulations, as long as the particle shape does not deviate too much from a sphere, and that e / = 0. To go further, the next step would be to find a way to directly measure the rolling friction coefficient of the particles (as is done by AFM measurements for μ s ), instead of deducing it by an image analysis process.Other difficulties encountered here concern the irregular shape of crushed particles, and the diversity of irregular shapes, which might invalidate the 2D approximation invoked here and make it harder to characterise the rolling resistance for irregular crushed particles.It is also true that different types of solid contact between the particles can exist in the crushed particles.In this paper, we have only considered a scenario where a particle rolls over a 'mirror' particle, but the physics is probably much more complex and it could be interesting to numerically study the distribution of the types of contact in such a suspension.Being aware of this, we think it could be interesting to compare the numerical results of Singh et al. (2020) with more regularly defined shapes such as cubic particles or other regular polygonal particles (for which e = 0).Small hard fibres (a ratio 2) could also be an interesting shape.One can note that the angularity explored through the crushed particles in the present study remains close to spheres (α j → 180 • ).Studying cubes or rectangular shapes may then be interesting to explore smaller internal angles domains (α j → 90 • ).It would have the second advantage of increasing the value of μ r .Indeed, Singh et al. (2020) have shown that the influence of μ r on φ m is very large when 3 × 10 −2 μ r 3. Cubes present this dual advantage of having a well-defined shape and an expected higher friction coefficient: μ r is expected to be between 0.1 (dynamic) and 0.5 (static), for which Singh et al. (2020) have predicted a much lower jamming fraction, 0.44 φ ∞ m 0.53.Thus, determining φ ∞ m would show whether static or dynamic is more important.A rectangular shape (a ratio 2) offers two very different side lengths and allows one to study further the influence of the angular dependency of μ r .the OpenCV cv2.adaptiveThreshold website: https://docs.opencv.org/2.4/index.html).A rough delimitation of each particle in the picture is thus detected.However, the pixels belonging to the interior of a particle, whose grey level can be similar to the background, can be incorrectly identified as not being part of the particle.As a result, the interior of particles is 'filled' (see OpenCV cv2.floodFill) in order to correct it.The projected particles and their well-defined contours (red pixels in figure 20) are then detected through a watershed segmentation process (Vincent & Soille 1991, see OpenCV cv2.watershed).Finally, the contour of each particle is approximated as an irregular convex envelope (see OpenCV cv2.approxPolyDP, coloured in blue in figure 20) for which the (x, y) coordinates of each vertex V i (black dot except G in figure 20) are known (0 i n s − 1, with n s the number of sides of the polygon).
Appendix C. Theoretical approach to determine the dynamic rolling friction coefficient Let us consider that the centre of mass G travels left to right over a horizontal distance which is given by x G = δx + δx (see figure 26).Under these conditions, the work required to displace the crushed particle over this given distance can be calculated as where y G is the vertical vector displacement of the centroid G for which the force F N exerts a rolling resistance, and its norm is equal to δy = a j − h ij .On the other hand, the work needed to displace a disc, characterised by a dynamic rolling friction coefficient μ r,d , over a distance equal to x G is where δy = a j − h ij , δx j = L ij /2 + e ij and δx j = L jk /2 + e jk (with e ij and e jk positive or negative).Note that, in the case for which the particle rolls around V j from the right to the left, μ r,d is defined as In our case, the polygons are irregular and, thus, the required work to roll around one vertex is not the same for all the vertices of a particle, as is shown qualitatively in figure 23 (the volume of each coloured peak is different).Unlike the static friction coefficient, μ r,s , for which we have considered the maximum resistant torque to rolling motion, we define  22 of a partial crushed particle, considered as an irregular convex polygon with a centre of mass G. Two sides of the polygon/crushed particle, V i V j and V j V k , are outlined in dark blue.Tangential and normal forces, F T and F N (not represented here to avoid overloading the schema), respectively, are applied on the centre of mass G of the particle as shown in figure 22 to create a rolling motion from left to right around V j .Under these conditions, V i is the other vertex (with V j ) considered initially in direct contact with the mirrored particle, whereas V k is that which will be in contact at the end of the rotation around V j .The initial (prime) location of each point is displayed with a black point and annotated ( ), whereas the final (second) position is displayed with a transparent dot and annotated ( ).Analogous to figure 22, α j corresponds to the internal angle of the polygon/particle at V j (α j = V i V J V k ), θ ij and θ jk are the angles between the segments h ij and GV j , and h jk and GV j , respectively.Here h ij and h jk are the heights of G from the segment V i V j and V j V k , respectively.The eccentricity for the sides V i V j and V j V k corresponds to the lengths |e ij |/d and |e jk |/d, respectively.The angle ψ (respectively, ψ ) is the dilatancy angle when the particle rolls from the left to the right (respectively, from the right to the left) around V j .
μ r,d by considering the whole particle.The sum of the work for rolling over all vertices W p | j (or W p | j ) corresponds to the total work needed to displace a given particle over a distance equal to its perimeter P p , and the value of the dynamic friction coefficient μ r,d associated with the given particle can then be determined by the equation ) where: (i) δy j = a j − h ij if a i < L 2 ij + h 2 ij ; (ii) else δy j = 0.
determined: μ ∞ s = 0.27 and m = 0.35, whereas Lobry et al. (2019) estimated L c = 20 nN based on the material properties (PS particles).More recently, Arshad et al. (2021) directly measured μ ∞ s = 0.18 by AFM measurements and determined L c = 33.2nN and m = 0.54 by fitting their experimental results obtained for PS particles in an aqueous liquid by (1.

Figure 5 .
Figure5.Size distribution of spheres (blue) and crushed particles (orange).For a crushed particle, the projected area, denoted A p , is measured by microscopic image analysis.The diameter d corresponds to the diameter of a disc having the same area as the projected crushed particle: d = 2 × A p /π.

974Figure 7 .
Figure 7. Variation of the measured relative steady viscosity difference between C PS40 and S PS40 , η C PS40 r − η S PS40 r , normalised by the measured relative steady viscosity for the suspension made of spherical PS particles, η S PS40 r

Figure 8 .
Figure 8. Variation of (a) the consistency factor K and (b) the shear-thinning index n with solid volume fraction φ for the suspensions S PS40 made of spherical PS particles (blue discs) and C PS40 made of crushed PS particles (orange squares), deduced from (3.1).

Figure 10 .
Figure 10.Variation of the jamming fraction φ m with shear stress Σ 12 for the suspensions S PS40 made of spheres PS (blue circle) and C PS40 made of crushed PS particles (orange square), deduced from figure 9 and (3.3) with α 0 = 0.85.The coloured area for each suspension is related to the possible range of α 0 : 0.65 α 0 1.

Figure 11 .
Figure 11.Example of transient viscosity response as a function of the accumulated strain γ during a shear reversal experiment for S PS40 (in blue) and C PS40 (in orange) with Σ 12 = 10 Pa and φ = 0.51.γ = 0 corresponds to the moment when the flow is reversed.Inset: Enlarged view to better visualise the minimum value of viscosity, η min reached during the transient when γ = γ min .

Figure 12 .
Figure12.Variation of the measured relative minimum value of the viscosity, η min /η 0 , with applied shear stress Σ 12 for (a) the suspension made of spherical PS and (b) the suspension made of crushed PS particles.The colour code labels the solid volume fraction φ of the suspension in the same way as in figure6.For each given φ, the experimental measurements (coloured symbols) for Σ 12 45 Pa are fitted by a power law (coloured straight line): Σ 12 = K min × γ ν min with γ = Σ 12 /η min .The parameters resulting from these fits are shown in figure13.

Figure 13 .
Figure 13.Variation of (a) the consistency factor, denoted K min , and (b) the shear-thinning index, denoted n min , with the solid volume fraction φ for the suspension made of spherical PS particles (blue discs) and the suspension made of crushed PS particles (orange squares), when η = η min and Σ 12 45 Pa.
φ.A closer examination shows that γ C PS40 min is nearly equal to γ S PS40 min at the highest volume fraction (φ = 0.51) whereas it is smaller than γ S PS40 min at smaller volume fractions with the largest difference occurring φ = 0.43, the lowest volume fraction studied.Furthermore, one can note that the experimental data are well-predicted by the simulations fromPeters et al. (2016), conducted on a non-Brownian suspension of frictional spheres, characterised by a combination of a sliding friction coefficient, 0 μ s 1, and a relative roughness height, h r /d = 10 −2 .

Figure 14 .
Figure 14.Characteristic strains γ min for the minimum viscosity (open symbols) and γ 0.5 for the partial recovery (filled symbol) as a function of volume fraction φ.The graph on the right is just an enlarged view of the left graph.The experimental measurements (Σ 12 = 10 Pa) from this study are represented for the suspensions S PS40 (blue discs) and C PS40 (orange squares).From the literature,Peters et al. (2016) have performed several numerical simulations to study the influence of the sliding friction coefficient μ s and the relative roughness h r /d on γ min and γ 0.5 .Some of their results are also plotted (black symbols) on the graph for h r /d = 10 −2 (empty symbols for γ min and filled symbols for γ 0.5 ): ( ) μ s = 0; (•) μ s = 0.5; (∇) μ s = 1.Some experimental measurements ( ) fromPine et al. (2005) of the critical strain γ c for which irreversibility occurs are also plotted, as well as the power law (−) resulting from their numerical simulations: γ c = 0.14φ −1.93 .

Figure 15
Figure 15.(a) Variation of the sliding friction coefficient μ s with Σ 12 (see (3.6)) as computed by Brizmer et al. (2007) and Lobry et al. (2019) (• • •) and the variation used to fit the experimental data φ m (Σ 12 ) (−−).The solid part of the curve corresponds to the experimentally studied range of Σ 12 .(b) Jamming volume fraction φ m as function of shear stress Σ 12 .The experimental measurements for S PS40 (blue discs) and C PS40 (orange squares) are fitted by the model obtained by combining (3.6) and (3.7),where μ ∞ s = 0.2, Σ c = 10 Pa, m = 0.5 and φ 0 m = 0.65.Here φ ∞ m and X p are free parameters.The best fit gives φ m = 0.555 and X p = 2.3 in the case of S PS40 .For C PS40 , we determined φ m = 0.536 and X p = 1.8.The two vertical dashed red straight lines on each graph delimit the range of Σ 12 explored experimentally.
and Lobry et al. (2019) determined φ m ≈ 0.56 and φ m ≈ 0.546 when μ s → +∞, respectively, whereas Le et al. (2023) obtained φ ∞ m ≈ 0.55 by studying experimentally the same suspension as S PS40 .Peters et al. (2016) found φ m ≈ 0.61 when μ s = 0.3 and Lobry et al. (2019) determined φ m 625.Regarding C PS40 , one can observe that (rolling mode) and having deduced the values of φ ∞ m for each type of suspension (see § 3.2.2), a value of the rolling friction coefficient μ φ ∞ m r for each suspension can be predicted from the rheological measurements: μ φ ∞ m r | S PS40 = 0.03 ± 0.02 and μ φ ∞ m r | C PS40 = 0.10 ± 0.01.Note that the uncertainty in μ φ ∞ m r for each suspension is due to the uncertainty in the value of φ ∞ m related to the possible range of α 0 (see (3.3)).

974Figure 17 .
Figure 17.Variation of the jamming volume fraction, φ m , as a function of the sliding friction coefficient, μ s .The coloured symbols are deduced from the experimental determination of φ m for S PS40 (blue circles) and C PS40 (orange squares), and μ s is determined from (3.6) with μ ∞ s = 0.2, Σ c = 10 Pa and m = 0.5.The coloured lines are displayed from the variable sliding friction model (see (3.6) and (3.7)) used to fit the experimental data in figure 15 for the suspensions S PS40 (blue solid line) and C PS40 (orange dashed line).The variation of φ m (μ s ) is compared for each suspension with the numerical results from Singh et al. (2020) (open black symbols) for μ r = 3 × 10 −2 (circles) and μ r = 1 × 10 −1 (squares).

Figure 18 .
Figure18.Variation of the shear viscosity, η r , against volume fraction, φ, for two different suspensions: the first made of the same spheres as in S PS40 (blue circles) and the second made of the same crushed particles as in C PS40 (orange squares).Here, the PS particles are dispersed in an aqueous solution (η 0 = 10 −3 Pa s).The two suspensions sheared at Σ 12 = 10 −2 Pa in a vane tool geometry are frictionless (μ s = 0).A fit of the experimental data by (3.3) (0.65 α 0 1) gives φ m ≈ 0.66 ± 0.01 for the suspension made of spheres (blue solid line) whereas φ m ≈ 0.66 ± 0.02 for the suspension made of the non-spherical particles (orange dashed line).

Figure 20 .
Figure 20.Examples of 2D approximation by an irregular convex polygon for particles composing the suspensions S PS40 (a,c) and C PS40(b,d).The sides of a polygon are coloured in blue whereas the vertices are marked as black dots.The centre of mass G of the polygon/particle is also located and orange straight lines connecting the centroid G to each vertex are displayed.The red pixels delimit the contour detected by the segmentation process.Scale: 'diameter' of particle ∼ 80 px.

≈
the number of segments per polygon and the length of the sides are comparable between the spheres from S PS40 and the 'reference' perfectly smooth spheres ( 35 and L ij /d S PS40 ∼ L ij /d ref ∼ 10 −1 ); (iii) the angles are nearly the same ( α j S PS40 ∼ α j ref ≈ 170 • ); the same global angle since 150 • α j 180 • for ∼88 % of α the same average length of polygon sides ( (L ij /d) C PS40 ∼ (L ij /d) S PS40 ∼ 0.1).Moreover, the mean normalised eccentricity (|e ij |/d) C PS40 remains globally small.For instance, approximately 60 % of the sides of polygon for crushed particles are characterised by a ratio (|e ij |/d) 5 × 10 −2 whereas it is 80 % for the spheres of S PS40 .Approximately 85 % of the ratios |e ij |/d are less than 10 −1 for crushed particles, whereas 95 % are less than 10 −1 for the spheres of S PS40 .974 A36-29 https://doi.org/10.1017/jfm.2023.
4.2.1.Theoretical approachStudies of granular media byWensrich & Katterfeld (2012),Wensrich, Katterfeld & Sugo (2014),Agarwal et al. (2021) andTripathi et al. (2021) have shown that an order of magnitude of μ r for usage in DEM simulations can often be obtained by measuring the ratio of the average contact eccentricity e to the projected particle diameter d: μ r ≈ e /d.This ratio is plotted in figure21(e) for the two studied suspensions in the present paper.We have measured |e ij |/d S PS40 ∼ 5 × 10 −2 and |e ij |/d C PS40 ∼ 10 −1 .Interestingly, one can observe that these values are in quite good agreement with the values previously predicted by the combination of the works of Singh et al. (2020) and the determination of φ m (Σ 12 ): μ φ ∞ m r | S PS40 ≈ 0.03 ± 0.02 and μ φ ∞ m r s for the two studied suspensions (recall that we determined μ φ ∞ m r | S PS40 ≈ 0.03 ± 0.02 and μ φ ∞ m r | C PS40 ≈ 0.10 ± 0.01), because the static rolling friction coefficient, μ r,s , should characterise the maximum rolling resistance.Figure 24(b) displays the measured distribution of the dynamic rolling friction coefficient, μ r,d , associated with the crushed particles of C PS40 (in orange) and the spheres of S PS40 (in blue).As expected, the values of μ r,d are smaller than the values of μ r,s for each studied suspension, and we still observe that μ S PS40 r,d is globally smaller and less distributed than μ C PS40 r,d .We found the following averaged values: μ r,d | S PS40 ≈ 0.02 and μ r,d | C PS40 ≈ 0.03.Here, one can note that μ r,d μ φ ∞ m r for each studied suspension.
coupled with the simulations of Singh et al. (2020) allowed us to obtain values of the apparent rolling friction coefficient for both types of suspension: μ φ ∞ m r = 0.03 ± 0.02 and μ φ ∞ m r = 0.10 ± 0.01 for spheres and crushed particles, respectively.The last part of the present paper has been focused on an experimental estimation of the rolling friction coefficient for both types of particles studied in the present paper.Faced with the difficulty of performing a direct experimental measurement for such 974 A36-37 https://doi.org/10.1017/jfm.2023.