Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Flow of wormlike micellar solutions around microfluidic cylinders with high aspect ratio and low blockage ratio

Simon J. Haward*a, Naoyuki Kitajimab, Kazumi Toda-Petersa, Tsutomu Takahashib and Amy Q. Shena
aMicro/Bio/Nanofluidics Unit, Okinawa Institute of Science and Technology, Onna, Okinawa 904-0495, Japan. E-mail:; Tel: +81 98 966 8289
bFaculty of Engineering, Nagaoka University of Technology, Nagaoka, Niigata 940-2188, Japan

Received 15th October 2018 , Accepted 12th December 2018

First published on 18th January 2019

We employ time-resolved flow velocimetry and birefringence imaging methods to study the flow of a well-characterized shear-banding wormlike micellar solution around a novel glass-fabricated microfluidic circular cylinder. In contrast with typical microfluidic cylinders, our geometry is characterized by a high aspect ratio α = H/W = 5 and a low blockage ratio β = 2r/W = 0.1, where H and W are the channel height and width, and the cylinder radius r = 20 μm. The small cylinder radius allows access up to very high Weissenberg numbers 1.9 ≤ Wi = λMU/r ≤ 3750 (where λM is the Maxwell relaxation time) while inertial effects remain entirely negligible (Reynolds number, Re < 10−4). At low Wi values, the flow remains steady and symmetric and a birefringent region (indicating micellar alignment and tensile stress) develops downstream of the cylinder. Above a critical value Wic ≈ 60 the flow transitions to a steady asymmetric state, characterized as a supercritical pitchfork bifurcation, in which the fluid takes a preferential path around one side of the cylinder. At a second critical value Wic2 ≈ 130, the flow becomes time-dependent, with a characteristic frequency f0 ≈ 1/λM. This initial transition to time dependence has characteristics of a subcritical Hopf bifurcation. Power spectra of the measured fluctuations become complex as Wi is increased further, showing a gradual slowing down of the dynamics and emergence of harmonics. A final transition at very high Wic3 corresponds to the re-emergence of a single peak in the power spectrum but at much higher frequency. We discuss this in terms of possible flow-induced breakage of micelles into shorter species with a faster relaxation time.

1 Introduction

In an aqueous solution above a critical micelle concentration, amphiphilic surfactant molecules spontaneously self-assemble into large aggregates known as micelles, organized so as to shield their hydrophobic tails from water.1 A wide variety of micelle morphologies are possible, depending on environmental conditions such as temperature, pH, concentration, salinity, and also on factors such as the surfactant packing parameter and flow conditions. For solutions of ionic surfactants such as cetyl pyridinium chloride (CPyCl), the combination with a strongly binding counterion (or co-surfactant) such as sodium salicylate (NaSal), leads to a transition from spherical micelles at low concentration through rodlike and eventually to giant wormlike micelles as the surfactant concentration is increased.2–4 Wormlike micelles (WLMs) are elongated, semi-flexible cylindrical aggregates, which similarly to polymer chains, impart elastic properties to their solutions, rendering the bulk fluid “viscoelastic”.5–7 However, unlike polymers, WLMs exist in a dynamic equilibrium, undergoing constant scission and recombination. This provides stress relaxation mechanisms unavailable to polymer chains, allows the possibility of morphological changes (e.g. micellar concatenation, or gel formation) to be caused by flow (referred to in general as flow-induced structures), and further means that bulk viscoelastic properties can recover subsequent to micellar degradation that can occur in strong flows. Their particular and highly tunable properties make wormlike micellar solutions extremely useful in a wide variety of industrial applications, as drag-reduction agents, in oilfield applications and also as rheology modifying additives in many household products.2–4,8–15 Given their widespread applications, gaining a thorough understanding of the behaviour of this imprortant class of fluids under different regimes of flow is of considerable importance.

Flow around a cylinder is considered as a benchmark system for studying the dynamics of viscoelastic fluids.16–25 The flow field around a cylinder confined between two walls provides an interesting mix of shearing and extensional kinematics. The leading axial stagnation point upstream of the cylinder causes a compressional flow and divergence of streamlines around the sides of the cylinder. There are regions of transient elongational flow as fluid is accelerated into the gaps between the cylinder and the side walls. In these gaps the shear rate can be high relative to the bulk. Downstream of the cylinder, in the locality of the trailing axial stagnation point, fluid is subjected to high extensional rates and residence times.26–28 Flow around a cylinder can be considered as a two-dimensional (2D) representation of flow around a sphere, and by extension to flow around a canonical “object”, which is clearly relevant to understanding many situations (e.g. flow separation, particle suspension or sedimentation, and flow through porous beds). However, flow around a cylinder is far more convenient to study experimentally than a sphere since (1) it can be held stationary in a channel and studied under steady-state flow conditions, (2) the flow rate is more readily controllable, and (3) the planar nature of the flow is easier to visualize. 2D numerical simulation of flow around a cylinder is also much simpler and less computationally expensive than simulating the three-dimensional (3D) flow around a sphere. As such, the flow around a cylinder provides a convenient approach to test numerical solutions to constitutive equations against experimental results, with a variety of flow kinematics incorporated in a single geometry and with significant real-world relevance.

Geometrically, a cylinder of radius r confined inside a rectangular channel of height H and width W can be characterized by two ratios of lengthscales: (1) the aspect ratio of the channel α = H/W, and (2) the blockage ratio β = 2r/W. The aspect ratio of the cylinder itself is given by AR = H/2r = α/β, see schematic diagram presented in Fig. 1.

image file: c8sm02099j-f1.tif
Fig. 1 Graphical depiction of the high-aspect ratio microfluidic cylinder device: (a) schematic drawing of the flow geometry showing the important characteristic dimensions and the coordinate system. Q represents the volumetric flow rate. (b) and (c) show micrographs of the flow geometry taken from top and side-on views, respectively. (d) Complete device in fused silica, assembled with stainless steel inlet and outlet connectors.

There are several previous experimental studies of wormlike micellar solutions flowing around macro-scale cylinders.29–31 Of particular interest, Moss and Rothstein30,31 used a combination of flow velocimetry, pressure drop measurements and flow-induced birefringence visualization to study the flow of two different wormlike micellar solutions around cylinders confined inside channels with H = W = 50 mm (α = 1) and with radii r = 2.5 and r = 5 mm (β = 0.1 and β = 0.2). The two micellar fluids were composed of 50 mM cetyltrimethyl ammonium bromide (CTAB) mixed with 50 mM NaSal, and 100 mM CPyCl mixed with 50 mM NaSal. For both fluids at low flow rates the flow around the cylinders was symmetric and steady. As the flow rate was increased, the appearance of birefringence (indicating micellar alignment and stress) around the cylinders and in an extended downstream wake was observed. For the CTAB/NaSal solution, the flow was found to become unstable above flow rates corresponding to a critical Weissenberg number Wic ≈ 4. Here, the Weissenberg number, which describes the relative importance of elastic to viscous forces in the flow, is defined by Wi = λMU/r, where U is the average velocity of the fluid in the channel, U/r describes a characteristic deformation rate around the cylinder and λM is the Maxwellian relaxation time of the fluid determined from its rheological response in linear oscillatory shear. The instability was manifested by asymmetry and motion of the stress birefringence in the downstream wake of the cylinder. Since inertial effects as characterized by the Reynolds number Re = ρUr/η (where ρ and η are the fluid density and viscosity, respectively) were low, and the flow destabilized above a critical Wi, the instability was considered as being elastic in origin. A growth of the stress followed by a sudden decay, considered to be due to micelle rupture, was reported, although the occurrence of rupture events was described as chaotic since no characteristic frequency for them was found. Such stress growth and micelle rupture in the downstream wake of the cylinder explains the unsteady settling of particles in wormlike micellar fluids,32–34 and confirms the importance of the strong extensional flow and high fluid strain in the wake for driving such phenomena. Dey et al. have recently shown how elastic instabilities caused by flow around a deformable elastic cylinder can feed back onto the cylinder itself, driving erratic motions of the structure, somewhat analogous to the vortex-induced-vibrations in high Re Newtonian flows around obstacles.35,36 We note that for Weissenberg numbers up to ≈10 (the maximum that could be achieved due to the scale of their experimental set up), Moss and Rothstein reported no elastic instability in the CPyCl/NaSal solution that they studied and the flow remained always steady and symmetric.30,31

Due to the small characteristic size scale (e.g. the cylinder radius), elastic forces can be considerably enhanced in micro-scale viscoelastic flows, in which extremely high Weissenberg numbers can be achieved while keeping the Reynolds number, and hence inertial forces, negligible. Microfluidic flow geometries thus provide ideal platforms for studying elastic instabilities and the non-linear dynamics arising from the feedback between flow-induced microstructural rearrangements and the flow field itself. Elastic instabilities are frequently manifested as asymmetries, where the fluid selects a preferential path through a geometrical configuration. A well studied example (with both viscoelastic polymer and WLM solutions) is the onset of asymmetric flow through the cross-slot geometry above a critical flow rate.37–44 Using the two-species Vasquez–Cook–McKinley (VCM) model for WLMs,45 that incorporates the possibility of chain halving and recombination into the constitutive equation, Kalb et al. have shown that micelle breakage rate has a significant impact on this phenomenon.43,44 WLMs that break more readily are less prone to instability.40,44

