# Elasto-inertial dissipation in particle-laden viscoelastic Taylor-Couette flow

Charles Carré, Masoud Moazzen, Tom Lacassagne, Seyed Amir Bahrani

► To cite this version:

Charles Carré, Masoud Moazzen, Tom Lacassagne, Seyed Amir Bahrani. Elasto-inertial dissipation in particle-laden viscoelastic Taylor-Couette flow. Journal of Fluid Mechanics, 2024, 997, pp.A19. 10.1017/jfm.2024.781. hal-04730882

HAL Id: hal-04730882

https://hal.science/hal-04730882v1

Submitted on 10 Oct 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Elasto-inertial dissipation in particle-laden viscoelastic Taylor-Couette flow

Charles Carré, Masoud Moazzen, Tom Lacassagne & S. Amir Bahrani*

October 10, 2024

## Abstract

Many natural and industrial processes involve the flow of fluids made of solid particles suspended in Non-Newtonian liquid matrices, which are challenging to control due to the fluid's non-linear rheology. In the present work, a Taylor-Couette canonical system is used to investigate the flow of dilute to semi-dilute suspensions of neutrally buoyant spherical particles in highly elastic base polymer solutions. Friction measurement synchronized with direct flow visualisation are combined to characterize the critical conditions for the onset of elasto-inertial instabilities, expected here as a direct transition to Elasto-Inertial Turbulence (EIT). Adding low particle volume fraction (≤ 2%, dilute regime) does not affect the nature of the primary transition and reduces the critical Weissenberg number for the onset of EIT, despite a significant decrease in the apparent fluid elasticity. However, for particle volume fractions ≥ 6% (semi-dilute regime), EIT is no longer observed in the explored Reynolds range, suggesting an apparent relaminarization with yet not further decrease in fluid elasticity. Instead, a new regime termed here Elasto-Inertial Dissipative (EID), was uncovered. It originates from particle-particle interactions altering particle-polymer interactions and occurring under elasto-inertial conditions comparable to that of EIT. Increasing particle volume fraction in the semi-dilute regime and, in so, the particle contribution to the overall viscosity, delays the onset of EID similarly to what was observed previously for EIT in lower elasticity fluids. After this onset, a decrease in the pseudo-Nusselt number observed with increasing inertia and particle-to-polymer concentration ratio confirms a particle-induced alteration of energy transfers in the flow.

## 1 Introduction

Suspensions of solid particles in a Newtonian or Non-Newtonian liquid matrix, are frequently encountered in nature and industrial applications (cement, toothpaste, 3D printing material, avalanches, blood flow, to name a few). When flowing at various scales, such fluids demonstrate a wide variety of non-linear dynamic behaviours, arising from inertia, complex rheology of the suspending liquid phase, fluid-particle interactions or particle-particle interactions. They induce a dynamic complexity of the flow, making it difficult to control its stability, mixing, and heat transfers efficiency as well as to predict overall process performance. Despite recent studies (Tanner, 2020; Zhang & Shaqfeh, 2023; Baroudi et al., 2023; Boulafentis et al., 2023, 2024), experimental data is still needed to better understand the hydrodynamic behaviour of complex suspensions and confront numerical models (Kang & Mirbod, 2021; Lin et al., 2023; Song et al., 2021; Kang et al., 2024; Lin et al., 2024).

In order to gather experimental data and deepen understanding on complex fluid flows, it is convenient to consider well-studied canonical flows. In the present work, the canonical Taylor-Couette flow (TCF) is used (Taylor, 1922; Andereck et al., 1986; Fardin et al., 2014). It consists in two concentric cylinders of radii $r_i$ and $r_e$ (with $r_e > r_i$, see figure 1) and height $h$. Here the case considered is that where only the inner cylinder is rotated at rotation speed $\Omega$ and the outer cylinder is fixed. The fluid entrapped in the gap $\delta = r_e - r_i$ is subjected to a shear rate $\dot{\gamma} = \frac{\Omega r_i}{\delta}$. This configuration has been widely used in the recent literature to characterize instabilities and transition to chaos in a wide range of complex fluids (Fardin et al., 2014; Boulafentis et al., 2023). In that context, it was for example found that the presence of non-colloidal neutrally buoyant particles in a viscous Newtonian base fluid alters the nature and onset of the primary bifurcation (Majji et al., 2018; Ramesh et al., 2019), as well as the friction dynamics, as evidenced by the evolution of the torque exerted on the inner cylinder $\tau$ as a function of the Reynolds number $\mathcal{R}$ (Dash et al., 2020; Moazzen et al., 2022; Kang & Mirbod, 2021).

$$ \mathcal{R} = \frac{\rho r_i \Omega \delta}{\mu} = \frac{\rho \dot{\gamma} \delta^2}{\mu} \qquad (1) $$

<footnote>*IMT Nord Europe, Institut Mines Télécom, Univ. Lille, Center for Energy and Environment, F-59000 Lille, France. amir.bahrani@imt-nord-europe.fr</footnote>
<figure><img src="image_2.png" /><figcaption>Figure 1: left) Taylor-Couette cell mounted on an MCR102 Anton Paar rheometer with geometrical parameters symbols and illustration of the imaging system (see details in Moazzen et al. (2022, 2023)) and quasi-steady inner cylinder ramp-up protocol. right) Geometrical cell parameters values.</figcaption></figure>

$\mathcal{R}$ is the non dimensional flow control parameter characterizing the onset of inertial instabilities (see 1) (with $\rho$ is the fluid density and $\mu$ is the suspension dynamic viscosity). On the other hand, particle-free viscoelastic fluids exhibit their own rich set of specific dynamics: low $\mathcal{R}$ transitions to modes that are not evidenced in Newtonian fluids (Boulafentis et al., 2023), transition to well known modes but triggered by elastic or elasto-inertial mechanisms (Larson et al., 1990), or ultimately transition to a form of chaos resembling turbulence at much lower (Dutcher & Muller, 2013; Boulafentis et al., 2023) or even negligible inertia (elastic turbulence, Groisman & Steinberg (1998, 2000, 2004)). The flow control parameter for such behaviours is the Weissenberg number $W = \lambda\dot{\gamma}$, with $\lambda$ the relaxation time of the viscoelastic solution. $W$ compares elastic to viscous forces. Additionally, it is worth defining the elastic number $\mathcal{E}$ (see 2), which quantifies the importance of the elastic behaviour in the given Taylor-Couette geometry.

$$ \mathcal{E} = \frac{W}{\mathcal{R}} = \frac{\lambda\mu}{\rho\delta^2} \qquad (2) $$

