# Start-up inertia as an origin for heterogeneous flow

## Abstract

For quite some time non-monotonic flow curve was thought to be a requirement for shear banded flows in complex fluids. Thus, in simple yield stress fluids shear banding was considered to be absent. Recent spatially resolved rheological experiments have found simple yield stress fluids to exhibit shear banded flow profiles. One proposed mechanism for the initiation of such transient shear banding process has been a small stress heterogeneity rising from the experimental device geometry. Here, using Computational Fluid Dynamics methods, we show that transient shear banding can be initialized even under homogeneous stress conditions by the fluid start-up inertia, and that such mechanism indeed is present in realistic experimental conditions.

## I Introduction

Soft Glassy Materials (SGMs) present a theoretical concept for fluids that can be solid-like or liquid-like at rest, but all
posses a microstructure with strongly interacting building blocks Divoux *et al.* [2016].
An extreme class of these are the yield stress fluids, which by virtue of their internal structure behave as solids under an imposed
mechanical stress below the yield point, yet flow like fluids once the stress is increased beyond this limit Møller *et al.* [2006]; Divoux *et al.* [2010]; Moller *et al.* [2009].
In addition to their extensive utilization in
commercial products as e.g. toothpastes, paints and process suspensions Coussot *et al.* [2002a],
yield stress fluids are also of profound theoretical interest due to their rich and complex rheological behavior.
Indeed, the peculiarities of flow response are also witnessed in shear banding, where the flowing fluid
exhibits a spatially banded structure, each band possessing a distinct viscosity associated with a unique shear
rate Manneville [2008].
This banding can appear as a true steady-state structure Olmsted [2008] or as transient shear bands
(TSBs), which, while possibly extremely long lived, eventually are replaced by the homogeneous steady-state Moorcroft *et al.* [2011].
State-of-the-art observations find transient shear banding of simple yield stress fluids very robustly in different experimental
geometries Divoux *et al.* [2010].
Furthermore, even though non-thixotropic, such fluids display time- and shear rate or stress dependent response to step shear,
where the relaxation time apparently follows a power-law of both applied stress and shear rate Divoux *et al.* [2011a].

Steady-state shear-banding has been described as a property arising from the non-monotonic nature
of the fluid’s intrinsic flow curve Olmsted [2008]; Fielding [2007]. This allows for a mechanically
unstable flow regime, where the negative slope in the constitutive curve permits the existence of a non-unique correspondence between
shear rates and shear stresses.
Thus, in this special regime the fluid may separate into two bands
of either identical shear stress (gradient banding) or identical shear rates (vorticity banding) Olmsted [2008].
Such a steady state scenario is unexpected for fluids known to have monotonic intrinsic flow curves Divoux *et al.* [2010].

Against this reasoning, transient shear banding has recently been experimentally observed Divoux *et al.* [2010, 2011b],
and theoretically predicted even in fluids possessing monotonic flow curves Adams *et al.* [2011].
Numerous reasons for this type of banding have been proposed.
In their studies, involving two models, the fluidity modelPicard *et al.* [2002] and spatially resolved version of the soft glassy rheology (SGR) model Fielding *et al.* [2009],
Moorcroft and co-workers Moorcroft *et al.* [2011] predicted transient shear banding resulting from a mechanical instability
in the start-up flows. This was clearly seen in the instantaneous (stress-time) constitutive curve as a stress overshoot and
subsequent negative slope of the curve as it regresses towards the steady state.
In a more recent work, the same authors derive a more rigorous expression for the onset of this mechanical instability in multiple
flow scenarios Moorcroft and Fielding [2014].

