Search

We explore theoretically the interplay between shear banding and edge fracture in complex fluids, by performing a detailed simulation study within two constitutive models: the Johnson-Segalman model and the Giesekus model. We consider separately parameter regimes in which the underlying constitutive curve is monotonic and non-monotonic, such that the bulk flow (in the absence of any edge effects) is respectively homogeneous and shear banded. Phase diagrams of the levels of edge disturbance and of bulk (or quasi-bulk) shear banding are mapped as a function of the surface tension of the fluid-air interface, the wetting angle where this interface meets the walls of the flow cell, and the imposed shear rate. In particular, we explore in more detail the basic result recently announced in Hemingway et al. (2018): that precursors to edge fracture can induce quasi-bulk shear banding. We also appraise analytical predictions that shear banding can induce edge fracture, Skorski and Olmsted (2011). Although a study of remarkable early insight, that study made some strong assumptions about the nature of the "base state", which we examine using direct numerical simulation. The basic prediction that shear banding can cause edge fracture remains valid, but with qualitatively modified phase boundaries., Comment: 14 pages, 6 figures

We perform a detailed theoretical study of the edge fracture instability, which commonly destabilises the fluid-air interface during strong shear flows of entangled polymeric fluids, leading to unreliable rheological measurements. By means of direct nonlinear simulations, we map out phase diagrams showing the degree of edge fracture in the plane of the surface tension of the fluid-air interface and the imposed shear rate, within the Giesekus and Johnson-Segalman models, for different values of the nonlinear constitutive parameters that determine the dependencies on shear rate of the shear and normal stresses. The threshold for the onset of edge fracture is shown to be relatively robust against variations in the wetting angle where the fluid-air interface meets the hard walls of the flow cell, whereas the nonlinear dynamics depend strongly on wetting angle. We perform a linear stability calculation to derive an exact analytical expression for the onset of edge fracture, expressed in terms of the shear-rate derivative of the second normal stress difference, the shear-rate derivative of the shear stress (sometimes called the tangent viscosity), the jump in shear stress across the interface between the fluid and the outside air, the surface tension of that interface, and the rheometer gap size. Full agreement between our analytical calculation and nonlinear simulations is demonstrated. We also elucidate in detail the mechanism of edge fracture, and finally suggest a new way in which it might be mitigated in experimental practice. Some of the results in this paper were first announced in an earlier letter. The present manuscript provides additional simulation results, calculational details of the linear stability analysis, and more detailed discussion of the significance and limitations of our findings., Comment: 18 pages, 12 figures; v2: updated to post-referee version

We study theoretically pressure driven planar channel flow of shear thinning viscoelastic fluids. Combining linear stability analysis and full nonlinear simulation, we study the instability of an initially one-dimensional base state to the growth of two-dimensional perturbations with wavevector in the flow direction. We do so within three widely used constitutive models: the microscopically motivated Rolie-Poly model, and the phenomenological Johnson-Segalman and White-Metzner models. In each model, we find instability when the degree of shear thinning exceeds some level characterised by the logarithmic slope of the flow curve at its shallowest point. Specifically, we find instability for $n

Despite decades of research, the question of whether solutions and melts of highly entangled polymers exhibit shear banding as their steady state response to a steadily imposed shear flow remains controversial. From a theoretical viewpoint, an important unanswered question is whether the underlying constitutive curve of shear stress $\sigma$ as a function of shear rate $\dot{\gamma}$ (for states of homogeneous shear) is monotonic, or has a region of negative slope, $d\sigma/d\dot{\gamma}<0$, which would trigger banding. Attempts to settle the question experimentally via velocimetry of the flow field inside the fluid are often confounded by an instability of the free surface where the sample meets the outside air, known as "edge fracture". Here we show by numerical simulation that in fact even only very modest edge disturbances - which are the precursor of full edge fracture but might well, in themselves, go unnoticed experimentally - can cause strong secondary flows in the form of shear bands that invade deep into the fluid bulk. Crucially, this is true even when the underlying constitutive curve is monotonically increasing, precluding true bulk shear banding in the absence of edge effects., Comment: 5 pages, 4 figures; v2: updated to post-referee version

We study theoretically the edge fracture instability in sheared complex fluids, by means of linear stability analysis and direct nonlinear simulations. We derive an exact analytical expression for the onset of edge fracture in terms of the shear-rate derivative of the fluid's second normal stress difference, the shear-rate derivative of the shear stress, the jump in shear stress across the interface between the fluid and the outside medium (usually air), the surface tension of that interface, and the rheometer gap size. We provide a full mechanistic understanding of the edge fracture instability, carefully validated against our simulations. These findings, which are robust with respect to choice of rheological constitutive model, also suggest a possible route to mitigating edge fracture, potentially allowing experimentalists to achieve and accurately measure stronger flows than hitherto., Comment: 5 pages, 4 figures (plus Supplementary Material); v2: updated to post-referee version, incl. minor update to Fig. 1

We reconcile two scaling laws that have been proposed in the literature for the slip length associated with a moving contact line in diffuse interface models, by demonstrating each to apply in a different regime of the ratio of the microscopic interfacial width $l$ and the macroscopic diffusive length $l_D= (M\eta)^{1/2}$, where $\eta$ is the fluid viscosity and $M$ the mobility governing intermolecular diffusion. For small $l_D/l$ we find a diffuse interface regime in which the slip length scales as $\xi \sim(l_Dl)^{1/2}$. For larger $l_D/l>1$ we find a sharp interface regime in which the slip length depends only on the diffusive length, $\xi \sim l_D \sim (M\eta)^{1/2}$, and therefore only on the macroscopic variables $\eta$ and $M$, independent of the microscopic interfacial width $l$. We also give evidence that modifying the microscopic interfacial terms in the model's free energy functional appears to affect the value of the slip length only the diffuse interface regime, consistent with the slip length depending only on macroscopic variables in the sharp interface regime. Finally, we demonstrate the dependence of the dynamic contact angle on the capillary number to be in excellent agreement with the theoretical prediction of \cite{Cox1986}, provided we allow the slip length to be rescaled by a dimensionless prefactor. This prefactor appears to converge to unity in the sharp interface limit, but is smaller in the diffuse interface limit. The excellent agreement of results obtained using three independent numerical methods, across several decades of the relevant dimensionless variables, demonstrates our findings to be free of numerical artifacts., Comment: 24 pages, 5 figures