In this work, emphasis is placed on fluids of moderate to high elasticity ($\mathcal{E} \gg 10^{-2}$, Dutcher & Muller (2013)) flowing at commensurate inertia, for which viscoelastic properties give rise to elasto-inertial transitions: instabilities that are either due to elastic mechanisms but arising at non negligible $\mathcal{R}$ values, or inertial but modified by elasticity. There, transition to a chaotic flow state labelled elasto-inertial turbulence (EIT) is observed (Groisman & Steinberg, 1996; Boulafentis et al., 2023; Lin et al., 2023). The onset and the pathway towards EIT, as well as its final structure, have been found to depend on various parameters, among which the nature, conformation, and concentration of the additives (most of the time polymers) conferring viscoelastic properties to the fluid. In particular, the dynamic properties of EIT differ from those of inertial turbulence both in term of friction scaling and spectral distribution (Moazzen et al., 2023; Boulafentis et al., 2023, 2024).

The aim of the present paper is to investigate the behaviour of fluids combining both complex features, particles and elasticity, and characterize elasto-inertial transitions in such systems. A first experimental study on this topic in a Taylor-Couette configuration (Lacassagne et al., 2021) indicated that low particles volume fraction may promote faster transition to EIT in terms of $\mathcal{R}$, whereas higher particle volume fractions (more than 5%) were found to delay transition to EIT. This former study employed visualisation to detect transitions and characterize EIT properties for a single moderately elastic base fluids ($\mathcal{E} = O(10^{-1})$, one polymer concentration) and four particles volume fractions. The present work sets at pushing the investigation further and unravel physical mechanisms by 1) coupling flow visualisation with torque measurements in order to get information on the friction dynamics during transitions and flow states, 2) extending the approach to more elastic base fluids ($\mathcal{E} = O(10^{1})$) for which direct transition from the base laminar state to fully chaotic EIT is evidenced (Groisman & Steinberg, 1996; Moazzen et al., 2023), and 3) covering a wider range of particle concentrations ($\Phi = 0, 2, 6, 10, 14\%$).

# 2 Materials and methods

## 2.1 Experimental setup

The experimental set-up is similar to the one used in (Moazzen et al., 2022, 2023), and its geometrical parameters are represented in figure 1.
With the outer cylinder kept at rest, the inner cylinder was rotated with angular velocity $\Omega$, controlled by the rheometer's motor. The maximum inner cylinder rotation frequency achievable bounded the maximum respective $R$ and $W$ experimentally accessible to values $R_{max}$ and $W_{max}$ reported in table 1, depending on each working fluid viscosity and elasticity. Experiments were performed using slow acceleration (or ramp-up) protocol (see figure 1). The difference maximum shear rate was kept constant at $\dot{\gamma}_{total} = 750 \text{ s}^{-1}$ and $\dot{\gamma}$ increased with steps of $2 \text{ s}^{-1}$. At each step, the shear rate was held constant for a period of $t = 30 \text{ s}$. This resulted in non-dimensional acceleration rates of order $(10^{-4})$ at most, well into the quasi steady assumption (Dutcher & Muller, 2013; Boulafentis et al., 2023).

Experiments were conducted with viscoelastic base fluids made of HPAM (Partially Hydrolised PolyAcrylamide, molecular weight $M_w = 15 - 20.10^6 \text{ g/mol}$) at concentration of $\tilde{C} = 300, 400, 500 \text{ ppm}$, dissolved in a solution of distilled water (31% in volume), glycerol (69% in volume) and sodium chloride (2% in mass), and loaded with mono-disperse spherical PMMA (Polyethyl methacrylate) particles of diameter of $d_p = 50 \text{ µm}$ at volume fractions $\Phi = 0, 2, 6, 10, 14 \%$. The particle Reynolds number $R_p = \frac{\rho\dot{\gamma}dp^2}{4\mu}$ was such that $R_p < 0.02$ for all fluids and experiments. The density of solid particles and viscoelastic base fluid were matched ($\rho = 1192 \text{ kg.m}^{-3}$), resulting in neutrally buoyant particle suspensions.

## 2.2 Rheological characterisation

The shear-viscosity of the working fluids was measured as a function of the shear rate $\dot{\gamma}$, using an MCR302 Anton Paar rheometer with a 50 mm diameter parallel-plate geometry (chosen over the conventional cone-plate one to limit inertial and centrifugal effects on particle loaded samples) at a 0.5 mm gap. Measurements confirmed that the working fluids could be considered as Boger fluids (James, 2009; Boger, 1977), their viscosity being shear-rate independent in the considered range (almost two order of magnitudes of $\dot{\gamma}$). This implies that the Reynolds number increased linearly as a function of the applied shear-rate in the Taylor-Couette setup. Averaged viscosity values are reported in figure 2 (left) and in table 1. As expected, the apparent viscosity of the fluid increased with increasing polymer concentration, but most importantly with particle volume fraction (Zhang & Shaqfeh, 2023).

The elasticity of the working fluids was measured using first normal stress differences $N_1$ (Lin et al., 2014; Schäfer et al., 2018; Moazzen et al., 2023), with the same parallel-plate geometry. For a Boger fluid, it is known that

$$N_1 = 2\mu_p\lambda_{N_1}\dot{\gamma}^2 \quad (3)$$

with $\mu_p$ the polymer viscosity contribution ($\mu_p = \mu - \mu_{sol}$, with $\mu_{sol} = 33.3 \text{ mPa.s the solvent viscosity}$), and $\lambda_{N_1}$ the relaxation time. $\lambda_{N_1}$ was estimated for all fluids inverting the previous equation and averaging on the range of $\dot{\gamma}$ over which the power 2 scaling is valid. The nominal elastic number of each base working fluid $\mathcal{E}^0$ was computed using the apparent overall suspension viscosity $\mu$ (reported in table 1) and the elastic time-scale of the corresponding particle-free fluid (each fluid at the same polymer concentration displays the same $\mathcal{E}^0$, reported in table 1, regardless of particle addition). Apparent elastic number values $\mathcal{E}$ were computed for each fluid using directly measured $\lambda_{N_1}$ and $\mu$. Additionally, measurements of apparent extensional time-scales $\lambda_e$ were performed using a slow retraction capillary thinning protocol (Campo-Deaño & Clasen, 2010) on a bespoke setup (analog to Boulafentis et al. (2024)). Thread thinning dynamics of a fluid filament created between two plates slowly separating was fitted in the exponential-thinning regime such that