Flow of WLMs around micro-scale cylinders may be a relevant model for understanding, e.g. transport of proppant particles used in fracturing fluids and for flow around elements of microporous structures (frequently modelled by microfluidic post arrays) and in which elastic instabilities have been shown to play an important role.46–49 However, a drawback of most microfluidic cylinder experiments to date is that, due to typically employed fabrication techniques (such as soft lithography or micromilling), their geometries have been restricted to low aspect ratios (α < 1) and to high blockage ratios (typically β ≥ 0.5).50–55 The low aspect ratio means the flow field is strongly inhomogeneous along the cylinder axial direction. Perhaps more significantly, the high blockage ratio results in the formation of narrow gaps between the sides of the cylinder and the channel walls. In this case, viscoelastic fluids respond most strongly to the squeezing flow through the gaps, resulting in the generation of elastic instabilities and recirculations upstream of the cylinder.51,52,54 The similarity of such upstream recirculations with those seen for viscoelastic flows into simple abrupt or hyperbolic contractions,56–58 demonstrates that cylinder geometries with high values of β are really only examining the effects of the contractions and not of the cylinder per se. The effects of the stronger extensional flows in the up- and downstream wake regions, are significantly suppressed as a result. It is well-recognized that to clearly observe the extensional flow in the cylinder wake regions, the blockage ratio should be much lower than β = 0.5,59 but this is difficult to achieve at the microfluidic scale using mainstream microfabrication methods.

In a recent work, we used selective laser-induced etching (SLE)60 to fabricate glass microfluidic cylinder devices (r = 20 μm) in a new geometrical regime featuring both a high aspect ratio α = 5 and a low blockage ratio β = 0.1.61 We used the devices to study the steady viscoelastic flow of dilute polymer solutions, clearly demonstrating the enhancement of elastic effects visible in the cylinder wakes.61 In the present work, we employ one of the same devices to examine the flow of a well-studied wormlike micellar system composed of CPyCl and NaSal. The flow field around the cylinder is characterized in detail using time-resolved micro-particle image velocimetry (μ-PIV) and we also make time-resolved flow birefringence measurements using a state-of-the-art high-speed polarization camera. As the Weissenberg number is increased for negligible Reynolds number ≪ 1, we observe a rich progression of flow instabilities manifested in the downstream wake, beginning with a steady flow asymmetry that develops spatio-temporal fluctuations at higher Wi and that are evident in both the flow field and the birefringence imaging. In contrast to the work of Moss and Rothstein using a CTAB/NaSal solution,30,31 spectral analysis reveals specific characteristic frequencies for the fluctuations, that vary with the Wi in a complex way. We discuss these effects in terms of micellar structural transitions caused by the flow.

2 Experimental methods

2.1 Microchannel design and fabrication

The microfluidic cylinder device used in this study is represented schematically in Fig. 1a and photographically in Fig. 1b–d. The device was fabricated by SLE in fused silica glass using a “LightFab” 3D printer (LightFab GmbH, Germany).62 The SLE fabrication is a two-step process: firstly, the shape of the channel is written within an initially solid block of fused silica by irradiation with a femtosecond laser, and secondly, the irradiated material is preferentially removed by ultrasonication at 80 °C in a wet chemical etchant (KOH). The volume occupied by the cylinder is simply not irradiated during the laser writing step and is left behind within the channel following the chemical etching step, resulting in a monolithic glass construction. The device fabrication is completed by simply bonding stainless steel tubing connectors to the inlet and outlet to the glass channel using 2-part epoxy, see Fig. 1(d). The channel has inner dimensions of W = 400 μm and H = 2 mm. This provides a relatively high aspect ratio α = 5, and hence a good approximation to a two dimensional flow, which was confirmed in our earlier study.61 The total channel length is 25 mm and the cylinder is located half way along it. The cylinder radius is r = 20 μm, which provides a low blockage ratio β = 0.1 and a cylinder aspect ratio AR = 50. Due to the high modulus of fused silica (≈75 GPa), the cylinder can be made so long and slender without deforming under an imposed flow. Even at the high magnifications used in the visualization experiments to be described in the following sections, and at the highest imposed flow rates, there was no observable deflection of the cylinder. As can be seen from Fig. 1b–d, the completely transparent glass construction provides useful optical access to both the xy and xz planes, which allows a detailed characterization of the flow.

2.2 Materials

The wormlike micellar test solution used in this study is composed of 100 mM CPyCl and 60 mM NaSal (both supplied by Sigma Aldrich) dissolved in deionised water. This is a well-studied surfactant/counterion system known to form giant wormlike micelles at the selected concentrations.3,6 The solution was prepared by adding weighed quantities of the CPyCl and NaSal powders to the appropriate volume of deionised water and stirring the mixture vigorously for three days. The fluid was then allowed to equilibrate at room temperature in dry, unlit conditions for ten days prior to experimentation.

The steady shear rheology of the test solution was measured at 24 °C (the same temperature as subsequent experiments around the microfluidic cylinder) using a DHR3 stress-controlled rotational rheometer (TA Instruments Inc.) fitted with a 40 mm diameter 1° stainless steel cone-and-plate fixture. In order to access higher shear rates, additional measurements were carried out using an m-VROC microfluidic slit rheometer (Rheosense Inc.) fitted with a 3 mm wide, 50 μm deep chip with flush-mounted 40 kPa full-scale pressure sensors.63 The results of these measurements are presented in Fig. 2a in the form of both viscosity η and shear stress σ as a function of the applied shear rate [small gamma, Greek, dot above]. The flow curve shows the typical features expected for this WLM fluid formulation.41,64,65 A constant viscosity (η ≈ 50 Pa s) is observed at low shear rates ([small gamma, Greek, dot above] ≲ 0.2 s−1) followed by a rapid transition to a stress-plateau region where the viscosity thins with a power-law index n ≈ 0. The stress-plateau, indicative of shear-banding,66,67 spans around 3 decades in shear rate, before the stress begins to increase with shear rate according to σ[small gamma, Greek, dot above]0.75 for [small gamma, Greek, dot above] ≳ 2000 s−1. The steady shear rheology is well described by the Carreau–Yasuda generalized Newtonian fluid (GNF) model,68 as shown by the solid lines in Fig. 2a:

image file: c8sm02099j-t1.tif(1)
where η is the infinite-shear-rate viscosity, η0 is the zero-shear-rate viscosity, [small gamma, Greek, dot above]* is the characteristic shear rate for the onset of shear-thinning, n is the power-law exponent in the shear-thinning region and a is a dimensionless fitting parameter that controls the rate of the transition between the constant viscosity and the shear-thinning regions. Values for all of these parameters are provided in Table 1. This GNF model accurately describes the shear-thinning behaviour, but makes no account for fluid viscoelasticity, therefore its applicability is restricted. Nevertheless, this simple model can be used to predict fully developed velocity profiles, and to provide values for the characteristic viscosity of the fluid at arbitrary shear rates applied in the microfluidic channel.

image file: c8sm02099j-f2.tif
Fig. 2 Rheological characterization of the 100/60 mM CPyCl/NaSal wormlike micellar solution at 24 °C. (a) Flow curve obtained in steady shear using cone-and-plate rotational rheometry and an m-VROC microfluidic slit rheometer. (b) Storage and loss moduli obtained by small amplitude oscillatory shear in the cone-and-plate geometry. Insert shows the Cole–Cole plot with fitted semi-circle of diameter 0.92.
Table 1 Summary of the rheological properties of the 100/60 mM CPyCl/NaSal WLM test solution at 24 °C
η0 [Pa s] η [mPa s] [small gamma, Greek, dot above]* [s−1] n a λM [s] G0 [Pa]
50 5 0.25 0 4 1.8 28.5

The storage and loss moduli, G′ and G′′ of the WLM test solution were determined by small angle oscillatory shear (SAOS) measurements made at 24 °C on the DHR3 rheometer with the 40 mm diameter 1° cone-and-plate oscillating over a range of angular frequency 0.01 ≤ ω ≤ 100 rad s−1 and with a 5% strain. Strain amplitude sweeps were carried out in order to ensure that a 5% strain was within the linear region over the full range of ω. The results of the SAOS measurements are presented in Fig. 2b, and are well-described over most of the frequency range using a simple single mode Maxwell model:

image file: c8sm02099j-t2.tif(2)
providing values for the Maxwell relaxation time λM = 1.8 s and plateau modulus G0 = 28.5 Pa, again in reasonable agreement with previous measurements.41,64,65 Following the methods described by Cates and coworkers,69–71 the characteristic micellar breakage time can be estimated as λbreak ≈ 0.17 s and the reptation time scale as λrepλM2/λbreak ≈ 19 s.