Transient shear banding in amorphous solids has been successfully modeled by Shi and Falk with
shear transformation zone (STZ) theory, in which the plastic deformations occur in specific zones, activated by an effective
temperature Shi and Falk [2005]. At present, Hinkle and Falk are extending this framework to study
transient shear banding in simple Yield Stress Fluids (YSF) Hinkle and Falk [2016]. These studies properly explain the movement and limited life-time
of the shear bands, however, they leave open detailed origin of the band initiation. To study this aspect,
a slightly different approach was proposed in Illa *et al.* [2013],
where a phenomenological, simple scalar -model is utilized. In this case, the local microstructure of a time-dependent fluid is
encoded in the parameter, which refers to the local immobilized volume fraction of the fluid. In the framework
of this model, transient shear banding can initiate due to the small stress gradients that are induced by the rotational Couette flow geometry,
often utilized in experimental studies. These stress variations
generate shear rate heterogeneities, which in turn result in the non-uniform disintegration of the underlying fluid structure,
observed as the formation of a TSB Illa *et al.* [2013]; Lehtinen *et al.* [2013]. However, in experimental geometries,
such as the cone-and-plate, where TSB is also found, the stress heterogeneity is considered to be negligible.

In this work, we consider the role of fluid inertia to initiate transient shear bands during the start-up, by solving
the Navier-Stokes equations for incompressible flow simultaneously to the evolution of a scalar field variable describing
the complex fluid structure. The purpose is to show the effect of accelerating flow to the transient shear band formation,
neglected earlier.
As reviewed above, to generate transient shear bands, some models pre-initialize by hand a non-homogeneous structure
profile Hinkle and Falk [2016]; Adams *et al.* [2011] and others rely on the geometrical stress heterogeneity Illa *et al.* [2013]; Adams *et al.* [2011].
While especially the latter mechanism compares reasonably to most experimental setups
(Couette, parallel plates), its relevance in the shear band initialization may be questioned in others (cone-and-plate).
Our results suggest a plausible mechanism initializing transient shear bands even under perfectly homogeneous
initial conditions due to the flow acceleration mechanism.

Utilizing a Computational Fluid Dynamics (CFD) approach, we apply the -model as in Ref. Illa *et al.* [2013]
in a homogeneous flow (planar Couette) scenario, where, as mentioned, any and all stress heterogeneities are due to inertial
effects alone.
We emphasize here, that, as will be explained later, the key concept in this model is the fact that the structure breakdown
is proportional to the shear rate, a fact based on advanced experiments Divoux *et al.* [2010].
Thus, the qualitative behavior shown here, is not influenced by this particular model, but is general for all the
models and real physical materials that have this property.
The paper is organized as follows: first, we detail the theoretical framework and the implementation. Then,
we proceed to comment on the relevant results obtained with our approach and finally, the paper finishes with concluding
remarks.

## Ii The model

As the aim here is to examine TSBs in a time-dependent fluid, a necessary
prerequisite for modeling transient shear-banding theoretically involves coupling a
structural model describing the internal structure of the fluid to the subsequent flow dynamics. Additionally, as the focus
is on providing a minimal example that presents qualitative evidence of transient shear-banding
driven by inertial effects, this structural model is chosen from a class of simple phenomenological models
based on the abstract scalar
structural parameter Coussot *et al.* [2002b, a]; Cheng and Evans [1965]; Cheng [2003]; Mujumdar *et al.* [2002]; Illa *et al.* [2013].
The specific -model addresses the internal structure of the fluid through the
structural parameter , describing the immobilized volume fraction present in the fluid Illa *et al.* [2013].
This model allows fine-tuning the fluidization exponent to match the experimental values of various complex fluids, while other similar scalar models (such as the -model by Coussot et al. Lehtinen *et al.* [2013]) inherently produce a fixed fluidization exponent, therefore
limiting their applicability for this purpose. Indeed, to the best of our knowledge, the -model is the simplest one possessing this crucial feature. Furthermore, the fluidization behavior of more elaborate models, such as the SGR-model Moorcroft and Fielding [2013]; Sollich *et al.* [1997], are presently not reported.
In the -model, the progression of the immobilized volume fraction, , in time can be due to shear (constructive and destructive) and
shear-independent motion of the structural elements (constructive), reflected in the model as two superimposed
kernels. The time-evolution equation for is Illa *et al.* [2013]

(1) |

