Using Acoustic Perturbations to Dynamically Tune Shear Thickening in Colloidal Suspensions

TL;DR

This paper demonstrates how acoustic perturbations via piezoelectric transducers can actively and dynamically control shear thickening in colloidal suspensions by disrupting force chains, offering a versatile method for flow regulation.

Contribution

It introduces a novel active control method using acoustic perturbations to tune shear thickening in colloidal suspensions, adaptable to various flow geometries.

Findings

Acoustic perturbations can dethicken colloidal suspensions by disrupting force chains.

The technique is easily adaptable to different flow geometries.

Active control of suspension viscosity is achievable through acoustic methods.

Abstract

Colloidal suspensions in industrial processes often exhibit shear thickening that is difficult to control actively. Here, we use piezoelectric transducers to apply acoustic perturbations to dynamically tune the suspension viscosity in the shear-thickening regime. We attribute the mechanism of dethickening to the disruption of shear-induced force chains via perturbations that are large relative to the particle roughness scale. The ease with which this technique can be adapted to various flow geometries makes it a powerful tool for actively controlling suspension flow properties and investigating system dynamics.

Full text

††preprint: APS/123-QED

The orders-of-magnitude increase in viscosity that arises under high shear makes dense suspensions ideal for numerous industrial applications including shock absorption, damping, soft-body armor, astronaut suits, and curved-surface polishing Brown and Jaeger (2014); Wagner and Brady (2009); Fischer et al. (2006); Lee et al. (2003); Li et al. (2015); Cwalina et al. (2016). The challenge in using such shear thickening fluids, however, is that this same increase in viscosity can lead to jamming and failure of pumping and mixing equipment driving the flows. The ability to manage these limitations of this important technological material remains challenging Brown and Jaeger (2014); Barnes (1989) because it requires actively tuning the suspension viscosity. Shear thickening viscosity previously has been tuned passively by changing the physical properties of the suspension constituents, such as the volume fraction ($\phi$) Cwalina and Wagner (2016); Mari et al. (2014); Denn et al. (2018), particle size Guy et al. (2015), particle shape Cwalina et al. (2017); Hsiao et al. (2017), roughness Hsu et al. (2018); Lootens et al. (2005), surface chemistry Yang et al. (2017), and solvent attributes Wagner and Brady (2009); Brown et al. (2010); James et al. (2018); Maranzano and Wagner (2001); all of which affect the formation of the force chains responsible for thickening Cates et al. (1998); Yang et al. (2017); Majmudar and Behringer (2005); Mari et al. (2014); Seto et al. (2013); Morris (2018); Bi et al. (2011); Comtet et al. (2017); Clavaud et al. (2017); Lin et al. (2015); Royer et al. (2016); Fernandez et al. (2013); Seto et al. (2013); Wyart and Cates (2014); Brown and Jaeger (2012); Guy et al. (2015); Thomas et al. (2018); Rathee et al. (2017). However, active tuning to change the flow properties on demand without changing the physical properties of the suspension constituents or without modifying the suspension has until recently remained largely unexplored.

Recently, it was shown that macroscopic boundary oscillations can be used to actively tune shear thickening in a dense suspension Lin et al. (2016). The dethickening mechanism entails disruption of the force chains through application of an oscillatory shear flow orthogonal to the primary flow direction. Further simulations explored the parameter space for active tuning and showed that this mechanism is robust and that such orthogonal mechanical perturbations can tune the suspension viscosity over a wide range of shear rates and volume fractions Ness et al. (2018). Unfortunately, using macroscopic boundary oscillations to introduce orthogonal perturbations is not practical for many applications.

Here, we determine whether externally applied acoustic perturbations can be used to actively tune the suspension viscosity in the shear thickening regime. The advantage of this approach is that acoustic perturbations can controllably manipulate particles Collins et al. (2018); Sehgal and Kirby (2017); Naseer et al. (2017); Lim et al. (2019); Antfolk et al. (2014); Destgeer et al. (2019); Muller et al. (2013); Ding et al. (2012) and can be applied via readily available piezoelectric transducers that are bonded to otherwise fixed surfaces Antfolk et al. (2014); Laurell et al. (2007); Lenshof et al. (2012). The key principle motivating our work is that nanoscale acoustic disturbances will locally perturb particles and break the force chains responsible for thickening (Fig. 1). From previous studies James et al. (2018); Hsiao et al. (2017); Denn and Morris (2014); Lin et al. (2016); Hsu et al. (2018), we hypothesize that, to break force chains and achieve dethickening, such perturbations should generate relative particle displacements greater than the particle surface roughness scale and on time scales faster than the turnover rate for the force chains. In our system, the acoustic perturbations can generate non-uniform particle displacements Bruus (2012) over the length of the force chains, and because these force chains are fragile Cates et al. (1998); Majmudar and Behringer (2005), there can be a weak link in the chain where two particles are displaced beyond the surface roughness, resulting in breaking of the chain. In addition, in dense suspensions where interparticle distance is small, the particle contacts may further be disrupted by the acoustic scattering effects from nearby particles Silva and Bruus (2014) and bulk flows due to acoustic dissipation Wiklund et al. (2012). Furthermore, a typical ultrasonic transducer can apply these acoustic perturbations at time scales that is orders of magnitude shorter than the force chain turnover time, which is typically $\sim 1/\dot{\gamma}$ and O($\sim$1s) in our system, where $\dot{\gamma}$ is the strain rate. Thus, we anticipate that acoustic perturbations can disrupt the force chains in a thickened suspension and actively tune the viscosity.

