Mechanical History Dependence in Carbon Black Suspensions for Flow Batteries: A Rheo-Impedance Study.

We studied the effects of shear and its history on suspensions of carbon black (CB) in lithium ion battery electrolyte via simultaneous rheometry and electrical impedance spectroscopy. Ketjen black (KB) suspensions showed shear thinning and rheopexy and exhibited a yield stress. Shear step experiments revealed a two time scale response. The immediate effect of decreasing the shear rate is an increase in both viscosity and electronic conductivity. In a much slower secondary response, both quantities change in the opposite direction, leading to a reversal of the initial change in the conductivity. Stepwise increases in the shear rate lead to similar responses in the opposite direction. This remarkable behavior is consistent with a picture in which agglomerating KB particles can stick directly on contact, forming open structures, and then slowly interpenetrate and densify. The fact that spherical CB particles show the opposite slow response suggests that the fractal structure of the KB primary units plays an important role. A theoretical scheme was used to analyze the shear and time-dependent viscosity and conductivity. Describing the agglomerates as effective hard spheres with a fractal architecture and using an effective medium approximation for the conductivity, we found the changes in the derived suspension structure to be in agreement with our qualitative mechanistic picture. This behavior of KB in flow has consequences for the properties of the gel network that is formed immediately after the cessation of shear: both the yield stress and the electronic conductivity increase with the previously applied shear rate. Our findings thus have clear implications for the operation and filling strategies of semisolid flow batteries.

Introduction

In light of climate change, recent years have seen a rapid adoption of renewable energy production.[1] Because of their inherently variable nature, renewables have placed considerable strain on the power grid which must match energy production to demand. A possible solution to this problem is to store energy in batteries.[2] Semisolid flow batteries (SSFBs), a recently developed configuration,[3−10] are considered to be especially promising for such applications. SSFBs use two fluid electrodes, an anolyte and a catholyte, in place of traditional solid electrodes. The use of fluid electrodes decouples the energy of a SSFB, which depends on the size of storage tanks, from its power, which depends on the size of the reactor. Additionally, SSFBs may allow easy lifecycle management through the modification or replacement of their fluid electrodes.

SSFB electrodes are mixtures of conductive nanoparticles (CNPs) and electrochemically active particles (EAPs) dispersed in an electrolyte solution. In most SSFBs, the EAPs intercalate and deintercalate lithium whereas the CNPs wire the EAPs to the current collectors. The continuous phase, a mixture of linear and cyclic carbonate solvents with a high concentration of dissolved lithium salt, provides an ion-conducting medium and a source of lithium ions. The CNPs used in SSFBs are typically superconductive carbon blacks such as Ketjen black (KB) or Timcal SuperP.[3,11] These carbon blacks are submicrometer-sized permanently fused aggregates of hollow spherical subunits.[12,13]

The self-assembly of the CNPs, which depends on their morphology, colloidal interactions, and shear conditions, is of crucial importance to both the electrical and mechanical performance of SSFBs. For colloidal carbon black (CB) units in SSFB media,[3,5,11,14] the van der Waals attractions should be dominant because their electrostatic interactions are strongly screened.[11,15] As a result, CB particles tend to form large cohesive structures. At rest (assuming a high enough concentration), they form a space-filling network that can conduct electrons and suspend EAPs against gravity through yield stress.[3,11] In flow (e.g., pumping and stirring), this network will be broken down into agglomerates. This leads to shear thinning and much lower electronic conductivity.[3,11] Henceforth in this article we differentiate a permanently fused primary aggregate from a reversibly flocculated cluster of these particles, which we will call an agglomerate.

Given the novelty of SSFBs, significant developments are still required to optimize their performance. Mechanical protocols should be a part of this optimization: fluids containing adhesive particles generally produce nonequilibrium structures that can depend on mechanical history.[16−19] In this article, we focus on how shear and its history influence the electrical and rheological properties of CB suspensions in an SSFB solvent.