$$\lambda_e = \frac{-t}{3 \ln \left( \frac{d}{d_{max}} \right)} \quad (4)$$

with $d(t)$ the filament diameter evolving with time ($t$) and $d_{max}$ the maximum (initial) diameter. All rheological data (shear and extension) were consolidated by averaging over at least three and up to six repeated tests. The evolution of apparent time-scale normalized by that of particle free samples for both normal stress and extensional methods is plotted as a function of $\Phi$ in figure 2.

## 2.3 Measurement techniques

### 2.3.1 Flow visualisation

The flow behaviour was probed using a coupling of direct visualisation an torque measurement on the inner cylinder (using the built-in torque sensor of the rheometer). Flow visualisation was achieved by seeding the flow with anisotropic reflective mica flakes, which enable the monitoring of the flow spatial structure (Boulafentis et al., 2023). The light reflected by the mica flakes was recorded using a camera at a frame rate of 70 images per second (see figure 1). From the reflected light signal along a vertical column of the test section, two types of
<figure><img src="image_3.png" /><figcaption>Figure 2: a) Constant shear viscosity (averaged on a 10 to 100 s⁻¹ shear-rate range) for all polymer and particle concentrations. An MCR302 rheometer with parallel-plate 50 mm geometry with a 0.5 mm gap was used. b) Elastic time-scale λ for all fluid samples, normalized by the particle-free measured value, obtained from first normal stress differences (empty markers, same rheometer) or extensional measurements (full markers - bespoke extensional rheometer setup).</figcaption></figure>

<table><caption>Table 1: Rheological viscoelastic suspensions properties and associated dimensionless numbers. W_c = λ_N1 × γ̇_c, R_c = ρδ²γ̇_c/μ and γ̇_c are the respective critical Wiessenberg number, Reynolds number, and shear-rate for the onset of EIT for Φ = 0 or 2%, or new elasto-inertial dissipative (EID) mechanisms (in italic, for Φ ≥ 6%). Superscript 0 indicates quantities computed using the base fluid relaxation time λ⁰ = λ_N1(Φ = 0%). Apparent relaxation time-scales λ = λ_N1(Φ, Ĉ) are used otherwise. The maximum particles Reynolds number is R_p^{max} = ργ̇_max dP² / (4μ).</caption><thead><tr><th>Ĉ [ppm]</th><th>Φ [%]</th><th>μ [Pa.s]</th><th>λ<sub>N₁</sub> [s]</th><th>λ<sub>e</sub> [s]</th><th>ε⁰</th><th>ε</th><th>W<sub>c</sub></th><th>W<sub>max</sub></th><th>R<sub>c</sub></th><th>R<sub>max</sub></th><th>R<sub>p</sub><sup>max</sup></th></tr></thead><tbody><tr><td rowspan="5">300</td><td>0</td><td>0.0508</td><td>0.1092</td><td>0.482</td><td></td><td>2.1</td><td>93.8</td><td>153.5</td><td>45.3</td><td>74.2</td><td>0.0190</td></tr><tr><td>2</td><td>0.0503</td><td>0.0332</td><td>0.122</td><td></td><td>0.6</td><td>26.5</td><td>46.7</td><td>42.5</td><td>74.9</td><td>0.0192</td></tr><tr><td>6</td><td>0.0484</td><td>0.0456</td><td>0.142</td><td>2.07</td><td>1.1</td><td>29.1</td><td>64.1</td><td>27.5</td><td>60.4</td><td>0.0155</td></tr><tr><td>10</td><td>0.0566</td><td>0.0217</td><td>0.121</td><td></td><td>0.5</td><td>13.0</td><td>30.5</td><td>28.4</td><td>66.6</td><td>0.0170</td></tr><tr><td>14</td><td>0.0782</td><td>0.0713</td><td>0.108</td><td></td><td>2.1</td><td>50.7</td><td>100.1</td><td>24.4</td><td>48.2</td><td>0.0123</td></tr><tr><td rowspan="5">400</td><td>0</td><td>0.0584</td><td>0.1971</td><td>0.795</td><td></td><td>4.3</td><td>197.3</td><td>276.9</td><td>46.0</td><td>64.6</td><td>0.0165</td></tr><tr><td>2</td><td>0.0599</td><td>0.1538</td><td>0.579</td><td></td><td>3.4</td><td>134.0</td><td>216.0</td><td>39.0</td><td>62.9</td><td>0.0161</td></tr><tr><td>6</td><td>0.0743</td><td>0.1647</td><td>0.485</td><td>4.28</td><td>4.6</td><td>145.1</td><td>231.4</td><td>31.8</td><td>50.7</td><td>0.0130</td></tr><tr><td>10</td><td>0.0756</td><td>0.1282</td><td>0.412</td><td></td><td>3.6</td><td>137.0</td><td>180.1</td><td>37.9</td><td>49.9</td><td>0.0128</td></tr><tr><td>14</td><td>0.0694</td><td>0.1169</td><td>0.324</td><td></td><td>3.0</td><td>99.0</td><td>164.2</td><td>32.7</td><td>54.3</td><td>0.0139</td></tr><tr><td rowspan="5">500</td><td>0</td><td>0.0680</td><td>0.3793</td><td>1.020</td><td></td><td>9.6</td><td>398.5</td><td>532.9</td><td>41.4</td><td>55.4</td><td>0.0142</td></tr><tr><td>2</td><td>0.0674</td><td>0.1199</td><td>0.486</td><td></td><td>3.0</td><td>147.1</td><td>168.5</td><td>48.8</td><td>55.9</td><td>0.0143</td></tr><tr><td>6</td><td>0.0768</td><td>0.1973</td><td>0.435</td><td>9.61</td><td>5.7</td><td>144.6</td><td>277.1</td><td>25.59</td><td>49.0</td><td>0.0126</td></tr><tr><td>10</td><td>0.0910</td><td>0.1971</td><td>0.485</td><td></td><td>6.7</td><td>170.1</td><td>277.0</td><td>25.45</td><td>41.4</td><td>0.0106</td></tr><tr><td>14</td><td>0.1006</td><td>0.2309</td><td>0.422</td><td></td><td>8.7</td><td>263.2</td><td>324.4</td><td>30.40</td><td>37.5</td><td>0.0096</td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>Figure 3: Coupling of flow visualisation with space-time diagrams (top row), frequency maps (middle row) where $f_{max} = \frac{2\pi}{\Omega_{max}}$ correspond to maximal inner cylinder frequency, and torque measurements (bottom row), applied to samples with $\tilde{C} = 300$ ppm, $\Phi = 2\%$ (left) and $\Phi = 6\%$ (right). Vertical dashed red lines indicate the critical $R_c$ value for the onset of EIT when this flow state is evidenced (left), or a new regime discussed hereinafter (right, also indicated by a red circle). Linear ridges visible in frequency maps (middle row) correspond to the rotation frequency of the inner cylinder detected by the method.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Figure 4: Detection of critical conditions (here $\dot{\gamma}_c$) for the transition to non-laminar friction regimes, for $\tilde{C} = 300$ ppm and $\Phi = 0, 6 \& 14\%$. Full lines curve are measured torque signals, dotted lines are laminar trends, and dashed lines are the norm of the difference between measured torque and laminar trend. All values are scaled by the maximum torque measured at the end of the run $\tau_{max}$. Black circles and vertical dashed lines represents detected critical $\dot{\gamma}$ where laminar trend and torque curve diverge. The colors scheme for particle concentrations is the same as in figure 2.</figcaption></figure>