The rheological response of viscoelastic fluids to extensional flows can be characterized from measurements of the necking rate of liquid bridges in the so-called “elasto-capillary” thinning regime.72,73 Such measurements are commonly performed using the commercially available capillary breakup extensional rheometer (CaBER, Thermo Haake), in which fluid is loaded between two circular end plates which are subsequently separated by a distance sufficient to initiate capillary-driven thinning. However, it has been shown that with CPyCl/NaSal WLM solutions the thinning dynamics are sensitive to the imposed experimental conditions (i.e. plate diameter, initial and final plate separation and also the rate of plate separation).74 Consequently, CaBER results for WLM solutions may not be entirely reliable, although experiments performed under typically employed conditions provide extensional relaxation times close to those obtained from linear oscillations in shear and also evidence a significant degree of strain hardening, with the extensional viscosity reaching ηE ≈ 10 to 100 × η0 at high strains.40,41,74

2.3 Flow control and dimensionless groups

The flow in the rectangular microchannel is driven at controlled volumetric flow rate Q using two Nemesys low-pressure syringe pumps (Cetoni GmbH), one infusing at the inlet and the second withdrawing at an equal and opposite rate from the outlet. The average flow velocity inside the channel is then U = Q/(WH). The dimensionless Reynolds number decribes the relative importance of inertial to viscous forces in a flow. For flow around the cylinder we define the Reynolds number as:
image file: c8sm02099j-t3.tif(3)
Here, since the cylinder is located half way between the channel side walls where there is negligible shear, rather than using a rate-dependent viscosity we simply take the zero-shear-rate viscosity η0 = 50 Pa s as being representative. The fluid density ρ is assumed equal to that of the aqueous solvent. For the maximum flow rates imposed in our experiments with the WLM solution, Re < 10−4.

A nominal deformation rate for the flow around a cylinder can be defined [small epsi, Greek, dot above] = U/r. On this basis, we define the dimensionless Weissenberg number for the flow, which describes the relative importance of elastic to viscous forces, as:

image file: c8sm02099j-t4.tif(4)
Our experiments span a range of 1.9 ≲ Wi ≲ 3750.

An elasticity number can be defined as El = Wi/Re, thus describing the relative importance of elasticity and inertia. The elasticity number is independent of the imposed flow conditions, depending only on the fluid properties and the geometry. Due to the high fluid viscosity and long relaxation time of the WLM solution and the microscopic cylinder radius, the elasticity number of our experiment is very high (El ≈ 2 × 108), indicating that elastic forces completely dominate over inertia.

2.4 Flow visualizations

2.4.1 Micro-particle image velocimetry (μPIV). Quantitative spatially-resolved flow velocimetry in the microfluidic cylinder geometry is performed using a volume illumination micro-particle image velocimetry (μ-PIV) system (TSI Inc., MN). The WLM solution is seeded with a low concentration (≈0.02 wt%) of fluorescent microparticles (1 μm diameter Fluoromax red, Thermo Scientific Inc.) with excitation/emission wavelengths of 542/612 nm. The flow geometry is placed on the stage of an inverted microscope (Nikon Eclipse Ti) and the plane of interest (either z = 0 or y = 0 plane) is brought into focus with a 5×, NA = 0.15 numerical aperture Nikon PlanFluor objective lens. The measurement width over which microparticles contribute to the determination of the velocity field is δm ≈ 110 μm,75 which is ≈0.05H and ≈0.25W. A dual-pulsed Nd:YLF laser with a wavelength of 527 nm and a time separation between pulses Δt illuminates the fluid, exciting fluorescence of the microparticles. Images of the fluorescing particles are captured in pairs in synchronicity with the pairs of laser pulses using a high speed 1280 × 800 pixel CMOS imaging sensor (Phantom MIRO) operating in frame-straddling mode. At each flow rate examined, the time Δt is set so that the maximum displacement of particles between the two images in each pair is around 8 pixels. For steady flows at lower Wi, 50 image pairs are captured and processed using an ensemble average PIV cross-correlation algorithm (TSI Insight 4G) and Nyquist criterion. For unsteady flows at moderate Wi, 400 images pairs are captured at a rate of 20 Hz, and processed both individually and by ensemble averaging. At the highest imposed Wi, 600 image pairs are captured at a rate of 60 Hz, and processed both individually and by ensemble averaging. The processing yields 2D velocity vectors on a square grid with spacing 25.6 × 25.6 μm. For data aquired in the z = 0 plane we obtain velocity components u and v in the x and y directions respectively, while data aquired in the y = 0 plane yields components u and w in the x and z directions, respectively. Subsequent image analysis, generation of contour plots and streamline traces is performed using the software Tecplot Focus (Tecplot Inc., WA).
2.4.2 Flow-induced birefringence (FIB). Optical retardation due to micellar alignment in the flowing test solution (normally referred to as flow induced birefringence, or FIB) is measured using a state-of-the-art high speed polarization camera (CRYSTA PI-1P, Photron Ltd, Japan), that combines a “micro-polarizer array” in front of a CMOS imaging sensor. The micro-polarizer array consists of 1024 × 1024 linear polarizing elements, each the size of one pixel of the 1024 × 1024 CMOS sensor. The polarizing elements are arranged in sets of 2 × 2 with orientations of 0°, 45°, 90°, and 135°. The microfluidic device is placed on the imaging stage of an inverted microscope (Nikon Eclipse Ti) and either the z = 0 or y = 0 plane is brought into focus with a 5× objective lens. White light is passed through a band-pass filter and a circular polarizer, before passing through the sample and being focussed through the micropolarizer array and onto the imaging sensor of the camera. By analysing the light intensity received at the individual pixels in each 2 × 2 grouping, a spatially-resolved and quantitative measurement of the retardation, R and orientation angle θ is obtained. The retardation is an extrinsic quantity that describes the total phase shift occurring as light polarized at θ and θ + 90° passes through the sample. The phase shift is due to the difference in the refractive indices (n1 and n2) along the two directions (i.e. due to the birefringence Δn = n1n2). The birefringence itself is an intrinsic quantity related to the retardation by Δn = R/[small script l], where [small script l] is the pathlength through the birefringent material.76,77 In this paper we report the measured values of the retardation and we do not make a conversion to birefringence by assuming a value for [small script l]. Images are captured at 125 Hz and have a spatial resolution of ≈3.6 μm per pixel.

3 Results and discussion

3.1 Preliminary control experiments

Before examining the flow of the WLM solution around the microfluidic cylinder, control experiments were made in order to check the nature of the Newtonian flow field around the cylinder. Fig. 3a shows the normalized Newtonian velocity magnitude field (|V|/|V|max, where image file: c8sm02099j-t5.tif) around the cylinder as measured in the z = 0 plane at a low Reynolds number Re = 0.007. The flow is clearly symmetric about y = 0 and about x = 0, as expected for Newtonian flow at low Re. Fig. 3b shows profiles of the normalized streamwise velocity (u/U) taken along the line y = 0 for a range of Re. The experimental data is compared with the prediction of a 3D finite volume numerical simulation performed using ANSYS Fluent. For this we simulated the flow in a section of channel 6.5 mm in length with the cylinder located 4.5 mm downstream of the inlet. A uniform velocity profile was applied at the inlet, a constant pressure boundary condition was applied at the outlet and a no-slip condition was applied at all solid boundaries in the flow. The simulation volume was divided into 5[thin space (1/6-em)]286[thin space (1/6-em)]526 cells with a minimum cell size of 0.0012 mm and the Navier–Stokes equations were solved using the semi-implicit method for pressure-linked equations (SIMPLE) algorithm. We ensured that flow was fully developed upstream of the cylinder and that flow profiles were independent of further mesh refinement. The good agreement between the numerical simulation and the experimental results provides confirmation that the microchannel is of good quality and that our experimental measurement protocols are sound.
image file: c8sm02099j-f3.tif
Fig. 3 Flow velocimetry in the z = 0 plane with Newtonian fluid around the microfluidic cylinder. (a) Normalized velocity magnitude |V|/|V|max, with superimposed streamlines, for Newtonian flow at Re = 0.0007. (b) Normalized velocity profiles of u/U taken along the line y = 0 at various values of Re. The experimental data is compared with a 3D finite volume numerical simulation at low Re performed using ANSYS Fluent (dashed black curve).

In addition, we measured the flow field in the channel upstream of the cylinder using both Newtonian and WLM fluid, see Fig. 4. At the low Re of our experiments, the Newtonian flow profiles through both the channel width and depth are independent of the flow rate and match the infinite series analytical solution given by Shah and London,61,78 which is shown by the dashed black curves. For Newtonian fluid, the fully-developed flow profile is parabolic across the channel width (Fig. 4a), but is more uniform through the channel depth (Fig. 4b), as expected for a channel of this high aspect ratio, α = 5.