Although the case of CB in an SSFB solvent is rather new,[3,11,15] we anticipate that some existing insights into the behavior of CB suspensions will also be applicable to our system. Suspensions of reversibly agglomerating colloids, including those of CB, are known to be strongly shear- and history-sensitive.[11,20−25] The storage modulus and yield stress of CB gels in oils have even been found to depend predictably on the preshear used to prepare them.[23,24] In flow, the viscosity of agglomerating colloidal suspensions typically decreases in time after an increase in shear stress or rate, a well-studied property known as thixotropy.[21,22] The opposite effect, a temporal increase in viscosity after an increase in shear stress or rate, known as rheopexy or antithixotropy, is much less common but has been observed in CB suspensions.[23,26−28] On a mechanistic level, rheopexy is thought to be caused by flow-induced flocculation[26] or by the ability of fractal structures to rearrange (when shear is lowered) into more densified agglomerates.[23,24]

The electrical impedance of CB suspensions in shear flow has been studied much less often, but also here a dependence on the shear rate was found.[11,14,29] The study of the suspension’s rheology and impedance in conjunction capitalizes on the tight connection of both behaviors to the (dynamic) microstructure of the agglomerates, which can be difficult to measure with optical or scattering techniques, especially for concentrated suspensions.

The present study makes use of a home-built rheo-impedance setup to characterize the influence of the mechanical history on both rheological and electrical properties. This approach addresses both the practical aspect of optimizing SSFB performance via mechanical protocols and the more fundamental aspect of understanding the underlying processes. Both stepwise changes in shear rate and prolonged shear are explored to identify transient responses. The cessation of shear is included as a special case, where the system is left with only Brownian forces to possibly reorganize its structure. Our analysis of the structural changes at the microscale is supported by a theoretical model in which the viscosity of effective hard-sphere agglomerates[30] is combined with an effective medium theory for the conductivity.[31]

Materials and Methods

Suspension Preparation

Ethylene carbonate (EC) and dimethyl carbonate (DMC) were obtained from Sigma-Aldrich (anhydrous, 99%+ purity). Binary mixtures of EC and DMC were 1:1 by mass. LiPF6 was purchased from Alfa Aesar (98% purity). Ketjen black EC 600JD powder (KB) was obtained from AkzoNobel (The Netherlands). Thermax N990 powder was donated by Cancarb. SEM samples were dispersed in acetone, dried on silica wafers, and then imaged. All other sample preparations and experiments were carried out in an MBraun argon-filled glovebox (O2, H2O below 5 ppm).

The two types of CB used in this study differ in morphology (see also Figure ): for N990, the unit particles are more or less compact spheres, whereas for KB, they are more fractal-like as a result of fusion between spherical subunits. In this article, we designate the KB units as aggregates, whereas clusters formed via reversible flocculation are called agglomerates. All suspensions of KB (of key interest for this work) were prepared at a concentration of 1% by weight (wf = 0.01). Because the KB aggregates are highly porous and their subunits are hollow shells, the occupied volume fraction in suspension is much higher than wf. A crude estimate can be obtained by multiplying wf by the ratio of the density of graphite (2267 kg/m3) to the tap density of KB (100 kg/m3), giving a volume fraction of 23%. Particles were first wetted by the pure (binary) solvent in polypropylene containers. After 8 h, LiPF6 salt was added via a concentrated solution to bring its concentration to 1 M (viscosity 4 mPa·s[32]). After hand shaking, at least 8 h was allowed to let the particles equilibrate. Samples were then homogenized by rotor stator mixing (Ultraturrax) at 15 000 rpm for 2 min and subsequently loaded into the rheometer.

Figure 8

Morphology of KB (left) and N990 (right). Particles are agglomerated by drying.

In exploratory experiments, we found similar trends as reported in this article. However, the quantitative behavior of KB suspensions in SSFB media appeared to be sensitive to preparation protocols. This could be related to the poor wetting by the solvents, the sensitivity of the media,[15] and the strong dependence of macroscopic properties on the volume fraction (S.I. Figure S1 and ref (11)). Thus, for KB we performed all reported measurements (except those in Figure ) on the same sample.

Figure 4

Intrinsic flow curve of agglomerates formed at different indicated shear rates (lines to guide the eye). The data point indicated by (x) represents the viscosity at 10 000 s–1 after prolonged shearing.

Rheoimpedance Measurements