where () is a kinetic constant for the shear growth (destruction), is the magnitude of the strain rate (shear rate) and both and describe the sensitivity of towards shearing. Furthermore, , and describe the growth of due to shear-independent effects. Again pursuing a minimal example, setting in Eq. (1) and defining the sample history by fixing the initial immobilized volume fraction allows for a minimal inspection of an internally structured fluid with a steady-state of , with the structure breakdown rate set by and proportional to the powers of the dimensionless shear rate . The proportionality gives the minimum volume fraction reached by a physical system. For instance in an aggregating colloidal suspension, this can be thought to describe the monomer volume fraction or in a microgel this could be the volume fraction at the maximally compressed state of the sponge like elements. We define . The parameter determines the system sensitivity to shear; a larger number means longer relaxation time.

The -parameter can be incorporated into the flow quantities, as done here, by the
empirical Krieger-Dougherty relation Illa *et al.* [2013]; Lehtinen *et al.* [2013]

(2) |

where is the viscosity of the fluid, describes the viscosity of the suspending
matrix and denotes the jamming volume fraction, which describes a completely jammed
configuration. This gives an additional constraint to the kinetic constants so that .
The quantities used in Eq. (2) are set to ,
(water) and Illa *et al.* [2013].

The final step in the model development is to include Eq. (2) in the Navier-Stokes equations that describe the flow field completely. The incompressible Navier-Stokes equations read

(3) |

and

(4) |

where Eq. (3) is incorporated in the latter equation, physically implying the conservation of momentum. In these equations, is the velocity field, is the pressure, is the density of the fluid and describe bodily forces (e.g. gravity) acting on the fluid (here, ). Additionally, is the kinematic viscosity of the fluid, readily obtained from Eq. (2) applying the relation . The inertial effects, generally omitted in rheological modeling, which usually deals with creeping flow conditions, are included in the convection and acceleration terms on the left-hand side of Eq. (4).

Transient shear-banding emerges as a property in the fluidization of an internally structured fluid experiencing spatial stress variations and subsequent shear rate inhomogeneities, leading to spatially non-uniform structural disintegration. As the viscosity is assumed a function of the structural integrity (here modeled by coupling Eqs. (2) and (1)), this is experimentally seen as viscosity bands across the measurement gap. Accordingly, to isolate any sources other than the inertial contributions for such variations, the planar Couette geometry shown in Fig. 1 is utilized in this work. The fluid, confined between the two parallel plates, is subjected to simple shear as the upper plate moves and the lower one remains stationary. In this scenario, the stress is uniform over the gap between the plates, ensuring that all stress variations are due to inertial effects alone.

This approach is implemented in a Computational Fluid Dynamics (CFD) environment.
The flexible OpenFOAM® software package Weller *et al.* [1998]; Chen *et al.* [2014] that
includes optimized numerical solvers, meshing, and visualization tools is applied and extended to account for the scenario presented here.
The simple shear geometry of Fig. 1 is meshed to a 2-D 6060 grid with periodic boundary conditions at the open
boundaries. The practical implementation extends the capabilities of nonNewtonianIcoFoam solver to allow for time-dependent viscosity.
The system of differential equations, is solved using the PCG (PCBiCG) linear solver with DIC (DILU) preconditioning and smoothing,
setting the absolute tolerance close to the machine precision of the 64-bit system ().
The kernel constants of the -model in Eq. 1 were set to , , and .
Additionally, the gap width was set to mm. For brevity, the following results refer to the normalized
gap (scaled by ) when necessary.
The qualitative results are independent of the selection of these parameters. However, varying the shifts the fluidization curves (see Fig. 4) in the horizontal direction. Changing the ratio to a smaller value improves the contrast of the shear band edge, i.e. makes the kink in the velocity profile steeper, and shifts the fluidization curves in the vertical direction having similar effect as increasing the (Fig. 4).

## Iii Results