image file: c8sm02099j-f4.tif
Fig. 4 Profiles of the normalized streamwise flow velocity u/U upstream of the cylinder (a) across the channel width and (b) through the channel depth. Experimental data for the WLM solution at various values of the nominal wall shear rate [small gamma, Greek, dot above]W is averaged over a 1 mm section of channel (−5 < x < −4 mm). Solid lines are numerical predictions computed assuming a no-slip boundary condition and rheology described by the Carreau–Yasuda GNF model. Dashed lines indicate the analytical solution for fully-developed Newtonian flow.

For the shear-banding WLM solution, the measured velocity profiles depend strongly upon the imposed flow rate. In Fig. 4 we present examples for a few different values of the nominal wall shear rate [small gamma, Greek, dot above]W = 6U/W, spanning the range of flow rates used in the experiments around the cylinder. At the lowest wall shear rate [small gamma, Greek, dot above]W = 0.31 s−1 (black squares), the flow profile across the channel width is significantly flattened compared with the Newtonian case, while through the depth the profile remains Newtonian-like in shape. As the wall shear rate is increased there is a general trend towards more pluglike flow through both the width and the depth directions. This is expected due to the localization of shear (or formation of shear bands) at the channel walls.41,65,79–81 We note the appearance of possible wall slip in some cases, however we cannot confirm this confidently from these results since the spacing of vectors precludes resolution arbitrarily close to the wall. Further, the solid lines (coloured according to the corresponding experimental data points) show the results of numerical predictions of the velocity profiles. These were computed using a time-marching scheme to numerically solve the transient Stokes equation, taking the Carreau–Yasuda GNF model as the constitutive equation for the shear-rate-dependent viscosity and assuming a no-slip boundary condition at the walls. In general the shape of the flow profile is predicted quite well by this simple model. We do not observe the occurrence of localized high velocity “jets” of fluid developing through the channel depth that have been reported by some authors for flows of WLM solutions in high aspect ratio rectangular channels.82,83

3.2 Wormlike micellar flow around the cylinder

For the flow of the shear-banding wormlike micellar solution around the cylinder, we observe a complex sequence of flow transitions as the Weissenberg number is increased. Examples of time averaged flow velocimetry and FIB imaging results, broadly illustrating the nature of the flow over the range of imposed Wi are presented in Fig. 5. At low Wi = 3.75 (Fig. 5a), the flow is steady and symmetric about y = 0. The flow field contrasts with that of a Newtonian fluid (shown in Fig. 3a) due to the shear localization at the channel walls and the more pluglike flow profile. A weak birefringent signal is detected in a wide region around the cylinder, with a relatively intense signal measured in the immediate downstream wake. As the Weissenberg number is incremented, the flow remains steady and symmetric, and the FIB signal becomes more intense with an extended wake developing downstream of the cylinder (e.g. Fig. 5b, Wi = 37.5). The region of high birefringence indicates localized micellar alignment and tensile stress and has a qualitatively similar form to the birefringence reported by Moss and Rothstein for the steady symmetric flow of a similar WLM fluid around a macroscale cylinder.31 As the Weissenberg number is incremented further (e.g. Fig. 5c, Wi = 93.8), a clear asymmetry becomes apparent in the flow field, with the fluid taking a preferential path around one side of the cylinder. In the example shown, the flow is faster for y > 0 than for y < 0, but it is important to note that this asymmetry has also been observed to develop in the opposite sense, with the direction being selected seemingly at random. As the Weissenberg number is increased to Wi = 150 and beyond (see Fig. 5d and e), the degree of flow asymmetry increases, with virtually all of the fluid passing the cylinder for y > 0 and the development of an effectively stagnant region on the opposite side of the cylinder. The asymmetry is also evident in the birefringence downstream of the cylinder, which is progressively shifted towards the channel side wall at y = −0.2 mm.
image file: c8sm02099j-f5.tif
Fig. 5 Time-averaged velocity fields (left column) and corresponding time-averaged retardation fields (right column) in the xy plane (focus at z = 0) for the wormlike micellar solution flowing around the cylinder at the Wi indicated (the Reynolds number is always negligible Re < 10−4). Unsteadiness is apparent in both the velocity and retardation fields for Wi ≳ 130 see accompanying movies provided in the ESI.

The asymmetry shown in Fig. 5c–e contrasts significantly with flow patterns reported previously for WLM solutions in microscale cylinder geometries with higher blockage ratios.53,54 Such asymmetry in the flow was also not reported by Moss and Rothstein for a similar CPyCl/NaSal WLM fluid in a macroscale cylinder device of similar blockage ratio, however their experiments were restricted to Wi ≲ 10 due to limitations conferred by the scale of their setup.31 Similar-looking flow asymmetries were recently reported for the flow of a CTAB-based WLM fluid around a flexible macroscale cylinder and were attributed to the off-center displacement of the cylinder in the flow.36 However in the present experiment the cylinder is considered effectively rigid since no deflection of the cylinder is detectable at the magnifications employed in the experiments. Indeed, our best estimates of the deflection of the cylinder at its midpoint Δxmax, based on the solution to the Euler–Bernouli equation for a circular cylinder with two pinned ends and estimates of the force per unit length on the cylinder, indicate ΔxmaxO(1 μm), which is on the same scale as the cylinder surface roughness.

At the highest Weissenberg numbers examined in our experiments (see e.g. Fig. 5f for Wi = 1875), we observe that the flow field starts to recover some symmetry. This re-symmetrization, though not complete, is accompanied by a significant increase in the magnitude of the birefringence downstream of the cylinder. It is important to note that for Wi ≳ 130 time-dependence is evident in both the flow field and the birefringence. Representative movies to demonstrate these spatio-temporal fluctuations at higher values of the Weissenberg number are included in the ESI. In the ESI, we also present a similar sequence of time-averaged velocity and retardation fields as shown in Fig. 5 but viewed in the orthogonal (x, z) plane.

3.2.1 Time-averaged flow-induced birefringence. For analysis of the flow-induced birefringence signals, we extract data from a location on the flow axis a distance of 5 radii downstream of the cylinder R|0.1mm,0mm. In Fig. 6 we present the time-averaged values [R with combining macron] as a function of the imposed Wi, where the error bars represent a standard deviation about the mean. We observe transitions through four broad regimes of flow as the Weissenberg number is increased. At lower Wi, there is a regime (I) in which the birefringence increases steadily and monotonically with increasing Wi, and the standard deviation is low (error bars are smaller than the data points). Regime I corresponds to steady symmetric flow around the cylinder. As the curve passes through an inflection point at Wic ≈ 60, regime II corresponds to the steady flow asymmetry. In regime II, the mean birefringence in the wake tends towards a local maximum value while the standard deviation remains low since the flow remains steady. Regime II ends abruptly as the flow becomes time-dependent in regime III at Wic2 ≈ 130. Here, error bars are large due to the significant fluctuations and the mean birefringence in the wake initially decreases, then increases with the imposed Wi. Finally, at Wic3 ≈ 1000 a third transition to regime IV is marked by a distinct jump in the mean birefringence and a reduction in the standard deviation of the fluctuations. Regime IV corresponds to the re-symmetrization of the flow as depicted in Fig. 5f.
image file: c8sm02099j-f6.tif
Fig. 6 Time-averaged retardation [R with combining macron] measured at the point (x, y) = (0.1 mm, 0 mm) as a function of the imposed Weissenberg number. Error bars represent the standard deviation of R over the same time period. The transitions between the various flow states are apparent, as described in the text.
3.2.2 Characterization of the steady flow asymmetry. In order to characterize the development of the steady flow asymmetry, we use time averaged velocity fields to compute an asymmetry parameter I based on the streamwise flow velocity measured on either side of the cylinder, halfway between the cylinder and the channel side wall:
image file: c8sm02099j-t6.tif(5)
where u1 = u|0mm,0.11mm and u2 = u|0mm,−0.11mm, as illustrated in Fig. 7a and b.

image file: c8sm02099j-f7.tif
Fig. 7 Determination of the asymmetry parameter I as a function of Wi.

For lower Weissenberg numbers u1 = u2 and hence I = 0. However, as the Weissenberg number exceeds a critical value Wic and the flow around the cylinder becomes unbalanced, u1u2 and hence I increases, see Fig. 7c. Close to the transition, the growth of the asymmetry parameter with Wi is well described by a simple Landau-type quartic potential minimized as:

Wi = Wic[gI2 + hI−1 + 1], (6)
with the critical Weissenberg number Wic = 61 and the growth rate coefficient g = 1. A small asymmetric term h = 5 × 10−4 is included in order to account for biases due to system imperfections. The fit with eqn (6) indicates that the transition to asymmetric flow is a supercritical pitchfork type birfurcation. As shown in Fig. 7c, we have performed several tests (with good agreement) including sweeps of increasing and decreasing Wi, which do not show hysteresis. The small value of h relative to g suggests the bifurcation is almost perfect and hence provides affirmation of the quality of our SLE-microfabricated device. As shown in Fig. 7c, the onset of time-dependent flow (i.e. for Wi ≳ 130) is marked by a significant increase in the scatter of data points.