Parallel rheological and electrical measurements on CB suspensions were performed on a stress-controlled Haake Rheostress RS600 rheometer with a home-built adaptation for electrical impedance spectroscopy (EIS) measurements (Figure ). A 60-mm-diameter parallel plate geometry was designed with both stainless steel plates also acting as electrodes. The rheometer rotates and measures torques on the upper geometry. To allow a frictionless low-noise electrical connection to the rotating plate, a mercury-based solvent trap was designed. A glass disk was used to electrically isolate the lower plate from the rheometer body. Because of its high thermal conductivity, glass also allows accurate temperature control (at 25 °C for all experiments). Argon was used to drive the rheometer air bearing.

Figure 1

Schematic of rheo-impedance setup. (Gray) Measuring geometry with stationary lower and rotating upper plates. (Light gray) Mercury solvent trap. (Red) Sample. (Blue) Isolating dielectric. (Brown) Polypropylene body. (White) Cover.

EIS was performed using a four-terminal configuration. The top plate was excited by a sinusoidal voltage, and the bottom was connected to the virtual ground of a transimpedance preamplifier (HF2TA, Zurich Instruments). A buffer preamplifier (HF2CA, Zurich Instruments) was used to measure the potential difference between the plates. An impedance spectroscope (HF2IS, Zurich Instruments) was used to extract the complex impedance from the current and voltage signals. AC frequency sweeps were performed from 10 to 100 MHz (total time 550 s) with a maximum applied voltage amplitude of ∼100 mV, which was within the linear response range of the samples. The frequency-dependent impedance of the empty measurement setup was calibrated out using the open short technique. Custom LabVIEW programs were used to synchronize rheological and electrical measurements.

Rheo-impedance experiments were performed with a plate–plate gap of 250 μm. Suspensions of CB have been reported to exhibit nonideal rheometric behaviors such as wall slip, shear banding,[23] and vorticity alignment of flocs,[33] especially during start up and in low-shear-rate regimes. We avoided these effects in shear flow experiments by employing relatively high shear rates, where the samples behave as low-viscosity liquids. To obtain a reproducible starting state (Figure ), the suspensions were presheared at 10 000 s–1 for 200 s, after which the shear rate was slowly ramped down to 316 s–1 and then held for 200 s. Step changes in shear rate were followed by a dwell time (τ) that was different per type of experiment.

Figure 2

Schematic of the general experimental protocol for preshearing and subsequent measurements of the electrical impedance and yield stress.

Measurements of the impedance during shear were made at a frequency of 0.5 Hz to allow real-time monitoring of the electronic resistance (Section 2.3). Measurements after the cessation of shear were made after a 200 s rest period. Impedance spectra were measured before the yield stress. In the measurement of the latter, the shear stress was ramped up at 62 s per stress decade while measuring the strain. The log(strain) versus log(stress) curve was fitted with two straight lines, and the stress at the intersection was taken to be the yield stress.[34]

Impedance Analysis

Because flow battery electrodes contain conductive carbon in an electrolyte solution, they are capable of both ionic and electronic conduction. To monitor the KB network electrically, it is necessary to separate these contributions. This is done using AC impedance spectroscopy. The KB suspensions have two ionic charge carriers (Li+, PF6–) and an electron-conducting carbon network. During measurements, the sample is in contact with electrodes that are ionically blocking but electronically reversible. Ignoring geometric capacitance, the response of this system can be simplified[35] to that of an ionic resistor and an ionic double layer capacitor, both in parallel to the electronic resistance of the carbon black network. Thus, the high-frequency real impedance gives us a parallel combination of both resistances whereas the low-frequency real impedance gives the electronic resistance. The data in Figure illustrate that the behavior of our KB suspension at rest can be analyzed with this approach. The very high frequencies (Figure insets) are affected by uncompensated parasitics.[36]

Figure 3

Nyquist with equivalent circuit (left) and Bode (right) plots of the impedance response of the KB suspension. The insets show a zoom around high-frequency data (>10 kHz). Filled symbols: KB at rest. Open symbols: KB under 1000 s–1 shear (only two frequencies measured). The 50 Hz points have been removed.