visualisation tools were constructed. Firstly, space-time diagrams were built by stacking intensity signal obtained for this vertical at various inner cylinder rotation speeds and thus $R_c$ values. The change in flow structure, due for example to a transition to EIT, (see figure 3) is clearly visible on such graphs: the laminar base flow, called Circular Couette Flow (CCF) is expressed by a uniform grey color while EIT regime's signature is a highly chaotic spatial signal. Secondly, temporal Fast Fourier Transform (FFT) of intensity signals were computed according to a protocol detailed in (Lacassagne et al., 2021; Moazzen et al., 2023), and resulting spectra were stacked vertically as a function of control parameters (frequency map). This allows to discard between stationary flow states with no specific spectral signatures (such as CCF) and flow state with characteristic frequencies or energy spectra, such as EIT, and thus detect the onset of the latter in an even more accurate way (see figure 3). In such figures, the straight darker lines correspond to the inner cylinder rotation frequency, which is captured by the method. Frequency axis are scaled by the maximal inner cylinder frequency $f_{max} = \frac{2\pi}{\Omega_{max}}$.

### 2.3.2 Torque measurement

Simultaneously, the rheometer allowed to measure the torque needed to rotate the inner cylinder (Moazzen et al., 2022, 2023; Dash et al., 2020), providing insight into the friction dynamics, and making it possible to evidence transitions to other friction regime, independently of any visual signature. This is further illustrated in figure 4.
<figure><img src="image_6.png" /><figcaption>Figure 5: Critical value $R_c$ for the onset of EIT as a function of the base fluid elastic number $\varepsilon^0$ (a) confronted to previous literature data (b). The horizontal dashed line in a) represents the critical Reynolds number $R_c^N$ for the onset of the primary instability in a Newtonian fluid in the same Couette geometry (Moazzen et al., 2022). Errors bars are derived from uncertainty propagation from viscosity (dominant) and detection (negligible) uncertainties. c) Critical $W_c$ for the onset of EIT (squares) compared to purely elastic instability threshold computed with the Pakdel-McKinely criterion adapted to the Taylor-Couette geometry (plus signs, see equation 5 proposed by Schäfer et al. (2018)). Empty squares are present work values for $\Phi = 0\%$ and full ones for $\Phi = 2\%$.</figcaption></figure>

In CCF (laminar) regime, the torque is expected to increase linearly with $R$ with a slope proportional to the fluid viscosity (Eckhardt et al., 2007; Dubrulle et al., 2005; Couette, 1890). Any departure of the torque measurements (figure 4 full lines) from this linear trend (figure 4 dotted lines), which becomes even more visible when plotting the difference between the laminar prediction and the torque measurement (figure 4 dashed lines), indicates a change in the friction regime.

For EIT, this corresponds to a sudden jump in measured torque (Moazzen et al., 2023; Martínez-Arias & Peixinho, 2017), followed by a gradual convergence towards a similar laminar-like torque-$R$ slope, yet shifted upwards (Moazzen et al. (2023), figure 3 bottom left subplot, figure 4 left plot). The coupling of all three indicators (space-time diagram, frequency map, torque dynamics) allows to detect the critical $R$ value for the onset of EIT when it arises, labelled $R_c$, as illustrated in figure 3 left column, with a high degree of accuracy (less that 0.5% of variation on $R_c$ obtained with the three methods, on all used fluids).

When no EIT is detected, no jump is reported. Yet, a change in trend is still observed. This takes the form of decrease in the slope of the measured torque-$R$ curve at a given critical $R$ value (circle in figure 3 bottom right plot) or shear rate (circles in figure 4 middle and right plots), corresponding to a continuous increase in the measurement-to-laminar prediction difference (figure 4 dashed lines).

It is remarkable that in pure Newtonian or particle loaded Newtonian fluids, primary and higher order bifurcations are always accompanied by secondary flows and lead to an increase in the torque-$R$ curve slope (Moazzen et al., 2022; Kang & Mirbod, 2021; Dash et al., 2020; Ramesh et al., 2019). Here, the milder slope suggests that additional energy dissipation may occur elsewhere in the flow reducing the friction dissipation at the inner cylinder as compared to the laminar trend, as will be discussed later.

# 3 Results and Discussions

Figure 3 illustrates the two types of results encountered. For particle free fluids the direct transition from CCF to EIT as depicted in (Groisman & Steinberg, 1996; Moazzen et al., 2023) and expected in this setup from our previous study (Moazzen et al., 2023) is indeed evidenced. A similar transition is also found in low particle volume fraction fluids (2%), as illustrated in figure 3 left column ($\tilde{C} = 300$ ppm, $\Phi = 2\%$). However, for all higher particle volume fractions, no such transition is observed in the investigated $R$ range, as illustrated in figure 3 right column ($\tilde{C} = 300$ ppm, $\Phi = 6\%$): the space-time diagram displays no clear sign of chaotic behaviour, the frequency maps no clear change in spectral signature, and the torque measurement no clear jump. The two types of behaviour will be discussed separately in what follows.
## 3.1 Particle free and dilute cases: EIT