We test this idea on a dense silica colloidal suspension undergoing controlled shear and simultaneous acoustic excitation. Our suspension consists of charge-stabilized silica particles, 2 $\mu$m in diameter, in dipropylene glycol at volume fractions of $\phi=0.53$ and $\phi=0.50$. Using atomic force microscopy, we measured the average particle surface roughness to be $\sim$2 nm, which is well below the estimated dipolar oscillation amplitude (Fig. 2). Our testing apparatus consists of a piezoelectric disk (APC International, Material 841) of diameter 21 mm and thickness 1.80 mm bonded via epoxy to an aluminum (6061-T6) bottom plate of diameter 19 mm and thickness 8.57 mm (Fig. 3a). The acoustic perturbations are generated by exciting the piezoelectric crystal in the thickness mode at a resonance frequency $f_{r}=$1.15 MHz. This mode applies perturbations in the gradient direction of the primary shear flow. The bottom plate thickness is optimized for maximum energy transfer to the suspension. The piezo-plate setup is integrated with an Anton-Paar MCR702 Rheometer, and a glass top plate is used to apply the primary shear flow and measure the shear viscosity. The suspension is confined between the two plates and the gap is set to 0.64 mm (See SI Section I for calibration). Using this setup, we quantify the effects that acoustic perturbations have on the shear-thickening behavior of our suspensions.

In our measurement protocol, we apply a steady shear to thicken the suspension. After the suspension reaches steady state, we add an amplitude-modulated (AM) acoustic perturbation. The voltage signal used to drive the piezoelectric element has the form $V=V_{0}[1+\sin{(2\pi f_{m}t+\Phi_{0})]}\sin{(2\pi f_{r}t)}$, where $f_{r}$ is the resonance frequency, and $f_{m}=0.2\text{ Hz}$ is the modulation frequency (Fig. S2). The phase $\Phi_{0}$ is set arbitrarily and the voltage $V_{0}$ is set at 2.5 V to obtain maximum peak-to-peak voltage ($V_{pp}$) of 10 V. This approach quantifies the dynamic, phase-sensitive, and power-dependent viscosity response of the suspension in a single measurement. The dynamic response probes the important time scales that govern the formation and breakup of force chain. The phase-sensitive response obtained from the controlled modulation eliminates noise in the temporal measurements of the viscosity. Finally, the power-dependent response quantifies the efficiency of this dethickening method. We perform these measurements for a range of strain rates ($\dot{\gamma}$) over which the fluid behavior varies from a Newtonian to a fully thickened state.

As hypothesized, the viscosity response of the thickened suspension depends sensitively on the acoustic perturbations (Fig. 3b and Fig. S4). For $\phi=0.53$ suspension sheared at strain rates corresponding to the thickened regime ($\dot{\gamma}=0.703\text{ s}^{-1}$), the instantaneous relative viscosity, $\eta_{r}$, oscillates when the acoustic perturbations are applied. The oscillations result from the amplitude modulation of the acoustic perturbations with the greatest decrease in viscosity arising from the largest perturbation amplitude. The minima of these viscosity oscillations are still above the Newtonian viscosity, which suggests that a higher acoustic power is required to break up all the force chains. In contrast, for strain rates corresponding to the unthickened Newtonian regime ($\dot{\gamma}=0.111\text{ s}^{-1}$) we observe no modulations in the viscosity. This difference in response is consistent with the proposed mechanism that the acoustic perturbations break the shear induced force chains responsible for thickening leaving the other suspension properties largely unchanged.

For strain rates corresponding to the transition regime between the Newtonian and fully thickened state ($\dot{\gamma}=0.596,0.352\text{ s}^{-1}$), the acoustic perturbations are sufficient to dethicken the suspension viscosity to the value in the Newtonian regime. Interestingly, the maximum viscosity during the time when AM perturbations are applied does not recover fully to the steady state value. We interpret this response to indicate that the AM frequency is too rapid for the force chains to fully form between successive oscillations at these strain rates. This picture is supported by the fact that the viscosity recovery time when the perturbations are turned off is much longer than the AM oscillation period.