Under strong shear, it is unlikely to have a continuous percolated network. (For low shear rates, dynamic percolation may be possible.[37]) However, a measurement at two frequencies on the same KB suspension, but now under high shear (open symbols of Figure ), shows that the low-frequency real impedance is finite whereas the imaginary impedance is close to zero. This implies that the sheared suspension has a finite electronic resistance.[35] Our measurements on other carbon black samples show that the overall impedance response under shear is similar in shape to that at rest (Figure S2 in S.I.), justifying our use of a single low frequency to probe the electronic resistance. Additionally, this is in qualitative agreement with earlier findings on carbon black filled polymers, which are known to conduct electrons via two mechanisms. (At high volume fractions, when the network is continuous, conduction is graphitic.[38] The resistance is still determined by constriction and ohmic contributions.[39]) When the network is broken, the electrons can tunnel across small gaps (<10 nm) modulated by thermal voltage fluctuations.[11,12,38,39] For larger gaps, however, the intrinsic electronic conductivity of the medium dominates.[39]

Results and Discussion

Intrinsic Flow Curves

We first consider the mechanical properties of our material in flow. To characterize KB, we measured a series of so-called “intrinsic” flow curves[23] where each of the curves corresponds to a unique preshear rate (PSR) that was maintained for 1000 s in order to reach a steady state. Data points within each curve were then obtained by stepping from the PSR to the shear rate of interest and recording the viscosity quickly. (As shown later, the viscosity almost immediately reached an extremum; it was this extremum that was recorded.) All intrinsic flow curves (Figure ) are strongly shear-thinning, as expected for agglomerating suspensions.[27,30,40−43]

Their dependence on the PSR reveals interesting information that is specific to our KB system. First, the three curves with the highest PSR coincide, indicating that a PSR of 1000 s–1 or more brings the material to a similar microstructural state. This suggests that the agglomerates are largely broken down under these conditions: generally for agglomerating systems a stronger shear flow causes more breakup. After prolonged (2000 s) shearing at the highest accessible shear rate of 10 000 s–1 (from an arbitrary starting state), the sample showed a similar viscosity to the corresponding points of the intrinsic flow curves for PSR ≥ 1000 s–1: this again suggests that the agglomerates are largely broken down (i.e., into primary aggregates) by prolonged shear above 1000 s–1. Second, it becomes clear that in order to reach a steady state even the highest preshears need to be maintained for a significant duration, i.e., much longer than the data acquisition time (of O(1 s)) of the intrinsic flow curve. This is shown by the curves with PSR ≤ 316 s–1, for which the viscosities are different from the curves at higher PSRs, even for measurement shear rates ≥1000 s–1. Third, the fact that the intrinsic flow curves with smaller PSRs have lower viscosities is remarkable. In earlier suggested mechanistic pictures,[30] where shear leads to both breakdown and the densification of agglomerates into smaller units, a lower PSRs should have resulted in higher viscosities, which is the opposite of what we observe. This points to a mechanism that can enhance the viscosity via agglomerate breakup in KB suspensions. The remarkable fact that the lowest possible viscosity (21 mPa·s) for our system was measured after the smallest preshear rate (100 s–1) also points in this direction.

Transient Behavior

To examine this intriguing behavior, we performed experiments in which the material’s response to shear rate steps was followed over time.[21] Both the viscosity and the electronic resistance (Section 2.3) were measured in parallel. All steps in the shear rate were taken from the same reference condition, for which we chose 316 s–1, based on Figure . Some typical results of these measurements are found in Figure , where both steps up (to 562 and 1000 s–1) and down (to 177 and 100 s–1) are illustrated.

Figure 5

Evolution of viscosity and low-frequency real impedance of the KB suspension after a shear rate step from 316 to 100 (black), 177 (red), 562 (green), and 1000 s–1 (blue). The inset is a close-up view of the response of the electronic resistance around the step.