For all fluids where EIT has been evidenced (all polymer concentrations at $\Phi = 0\%$ and $\Phi = 2\%$), $R_c$ values for the onset of EIT are represented in figure 5 a) as a function of the base fluid $\mathcal{E}^0$ value. The $\mathcal{E}^0$ values are much higher than the ones reported in previous literature for elasto-inertial TCF of particle loaded viscoelastic fluids (Lacassagne et al., 2021) and also significantly higher than our previous work (Moazzen et al., 2023) (mostly owing to the difference in working fluids: increased solvent viscosity, high polymer molecular weight, high polymer concentrations - see figure 5 b)). A reduction in $R_c$ with increasing $\mathcal{E}^0$ could have been expected from previous studies (Moazzen et al., 2023; Groisman & Steinberg, 1996). It here appears that in the particle-free ($\Phi = 0\%$) case, increasing $\mathcal{E}^0$ (through an increase in polymer concentration) only leads to a weak reduction of $R_c$ (figure 5, full markers). From figure 5 c), it can clearly be seen that the onset of EIT (square symbols) occurs at critical $W_c$ values much higher than what would have been expected from the Pakdel-McKinley criterion (plus symbols) for the onset of purely elastic instability in curves streamline flows (Schäfer et al., 2018; Pakdel & McKinley, 1996). This criterion is expressed as:

$$W_c^{\text{PMcK}} = \frac{M_{\text{crit}}}{\sqrt{2(1-\beta)}} \times \frac{1}{\sqrt{\epsilon}} \qquad (5)$$

where $M_{\text{crit}} = 1.8$, $\beta = \frac{\mu_s}{\mu}$ with $\mu_s$ is the solvent viscosity and $\epsilon = \frac{\delta}{r_i}$ (Pakdel & McKinley, 1996). Thus, the chaotic behaviour is confirmed to be elasto-inertial (and not purely elastic) in nature, requiring both non negligible inertia and elasticity to develop, even at such high $\mathcal{E}^0$ values.

When considering dilute ($\Phi = 0 \& 2\%$ here) particle loaded fluids for which EIT is observed in the experimental range, it appears that this monotonic decrease in $R_c$ with $\mathcal{E}$ is no longer observed, confirming that particle do mitigate the transition to EIT (figure 5, empty markers), as expected from (Lacassagne et al., 2021). The $R_c$ values are essentially identical for the $\Phi = 0\%$ and $\Phi = 2\%$ cases for all three polymer concentrations, with possibly a slight decrease of $R_c$ upon particle addition for the two lower $\mathcal{E}$ values ($\tilde{C} = 300$ and $\tilde{C} = 400$ ppm fluids), but a slightly delayed transition at $\Phi = 2\%$ in the higher $\mathcal{E}$ case. The 2% particle addition yet causes a significant decrease of the critical $W$ values, as illustrated in figure 6 c), (still much higher than the purely elastic instability threshold) which matches the trend observed when revisiting data from (Lacassagne et al., 2021). It thus implies that in the dilute $\Phi$ limit, the effect of particle addition is mostly to promote elasto-inertial instabilities by reducing the amount of elasticity ($W$) needed at a given inertia ($R$) for EIT to arise.

## 3.2 Absence of EIT in the semi-dilute regime

A striking result though, is that for particle volume fractions higher than 2%, in what can be referred to as a semi-dilute particle concentration regime (Moazzen et al., 2022), no transition to EIT is observed in the investigated $R$ range (see 3 right column). The rheometer limit is reached when $\dot{\gamma} = 1405 \text{ s}^{-1}$ for each samples, which, since viscosity increase with particles concentration, leads to a decrease in the maximum attainable value $R_{max}$ and $W_{max}$. EIT not being observed here thus implies that it has been either delayed beyond $R_{max}$ or $W_{max}$, or completely suppressed, in such semi-dilute suspensions. The question naturally arising is then: Could the effect of particle addition simply be to globally reduce the fluid's elasticity, resulting in apparent $\mathcal{E}$ values low enough for elasto-inertial instabilities no longer to be expected? At this point, it is thus worth focusing on the effect of particle addition on the apparent viscoelastic time-scales of particle loaded fluids. Those are measured by two protocols, as detailed in section 2.3: first normal stress differences and elongational experiments. Results are reported in figure 2 b) (scaled by the $\Phi = 0\%$ value, $\lambda^0$) and in table 1. As a reminder: the apparent elastic numbers $\mathcal{E}$ are computed for each fluid from $\lambda_{N1}$ values and also reported in table 1. It appears that 2% particle addition leads to a steep decrease of $\lambda$ and thus of $\mathcal{E}$. Although the elastic time scale keeps a weakly decreasing trend with particle addition beyond 2% (confirmed by both measurement techniques) the simultaneous increase in fluid viscosity makes it so that $\mathcal{E}$ remains in the same range on the whole set of $\Phi$ investigated. It means that all fluids are still in the high elasticity regime ($\mathcal{E} \sim O(10^0)$), where direct transition from CCF to EIT would be expected (Groisman & Steinberg, 1996; Moazzen et al., 2023). Thus, the suppression or delay of EIT upon particle addition beyond 6% can not solely be explained by an apparent reduction of the fluid overall elasticity. Moreover, should that be the case, the difference in behaviour would be expected between the 0% and 2% cases for which the decrease in $\mathcal{E}$ is the more evident, and not between the 2% and 6% cases for which $\mathcal{E}$ stay equivalent or even increase. The change in behaviour between 2% and 6% has to be understood in light of local flow dynamics and fluid particle interactions, and in that context, torque measurements provide insightful information.

Figure 6 a) displays the evolution of the torque measured as a function of $W = \lambda_{N1}\dot{\gamma}$ for the example of the $\tilde{C} = 300$ ppm ($\mathcal{E}^0 = 2.07$) fluid in one dilute case where EIT is observed (2%, $\mathcal{E} = 0.62$) and two semi-dilute cases for which it is not reported (6% - $\mathcal{E} = 1.06$ and 10% - $\mathcal{E} = 0.46$). As a confirmation of the above statement on the elasticity levels remaining high, it can be seen that the $W$ range covered is similar in all three cases. The transition
<figure><img src="image_7.png" /><figcaption>Figure 6: a) Torque measured on the inner cylinder as a function of $W = \dot{\gamma}\lambda_{N_1}$. Vertical dashed lines indicate the $W_c$ values for trend change. Dashed black linear trends correspond to the laminar flow state. b1) and b2) are $R_c$ as a function of $\Phi$, and $W_c$ as a function of $\mathcal{E}$, respectively, with squares reporting the onset of EIT and crosses the onset of EID, at approximately constant $R_c$ values for each regime, indicated by dashed lines. c) Values of $W_c$ (same symbols) scaled by the zero-particle value $W_c(\Phi = 0\%)$, as a function of the particle-related viscosity increase ($\mu_{\Phi}/\mu$). The grey area correspond to EIT data points and the vertical dashed black line a $\mu_{\Phi}/\mu = 0.1$ limit from dilute to semi-dilute concentration regime.</figcaption></figure>