The Navier-Stokes equation is characterized (at steady state) by the Reynolds number, which is represented here by the shear rate. In addition in the relaxation of the model depends on the shear rate exponent and the initial condition . Thus, this set of parameters define the TSB and is studied in what follows.
In Fig. 2 (, ), the linear velocity profile
is displayed in the homogeneous, planar Couette start-up flow scenario. In experiments, these 1-D pictures
are accessed by e.g. ultrasonic speckle velocimetry Divoux *et al.* [2010] and, in general, are used to identify shear
localization effects.
Observing the evolution of the velocity in the flow direction allows to monitor the shear localization here, as well.
As demonstrated in the example scenario of Fig. 2, shear banded flow profile develops in the course of a few seconds from
the on-set of the flow. Here we emphasize that no external initial disturbance to the model is required and the shear localization is induced by the stress inhomogeneity due to the initial acceleration of the fluid, which
interacts with the viscosity evolution as discussed in detail for the creeping flow case in Ref. Illa *et al.* [2013].

With these model parameters, the TSB persist for up until s where complete fluidization and the steady-state
flow profile is reached. This certainly is a time-scale comparable to the experiments, where long lasting transient
banding was observed Divoux *et al.* [2010].
Note that since the mechanism here is based on the flow model interaction with the time-dependent shear thinning of the fluid, any
model having this property will potentially experience shear banding induced by the same effect. Experimentally this behavior is
present in systems showing microstructural disintegration due to shear, where the time-dependence is a power-law (typically of exponents
in the range of 2-4) of the shear rate or shear stress Divoux *et al.* [2010].
Conversely, if the structural decomposition is linear (power-law exponent 1.0),
the inertial effects and the subsequent stress
variations are incapable of generating a TSB. This is illustrated in Fig. 3, where the model parameters are identical to those in Fig. 2
aside from the exponent of the -model ().
Here, due to the linear dependence between and , the structural breakdown is uniform and
the shear localization does not persist up to the time scales observed in Fig. 2. Instead, the fluidization time is simply established
by the viscous time scale ( denoting the gap width), as for a simple (Newtonian) liquid. Judging by this expression,
the high viscosity of the jammed initial configuration should yield an extremely low fluidization time. Indeed, as seen in Fig. 3, the initial
fluidization occurs within seconds from the moment the flow commences.

On the other hand, in the model, shear localization occurs only on certain conditions:
the formation of a TSB is largely dependent on the model parameters as in our earlier work Illa *et al.* [2013],
where both the values of sample history variable and the kinetic exponent are able to predict the transient shear banding.
However, here, the imposed shear rate also influences this process,
since it is this quantity that provides the necessary stress heterogeneity in the beginning of the flow.
This point is illustrated in Fig. 4, where the fluidization time (here, the time the TSB persists) is represented as a function
of the applied shear rate . For each ,
was determined from the linear velocity profile (as shown in Fig. 2) as follows:
two separate linear curves were fitted to the profile, one spanning the area
and the other .

Once the slopes ( and ) of these curves reached a certain threshold (here: ), the profile was determined to having reached the homogeneous flow state. The results are displayed in Fig. 4, where it is observed that a sufficient value of (and accordingly, sufficient stress heterogeneity) is required for transient shear banding, otherwise an uniform breakdown of the underlying structure is witnessed (). Once this is reached, the time required for complete fluidization decays rapidly. Moreover, if the structural decomposition is linear (), the breakdown is uniform regardless of . Note, that the discrete ramp in the fluidization times should not be interpreted literally as the transformation from strong shear localization (TSB) to homogeneous flow is of continuous nature. However, to quantify the fluidization times, we have resorted to the procedure introduced above, which essentially imposes a cut-off criterion (and the discrete jump) to the results.

Additionally, a complementing point of view on TSB formation is provided by various temporal plots. In Fig. 5, the apparent shear stress , measured at the shearing plane, is displayed as a function of time for various applied shear rates . As the flow progresses with increasing time, the microstructure of the fluid disintegrates and the local jammed volume fraction at the shearing plane decreases, leading to a monotonic decrease of the effective viscosity as well, as the two are linked by Eq. (2). This, in turn, leads to a monotonic decrease in the shear stress, as observed in the figure. Furthermore, the time at which the stresses reach their steady-state values seems to concur well with the fluidization time observed in Fig. 4 as one would expect.