Both the viscosity and the electronic resistance show a two-staged response. The initial response, a large increase in viscosity and a small decrease in electronic resistance for a step down in rate (transition from region I to region II) is fast, with a time scale of O(1 s). It is not resolvable, considering the short time that the rheometer needs to step the shear rate and the low frequency of the impedance measurement. For the steps to higher shear rates, the opposite trends were observed. The secondary response (region II) is significantly slower O (100–1000 s). The slow changes in the viscosity and resistance appear to have similar time scales, suggesting that they both probe changes in the microstructure of the KB.

Remarkably, for both signals, the secondary response opposes the direction of the initial response. For the viscosity, this amounts to a relatively small correction (rheopexy, shear thinning), but for the electronic resistance, the secondary effect is strong and causes a reversal of the overall effect. The shape of the decays of the viscosity and electronic resistance is well described by stretched exponentials; such behavior has been found for the viscosity of thixotropic systems.[21] The opposing change in the second stage is found irrespective of whether the shear rate is increased or decreased. Yet the direction in which the shear rate is changed still determines the sign of all changes. This is consistent with a picture in which all structural transitions can be reversed via the shear rate.

Mechanistic Picture

We interpret our findings in Figures and 5 with a picture in which KB reversibly agglomerates via two distinct mechanisms (each with its own time scale and influence on the microstructure). Both the viscosity and electronic resistance depend on the concentration and morphology of the agglomerates, but in different ways.

Flow curves of weakly agglomerating suspensions have been successfully described by modeling the agglomerates along with their immobilized solvent as effective hard spheres.[40,44] Shear thinning is then explained via a lowering of the effective hard sphere volume fraction. As the shear rate is increased, higher shear stress causes the agglomerates to break down into smaller structures. Because of the fractal build up, the latter then occupy less total volume (Figure ). Although (strictly speaking) perfectly fractal scaling is rarely observed, it is has often been found that the structure of agglomerates can be fairly well described with this concept.[30,40,43,44]

Figure 6

Cartoon of fractal agglomeration and imperfect percolation (left) in shear flow.

To understand the (changes in) electronic conductivity of the suspension under flow, we need to take into account the formation of (transient) pathways with missing links (Figure , left). As such gaps dominate the resistance (Section 2.3), the macroscopic electronic resistance of the sheared suspension should be determined by the concentration of agglomerates and their typical distance of closest approach.

We now proceed to a qualitative interpretation of Figure , focusing first on the instant process. The jumps in viscosity and electronic resistance may be explained by agglomeration due to lower shear forces when the rate is stepped down. The immediate sticking between colliding fractal entities leads to more open agglomerates that enclose (and hence immobilize) more solvent, causing the hydrodynamic volume fraction and thus the suspension viscosity to increase. (Note: if the shear step were small, then a binary collision would lead to an unstable agglomerate, which would be quickly eroded.) The same agglomerate restructuring also changes the ensemble of interrupted conductive pathways. The formation of fewer but larger agglomerates after a step down in shear rate diminishes both the average gap between agglomerates and the number of gaps an electron needs to traverse along a conductive path. Both of these effects should lower the electronic resistance. However, they are counteracted by a decrease in the number of electronic pathways. This could explain why the overall decrease in resistance is small.

As shear is continued after the step down in rate, the population of agglomerates is rejuvenated continually: bonds between KB units are incessantly broken and formed by shear. This provides a means for the agglomerates to overcome kinetic barriers and discover more stable states. As the unit particles (and their agglomerates) are (roughly) fractal in nature, the randomness of collisions will occasionally cause two particles to interpenetrate (Figure ). Because of the enhanced contact area and generally lower aspect ratio, these agglomerates will be more stable (compared to those formed by peripheral contacts). However, if the shear stress gets high enough, they can be torn apart again, making the process reversible. Clearly, a consolidated interpenetration should lead to an overall densification of agglomerates, leading to a lower hydrodynamic volume fraction and hence a reduced viscosity. It also follows from mass conservation that the average gap between agglomerates will increase in this scenario. Additionally, the number of conduction pathways will decrease. Both effects should lead to the observed increase in the electronic resistance (after a stepdown in shear rate).

Figure 7

Cartoon of interpenetrating agglomeration.