to EIT in the dilute 2% case is evidenced by a clear jump in torque values (Moazzen et al., 2023; Martínez-Arias & Peixinho, 2017). No such behaviour is observed for the two other cases. However, as mentioned in section 2.3, close inspection of the torque curves reveals a change in slope at $W$ values (reported in table 1) close to the critical one for the onset of EIT in the 2% case, where the torque-$R$ curve is found to depart from its laminar, CCF trend (see figure 3 bottom right plot, and 4 middle and right subplots). Thus, although there is neither visual nor dynamic (torque) evidence of EIT, the torque measurement suggests that there is indeed a modification of the dynamics, in terms of friction. The latter observation, supports the assumption that EIT is not only delayed, but rather replaced by another mechanism qualified here as Elasto-Inertial Dissipative (EID), which, to the best of the authors knowledge, has never been reported before. EID is not associated to observable flow structures and thus contributes to an apparent remanuarization of the flow by suppression of the large scale chaotic patterns, but still corresponds to a different friction dynamics than the initial laminar CCF. Indeed, the fact that the torque increase with increasing inertia is milder than the laminar prediction supports the hypothesis that additional energy dissipation is occurring in the flow, reducing friction at the inner cylinder, hence the term "dissipative regime". It is expected that particle-particle and particle-polymer interactions are at the root of such mechanism, considering that EIT can only be observed in the semi-dilute range ($\Phi > 6\%$). This assumption will be discussed further. It should here be recalled that EIT is instead characterised by a torque increase (jump) and then convergence towards a linear trend with similar slope as that found in CCF (see figure 3 bottom left) corresponding to a saturation in the EIT dynamics as reported in Moazzen et al. (2023), figure 6). This is thus a first fundamental difference between EIT and EID.

## 3.3 Onset of EID and EIT

Further information can be gained on this new transition by plotting the evolution of all $R_c$ and $W_c$ values (for the onset of EIT with squares or EID with crosses) as a function of $\Phi$ and $\mathcal{E}$ on figure 6 b1) and b2) respectively. This shows that $R_c$ is both quasi $\Phi$-independent and $\mathcal{E}$-independent, although slightly lower for the onset of EID than for EIT, suggesting that EID mechanism is actually occurring sooner (in terms of $R$ during the ramp-up protocol) than EIT in semi-dilute particle loaded fluids, replacing it under comparable elasto-inertial conditions. Figure 6 c) compiles the previous results by plotting all $W_c$ (EIT or EID) normalized by the $\Phi = 0\%$ pure viscoelastic fluid value, as a function of the relative particle contribution to the overall viscosity $\mu_{\Phi}/\mu = (\mu - \mu(\Phi = 0))/\mu$. The graph splits up into two distinct regions. On the left hand side, the viscosity increase upon particle addition remains negligible ($\mu_{\Phi}/\mu < 0.1$), and EIT is promoted, as discussed earlier, with $W_c < W_c(\Phi = 0)$. On the right hand side, particles significantly increase the overall viscosity ($\mu_{\Phi}/\mu > 0.1$) and, in the present work, EIT is replaced by EID. In this second range, it appears that the critical $W$ values for the onset of EID increase with increasing $\Phi$ (and thus increasing $\mu_{\Phi}/\mu$). The effect of particle addition in the semi-dilute regime can thus be understood as
follows: firstly, the significant increase in overall viscosity µ leads to a stronger viscous dissipation mechanism, as illustrated by the steeper torque slopes in the laminar regime in figure 6 a) and figure in 4. Further increasing the inner cylinder rotational speed leads to an increase in W and R until the required elasto-inertial conditions for the onset of EID are reached, and in particular the constant critical R number. However, this requires higher ˙γ and thus higher W in the most particle loaded fluids, where the apparent increased viscosity promotes viscous dissipation and delays the onset of EID. Interestingly in Lacassagne et al. (2021), the same trend in the same µ$_{Φ}$/µ > 0.1 range was observed for the delay of the EIT onset in non-dilute particle suspensions (see figure 6 c), diamonds).

Connection can be made with the recent work of Zhang & Shaqfeh (2023), who used numerical simulations to show that under conditions of β = µs µ 0.5 and W ≫ 1 similar to that of the present work, particle-particle interactions significantly affect the viscoelastic stress hardening mechanisms known to arise in such suspensions by altering the particle-induced fluid stresses, leading to an increased stresslet. Interestingly, this effect was evidenced for Φ values of 2.5% or 5% quite close to the critical concentration value reported here for the EIT-EID shift, suggesting a possible connection between this semi-dilute mechanism and the observation of EID.

3.4 Pseudo-Nusselt number in EID

Additional insight into the EID regime developing after the transition can here only be brought by the torque information, as no flow perturbation is observed from the outer-cylinder perspective of flow visualisation. It is yet worth mentioning that the direct visualisation method employed here only gives access to flow properties visible from the flow periphery and the outer cylinder, while torque measurements probe the flow behaviour closer to the inner cylinder. This suggests that the origin of EID could be spatially localized either in the bulk of the flow, or close to the inner cylinder with yet no associated large scale secondary flow structures, since they cause no apparent outer flow disturbance and no increase in apparent torque on the inner cylinder as would many well known inertial, elastic, or elasto-inertial Taylor-Couette instabilities (Moazzen et al., 2022, 2023).

A relevant parameter to infer radial energy transfer properties from torque measurement is the the pseudo- Nusselt number N as defined in Eckhardt et al. (2000) and used for example in Moazzen et al. (2022, 2023):

$$
\mathcal{N} = \frac{J}{J_{lam}} = \frac{\frac{\tau}{2\pi\rho h}}{2\nu r_o^2 r_i^2 \frac{\Omega}{r_o^2 - r_i^2}} \quad (6)
$$

This quantity represents the flow capacity to transfer momentum in the radial direction, J being the momentum transport term and J$_{lam}$ its value in the laminar CCF state.

For all experiments in the EID regime, the critical pseudo-Nusselt value N$_{c}$ = N ( ˙ γ$_{c}$ ) is defined at the onset of EID, and the quantity Δ N = N − N$_{c}$ is plotted as a function of Δ R = R − R$_{c}$ . The example of the ˜ C = 300 ppm base fluid is illustrated in figure 7 a. This represents the variations of radial momentum transfer from the laminar reference within the EID regime, upon increasing inertia.