The root cause of the flow asymmetry remains uncertain, but one possibility is that it develops from an initially random sideways fluctuation of the highly-stressed downstream birefringent wake. Due to the enhanced extensional viscosity in the wake, we propose that such a fluctuation would cause an (initially small) imbalance between the flow resistance around either side of the cylinder, forcing the fluid to prefer the least resistant path. If this occurs while the shear rate on either side of the cylinder is in the shear-thinning region of the flow curve, where n ≈ 0, the resulting mismatch between the fluid viscosity on either side of the cylinder would consolidate the imbalance in flow resistance and maintain the asymmetry. Further increases in Wi would then cause the instability to grow as the viscosity mismatch increases. At the critical Weissenberg number Wic = 61 the average shear rate in the fluid either side of the cylinder is [small gamma, Greek, dot above]side ≈17 s−1, which is firmly in the stress-plateau region (see Fig. 2a), providing some support to our conjecture. In principle, this hypothesis could be tested using devices with even lower blockage ratios, which would give access to the same range of Wi, allowing the downstream birefringence to develop, while maintaining [small gamma, Greek, dot above]side < [small gamma, Greek, dot above]*. However, calculations indicate that given the same WLM test fluid as used here, a blockage ratio β < 0.002 would be required to achieve such a scenario, which would present a significant technical challenge in terms of device fabrication. Insight could also be gained by comparing the responses of non-shear thinning viscoelastic (Boger) fluids with shear thinning but relatively inelastic fluids such as xanthan gum solutions, for example. The importance of micelle breakage and reformation could perhaps be understood by contrasting the responses of two structurally similar polymers such as ethyl hydroxy–ethyl cellulose (EHEC), which does not self-assemble in solution, and its hydrophobically modified analogue (hmEHEC), which is micelle forming.84

3.2.3 Characterization of time-dependent flows. To illustrate graphically the onset of time-dependence in flow as recorded in the velocimetry measurements, Fig. 8a–d show space time diagrams composed of profiles of the x-component of the velocity (u) along the line x = 0 sampled at 20 Hz over a 20 s interval. As a control case, Fig. 8a shows the behaviour of a Newtonian fluid flowing at Q = 0.08 mL per minute (corresponding to a Reynolds number Re ≈ 6 × 10−4). Clearly the flow is highly symmetric to either side of the cylinder and (apart from random noise) does not fluctuate in time. Fig. 8b shows the response of the CPyCl/NaSal solution at a fairly modest Wi = 37.5, also indicating an essentially steady and symmetric flow around the cylinder. At a higher Wi = 93.8 (Fig. 8c), the flow has become significantly asymmetric, with clearly higher velocities recorded above than below the cylinder, however temporal variations in the flow appear to be of a similar magnitude as for the two previous cases shown in Fig. 8a and b. Finally, as the Weissenberg number is increased to Wi = 150 (Fig. 8d), the flow remains asymmetric and a clearly periodic variation in the flow velocity is evident on both sides of the cylinder. The data in Fig. 8a–d is used to evaluate a time-resolved asymmetry parameter I(t), shown in Fig. 8e. For the Newtonian fluid and the CPyCl/NaSal solution under steady conditions, I has an essentially constant value over the 20 s time interval. However, at Wi = 150, the CPyCl/NaSal solution displays a time-dependent value of I, based around a value of I ≈ 0.8 and with periodic spikes of higher asymmetry where I can be greater than 1 (implying that fluid momentarily flows backwards on one side of the cylinder). Fig. 8f shows the power spectral density (PSD) as a function of the frequency (f) obtained from the Fourier transform of the signals presented in Fig. 8e. For the Newtonian fluid and the CPyCl/NaSal solution at the two lower values of Wi, there are no peaks in the PSD confirming that the flow is indeed steady in all those cases. In contrast, at Wi = 150 a peak is observed in the PSD at a fundamental frequency f0 ≈ 0.6 Hz. We note that this is close to the reciprocal of the relaxation time of the WLM solution, f0 ≈ 1/λM, and the agreement seems unlikely to be coincidental. In fact, during preliminary studies, time-resolved velocimetry experiments were carried out at a slightly higher ambient temperature of ≈25.5 °C and resulted in a fundamental frequency f0 ≈ 0.8 Hz, at a similar Wi ≈ 143. SAOS measurements carried out at this same temperature yielded a value for the relaxation time λM ≈ 1.2 s, again showing f0 ≈ 1/λM (data not shown).
image file: c8sm02099j-f8.tif
Fig. 8 Analysis of the time-dependent asymmetry parameter at lower Weissenberg numbers: (a–d) show space-time diagrams composed of the normalized u-component of the velocity field measured along the line x = 0 at a rate of 20 Hz and over a time period of 20 s. (a) Newtonian fluid, (b) CPyCl/NaSal at Wi = 37.5, (c) CPyCl/NaSal at Wi = 93.8, (d) CPyCl/NaSal at Wi = 150. Note that experiments (a) and (d) are carried at equal volumetric flow rates Q = 0.08 mL per minute. (e) Time-dependent asymmetry parameter I for the four cases depicted in (a–d), and (f) power spectra of the fluctuations shown in (e).

Time dependence in the birefringence measurements also show characteristic frequencies similar to those obtained from the μ-PIV. As illustrated in Fig. 9a and b, we extract a time-resolved retardation signal R(t) from the location (x, y) = (0.1 mm, 0 mm) and subtract the time-averaged value [R with combining macron]. The signal appears highly sinusoidal, though we note a sudden cut off on some of the peaks that is reminiscent of the slow build up and subsequent rupture of stress, as reported previously for related kinds of flows.23,30,31,35,36,85 This suggests that a periodic breakage and reformation of micelles in the highly stressed birefringent region plays an important role in driving the flow to become time dependent. The birefringent signal shown in Fig. 9b results in a power spectrum with a clear single peak at a frequency f0 ≈ 0.5 Hz (Fig. 9c), close to that determined from the velocimetry data at the same Wi = 150 near the onset of time dependence.

image file: c8sm02099j-f9.tif
Fig. 9 (a) Space-time diagram depicting the retardation measured over a 20 s interval 5 cylinder radii downstream (i.e. along the line x = 0.1 mm) at a Weissenberg number Wi = 150. (b) Time-resolved retardation signal R[R with combining macron] extracted along the dashed line shown in (a). (c) Power spectrum of the retardation signal, showing a single characteristic frequency f0 ≈ 0.5 Hz.

As the Weissenberg number is increased, we observe changes in the form of the temporal fluctuations in I and their power spectra (see Fig. 10a–c at Wi = 188, 375 and 450, respectively). In contrast to the time-dependence exhibited in Fig. 8e, at progressively higher Wi the asymmetry parameter appears to settle towards a value I ≈ 1 with periodic downward spikes to increasingly lower asymmetry. In the corresponding PSDs, we observe a gradual shift of the fundamental frequency to a lower value f0 ≈ 0.3 Hz (i.e. a kind of gradual period doubling transition) and the emergence of higher harmonics.

image file: c8sm02099j-f10.tif
Fig. 10 Analysis of the time-dependent asymmetry parameter for flow of CPyCl/NaSal solution at moderate Weissenberg numbers: (a) Wi = 188, (b) Wi = 300, (c) Wi = 450 (left) time-dependent asymmetry parameter I, and (right) power spectra of the fluctuations shown to the left. For comparison, the fluctuations and power spectra of a Newtonian fluid under equivalent imposed conditions are also shown.

With increasing Wi, changes are also evident in the fluctuations observed in the retardation signal (see Fig. 11a for Wi = 938). Here we observe a generally rather low level of birefringence behind the cylinder, with periodic spikes to higher intensities. The intense peaks in the birefringence correspond approximately to the spikes of low asymmetry seen in the I(t) traces (such as in Fig. 10c, for example). As the asymmetry returns to a value I ≈ 1 subsequent to a downward spike, the birefringence slowly relaxes back to a low intensity. Similar to the I(t) traces in Fig. 10, the time-dependent birefringence signal at Wi = 938 also results in a complex PSD exhibiting harmonics, as illustrated in Fig. 11b.

image file: c8sm02099j-f11.tif
Fig. 11 (a) Time-resolved retardation signals extracted at (x, y) = (0.1 mm, 0 mm) and, (b) corresponding power spectra of the fluctuations for Wi = 938.