The dependence of the interpenetrating agglomeration on shear rate can be understood as follows: the stability of interpenetrating KB units will depend on how deeply they are lodged within the agglomerate because only the peripheral zone will be exposed to local shear flow. This explains why a complete reversal of the process can only be achieved at high shear rates: here, the shear stresses are high and the agglomerates are small, thus exposing all KB units to shear. (This is also why we defined our samples by preshearing them at 10000 s–1.)

In our proposed mechanism of interpenetrating agglomeration, the shape of the KB units plays an important role. To test this hypothesis, we performed an experiment comparable to that in Figure but now with spherical unit particles (N990). This system was also strongly shear-thinning, indicating reversible agglomeration. A weight fraction of 30% was chosen to achieve a high shear viscosity similar to that of KB. Assuming that the hydrodynamic volume fractions of both suspensions are comparable when they are broken down to the unit particles, the only relevant difference should be the shape of the unit particle (Figure ).

The electrical resistance of the N990 suspensions turned out to be too high to be accessible at experimental frequencies (<10 mHz). As shown in Figure , similar to KB the N990 suspension also exhibits an immediate jump in viscosity when the shear rate is stepped down, consistent with instant open agglomeration. However, when shear is continued, the viscosity slightly increases (thixotropy, shear thinning), in contrast to KB where it decreases. This supports our explanation that the fractal morphology of KB is responsible for the occurrence of interpenetrating agglomeration.

Figure 9

Evolution of the viscosity of an N990 suspension after a shear rate step from 316 s–1. The colors correspond to shear rate steps described in Figure . Viscosity axis breaks and their scales are not uniform. One tick corresponds to 1 mPa·s.

Modeling

To further examine the mechanistic explanations proposed so far (as sketched in Figures and 7), we combined previously developed concepts for describing the viscosity of agglomerating systems[30] and the conductivity of filled polymers[39] into a simplistic theoretical scheme. A complete description of this scheme is given in the Appendix. The scheme uses the measured viscosity and electronic resistance as inputs, and the collision radius and fractal dimension of the agglomerates (along with other derived geometrical properties) are output. The Krieger–Dougherty equation[45] is used to calculate the hydrodynamic volume fraction from the relative viscosity, while part of the Potanin model[30,46] is used to set a relation between the agglomerate radius and the fractal dimension. An effective medium approach for the electronic resistance is used to construct an additional relation between the agglomerate radius and the fractal dimension. In this approach, agglomerates are assumed to occupy centers in a lattice. The resistivity of the suspension then depends on the resistance between adjacent occupied lattice sites and depends asymptotically on the fraction of occupied sites. Because the agglomerates have a fractal architecture, the resistance of a conductive contact depends not only on proximity but also on the local (sub)unit particle density in the outer shells of the agglomerate.

The above relations are solved numerically for each measured viscosity and electronic resistance combination to yield a solution (fractal dimension, collision radius). Several other structure parameters can then be derived. Application of the scheme to the measurements of Figure leads to the data shown in Figure . We observe that the trends used in our qualitative interpretation are reproduced by the model. For steps down in shear rate, the fast restructuring response leads to larger and less dense agglomerates while reducing both their concentration and spacing. During the slow interpenetrating agglomeration process, densification leads to a strong increase in the interagglomerate distance (and hence electronic resistance). For steps up in shear rate, the opposite trends occur, in line with expectations.

Figure 10

Analysis of the experiments of Figure with our theoretical model. Shear rate steps are made from 316 to 100 (black), 177 (red), 562 (green), and 1000 s–1 (blue).

Cessation of Shear

We revisit the experiment shown in Figure , now focusing on the last stage (region III) in which the flow is terminated after 1000 s of maintained shear at a constant rate. The (still monitored) electronic resistance shows an immediate sharp drop, indicating the formation of a space-filling network. After 200 s, the resistance reaches a steady state, justifying measurements of the full impedance spectra (Figure S3 in the S.I.) and the yield stress. Remarkably, both the low-frequency real impedance and the yield stress show a systematic dependence on the previously applied shear rate, with higher preshear leading to a higher yield stress and a lower electronic resistance (Figure ). This effect is of strong significance to SSFBs where a high yield stress and low electronic resistance are both required for efficient operation.

Figure 11