Another temporal plot is illustrated in Fig. 6, where the shear band edge is tracked over the normalized gap () as a function
of time . In this plot, the value corresponds to a scenario where the flow is completely localized in the region near the shearing
plane and where the TSB vanishes. For consistency, the shear band edge was determined by the linear fitting scheme
introduced earlier (see the commentary on Fig. 4). The edge location was determined as the intersection point of the two linear curves
as long as the steady-state condition ()) was not met. In the framework of the -model, the shear band edge
travels at a varying (non-constant) velocity, dependent on the applied shear rate . In addition, the fluidization process is gradual up to a critical point, at which
the fluid suddenly yields completely. This is observed in Fig. 6 as the discrete jump of the edge location to the value .
The exact location of this jump is dependent on the applied shear . The behavior witnessed here should be compared with the corresponding earlier results Illa *et al.* [2013],
where the shear band edge also proceeds in a highly non-linear manner. However, the critical point is not observed, and the edge travels independent of
the applied shear rate for all . The reason for these differences lies in the fundamentally different source of the
stress gradient, that serves as the catalyst for transient shear banding. In Ref. Illa *et al.* [2013], this stress variation is induced by the geometry and is essentially time-independent.
Here, the stress variation is a (time-dependent) result of the finite fluid inertia, which is directly influenced
by . This leads to a weaker localization of the fluid flow at low and strong localization at high . Indeed, the blue
curve in Fig. 6 (s) bears a strong resemblance to the results in Ref. Illa *et al.* [2013], while lower values of
yield a more convoluted profile: the velocity of the shear band edge increases at first, but then decreases and finally diverges as the
critical point is reached. Thus, rather than reaching the static plate, the shear band disperses suddenly as the stress gradient has diminished over time and is unable to sustain
the banded profile.

Finally, Fig. 7 provides a semilogarithmic “phase” diagram, in which the regions of localized and uniform shear are displayed.
The boundary indicates the minimum triggering value of ,
at each required to yield significant heterogeneous structure relaxation, triggered by the start-up stress overshoot due to
the fluid inertia.
Above this value, a TSB forms, and below it the coupling between the structure relaxation and the local shear heterogeneity
during the start-up inertia is not strong enough to trigger the TSB. At this regime the
structure relaxes homogeneously, manifested by uniform shear. The boundary between these two regimes moves to lower shear rates, as the
relaxation exponent is increased.
As Ref. Illa *et al.* [2013] explains, in the framework of this model, a kinetic exponent on the order of would be required
to match the duration of the TSB to experiments. Unfortunately, using here an exponent as high as this would not be permitted
due to numerical and computational limitations (the number of required iterations for convergence).
Extrapolating the monotonically decreasing trend based on the computed values of Fig. 6 to hints that
the minimum value for should be well below the typical (minimum) experimentally used
shear rates (). Therefore, we are inclined to argue based on the model,
TSB formation and the subsequent localized flow dynamics due to start-up inertia only should be recovered trivially for all
real-world complex fluids. This has the important implication for the experiments that TSB should appear in all simple yield stress fluids with such
fluidization exponents irrespective of the measuring geometry.

## Iv Conclusions

Earlier theoretical work proposes a scenario, where the transient shear banding is initialized by gradients in
stress Illa *et al.* [2013], or inhomogeneity of the sample Adams *et al.* [2011]; Hinkle and Falk [2016]; Moorcroft and Fielding [2013]. The present study brings a new mechanism, which shows
that transient shear banding can result even due to the start-up stress gradients caused by accelerating flow. As we show here,
depending on the linearity of the fluids’ viscosity response, and its initial state, this effect can be either negligible or very significant.
When the initial viscosity of the fluid is small, and/or the time-dependence of the fluid is close to linear
the stress gradient is not able to induce any inhomogeneities during the start-up flow.
However, according to our simulations, when the fluid at rest is in a high viscosity state
and shows time-dependence under constant shear rate with a large fluidization exponent
the start-up stress inhomogeneities initialize long lasting transient shear bands at all
experimentally accessible shear rates.
This occurs despite the fact that the influence of the Navier-Stokes instability is multiple orders of magnitude faster compared to the actual TSB. However, it creates a small but non-negligible inhomogeneity in the structure, which is amplified over time into a complete TSB.
An example of such material is the carbopol, a typical example of so called simple yield stress fluids, for which
the fluidization exponent is larger than 2 Divoux *et al.* [2010] and the viscosity diverges at rest.
Our findings give a reasonable explanation for the initialization of transient shear banding even in experimental geometries
exhibiting negligible stress heterogeneities, such as the cone and plate one.
This suggests that stress heterogeneities should be anticipated in all
strongly shear thinning fluids during the start-up phase of rheological experiments.