It is worth further clarifying the connection between the time-resolved asymmetry parameter and retardation signals under flow conditions approaching the onset of regime IV. The dynamics of these flows seem to result from the periodic accumulation and rinsing of the almost stagnant fluid located on one side of the cylinder (to negative y as depicted in the figures in this paper). As mentioned, the spike in birefringence correlates to a sudden reduction in the flow asymmetry parameter, which occurs as the stagnant fluid is suddenly washed downstream. Following this spike, the region of stagnant fluid re-establishes at negative y, once again forcing nearly all of the fluid to flow around the opposite side of the cylinder (i.e. to positive y). As a result, the geometry is effectively transformed into a 2[thin space (1/6-em)]:[thin space (1/6-em)]1 asymmetric contraction. Rather than being due to the strong stagnation point extensional flow in the cylinder wake, the micellar alignment and stress in the fluid is now largely generated by the relatively weak squeezing flow between the cylinder and the channel side wall at y = 0.2 mm. Downstream of the cylinder, the birefringence is only maintained by the shear along the streamline dividing the stagnant and the flowing regions. Consequently, the downstream birefringence relaxes over time. Periodically (presumably due to a build up of pressure upstream), the stagnant accumulated material at negative y is swept away again and travels downstream; the asymmetry parameter displays a negative spike and the birefringence intensity increases dramatically as the fluid is (momentarily) subjected to the stagnation point.

The appearance of a single characteristic peak in the power spectrum above Wic2 is potentially characteristic of a Hopf bifurcation, which have been observed previously in complex viscoelastic flows of polymer solutions through axisymmetric contractions.86,87 In an attempt to confirm this we examine the standard deviations of the fluctuations in the asymmetry parameter (sI), which are proportional to the amplitude of the fluctuations. In Fig. 12, we show sI2 as a function of the imposed Wi. At low Wi, fluctuations are small and sI2 is close to zero. At Wi = Wic2, the flow becomes time dependent and the amplitude of the fluctuations becomes suddenly large. Close to the transition, sI2 decays linearly with increasing Wi, which is characteristic of a subcritical (or backwards) Hopf bifurcation.88 This contrasts with earlier studies on the non-linear dynamics of polymer solutions in complex flows, for which transitions to time-dependent flows were characterized as supercritical Hopf bifurcations.86,87 The increase in sI2 observed in Fig. 12 for Wi ≳ 300 likely results from the onset of more complex dynamics, i.e. the presence of harmonics, as discussed above.

image file: c8sm02099j-f12.tif
Fig. 12 Square of the standard deviation of the fluctuations observed in the asymmetry parameter sI2, plotted as a function of the Weissenberg number. Close to the transition to time dependence above Wic2 ≈ 130, sI2 decays linearly with Wi, characteristic of a subcritical (backwards) Hopf bifurcation.

The progression of behaviour so far in the present experiments for 1.9 < Wi ≲ 1000 has quite strong similarities with earlier studies of the nonlinear dynamics of viscoelastic fluids in complex geometries. In particular, McKinley et al. (1991) studied a non-shear thinning polymeric Boger fluid composed of polyisobutylene in polybutene flowing through an abrupt axisymmetric contraction and used laser-Doppler velocimetry to characterize the time dependent fluctuations of the resulting corner vortices.87 The initial onset of time dependence (for Wi ≳ 1.5) was determined to be a supercritical Hopf bifurcation. Subsequently, at higher Wi period-doubling and quasi-periodic behaviour was reported. The sequence appears to be a classical route to chaos, and indeed at higher values of the Weissenberg number Wi ≳ 4, McKinley et al. reported an ultimate transition to a state of aperiodic fluctuations. Such a state would now most likely be referred to as ‘elastic turbulence’, which can be loosely described as self-sustained broadband spatio-temporally fluctuating motions driven by elasticity in the absence of inertia.89–92 Aperiodic fluctuations have been reported for strong extensional flows of viscoelastic WLM solutions in cross-slot devices at high Wi (including for similar CPyCl/NaSal solutions as used here).38,41 Elastic turbulence in WLM solutions has also been reported in shear-banded Taylor–Couette flows, due to the onset of an interfacial instabilitiy between the low and high shear rate bands.93–95

Given the aforementioned discussion, it might then be reasonable to expect a similar aperiodic or turbulent-like behaviour to emerge in the present system as the Weissenberg number is increased to even larger values into the regime IV beyond Wic3 ≈ 1000. However, this is not the case.

As shown in Fig. 13 for Wi = 1875 (i.e. within regime IV), a rather clear single peak arises at a high frequency f0 ≈ 11 Hz in the power spectrum of birefringence fluctuations.

image file: c8sm02099j-f13.tif
Fig. 13 (a) Time-resolved retardation signals extracted at (x, y) = (0.1 mm, 0 mm) and, (b) corresponding power spectra for Wi = 1875.

Fig. 14a and b show the time series and corresponding PSD for fluctuations in I in regime IV at a Weissenberg number of Wi = 2813. Here we very clearly observe a single sharp peak in the power spectrum at a frequency f0 ≈ 15 Hz (Fig. 14b). The Newtonian response under equivalent imposed flow conditions, shown for comparison, confirms that this frequency response is material and not the result of instrument noise. The transition into this fourth flow state occurs as the blockage of stagnant material to one side of the cylinder becomes permanently swept away; presumably, the pressure upstream of the cylinder reaches a level such that the presence of the blockage becomes unsustainable. With fluid now passing around both sides of the cylinder, a semblance of symmetry is recovered in the flow field. Fluid elements become once again subjected to the strong extensional flow in the downstream wake of the cylinder, leading to the observed jump in the time-averaged stress birefringence in the wake (shown in Fig. 6). We suggest that the fast temporal dynamics may emerge due to the consequent breakage of micelles into shorter species with a shorter relaxation time on the order λ* ∼ 1/15 s. Such a scenario would have a close consistency with the recent modeling work from Kalb et al., who used the two-species VCM model to study the effect of micelle breakage on asymmetries arising in flows of WLM solutions in the cross-slots geometry.38,40,41,43,44 Micellar solutions composed of long unbroken micelles displayed strong flow asymmetries above a certain Weissenberg number. The onset of asymmetry effectively weakened the streamwise velocity gradients, reducing the degree of micellar alignment and hence stress birefringence. In contrast, for solutions in which the WLMs broke easily into shorter species the flow remained symmetric, with consequently strong streamwise velocity gradients being maintained and hence a higher level of micellar orientation. The breakage of micelles in the cylinder wake, proposed to explain the present results, could potentially be confirmed directly by small-angle neutron scattering experiments.

image file: c8sm02099j-f14.tif
Fig. 14 Analysis of the time-dependent asymmetry parameter for flow of CPyCl/NaSal solution at Wi = 2813: (a) time-dependent asymmetry parameter I, and (b) power spectra of the fluctuations shown in (a). For comparison, the fluctuations and PSD for a Newtonian fluid under equivalent imposed conditions are also shown.

4 Discussion and conclusions

In this work we have examined the flow of a well-known shear-banding wormlike micellar fluid around a novel type of microfluidic cylinder geometry with a high aspect ratio and a low blockage ratio. Using selective laser-induced etching to fabricate the geometry in fused silica allows the cylinder to be made with a small radius and high aspect ratio while remaining effectively rigid and undeformed under even quite high imposed flow rates. The small radius thus allows us to study very high Weissenberg number flows (up to Wi = 3750) around the cylinder (and also means the Reynolds number remains negligible). In studies with a similar fluid in a similar (but macroscale) geometry, Weissenberg numbers were restricted to Wi ≲ 10, and the flow was found to always remain steady and symmetric.31 By contrast, in our experiments the flow remains steady and symmetric at the lower range of Weissenberg numbers examined Wi ≲ 60, but we observe an extremely rich sequence of instabilities as Wi is increased beyond this value.

A first instability is a supercritical pitchfork bifurcation at a critical Weissenberg number Wic ≈ 60, which is characterized by a steady flow asymmetry where more fluid passes around one side of the cylinder than the other. We suggest that this asymmetry may develop as a result of a small random fluctuation in the position of the initially symmetric birefringent wake downstream of the cylinder. A small random asymmetry in the location of the birefringence could lead to an imbalance in the flow resistance (and hence flow rate) at each side of the cylinder. Due to the difference in shear rate at either side of the cylinder and the strong shear thinning of the fluid viscosity, any small imbalance in the flow resistance could become compounded, hence leading to growth of the instability. As the asymmetry intensifies with increasing Wi an effectively stagnant blockage develops on one side of the cylinder while essentially all of the fluid is diverted around the opposite side of the cylinder. The flow geometry is thus effectively transformed from flow around a cylinder to flow through a 2[thin space (1/6-em)]:[thin space (1/6-em)]1 asymmetric contraction.

Above a second critical Weissenberg number Wic2 ≈ 130, the flow becomes time dependent, displaying a single characteristic frequency f0 in the power spectra of both the flow field and stress birefringence fluctuations which appears to be robustly connected to the single mode Maxwell relaxation time of the fluid: f0 ≈ 1/λM. By examining the standard deviations of the fluctuations close to this transition, we characterized it as possibly being a backwards Hopf bifurcation. The time-dependence appears to be controlled by the periodic accumulation and washing away of the stagnant fluid blocking the flow around one side of the cylinder.