Rest impedance and yield stress of KB suspensions after preshearing for 1000 s at indicated rates. (Lines are drawn to guide the eye.) Figure was used as the preshear protocol.

To explain these trends, we propose that the agglomerate structures as defined by the shear history are the building blocks of the solid network that forms after shear cessation. In other words, the agglomerates remain intact up to the point where they connect with other agglomerates to form a space-filling network. This freezing-in effect is plausible, as agglomerates can still undergo Brownian motion[47] while thermal forces are too weak to break bonds between the KB units.

The qualitative dependences of the electronic resistance and yield stress on the preshear rate can be rationalized using the same shear (history) induced changes in the agglomerate structure as discussed earlier. High preshears result in small agglomerates, which assemble into a highly branched network when shear is stopped. Conversely, preshearing at low rates leads to networks with fewer and thicker branches because the agglomerates were larger in size and fewer in number. The contact areas between the agglomerates that form the low preshear gel should be similar (peripheral contacts) and thus no more effective at transferring charge and mechanical force than those of the high preshear gel. However, in the latter material, the number density of branches will be higher, leading to a higher yield stress and a lower electronic resistance.

To further examine the role of the secondary response (region II in Figure ) in light of the results of Figure , we explored the effect of shear cessation at different times (Figure left). As before, a 200 s rest time was given between the cessation of shear and the measurement of the impedance spectrum (complete spectrum in S.I. Figure S4). The yield stress was subsequently measured. As shown in Figure right, both the final electronic resistance and yield stress were found to evolve in an S shape from their starting value. (In the case of the step to 100 s–1, steady state is not completely reached even after 1000 s.) Remarkably, their values for just 1 s of shear after the step are fairly similar to those without the step. This indicates that the structures of the gelled networks must also be very similar, despite the instant agglomeration (viscosity and resistance jumps) that occurred during the 1 s step. This suggests that the instant agglomeration process and gel formation on shear cessation have similar mechanisms. Additionally, it shows that the process responsible for the secondary response (interpenetration) is almost exclusively responsible for the change in rest network properties.

Figure 12

(Left) Evolution of viscosity and low-frequency real impedance after a shear rate step from 316 to 1000 (up) or 100 s–1 (down) with shear cessation at different (indicated) times. (Right) Rest network resistance and yield stress measured after shear cessation at the indicated time. (Lines are drawn to guide the eye.)

Conclusions

The influence of shear and its history on the rheology and electrical impedance of Ketjen black suspensions in SSFB solvent was studied with two objectives: to understand the microscopic origin of the shear-history-dependent behaviors and to identify the implications for the performance of SSFBs. Both aspects revolve around the microstructure of the agglomerated suspension and how it adapts to mechanical conditions. Given the difficulty of measuring relevant structural properties directly via microscopy or scattering, our simultaneous study of the rheology and electrical impedance provided a good alternative method.

The fast and slow transient responses to a change in shear rate clearly indicate that KB agglomerates exhibit two buildup (breakup) mechanisms after decreasing (increasing) the shear rate. The fast response was attributed to open agglomeration, and the slower response was ascribed to an interpenetration of the fractal-like KB units. The combination of these processes leads to a shear- and time-dependent agglomerate structure that gets frozen in on cessation of flow, with a higher preshear leading to a stronger gel with a lower electronic resistance, both critical to SSFB performance.

An implication of our findings is that for SSFBs operated in the absence of flow, a high preshear (through fast pumping or stirring) would be the best for minimizing the electronic resistance and maximizing the yield stress. For SSFBs operated in continuous flow, a high flow rate would minimize the electronic resistance. We expect that our findings will be of relevance not just to SSFBs but also to other systems that use fractal carbon blacks such as flow supercapacitors. Additionally, we foresee electrochemical cycling to play an important role in the rheo-impedance of fluid electrodes (and vice versa) and are currently investigating the same.

Table 1

Parameter Estimates Used in the Model

parameter	value
η0	4 mPa·s	ref [32]
a	5 nm	ref [12]
ϕmhyd	0.71	high-shear limit for hard spheres
β	0.72	ref [49]
z	12
df primary aggregate	1.8	ref [13]
Rh primary aggregate	90 nm	DLS (S.I. Figure S5)