Transient Shear Banding in Time-dependent Fluids
We study the dynamics of shear-band formation and evolution using a simple rheological model. The description couples the local structure and viscosity to the applied shear stress. We consider in detail the Couette geometry, where the model is solved iteratively with the Navier-Stokes equation to obtain the time-evolution of the local velocity and viscosity fields. It is found that the underlying reason for dynamic effects is the non-homogeneous shear distribution, which is amplified due to a positive feedback between the flow field and the viscosity response of the shear thinning fluid. This offers a simple explanation for the recent observations of Transient Shear Banding in time-dependent fluids. Extensions to more complicated rheological systems are considered.
pacs:83.60.Rs, 83.60.Wc, 83.50.Ax
The property that characterizes complex fluids is their non-trivial rheology, shear rate – stress relation. They are generally further categorized into shear thinning or shear thickening fluids. Both cases are additionally complicated by time dependence. Due to the stress-shear interaction, already small perturbations in the local stress can result in a positive feedback with the flow promoting shear instabilities in each case Head et al. (2002); Fielding and Olmsted (2004). The understanding of complex fluids is of enormous importance for many practical applications Braun and Rosen (1999) and the theory touches on many branches of physics. Recent advances allow to follow the suspension local velocity during a standard rheological experiment Bonn et al. (2008); Manneville (2008). Quantifying the local flow field simultaneously with rheological measurements gives the possibility to measure both the intrinsic and apparent rheology. This has lead to the discovery that a heterogeneous shear distribution in samples during such tests is ubiquitous. Shear banding Ovarlez et al. (2009) has been observed in many systems composed of substantially different building blocks, such as colloidal glasses, wormlike micelles, foams and granular matter Schall and van Hecke (2010). The current viewpoint, both phenomenologically and theoretically, is that a non-monotonic intrinsic flow curve, is what is common to most of these materials Ovarlez et al. (2009); Porte et al. (1997), but also other mechanisms have been suggested Besseling et al. (2010).
A branch of complex fluids are the simple yield stress fluids Ragouilliaux et al. (2007). These materials do not show aging phenomena (thixotropy). Therefore they are expected to have a monotonic intrinsic flow curve, and a steady state without shear bands Fall et al. (2010). However, recent experiments Divoux et al. (2010) display shear banding during startup flows in a rotational rheometer indicating time-dependent behavior. These, so called transient shear bands, can be very long lasting, but eventually vanish with a homogeneous steady state. The transient shear banding phenomenon tests our fundamental understanding of non-Newtonian fluids, and is also important for industrial processes and simply for understanding usual rheological measurements. A particular feature of the transient shear banding is that it appears to exhibit scaling familiar from critical phenomena: the time it takes for the transient to disappear (fluidization time ) is a power-law function of the shear rate or applied stress Divoux et al. (2010).
For the reasons leading to transient shear banding, see the summary by Adams et al. Adams et al. (2011), three main candidates are offered: i) with certain parameters, the dynamical equations are unstable amplifying small perturbations, which slowly quench towards the homogeneous steady state Fielding and Olmsted (2003). ii) The fluid flow curve is time dependent, and can, at different times, have non-monotonic shape Hayes et al. (2010). iii) Elastic stress overshoots cause instability in the flow Marrucci and Grizzuti (1983). Some theoretical models appear to produce transient shear banding including shear transformation zone theories Manning et al. (2009), a modified soft glassy rheology model Moorcroft et al. (2011), a simplified fluidity model Moorcroft et al. (2011), and a mesoscopic model of plasticity Jagla (2010). Such models of transient shear banding share the property of time-dependent reduction of the local stress under shear Manning et al. (2009); Moorcroft et al. (2011); Jagla (2010) as explained in Ref. Moorcroft and Fielding (2012). Literature reports experimental evidence of time-dependent rheology in carbopol gels Divoux et al. (2010) especially at small shear rates, as well as in other simple yield stress fluids Møller et al. (2009), appearing as slight hysteresis in the flow curves. Further details of the time-dependent flow curve hysteresis related to carbopol gels is reported in Ref. Divoux et al. (2011). Stronger hysteresis is inherent to thixotropic fluids Møller et al. (2009). Indeed, transient shear banding has been recently found also there well above the critical shear rates Martin and Hu (2012) indicating that transient shear banding could be a general property of soft glassy materials.
Motivated by these findings, we consider here a structural model to find the main ingredients of transient shear banding. It is for a simple time-dependent Newtonian fluid, on purpose neglecting other complications present in yield stress fluids and thixotropic fluids, such as elastic and yield stresses and critical shear rates. Spatial resolution, necessary for the transient shear banding, is obtained when the rheological model is coupled to the Navier-Stokes equation for laminar flow in a concentric cylinder Couette device. The transient shear banding here occurs due to the shear rate - viscosity coupling in the Navier-Stokes solution. We will demonstrate that such transients initiate since the shear stress exhibits a finite gradient in any rotational geometry. This gradient is amplified by the interaction with the shear thinning fluid. In the model, the timescales associated with the relaxation of the structure depend on the shear rate. A higher shear rate implies faster dynamics creating, with the coupled Navier-Stokes equations, an amplifying effect during the startup flow.
Model. — Most microscopic descriptions of the structure of a complex fluid are purely phenomenological, and based on a kinetic relaxation equation. The simplest -models stipulate an evolution equation for a structural parameter, Bonn et al. (2008); Mewis and Wagner (2009). This describes the internal order, such as the state of aggregation in colloids or the alignment of particles with shear Bonn et al. (2008). Usually, shear and/or temperature influences the temporal evolution of Mewis and Wagner (2009). The macroscopic rheology is obtained coupling to a constitutive equation. Variants of -models can be used to describe the flow curves of thixotropic, simple yield stress fluids, and shear thinning fluids Møller et al. (2006); Bonn et al. (2008). They generally assume homogeneous flow, and thus are not applicable to shear bands.
In aggregating suspensions and microgels a portion of the liquid is trapped due to the presence of solid structures Heath et al. (2006); Piau (2007). This is infact also the interpretation shared in the literature for carbopol microstructure, which is observed to be formed of elastic sponge-like elements Piau (2007); Oppong et al. (2006). Therefore, instead of using an abstract structure parameter (with no direct physical interpretation) it is more sensible to select the immobilized volume fraction , which describes the jammed fluid, as it has been done in approaches using the Population Balance Equations Heath et al. (2006). The simplest time-dependent rheology follows by including a temporal relaxation of the volume fraction to a steady state. Unlike the relaxation rate, the steady state volume fraction is independent of the applied shear rate corresponding to a Newtonian fluid. Such a model can be thought as a Taylor-expansion of the dynamical equation for the volume around the steady-state. The resulting kinetic equation is given
where () is the kinetic constant for the shear growth (destruction), is the magnitude of the shear rate, and (set to unity) both relate to the volume fraction sensitivity to shearing. , , and describe the growth of jammed volume fraction due to the shear independent motion of the structure elements. The special case of , presents a simple yield stress fluid, gives a non-monotonic flow curve indicating thixotropy, and produces shear thinning behavior. In what follows, we fix the parameters and to 0.665 and 1.0, respectively. Since we are concentrating on a minimum model showing transient shear banding, in the following we set , and parametrize the initial volume fraction instead of specifying the sample history and the associated parameters. It is well known that in any practical experiment the initial state depends on the sample history, and the shear independent structure dynamics. Therefore special attention has to be payed to the measuring protocol. The shear independent terms dominate the structure evolution at small shear rates making the flow curve non-Newtonian, but have negligible influence in the transient shear banding regime. There, it is reasonable to assume that both, the structure growth and destruction terms are determined by the shear rate justifying such approximation. The exact form of the rate kernels are presently unclear and here we choose a simple power-law dependence. The scaling factors and make the equation unitless, and could be incorporated to the kernels , , and equally well.
We use the well-known Krieger-Dougherty constitutive equation Krieger and Dougherty (1959), that reads
where is the jamming volume fraction, the liquid viscosity, and gives the scaling of the viscosity. From now on these quantities are set to (random sphere packing), mPas (water), and Quemada (2008). These values relate to the initial conditions at which the transient shear banding appears, but are irrelevant e.g. for the scaling relations of the fluidization times.
The steady state volume fraction is . As a consequence of leaving the shear independent terms out of the kinetic equation, the resulting steady state volume fraction is independent of the shear rate. As mentioned before, these shear independent processes are the ones that drive the system out of steady state implied by Eq. (1) when not sheared. This effect is incorporated to the model by simply initializing the system to the desired volume fraction before starting the structure evolution. For shear thinning fluids having no yield stress, the present approximation is valid for all shear rates. Fluids having a yield stress, are properly described in the range where shear localization does not appear i. e. the fluidization is complete. Besides fixing the steady state, and determine the rate of relaxation. Plugging to the constitutive equation and applying the Newtonian assumption , gives the Newtonian steady state flow curve in Fig. 1. The stress shows time and shear rate dependent exponential relaxation (the Inset), when started from , since the structure evolution follows Eq. (1).
For the spatial resolution we couple the model with a continuum description of the fluid flow in a Couette rheometer geometry Roussel et al. (2004). The incompressible Navier-Stokes equation for laminar flow in that has the analytical one-dimensional (radial) solution Phillips et al. (1992)
where is the radial distance from the cell center, is the local shear rate, is the angular velocity of the inner (outer) cylinder, and is the radius of the inner (outer) cylinder. Eq. (3) implies that only the relative angular velocity of the two cylinders matters. In order to mimic the experiments of Ref. Divoux et al. (2010) only the inner cylinder is rotated and the radii are set to mm and mm. The Eq. (3) implies a spatial dependence for the volume fraction . To this end, a radial discretization is applied with a uniform grid with N sampling points, at each of which a separate is evolved. Thus, the time evolution of the local flow field is obtained when the Equations (1) and (3) are iteratively solved using a forward Euler algorithm under a constant global shear rate. Solving the structure and flow field evolutions this way assumes that the time-scales of the inertial effects and the structure kinetics is very different. The flow field quickly adapts the changes caused by the slowly varying time-dependent viscosity.
The global shear rate is defined as the radial average shear rate for a Newtonian fluid in the Couette geometry. This reads
Results. — Varying the initial volume fraction and the kinetic exponent uncovers three different startup flow scenarios as illustrated in the schematic phase diagram shown in Fig. 2. An almost homogeneous relaxation of the volume fraction profile arises at initial volume fractions above the steady state value and low kinetic exponents (Fig. 2(a)). Increasing either makes the transient shear bands to appear (Fig. 2(b)). At initial volume fractions below the steady state, there is no transient shear bands (Fig. 2(c)). The boundaries of that phase (grey area) are qualitative, as the transition from the homogeneous relaxation Fig. 2(a) to the transient shear banding Fig. 2(b) occurs smoothly.
For a given value of , the initial value is chosen in order to be in the transient shear banding region (grey area in Fig. 2). Thus, the simulation starts close to jamming, at . Fig. 3(b) shows the local shear rates corresponding to the evolution of the volume fraction plotted in Fig. 3(a), from the numerical solution of Eqs. (2-3). The shear rate profiles show the development of a transient shear bands having two clearly distinct bands evolving towards a homogeneous flow. Comparing the velocity profiles, Fig 3(c), with the carpobol gel experiments Divoux et al. (2010) shows close similarities, even if the model here exhibits no yield stress.
The origin of transient shear banding is a self-feeding mechanism, which can be understood considering the Eqs. (1-3). When the fluid starts to evolve at a homogeneous high volume fraction, the decreases fastest at the regions, where has the largest value (Eq. (1)): close to the inner cylinder of the device. The accelerated decrease of (Eq. (2)) due to the faster relaxation rate, further increases the (Eq. 3) at the same location. This accelerates the decrease of the (Eq. (1)) at the same position creating a self-amplifying mechanism for the growth of transient shear banding. Since the steady state viscosity is constant, the will decrease elsewhere in the device, only with slower rate due to lower . Finally a homogeneous steady state profile is reached. If the feedback either from to (Eqs. (2) and (3)), or from to (Eq. (1)) is not strong enough, the transient shear banding does not appear. The intensity of the feedback is adjusted by the exponent and the derivative of the Eq. (2) at the corresponding . If at startup is below the steady state one, such a feedback loop does not exist; the growth of the is fastest at the high shear rates, thus promoting the increase of the and the decrease of the at high shear rate regions. This reduces the shear rate differences in the gap.
It has been observed experimentally, that in simple yield stress fluids the fluidization time decays with the global shear rate following a power-law. Similar behaviour is observed here (Fig. 4).
Here, the power-law decay follows , with fixed by/to the exponent in Eq. (1). The other parameters or even the form of the viscosity function (Eq. (2)) do not influence the power-law slope. The rest of the model parameters and the gap width simply change the vertical position of the resulting power-laws.
where represents the geometry effects, and in particular is crucial for creating the transient shear bands. This form reveals a natural way to rescale the time: . The measures dependent on the angular velocity can be rescaled with the global shear rate.
Other ways to illustrate transient shear banding, besides the velocity profiles, are to plot the temporal evolution of the local shear rate at the gap edges, and the band width Divoux et al. (2010). These quantities are plotted with the same rescaling of time and shear rate in Fig. 5 (a) and (b), respectively. The band edge is estimated as the position, where . The Fig. 5 shows that there is a short induction period at small times, during which the shear rate localizes near the rotor. During the relaxation period, the shear rate decreases towards the steady state value as the transient shear bands vanish.
Divoux et al. Divoux et al. (2011) demonstrated that the fluidization time exponents in stress and shear rate controlled experiments relate to the steady state power-law exponent. Here, this equals unity. As one would expect, the fluidization for the stress controlled case follows a power-law as does the shear rate controlled one . A structural kinetics model producing Herschel-Bulkley flow curve with a control of the exponent can be speculated to change also the the fluidization exponents towards the experimental ones Divoux et al. (2011). To the best of our knowledge, such a model does not exist. Studying the properties of more advanced structural kinetics equations on the fluidization is therefore left to future studies.
Conclusions. — We have studied transient shear banding, or fluidization of a simple time dependent fluid. We constructed a minimal model which lacks further complications such as the (visco)elasticity of the structure, a physically-motivated yield stress, the normal stresses suggested very recently Cheddadi et al. (2012) to play a role in steady state shear banding Møller et al. (2006), or elastic stress overshoots Adams et al. (2011). Our analysis indicates that such conditions are not required for transient shear banding. On the contrary, transient shear banding should be a general feature of complex fluids like colloidal suspensions and microgels.
The transient shear bands during the fluidization was found to originate from the initial shear inhomogeneity here arising from the Couette geometry, which is amplified by the positive feedback coming from the shear dependent relaxation rates of the fluid. This mechanism should be present in all similar scenarios of fluidization: studying the couplings built into the model implies, that the same effect should occur in all practical measuring geometries for shear thinning complex fluids. This is since even the smallest stress gradient is enough to trigger the transient shear bands due to the shear dependent response of the fluid. Experimental observations of stress signatures associated to transient shear banding found in a cone-and-plate geometry, where the stress gradient is extremely small, also support this finding Divoux et al. (2012). Furthermore, in Eq. 3 we neglected the inertial terms in the Navier-Stokes equation, which could play a role in the start-up flow in such situations (referring to cone-and-plate). The self-amplifying mechanism between the flow and the fluid structure would be the same, but the origin of the shear inhomogeneity would come from the inertial terms rather than the flow geometry, which was the case in the circular Couette studied here.
The kinetic exponent related to the dependence of on , is connected to the fluidization exponent, in both the stress and shear controlled cases. This indicates that studying the fluidization experimentally gives detailed information of the timescales of the internal relaxation processes in time-dependent fluids, which could be utilized to build proper structural models based on the experimental fluidization data. Such models could be established for instance around rheological models describing the volume fraction using Population Balances, which take into account the particle size distribution and concentration as rheological parameters Heath et al. (2006). Presently, work is devoded along these lines of research.
We are grateful to S. Manneville and T. Divoux for fruitful discussions. This work has been supported by the Effnet program in the Finnish Forest Cluster Ltd, and EU Framework 7 program SUNPAP. Also, the support from the Academy of Finland through the COMP center of excellence and the project number 140268 and within the framework of the International Doctoral Programme in Bioproducts Technology (PaPSaT) are acknowledged.
- D. A. Head, A. Ajdari, and M. E. Cates, Europhys. Lett 57, 120 (2002).
- S. M. Fielding and P. D. Olmsted, Phys. Rev. Lett. 92, 084502 (2004).
- D. B. Braun and M. R. Rosen, Rheology Modifiers Handbook - Practical Use and Application (William Andrew Publishing, 1999), first edition ed., ISBN 0815514417.
- D. Bonn, S. Rodts, M. Groenink, S. Rafa, N. Shahidzadeh-Bonn, and P. Coussot, Annu. Rev. Fluid Mech. 40, 209 (2008).
- S. Manneville, Rheologica Acta 47, 301 (2008).
- G. Ovarlez, S. Rodts, X. Chateau, and P. Coussot, Rheol. Acta 48, 831 (2009).
- P. Schall and M. van Hecke, Annu. Rev. Fluid Mech. 42, 67 (2010).
- G. Porte, J.-F. Berret, and J. L. Harden, J. Phys. II 7, 459 (1997).
- R. Besseling, L. Isa, P. Ballesta, G. Petekidis, M. E. Cates, and W. C. K. Poon, Phys. Rev. Lett. 105, 268301 (2010).
- A. Ragouilliaux, G., N. Shahidzadeh-Bonn, B. Herzhaft, T. Palermo, and P. Coussot, Phys. Rev. E 76, 051408 (2007).
- A. Fall, J. Paredes, and D. Bonn, Phys. Rev. Lett. 105, 225502 (2010).
- T. Divoux, D. Tamarii, C. Barentin, and S. Manneville, Phys. Rev. Lett 104, 208301 (2010).
- J. M. Adams, S. M. Fielding, and P. D. Olmsted, J. Rheol. 55, 1007 (2011).
- S. M. Fielding and P. D. Olmsted, Phys. Rev. E 68, 036313 (2003).
- K. A. Hayes, M. R. Buckley, H. Qi, I. Cohen, and L. A. Archer, Macromol. 43, 4412 (2010).
- G. Marrucci and N. Grizzuti, J. Rheol. 27, 433 (1983).
- M. Manning, E. G. Daub, J. S. Langer, and J. M. Carlson, Phys. Rev. E 79, 016110 (2009).
- R. L. Moorcroft, M. E. Cates, and S. M. Fielding, Phys. Rev. Lett. 106, 055502 (2011).
- E. A. Jagla, J. Stat. Mech.: Theory Exp. 2010, P06006 (2010).
- R. L. Moorcroft and S. M. Fielding, arXiv: 1201.6259v1 (2012).
- P. Møller, A. Fall, V. Chikkadi, D. Derks, and D. Bonn, Phil. Trans. R. Soc. A 367, 5139 (2009).
- T. Divoux, C. Barentin, and S. Manneville, Soft Matter 7, 8409 (2011), URL http://dx.doi.org/10.1039/C1SM05607G.
- J. D. Martin and T. Hu, Soft Matter 8, 6940 (2012).
- J. Mewis and N. Wagner, Adv. Colloid Interface Sci. 147-148, 214–227 (2009).
- P. C. F. Møller, J. Mewis, and D. Bonn, Soft Matter 2, 274 (2006).
- A. R. Heath, P. A. Bahri, P. D. Fawell, and J. B. Farrow, A.I.Ch.E. J. 52, 1641 (2006).
- J. M. Piau, J. Non-Newtonian Fluid Mech. 144, 1 (2007).
- F. K. Oppong, L. Rubatat, A. E. B. B. J. Frisken, and J. R. de Bruyn, Phys. Rev. E 73, 041405 (2006).
- I. Krieger and T. Dougherty, Trans. Soc. Rheol. 3, 137 (1959).
- D. Quemada, Appl. Rheol. 18, 53298 (2008).
- N. Roussel, R. L. Roy, and P. Coussot, J. Non-Newtonian Fluid Mech. 117, 85 (2004).
- R. J. Phillips, R. C. Armstrong, R. A. Brown, A. L. Graham, and J. R. Abbott, Phys. Fluids A 4, 30 (1992).
- P. Estellé, C. Lanos, and A. Perrot, Journal of Non-Newtonian Fluid Mechanics 154, 31 (2008), URL http://www.sciencedirect.com/science/article/pii/S03770257080%˙0013X.
- L. Bécu, P. Grondin, A. Colin, and S. Manneville, Colloids and Surfaces A: Physicochemical and Engineering Aspects 263, 146 (2004), ISSN 09277757, URL http://linkinghub.elsevier.com/retrieve/pii/S0927775704010064%˙.
- I. Cheddadi, P. Saramito, and F. Graner, J. Rheol. 56, 213 (2012).
- T. Divoux, D. Tamarii, C. Barentin, S. Teitel, and S. Manneville, Soft Matter 8, 4151 (2012).