## V Acknowledgements

This work was supported by the Academy of Finland through the COMP center of excellence and the project number 278367. The simulations were performed using the computer resources within the Aalto University School of Science \unichar8220Science-IT\unichar8221 project.

### References

- T. Divoux, A. Fardin, Marc, S. Manneville, and S. Lerouge, Annu. Rev. Fluid Mech. 48, 81 (2016).
- P. C. Møller, J. Mewis, and D. Bonn, Soft Matter 2, 274 (2006).
- T. Divoux, D. Tamarii, C. Barentin, and S. Manneville, Phys. Rev. Lett. 104, 208301 (2010).
- P. Moller, A. Fall, V. Chikkadi, D. Derks, and D. Bonn, Philos. Trans. A Math. Phys. Eng. Sci. 367, 5139 (2009).
- P. Coussot, Q. D. Nguyen, H. Huynh, and D. Bonn, Phys. Rev. Lett. 88, 175501 (2002a).
- S. Manneville, Rheologica Acta 47, 301 (2008).
- P. D. Olmsted, Rheologica Acta 47, 283 (2008).
- R. L. Moorcroft, M. E. Cates, and S. M. Fielding, Phys. Rev. Lett. 106, 055502 (2011).
- T. Divoux, C. Barentin, and S. Manneville, Soft Matter 7, 8409 (2011a).
- S. M. Fielding, Soft Matter 3, 1262 (2007).
- T. Divoux, C. Barentin, and S. Manneville, Soft Matter 7, 9335 (2011b).
- J. M. Adams, S. M. Fielding, and P. D. Olmsted, J. Rheol. 55, 1007 (2011).
- G. Picard, A. Ajdari, L. Bocquet, and F. M. C. Lequeux, Phys. Rev. E 66, 051501 (2002).
- S. M. Fielding, M. E. Cates, and P. Sollich, Soft Matter 5, 2378 (2009).
- R. L. Moorcroft and S. M. Fielding, J. Rheol. 58, 103 (2014).
- Y. Shi and M. L. Falk, Phys. Rev. Lett. 95, 095502 (2005).
- A. R. Hinkle and M. L. Falk, J. Rheol. 60, 873 (2016).
- X. Illa, A. Puisto, A. Lehtinen, M. Mohtaschemi, and M. J. Alava, Phys. Rev. E 87, 022307 (2013).
- A. Lehtinen, A. Puisto, X. Illa, M. Mohtaschemi, and M. J. Alava, Soft Matter 9, 8041 (2013).
- P. Coussot, Q. Nguyen, H. Huynh, and D. Bonn, J. of Rheol. 46, 573 (2002b).
- D. C. Cheng and F. Evans, Br. J. Appl. Phys. 16, 1599 (1965).
- D.-H. Cheng, Rheologica acta 42, 372 (2003).
- A. Mujumdar, A. N. Beris, and A. B. Metzner, J. Non-Newtonian Fluid Mech. 102, 157 (2002).
- R. L. Moorcroft and S. M. Fielding, Phys. Rev. Lett. 110, 086001 (2013).
- P. Sollich, F. Lequeux, P. Hébraud, and M. E. Cates, Phys. Rev. Lett. 78, 2020 (1997).
- H. G. Weller, G. Tabor, H. Jasak, and C. Fureby, Comput. Phys. 12, 620 (1998).
- G. Chen, Q. Xiong, P. J. Morris, E. G. Paterson, A. Sergeev, and Y.-C. Wang, Notices Amer. Math. Soc. 61 (2014).