In most particle-loaded Newtonian fluids, this pseudo-Nusselt number is increased compared to the laminar case by secondary flows at the primary bifurcation, showing an improved radial momentum transfer ability (Moazzen et al., 2022; Kang & Mirbod, 2021; Dash et al., 2020; Ramesh et al., 2019). In EIT, one would also report, Δ N > 0 and a Δ N increase with increasing Δ R (not shown here for the sake of plots clarity) consistently with the observations of Moazzen et al. (2023). However from Φ = 6%, this quantity decreases with increasing Δ R . showing an opposite N behavior between EIT and EID, and supporting the fact that momentum radial transfer is reduced with EID. Figure 7 a also suggests that increasing Φ leads to an even steeper decrease of Δ N .

To investigate this point, the slope ∂ Δ N ∂ Δ R (Δ R ) is locally computed for all fluids (those represented in figure 7 a and all the 400 and 500 ppm cases) and average values over the full Δ R range are represented in figure 7 b) (with the error bars obtained from the standard deviation on the same range), as a function of the particle- to-polymer concentration ratio Φ / ˜ C . A good collapse on a clear monotonic trend is observed positioning the particle-to-polymer concentration ratio as a key parameter of the EID behaviour. This is yet another indication of the possible connection to the mechanisms reported by Zhang & Shaqfeh (2023). The more particles and subsequent particle-particle interactions for a given number of polymer chains, the larger the stresslet increase and the stress hardening effects on a limited number of polymer chains, and thus the stronger the resulting local energy dissipation, reducing the radial momentum transfer towards the outer cylinder.

4 Conclusion

This study aimed at describing the influence of particle addition on the transition to EIT in high elasticity base fluids. Spherical particles were added in concentration regimes ranging from dilute to semi-dilute, in base fluids with high E values for which direct transition from CCF to EIT upon R increase would be expected. Transitions
<figure><img src="image_8.png" /><figcaption>Figure 7: a) $\Delta\mathcal{N} = \mathcal{N} - \mathcal{N}_c$ as a function of $\Delta\mathcal{R} = \mathcal{R} - \mathcal{R}_c$ in the EID regime, for a 300 ppm fluids with $\Phi = 6, 10, 14\%$. b) Mean of the slope calculated as $\frac{\partial\Delta\mathcal{N}}{\partial\Delta\mathcal{R}}$ for all experiments where EID was reported. Markers and colors are similar to the ones used in figures 2, 5, 6</figcaption></figure>

were monitored using a coupling of flow visualisation and inner-cylinder torque measurements. Transition to EIT in particle-free Boger fluids was indeed retrieved, at lower $\mathcal{R}_c$ and milder $\mathcal{E}^0$ dependency than in previous studies (Moazzen et al., 2023; Groisman & Steinberg, 1996), but still with a clear elasto-inertial nature. The addition particles at a low volume fraction ($\Phi = 2\%$) led to an earlier onset of EIT in terms of $W_c$ (but not obvious in terms of $\mathcal{R}_c$). However, the most striking result was that particle addition at higher volume fractions ($\Phi \ge 6\%$, semi-dilute regime) led to EIT being no longer observed in the $\mathcal{R}$ and $W$ range investigated, suggesting an apparent relaminarization of the flow with either a delay or a replacement of EIT, which could not be ascribed to a reduction in apparent fluid elasticity. Instead, a new dissipative mechanism, termed here Elasto-Inertial Dissipative (EID), was identified from the torque curves and friction properties. The critical conditions for the onset of EID also seem elasto-inertial and quite similar to that for the onset of EIT in pure fluids or dilute suspensions. It appears that the relative increase of fluid viscosity upon particle addition is a key parameter in the onset of EIT or EID, as confirmed by compiling results from the present work and Lacassagne et al. (2021): the dilute regime at $\Phi \le 5\%$ corresponds to a minor increase in viscous dissipation ($\mu_{\Phi}/\mu \le 0.1$) and earlier onset of EIT, while increasing $\Phi$ in the semi-dilute regime sees significant increase in viscous dissipation upon particle addition ($\mu_{\Phi}/\mu > 0.1$) and a delay in EIT or delayed onset of EID. The evolution of the pseudo-Nusselt number in the EID regime shows that particle addition reduces radial momentum transfers, commensurate with the $\Phi/\tilde{C}$ ratio. In summary, EID is thus expected to be 1) a local mechanism, 2) elasto-inertial in that it requires non negligible values of both $\mathcal{R}$ and $W$, and 3) driven by polymer-particle interactions as it arises only in the semi-dilute regime. The accurate determination of the critical conditions for the onset of EID ($\Phi, W, \mathcal{R}$ and geometrical parameters) could be the object of a future study. More generally, future studies should aim at providing local quantitative information on flow velocity fields (Boulafentis et al., 2024) and possible particle migration phenomena. This would allow to investigate local particle-polymer interaction mechanisms and their effect on polymeric stress field in a semi-dilute context (Zhang & Shaqfeh, 2023), which would be of great interest to unravel the mechanisms behind delayed EIT or EID in the Taylor-Couette flow of particle-loaded viscoelastic fluids.

## Acknowledgements

This work was partially funded by the French National Research Agency (ANR-22-CE05-0017-01). Financial support from Region Hauts de France is also gratefully acknowledged.

## Declaration of interests

The authors declare no conflict of interest.
References

ANDERECK, C. D., LIU, S. S. & SWINNEY, H. L. 1986 Flow regimes in a circular Couette system with independently rotating cylinders. *J. Fluid Mech.* **164**, 155–183.

BAROUDI, L., MAJI, M. V., PELUSO, S. & MORRIS, J. F. 2023 Taylor-Couette flow of hard-sphere suspensions: overview of current understanding. *Philos. Trans. R. Soc. A* **381**.

BOGER, D. V. 1977 Demonstration of upper and lower Newtonian fluid behaviour in a pseudoplastic fluid. *Nature* **265**.

BOULAFENTIS, T., LACASSAGNE, T., CAGNEY, N. & BALABANI, S. 2023 Experimental insights into elasto-inertial transitions in Taylor–Couette flows. *Philos. Trans. R. Soc. A* **381** (2243), 20220131.

BOULAFENTIS, T., LACASSAGNE, T., CAGNEY, N. & BALABANI, S. 2024 Coherent structures of elasto-inertial instabilities in Taylor–Couette flows. *J. Fluid Mech.* **986**, A27.