We extract the magnitude of acoustic dethickening as a function of strain rates using a phase-sensitive analysis of the instantaneous viscosity response curves (Fig. 4, see SI section IV for details). We observe that the application of the acoustic perturbations decreases the viscosity substantially in the regime where the suspension thickens. This response is sensitive to the strain rate, with the largest decrease occurring in the transition regime (Fig. 4a,b). We quantify this response by plotting the $\%$Reduction in viscosity versus $\tilde{\dot{\gamma}}$, the strain rate normalized by the strain rate at the onset of thickening (Fig. 4c). We find negligible decrease in the viscosity in the Newtonian regime ($\tilde{\dot{\gamma}}<1$), in which the force chains are mostly absent. We find the highest reduction in the transition region ($1<\tilde{\dot{\gamma}}<2$), in which the applied acoustic perturbations are sufficient to break up the majority of the force chains. This decrease in viscosity to nearly the Newtonian value effectively shifts the onset strain rate for thickening. Finally, we find that the $\%$Reduction decreases and plateaus in the fully thickened regime ($2<\tilde{\dot{\gamma}}$). This plateau is consistent with literature predictions that the force chain network saturates in the thickened regime Mari et al. (2014) and with the idea that the acoustic perturbations are only breaking up a fraction of this network at this power. The trends at each volume fraction are similar, but dethickening is lower in $\phi=0.50$ suspensions than in $\phi=0.53$ suspensions. This difference in effect magnitude is consistent with the results from boundary oscillation simulations Ness et al. (2018). Collectively, these data suggest that a further increase in power would shift the onset of the thickening to higher strain rates or stress scales Wyart and Cates (2014), and increase the dethickening in the thickened regime.

Using the AM protocol, we determine the power-dependent response of the suspension and observe a logarithmic decrease in the viscosity above a threshold power (Fig. 5). We define a normalized viscosity $\bar{\eta}=(\eta-\eta_{min})/(\eta_{max}-\eta_{min})$, where $\eta_{max}$ and $\eta_{min}$ are the maximum and minimum viscosities at given $\dot{\gamma}$ during a modulation period. We plot $\bar{\eta}$ versus power for strain rates in which the turnover time of force chains is substantially faster than the time period of the AM signal (e.g. blue curve in figure 3b). Strikingly, we find a threshold-like behaviour at very low powers (inset Fig. 5), where the viscosity decrease is very slow and negligible up to powers $\sim$5 mW (corresponding to $V_{pp}=1.4$ V). This trend is consistent with our mechanistic hypothesis that a minimum particle displacement is required to break force chains. Beyond this threshold, the viscosity initially decreases rapidly at lower powers up to $\sim$ 0.08 W (orange curve in Fig. 5), which suggests that a significant number of force chains are broken even at a small input power. At higher power, the viscosity decreases more slowly. This non-linear evolution suggests that although a higher power is required to completely eliminate the thickening effects of force chains, a large fraction of dethickening can be achieved by just small input power.

As the power is ramped down (green curve), we observe a hysteretic response indicating that while a larger power is needed to break up the force chains, a significantly lower power is needed to maintain this disruption during this period. These results suggest that the suspension retains some memory of the microstructure within each AM cycle even though the modulation time period is large enough for the force chains to fully recover. Our findings are consistent with the simulations of orthogonal boundary oscillations that indicate a dethickened viscosity can be maintained using pulsed perturbations Ness et al. (2018). Our ability to precisely measure these hysteresis loops shows that this amplitude modulation technique can determine the relevant time scales to form and break force chains. This understanding can then be leveraged to optimize strategies for achieving dethickening with minimal amount of power.

For the current implementation of these acoustic modulations, the power required to dethicken the suspension is over an order of magnitude higher than the power required for shearing the thickened suspension. Further studies, aimed at improving the coupling of acoustic energy to the suspension and applying perturbation intermittently, may change this balance. Moreover, the particle displacements that are driving these changes in the viscosity maybe complex in the dense suspensions because of the acoustic scattering from nearby particles, bulk acoustic flows, and steric constraints. Measurements or simulations of the exact particle displacements are required to confirm their hypothesized mechanistic link with the particle roughness. This understanding of the suspension microstructure upon acoustic excitation and development of strategies for enhancing disruptions should also enable more efficient control of the suspension viscosity.

Even in its current form, however, our method has paramount advantages in the applications where the goal is to increase the flow rate, unclog a system, or control the viscosity, without energy concerns. Such applications include high-throughput processing of dense suspensions, avoiding jamming in narrow conduits, 3D printing, and designing of smart materials. In each of these cases, the perturbations can be applied by simply bonding a piezoelectric element to a fixed surface, which makes this method easy to integrate with the existing practical systems without modifying their geometry. Furthermore, acoustic perturbations can tune thickening in these applications in multiple modes, thus providing flexibility in implementation (See SI Section VII for data on acoustic perturbations in the vorticity direction). Overall, our method has laid a strong foundation to robustly design smart transport systems that handle shear-thickening fluids.

We thank Anton Paar for use of the MCR 702 rheometer through their VIP academic research program, Abhishek Shetty for help with the rheometer, the Kirby group and the Cohen group for their insightful suggestions, Prof. Wolfgang Sachse at Cornell for lending equipment and useful discussions, and the Amit Lal group at Cornell for providing the impedance analyzer. This work is supported by NSF CBET award numbers 1804963, 1232666, and 1509308, partially supported by the Center on the Physics of Cancer Metabolism through Award Number 1U54CA210184-01, and performed in part at the Cornell NanoScale Facility, an NNCI member supported by NSF Grant NNCI-1542081. PS and MR contributed equally to this work.