As the Weissenberg number is increased beyond Wi ≈ 300, we observe a progressive slowing down of the dynamics, and the emergence of harmonics in the power spectra of the fluctuations. The progression of instabilities appears to be a fairly standard route towards aperiodic behaviour (or ‘elastic turbulence’), however this progression is abruptly halted at a third transition occuring at Wic3 ≈ 1000. Beyond Wic3 the stagnant blockage is permanently washed away, the flow field regains a degree of symmetry and there is a marked increase in the flow birefringence in the wake of the cylinder. There is also a significant increase in the frequency of the fluctuations, with peaks arising in the power spectra at around 10 to 15 Hz. Based on recent numerical simulations using the 2-species VCM model to study the effect of micelle breakage on flow instabilities in the cross-slot geometry,43,44 these observations seem consistent with the possible breakage of WLM chains into shorter sections with shorter relaxation times.

While somewhat similar sequences of nonlinear transitions have been reported for flows of viscoelastic polymer solutions in complex flow geometries, we believe this is the first such study revealing related phenomena in more complex WLM solutions, which have the ability to break and reform dynamically in flow. Fruitful insights into the importance of self-assembly in driving such transitions could be extracted from a comparative study of structurally similar polymers which either remain isolated or self-assemble when in solution. An example of such a system could be EHEC (non-self-assembling) and its hydrophobically modified analogue hmEHEC, which forms micelles. The quantitative results reported here will be of fundamental relevance to optimizing processes involving flows of viscoelastic surfactant solutions, for example in oilfield operations, and provide insight into the anomalous motions of particles settling in WLM fluids. In addition they will provide valuable data for benchmarking numerical simulations of viscoelastic constitutive equations designed to describe WLM rheology.

Conflicts of interest

There are no conflicts to declare.