CAMPO-DEAÑO, L. & CLASEN, C. 2010 The slow retraction method (SRM) for the determination of ultra-short relaxation times in capillary breakup extensional rheometry experiments. *J. Non-Newton. Fluid Mech.* **165** (23), 1688–1699.

COUETTE, M. 1890 Distinction de deux régimes dans le mouvement des fluides. *J. Phys. Theor. Appl.* **9**, 414–424.

DASH, A., ANANTHARAMAN, A. & POELMA, C. 2020 Particle-laden Taylor–Couette flows: higher-order transitions and evidence for azimuthally localized wavy vortices. *J. Fluid Mech.* **903**, A20.

DUBRULLE, B., DAUCHOT, O., DAVIAUD, F., LONGARETTI, P.-Y., RICHARD, D. & ZAHN, J.-P. 2005 Stability and turbulent transport in Taylor–Couette flow from analysis of experimental data. *Physics of Fluids* **17**, 095103.

DUTCHER, C. S. & MULLER, S. J. 2013 Effects of moderate elasticity on the stability of co- and counter-rotating Taylor–Couette flows. *J. Rheol.* **57** (3), 791–812.

ECKHARDT, B., GROSSMANN, S. & LOHSE, D. 2000 Scaling of global momentum transport in Taylor–Couette and pipe flow. *Eur. Phys. J. B* **18**, 541–544.

ECKHARDT, B., GROSSMANN, S. & LOHSE, D. 2007 Torque scaling in turbulent Taylor–Couette flow between independently rotating cylinders. *J. Fluid Mech.* **581**, 221–250.

FARDIN, M. A., PERGE, C. & TABERLET, N. 2014 “The hydrogen atom of fluid dynamics” – introduction to the Taylor–Couette flow for soft matter scientists. *Soft Matter* **10** (20), 3523.

GROISMAN, A. & STEINBERG, V. 1996 Couette–Taylor flow in a dilute polymer solution. *Phys. Rev. Lett.* **77**.

GROISMAN, A. & STEINBERG, V. 1998 Mechanism of elastic instability in Couette flow of polymer solutions: Experiment. *Phys. Fluids* **10** (10), 2451–2463.

GROISMAN, A. & STEINBERG, V. 2000 Elastic turbulence in a polymer solution flow. *Nature* **405**.

GROISMAN, A. & STEINBERG, V., 2004 Elastic turbulence in curvilinear flows of polymer solutions. *New J. Phys.* **6**.

JAMES, D. F. 2009 Boger Fluids. *Annu. Rev. Fluid Mech.* **41** (1), 129–142.

KANG, C. & MIRBOD, P. 2021 Flow instability and transitions in Taylor–Couette flow of a semidilute non-colloidal suspension. *J. Fluid Mech.* **916**, A12.

KANG, C., SCHATZ, M. F & MIRBOD, P. 2024 Hysteresis and ribbons in Taylor-Couette flow of a semidilute non-colloidal suspension. *Phys. Rev. Fluids* **9** (2), 023901.

LACASSAGNE, T., BOULAFENTIS, T., CAGNEY, N. & BALABANI, S. 2021 Modulation of elasto-inertial transitions in Taylor–Couette flow by small particles. *J. Fluid Mech.* **929**.

LARSON, R. G., SHAQFEH, E. S. G. & MULLER, S. J. 1990 A purely elastic instability in Taylor–Couette flow. *J. Fluid Mech.* **218**, 573.

LIN, F., SONG, J., LIU, N., WAN, Z., LU, X.-Y. & KHOMAMI, B. 2024 Maximum drag enhancement asymptote in turbulent Taylor–Couette flow of dilute polymeric solutions. *Journal of Non-Newtonian Fluid Mechanics* **323**, 105172.
LIN, F., SONG, J., ZHAO, Z., LIU, N., LU, X.-Y. & KHOMAMI, B. 2023 A novel transition route to elastically dominated turbulence in viscoelastic Taylor–Couette flow. *J. Non-Newton. Fluid Mech.* **312**, 104968.

LIN, Y., PHAN-THIEN, N. & CHEONG KHOO, B. 2014 Normal stress differences behavior of polymeric particle suspension in shear flow. *J. Rheol.* **58** (4).

MAJJI, M. V., BANERJEE, S. & MORRIS, J. F. 2018 Inertial flow transitions of a suspension in Taylor–Couette geometry. *J. Fluid Mech.* **835**, 936–969.

MARTÍNEZ-ARIAS, B. & PEIXINHO, J. 2017 Torque in Taylor–Couette flow of viscoelastic polymer solutions. *J. Non-Newton. Fluid Mech.* **247**, 221–228.

MOAZZEN, M., LACASSAGNE, T., THOMY, V. & BAHRANI, S. A. 2022 Torque scaling at primary and secondary bifurcations in a Taylor–Couette flow of suspensions. *J. Fluid Mech.* **937**.

MOAZZEN, M., LACASSAGNE, T., THOMY, V. & BAHRANI, S. A. 2023 Friction dynamics of elasto-inertial turbulence in Taylor–Couette flow of viscoelastic fluids. *Philos. Trans. R. Soc. A* **381** (2246), 20220300.

PAKDEL, P. & MCKINLEY, G. H. 1996 Elastic instability and curved streamlines. *Phys. Rev. Lett.* **77** (12), 2459–2462.

RAMESH, P., BHARADWAJ, S. & ALAM, M. 2019 Suspension Taylor–Couette flow: co-existence of stationary and travelling waves, and the characteristics of Taylor vortices and spirals. *J. Fluid Mech.* **870**, 901–940.

SCHÄFER, C., MOROZOV, A. & WAGNER, C. 2018 Geometric scaling of elastic instabilities in the Taylor–Couette geometry: A theoretical, experimental and numerical study. *J. Non-Newton. Fluid Mech.* **259**, 78–90.

SONG, J., WAN, Z.-H., LIU, N., LU, X.-Y. & KHOMAMI, B. 2021 A reverse transition route from inertial to elasticity-dominated turbulence in viscoelastic Taylor–Couette flow. *J. Fluid Mech.* **927**, A10.

TANNER, R. I. 2020 Computation and experiment in non- colloidal suspension rheology. *J. Non-Newton. Fluid Mech.* **281**, 104282.

TAYLOR 1922 Stability of a viscous liquid contained between two rotating cylinders. *The Royal Society* pp. 289 – 343.

ZHANG, A. & SHAQFEH, E. S. G. 2023 Rheology of non-Brownian particle suspensions in viscoelastic solutions. Part 1: Effect of the polymer concentration. *J. Rheol.* **67** (2), 499–516.