SJH, KT-P and AQS gratefully acknowledge the support of the Okinawa Institute of Science and Technology Graduate University (OIST) with subsidy funding from the Cabinet Office, Government of Japan. SJH and AQS also acknowledge funding from the Japan Society for the Promotion of Science (Grants-in-Aid for Scientific Research (C), Grant No. 18K03958 and 17K06173, respectively). SJH is grateful to TA Instruments for the donation of a DHR3 rheometer through their Distinguished Young Rheologist Award. We thank Mr San To Chan (OIST) for performing Newtonian flow simulations using ANSYS Fluent and we thank Photron Ltd for the loan of a CRYSTA PI-1P high-speed polarization camera.


  1. J. N. Israelachvili, Intermolecular and Surface Forces: With Applications to Colloidal and Biological Systems, Academic Press, London, 1985 Search PubMed.
  2. H. Hoffmann, H. Löbl, H. Rehage and I. Wunderlich, Tenside Deterg., 1985, 22, 290–298 CAS.
  3. H. Rehage and H. Hoffmann, J. Phys. Chem., 1988, 92, 4712–4719 CrossRef CAS.
  4. H. Hoffmann, Structure and Flow in Surfactant Solutions, 1994, ch. 1, pp. 2–31 Search PubMed.
  5. F. Nettesheim and N. J. Wagner, Langmuir, 2007, 23, 5267–5269 CrossRef CAS PubMed.
  6. H. Rehage and H. Hoffmann, Mol. Phys., 1991, 74, 933–973 CrossRef CAS.
  7. T. Shikata, S. J. Dahman and D. S. Pearson, Langmuir, 1994, 10, 3470–3476 CrossRef CAS.
  8. S. J. Candau and R. Oda, Colloids Surf., A, 2001, 183, 5–14 CrossRef.
  9. T. M. Clausen, P. K. Vinson, J. R. Minter, H. T. Davis, Y. Talmon and W. G. Miller, J. Phys. Chem., 1992, 96, 474–484 CrossRef CAS.
  10. T. S. Davies, A. M. Ketner and S. R. Raghavan, J. Am. Chem. Soc., 2006, 128, 6669–6675 CrossRef CAS PubMed.
  11. S. R. Raghavan, Langmuir, 2009, 25, 8382–8385 CrossRef CAS PubMed.
  12. K. Trickett and J. Eastoe, Adv. Colloid Interface Sci., 2008, 144, 66–74 CrossRef CAS PubMed.
  13. J. Yang, Curr. Opin. Colloid Interface Sci., 2002, 7, 276–281 CrossRef CAS.
  14. S. Kefi, J. Lee, T. Pope, P. Sullivan, E. Nelson, A. Hernandez, T. Olsen, M. Parlar, B. Powers, A. Roy, A. Wilson and A. Twynam, Oilfield Rev., 2005, 16, 10–23 Search PubMed.
  15. S. Ezrahi, E. Tuval and A. Aserin, Adv. Colloid Interface Sci., 2006, 128, 77–102 CrossRef PubMed.
  16. B. Kalpakci, E. E. Klaus, J. L. Duda and R. Nagrajan, Soc. Pet. Eng. J., 1981, 21, 709–720 CrossRef CAS.
  17. J. Holweg, P. O. Brunn and F. Durst, Proceedings, 4th Europ. Symp. on Enhanced Oil Recovery, DGMK, Hamburg, 1987, pp. 1007–1018 Search PubMed.
  18. J. Holweg, P. O. Brunn and F. Durst, in Progress and Trends in Rheology II, ed. H. Giesekus and M. Hibberd, Steinkopff, 1988, pp. 195–197 Search PubMed.
  19. P. O. Brunn and J. Holweg, J. Non-Newtonian Fluid Mech., 1988, 30, 317–324 CrossRef CAS.
  20. E. Ruckenstein, P. O. Brunn and J. Holweg, Langmuir, 1988, 4, 350–354 CrossRef CAS.
  21. P. O. Brunn and J. Holweg, The Flow of Surfactant Solutions Through Porous Media: Universal Laws, Springer, Netherlands, 1990, pp. 78–80 Search PubMed.
  22. J. Vorwerk and P. O. Brunn, J. Non-Newtonian Fluid Mech., 1994, 51, 79–95 CrossRef CAS.
  23. J. P. Rothstein, Rheol. Rev., 2008, 6, 1–46 Search PubMed.
  24. M. F. Torres, J. M. Gonzaléz, M. R. Rojas, A. J. Müller, A. E. Sáez, D. Löf and K. Schillén, J. Colloid Interface Sci., 2007, 307, 221–228 CrossRef CAS PubMed.
  25. M. R. Rojas, A. J. Müller and A. E. Sáez, J. Colloid Interface Sci., 2008, 326, 221–226 CrossRef CAS PubMed.
  26. R. Cressely and R. Hocquart, Optica Acta, 1980, 27, 699–711 CrossRef CAS.
  27. M. D. Chilcott and J. M. Rallison, J. Non-Newtonian Fluid Mech., 1988, 29, 381–432 CrossRef CAS.
  28. O. G. Harlen, J. M. Rallison and M. D. Chilcott, J. Non-Newtonian Fluid Mech., 1990, 34, 319–349 CrossRef CAS.
  29. J. R. Gladden and A. Belmonte, Phys. Rev. Lett., 2007, 98, 224501 CrossRef PubMed.
  30. G. R. Moss and J. P. Rothstein, J. Non-Newtonian Fluid Mech., 2010, 165, 1–13 CrossRef CAS.
  31. G. R. Moss and J. P. Rothstein, J. Non-Newtonian Fluid Mech., 2010, 165, 1505–1515 CrossRef CAS.
  32. A. Jayaraman and A. Belmonte, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2003, 67, 065301 CrossRef PubMed.
  33. S. Chen and J. P. Rothstein, J. Non-Newtonian Fluid Mech., 2004, 116, 205–234 CrossRef CAS.
  34. H. Mohammadigoushki and S. J. Muller, J. Rheol., 2016, 60, 587–601 CrossRef CAS.
  35. A. A. Dey, Y. Modarres-Sadeghi and J. P. Rothstein, J. Fluid Mech., 2017, 813, R5 CrossRef.
  36. A. A. Dey, Y. Modarres-Sadeghi and J. P. Rothstein, Phys. Rev. Fluids, 2018, 3, 063301 CrossRef.
  37. P. E. Arratia, C. C. Thomas, J. Diorio and J. P. Gollub, Phys. Rev. Lett., 2006, 96, 144502 CrossRef CAS PubMed.
  38. J. A. Pathak and S. D. Hudson, Macromolecules, 2006, 39, 8782–8792 CrossRef CAS.
  39. R. J. Poole, M. A. Alves and P. J. Oliveira, Phys. Rev. Lett., 2007, 99, 164503 CrossRef CAS PubMed.
  40. S. J. Haward and G. H. McKinley, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2012, 85, 031502 CrossRef PubMed.
  41. S. J. Haward, T. J. Ober, M. S. N. Oliveira, M. A. Alves and G. H. McKinley, Soft Matter, 2012, 8, 536–555 RSC.
  42. N. Dubash, P. Cheung and A. Q. Shen, Soft Matter, 2012, 8, 5847–5856 RSC.
  43. A. Kalb, L. A. Villasmil-Urdaneta and M. Cromer, Phys. Rev. Fluids, 2017, 2, 071301 CrossRef.
  44. A. Kalb, L. A. Villasmil-Urdaneta and M. Cromer, J. Non-Newtonian Fluid Mech., 2018, 262, 79–91 CrossRef CAS.
  45. P. A. Vasquez, G. H. McKinley and P. L. Cook, J. Non-Newtonian Fluid Mech., 2007, 144, 122–139 CrossRef CAS.
  46. S. De, J. van der Schaaf, J. A. M. Deen, N. G. Kuipers, E. A. J. F. Peters and J. T. Padding, Phys. Fluids, 2017, 29, 113102 CrossRef.
  47. D. Kawale, E. Marques, P. L. J. Zitha, M. T. Kreutzer, W. R. Rossen and P. E. Boukany, Soft Matter, 2017, 13, 765–775 RSC.
  48. D. Kawale, G. Bouwman, S. Sachdev, P. L. J. Zitha, M. T. Kreutzer, W. R. Rossen and P. E. Boukany, Soft Matter, 2017, 13, 8745–8755 RSC.
  49. S. De, P. Krishnan, J. van der Schaaf, J. A. M. Kuipers, E. A. J. F. Peters and J. T. Padding, J. Colloid Interface Sci., 2018, 510, 262–271 CrossRef CAS PubMed.
  50. N. François, D. Lasne, Y. Amarouchene, B. Lounis and H. Kellay, Phys. Rev. Lett., 2008, 100, 018302 CrossRef PubMed.
  51. S. Kenney, K. Poper, G. Chapagain and G. Christopher, Rheol. Acta, 2013, 52, 485–497 CrossRef CAS.
  52. X. Shi, S. Kenney, G. Chapagain and G. F. Christopher, Rheol. Acta, 2015, 54, 805–815 CrossRef CAS.
  53. C.-L. Sun and H.-Y. Huang, Biomicrofluidics, 2016, 10, 011903 CrossRef PubMed.
  54. Y. Zhao, A. Q. Shen and S. J. Haward, Soft Matter, 2016, 12, 8666–8681 RSC.
  55. K. P. Nolan, A. Agarwal, S. Lei and R. Shields, Microfluid. Nanofluid., 2016, 20, 101 CrossRef.
  56. L. E. Rodd, T. P. Scott, D. V. Boger, J. J. Cooper-White and G. H. McKinley, J. Non-Newtonian Fluid Mech., 2005, 129, 1–22 CrossRef CAS.
  57. L. E. Rodd, J. J. Cooper-White, D. V. Boger and G. H. McKinley, J. Non-Newtonian Fluid Mech., 2007, 143, 170–191 CrossRef CAS.
  58. T. J. Ober, S. J. Haward, C. J. Pipe, J. Soulages and G. H. McKinley, Rheol. Acta, 2013, 52, 529–546 CrossRef CAS.
  59. G. H. McKinley, R. C. Armstrong and R. A. Brown, Philos. Trans. R. Soc., A, 1993, 344, 265–304 CrossRef.
  60. J. Gottmann, M. Hermans and J. Ortmann, Phys. Procedia, 2012, 39, 534–541 CrossRef.
  61. S. J. Haward, K. Toda-Peters and A. Q. Shen, J. Non-Newtonian Fluid Mech., 2018, 254, 23–35 CrossRef CAS.
  62. G. Meineke, M. Hermans, J. Klos, A. Lenenbach and R. Noll, Lab Chip, 2016, 16, 820–828 RSC.
  63. C. J. Pipe, T. S. Majmudar and G. H. McKinley, Rheol. Acta, 2008, 47, 621–642 CrossRef CAS.
  64. J. Y. Lee, G. G. Fuller, N. E. Hudson and X. F. Yuan, J. Rheol., 2005, 49, 537–550 CrossRef CAS.
  65. T. J. Ober, J. Soulages and G. H. McKinley, J. Rheol., 2011, 55, 1127–1159 CrossRef CAS.
  66. R. W. Mair and P. T. Callaghan, Europhys. Lett., 1996, 36, 719–724 CrossRef CAS.
  67. S. M. Fielding, J. Rheol., 2016, 60, 821–834 CrossRef CAS.
  68. R. B. Bird, R. C. Armstrong and O. Hassager, Dynamics of Polymeric Liquids, John Wiley and Sons, New York, 1987 Search PubMed.
  69. M. E. Cates, Macromolecules, 1987, 20, 2289–2296 CrossRef CAS.
  70. M. S. Turner and M. E. Cates, Langmuir, 1991, 7, 1590–1594 CrossRef CAS.
  71. F. Kern, P. Lemarechal, S. J. Candau and M. E. Cates, Langmuir, 1992, 8, 837–840 CrossRef.
  72. V. M. Entov and E. J. Hinch, J. Non-Newtonian Fluid Mech., 1997, 72, 31–54 CrossRef CAS.
  73. S. L. Anna and G. H. McKinley, J. Rheol., 2001, 45, 115–138 CrossRef CAS.
  74. N. J. Kim, C. J. Pipe, K. H. Ahn, S. J. Lee and G. H. McKinley, Korea-Aust. Rheol. J., 2010, 22, 31–41 Search PubMed.
  75. C. D. Meinhart, S. T. Wereley and M. H. B. Gray, Meas. Sci. Technol., 2000, 11, 809–814 CrossRef CAS.
  76. G. G. Fuller, Optical Rheometry of Complex Fluids, Oxford University Press, New York, 1995 Search PubMed.
  77. J. A. Odell, in Handbook of Experimental Fluid Mechanics, ed. C. Tropea, A. L. Yarin and J. F. Foss, Springer-Verlag, Heidelberg, 2007, pp. 724–732 Search PubMed.
  78. R. K. Shah and A. L. London, Laminar flow forced convection in ducts: A source book for compact heat exchanger analytical data, Academic Press, New York, 1978 Search PubMed.
  79. C. Masselon, J. B. Salmon and A. Colin, Phys. Rev. Lett., 2008, 10, 038301 CrossRef PubMed.
  80. C. Masselon, A. Colin and P. D. Olmsted, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2010, 81, 021502 CrossRef PubMed.
  81. V. Lutz-Bueno, R. Pasquino, S. J. Haward, A. Q. Shen and P. Fischer, J. Rheol., 2017, 61, 769–783 CrossRef CAS.
  82. S. J. Haward, F. J. Galindo-Rosales, P. Ballesta and M. A. Alves, Appl. Phys. Lett., 2014, 104, 124101 CrossRef.
  83. P. F. Salipante, C. A. E. Little and S. D. Hudson, Phys. Rev. Fluids, 2017, 2, 033302 CrossRef PubMed.
  84. V. Sharma, S. J. Haward, J. Serdy, B. Keshavarz, A. Soderlund, P. Threlfall-Holmes and G. H. McKinley, Soft Matter, 2015, 11, 3251–3270 RSC.
  85. J. P. Rothstein, J. Rheol., 2003, 47, 1227–1247 CrossRef CAS.
  86. J. V. Lawler, S. J. Muller, R. A. Brown and R. C. Armstrong, J. Non-Newtonian Fluid Mech., 1986, 20, 51–92 CrossRef CAS.
  87. G. H. McKinley, W. P. Raiford, R. A. Brown and R. C. Armstrong, J. Fluid Mech., 1991, 223, 411–456 CrossRef CAS.
  88. T. Kalmár-Nagy, G. Stépán and F. C. Moon, Nonlinear Dyn., 2001, 26, 121–142 CrossRef.
  89. A. Groisman and V. Steinberg, Nature, 2000, 405, 53–55 CrossRef CAS PubMed.
  90. R. G. Larson, Nature, 2000, 405, 27–28 CrossRef CAS PubMed.
  91. A. Groisman and V. Steinberg, Nature, 2001, 410, 905–908 CrossRef CAS PubMed.
  92. P. C. Sousa, F. T. Pinho and M. A. Alves, Soft Matter, 2018, 14, 1344–1354 RSC.
  93. M. A. Fardin, D. Lopez, J. Croso, G. Gregoire, O. Cardoso, G. H. McKinley and S. Lerouge, Phys. Rev. Lett., 2010, 104, 178303 CrossRef CAS PubMed.
  94. M. A. Fardin, T. J. Ober, V. Grenard, T. Divoux, S. Manneville, G. H. McKinley and S. Lerouge, Soft Matter, 2012, 8, 10072–10089 RSC.
  95. M. A. Fardin and S. Lerouge, Eur. Phys. J. E: Soft Matter Biol. Phys., 2012, 35, 91 CrossRef PubMed.


Electronic supplementary information (ESI) available: Movies showing time-resolved flow and retardation fields captured for the wormlike micellar test fluid around the cylinder at various Wi. See DOI: 10.1039/c8sm02099j

This journal is © The Royal Society of Chemistry 2019