Spatiotemporal flow instabilities of wormlike micellar solutions in rectangular microchannels

Spatiotemporal flow instabilities of wormlike micellar solutions in rectangular microchannels S. J. Haward, F. J. Galindo-Rosales, P. Ballesta, and M. A. Alves Citation: Applied Physics Letters 104, 124101 (2014); doi: 10.1063/1.4869476 View online: http://dx.doi.org/10.1063/1.4869476 View Table of Contents: http://scitation.aip.org/content/aip/journal/apl/104/12?ver=pdfcov Published by the AIP Publishing Articles you may be interested in Shear banding in time-dependent flows of polymers and wormlike micelles J. Rheol. 58, 103 (2014); 10.1122/1.4842155 Role of shear-thinning on the dynamics of rinsing flow by an impinging jet Phys. Fluids 24, 093102 (2012); 10.1063/1.4752765 Passive flow control in microdevices using thermally responsive polymer solutions Phys. Fluids 18, 053103 (2006); 10.1063/1.2204077 Investigation of shear-banding structure in wormlike micellar solution by point-wise flow-induced birefringence measurements J. Rheol. 49, 537 (2005); 10.1122/1.1849179 Transition to shear banding in pipe and Couette flow of wormlike micellar solutions J. Rheol. 43, 897 (1999); 10.1122/1.551008 This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 193.136.33.132 On: Wed, 21 May 2014 01:53:10 Spatiotemporal flow instabilities of wormlike micellar solutions in rectangular microchannels S. J. Haward, a) F. J. Galindo-Rosales, P. Ballesta, and M. A. Alves Departamento de Engenharia Qu ımica, Faculdade de Engenharia da Universidade do Porto, Rua Dr Roberto Frias s/n, 4200-465 Porto, Portugal (Received 5 December 2013; accepted 11 March 2014; published online 25 March 2014) Flow velocimetry measurements are made on a non-shear-banding wormlike micellar solution within high-aspect-ratio rectilinear microchannels over a wide range of imposed steady flow rates. At the lowest and highest flow rates tested, Newtonian-like velocity profiles are measured. However, at intermediate flow rates the velocity field never stabilizes on the timescale of the experiments (up to several hours). Here, spatiotemporally dependent “jets” of high velocity fluid are observed to fluctuate within regions of essentially stagnant fluid. The reason for this flow instability remains undetermined, but it has significant consequences for many industrial applications and also for microfluidic rheometry of complex fluids. V C2014 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License. [http://dx.doi.org/10.1063/1.4869476] Wormlike micelles are formed by the self-assembly of surfactant molecules in solution. 1,2 They exist in a dynamic equilibrium with their solvent and for this reason are sometimes referred to as “living polymers.” The dynamic breakage and reformation of the micelles confers unique viscoelastic properties to the fluids, which have made them extremely useful in a wide range of applications including inkjet printing, enhanced oil recovery, district heating systems, and as rheology modifiers in countless consumer products. 3–5 In semidilute to concentrated wormlike micellar systems, the rate of micelle breaking and reformation can dominate the characteristic timescale of the fluid (the fastbreaking limit) and in such cases the resulting linear rheological response of the fluid obeys almost ideal Maxwellian behavior described by a single-mode Maxwell-Debye relaxation mechanism. 6,7 For this reason, semidilute wormlike micellar solutions are often considered to be simple model viscoelastic fluids. However, their non-linear rheological behavior (which is far more relevant to any real-world process or application) is anything but simple. Here, a range of poorly understood phenomena are encountered, including shear-banding, 8–12 flow-induced structure formation, 13,14 and elastic and interfacial instabilities. 15–21 In this Letter, we consider a semidilute wormlike micellar solution formed from 30 mM of the surfactant cetyltrimethyl ammonium bromide (CTAB) in the presence of 240 mM sodium salicylate (NaSal) flowing in a simple highaspect ratio rectangular rectilinear microchannel. The fluid is seeded with 1 lm diameter fluorescent particles (volume fraction /2104) to enable micro-particle image velocimetry (l-PIV) measurements, described below. Measurements of the steady shear rheology in a cone-andplate rotational rheometer at 25 C show the fluid is strongly shear-thinning, but does not display the clear plateau of constant stress that is characteristic of shear-banding micellar fluids (Fig. 1). 8,9 This rheological response is typical of fluids of the particular formulation used here. 22,23 The steady rheology is well described by the Carreau generalized Newtonian fluid (GNF) model gð_ cÞ¼g1þg0g1 ðÞ 1þ_ c=_ c ðÞ 2 hi n1 ðÞ =2 ;(1) where g1¼0:07 Pa s is the infinite-shear-rate viscosity, g0¼3:3 Pa s is the zero-shear-rate viscosity, _ cis the applied shear rate, _ c¼0:95 s1is the characteristic shear rate for the onset of shear-thinning, and n¼0:28 is the “power-law exponent” in the shear-thinning region. The microchannel employed in most of the experiments has a depth d¼275 65lm in the z-direction and width w¼2.2 mm in the y-direction, providing a high aspect ratio a¼w=d¼8 and a hydraulic diameter Dh¼2wd=ðwþdÞ ¼489lm. The total length of the channel in the direction of flow (i.e., the x-direction) is L¼6 cm, therefore L=Dh123. We define the origin of coordinates as being at the center of the channel cross-section and at the location of the channel inlet. Two of the side-walls are made from precision 250 lm-thick steel shim while the remaining two walls are formed from standard glass microscope slides. Bonds between the glass and steel are made using a solventresistant epoxy resin with a layer thickness of approximately 12–15 lm on each side, top and bottom. The flow inside the channel is driven at controlled volume flow rates 0:00022 lls 1Q1lls 1using a high precision neMESYS syringe pump and Hamilton Gastight syringes. To ensure essentially pulsation-free injection of fluid into the channel, syringe volumes are selected such that required flow rate is achieved by driving the pump at no less than 40the minimum specified rate. Typically, the pumping rate used corresponds to at least 100 increments-per-second of the pump stepper motor, or more than 2 orders of magnitude greater than the inverse relaxation time of the CTAB solution, 1=kffi_ c¼0:95 s1. At such high discretization, no feedback between the fluid drive and response is to be a) Author to whom correspondence should be addressed. Electronic mail: [email protected]. 0003-6951/2014/104(12)/124101/5 V CAuthor(s) 2014104, 124101-1 APPLIED PHYSICS LETTERS 104, 124101 (2014) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 193.136.33.132 On: Wed, 21 May 2014 01:53:10 expected. The average flow velocity inside the channel is given by U¼Q=wd and the nominal wall shear rate can be estimated as _ cwall 6U=d. The range of wall shear rates examined in the microchannel experiments is thus 0:01 s1_ cwall 40 s1, closely corresponding to the range of shear rates covered in the rotational rheometer (Fig. 1). The Reynolds number in the channel is defined as Re ¼qUDh=gð_ cÞ, where q1000 kg m3is the fluid density and the shear-rate-dependent viscosity, gð_ cÞ, is estimated using the nominal wall shear rate in Eq. (1). In all the experiments presented here Re <3103, thus inertial effects can be disregarded. The fully-developed flow of the CTAB solution through the rectangular channel is numerically simulated in Matlab (Mathworks, Inc.) by solving the transient Stokes equation using a time marching scheme until convergence is obtained. The fitted Carreau GNF model (Eq. (1)) is used as the constitutive equation for the shear-rate-dependent fluid viscosity. Note that this model accounts only for the shear-thinning and not for any memory effects that might be expected for a micellar solution. For l-PIV measurements, the flow channel is placed on the imaging stage of an inverted microscope (Olympus IX83) and a 4NA ¼0.13 objective lens is focused at the center-plane of the channel depth (i.e., at z¼0). A diodepumped pulsed 526.5 nm neodymium-doped yttrium lithium fluoride (Nd:YLF) laser (Litron Lasers LDY300) excites the fluorescent tracer particles in the fluid (Nile Red, Molecular Probes, Invitrogen Ex/Em ¼535 nm/575 nm). An epifluorescent filter cube is used to prevent the laser light from reaching the sensor of a high-speed CMOS camera (Phantom Miro M340, Vision Research), which thus only receives the light emitted by the fluorescent particles. An additional 1.9 lens is used to ensure the whole CMOS sensor is illuminated. This system allows instantaneous imaging across the entire channel width. Use of higher magnification, higher NA objectives would increase accuracy, but would require scanning of the focal position across the channel width in order to obtain the complete velocity profile. This would not only slow down the data acquisition, which is not desirable in a potentially time-varying flow field, but would also complicate the subsequent image and data analysis unnecessarily. Images are captured continuously at rates such that the average frame-by-frame particle displacement is around 3 pixels, which is optimal for subsequent PIV analysis. The analysis is performed using an adaptive cross-correlation PIV algorithm with interrogation areas that reduce from 64 64 to 16 16 pixels, implemented in Dynamic Studio software (Dantec Dynamics). The maximum full-frame acquisition rate of the Phantom Miro CMOS camera is 800 fps (framesper-second), while the actual acquisition rates used for the experimental measurements are 600 fps at the highest flow rates and proportionally slower at lower imposed flow rates. The adaptive PIV algorithm itself is capable of resolving flow velocity vectors for frame-to-frame pixel displacements of between 0.1 and 20 pixels. Initial control experiments involve seeding a Newtonian fluid (water) with fluorescent tracer particles and making l-PIV measurements for Newtonian flow at low Reynolds numbers. In this case, the streamwise velocity profiles measured across the channel width and at the center-plane for various imposed flow rates are selfsimilar and closely conform to the analytically computed profile (Fig. 2), thus confirming the suitable accuracy of the measurement technique. Depending upon the imposed flow rate, subsequent experiments with the shear-thinning viscoelastic CTAB solution reveal a sequence of entirely unpredictable spatiotemporally varying velocity profiles (Fig. 3). At very low flow rates, e.g., Q¼0:00022lls 1(corresponding to a nominal wall shear rate _ cwall 0:01 s1, Fig. 3(a)), the velocity profile across the width of the microchannel is similar to the Newtonian case, is well-described by the Carreau GNF model and remains essentially invariant over long time periods (>2 h of observation in this case). However, at intermediate flow rates, the velocity profiles across the channel become quite atypical (Figs. 3(b) and 3(c)). Here, as shown in the example given in Fig. 3(b) ðQ¼0:0555 lls 1;_ cwall 2s 1Þ, in some regions of the channel the fluid can be almost stagnant, with essentially all of the volume being transported along narrow “jet-like” sections. The peak flow velocity within the jet can reach values FIG. 1. Steady flow curves of viscosity (g) and shear stress (r) as a function of the shear rate (_ c) measured for the CTAB test fluid at 25 C and fitted with the Carreau model (Eq. (1), solid lines). FIG. 2. Results of l-PIV measurements made at the center-plane of the rectangular microchannel with a Newtonian fluid (water) at low Re, compared with the analytical solution obtained for fully developed laminar Newtonian flow. 124101-2 Haward et al. Appl. Phys. Lett. 104, 124101 (2014) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 193.136.33.132 On: Wed, 21 May 2014 01:53:10 of up to vx=U10. The position of the jet and the width of the jet region appear to be metastable, and fluctuate over time periods on the order of minutes, which is much longer than the relaxation time of the fluid (k1=_ c1:1 s). Experiments have been conducted in which the microchannel is translated along its entire length between the inlet (x¼0 mm) and the outlet ðx60 mmÞ. These high velocity jets can persist for the entire length of the microchannel, only varying slightly in width and position. At this volume flow rate, the average flow velocity in the channel is U0:1mms 1, so the nominal transit time for fluid along the length of the channel is 10 min. Fig. 3(c) shows velocity profiles obtained across the channel at various times for a flow rate of Q¼0:139lls 1ð_ cwall 5s 1Þ. In this particular example, at early times (e.g., 5 min after startup of flow) the profile appears almost plug-like, but with a superimposed high velocity region to the right-hand side of the channel centerline. After 20 min of continuous flow, the flow field looks as though it might be developing towards the expected plug-like profile predicted by the Carreau model. Subsequently, however, an extremely unusual velocity distribution develops in which two high-velocity fluid jets are observed close to the channel side walls and the lowest velocity is measured along the channel centerline. As previously, such velocity distributions appear to be metastable and can remain essentially unaltered for periods spanning many minutes of observation. At the highest flow rate tested ðQ¼1:0lls 1;_ cwall 36 s1Þ, initial fluctuations are observed before the flow velocity distribution approaches the prediction of the Carreau model and subsequently tends to stabilize (Fig. 3(d)). In order to provide confirmation that the jets observed at intermediate flow rates are not due to any pulsation from the syringe pump, experiments have been performed using a range of syringes. For example, at the volume flow rate of Q¼0:139 lls 1(_ cwall 5s 1, Fig. 3(c)) syringes of 100 ll and 2500 ll have been utilized. At this flow rate, these two syringes require the pump to increment at approximately 2400 and 100 steps-per-second, respectively. Use of different syringes to drive the flow results in no qualitative difference in the observation of the jets of fluid within the flow channel, indicating that pulsation of the syringe pump is not a causal factor. Fig. 4shows a space-time diagram constructed from velocity profiles made across the rectilinear microchannel at a single x-location approximately 40 mm downstream from the channel inlet (x=w18) over a 150 min time period for continuous flow of the CTAB solution at Q¼0:0555 lls 1ð_ cwall 2s 1Þ. Measurements commence approximately 10 min after the startup of flow and initially appear to show an almost uniform plug-like velocity profile across the channel. However, after approximately 15 min more of observation, the velocity distribution becomes asymmetric, with a velocity peak growing on the left-hand side and a corresponding trough developing to the right-hand side of the channel. Over the subsequent 10–15 min period, the peak velocity shifts towards the right-hand side of the channel and FIG. 3. Velocity profiles measured at the center-plane of the rectangular microchannel with the 30 mM CTAB 240 mM NaSal solution at 25 C at various flow rates: (a) Q¼0:00022 lls 1;_ cwall 0:01 s1; (b) Q¼ 0:0555 lls 1;_ cwall 2s 1; (c) Q¼ 0:139 lls 1;_ cwall 5s 1; and (d) Q¼1:0lls 1;_ cwall 36 s1. In each case the solid line represents the numerical prediction of the velocity profile determined using the fit of the Carreau model to the steady shear rheology (see Eq. (1), Fig. 1). The times quoted in the figure legends are the approximate times from the startup of flow at each respective flow rate. Profiles are obtained from approximately 40 mm downstream of the channel inlet. FIG. 4. Space-time diagram formed from l-PIV measurements made across the rectangular microchannel at a single x-location ðx40 mmÞover a 2.5 h time period. The volume flow rate is Q¼0:0555 lls 1, providing a nominal wall shear rate _ cwall 2s 1. 124101-3 Haward et al. Appl. Phys. Lett. 104, 124101 (2014) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 193.136.33.132 On: Wed, 21 May 2014 01:53:10 develops into a very high velocity jet centered at y=w0:2 with a maximum velocity of vx=U9. Subsequently, the jet shifts first left then right and also gradually reduces in intensity. Finally, after around 70 min a metastable state is achieved in which the flow is concentrated along the righthand side of the channel with the fluid at the left-hand side being almost stagnant; this state persists with only minor alteration for the remaining 80 min of observation. Fig. 5shows an example of a velocity field in the z¼0 plane of the microchannel for a case in which a significant spatial variation in velocity is observed along the x-direction. In this example, the flow rate is again Q¼0:0555 lls 1ð_ cwall 2s 1Þand the flow has been steadily imposed for 50 min. As fluid flows from the top towards the bottom of the field of view, the flow becomes focused into the typically observed high-speed jet of the form reported in Fig. 2(b). Flow profiles of the type shown here have not been reported previously. Ober et al. 22 conducted experiments with a CTAB solution of the same formulation as used here and in a microchannel of comparable dimensions, (w¼1 mm, d¼130 lm). However, measurements were only made of vxðzÞfor one value of y(0:2y=w0:3), so velocity distributions across the wide channel dimension (i.e., vxðyÞ) were not ascertained. The reported velocity profiles of vxðzÞevolved from parabolic at low flow rates towards a more plug-like form at higher flow rates, in qualitative agreement with numerical predictions 22 and with the results of previous authors with strongly shear-thinning and shear-banding fluids. 24,25 Use of a second rectangular microchannel of dimensions w¼1.50 mm, d¼250 lm, oriented vertically, permits measurements of vxðzÞfor the CTAB solution that can be qualitatively compared to the data of Ober et al. 22 For these experiments, a 20NA ¼0.4 objective is used for the l-PIV measurements (measurement depth 14lm, or 0:01wÞ. 26 The measurement plane is scanned through yin 20 lm increments between the midplane and a side wall, and the data from each plane are stacked to provide the full velocity distribution vxðy;zÞover half of the channel cross-section, as shown in Fig. 6for two different flow rates. The process of acquiring the full velocity fields in this manner takes up to 1 h, so there is most certainly some smearing of the data due to the time delay between acquisition at different measurement planes. However, it is clear that as the flow rate is increased the velocity distribution over the channel crosssection becomes much more uniform (i.e., plug-like), as expected for this shear-thinning fluid. The startling result is the clear presence of the high velocity jet of fluid observed near the side wall of the channel at the lower flow rate of Q¼0:069lls 1(Fig. 6(a)). The nominal wall shear rate in this case is _ cwall 4:4s 1, which is of very similar magnitude to the shear rates at which these anomalous fluid jets were observed in the previous microchannel (Figs. 3(b) and 3(c)). Obviously, by only making measurements in a single plane (at the position approximately marked by the dashed lines on Fig. 6) Ober et al. could not have observed such phenomena in their experiments. 22 Many questions remain about the origin of the spatiotemporally dependent velocity profiles reported here. It is clear that over a wide range of imposed flow rates, the stress distribution across the width of the channel diverges widely from the predictions of the shear-thinning model, therefore pointing towards an elastic cause for the instability. It is possible that this pressure difference is compensated by normal stresses due to fluid elasticity. On the other hand, the instability is only observed at wall shear rates corresponding to the power-law region of the flow curve (within the lowshear-rate Newtonian plateau of the flow curve the instability is not observed, and approaching the high-shear-rate plateau it is significantly suppressed); this suggests shear-thinning being particularly important as the cause. However, the fact that the GNF model does not reproduce this anomalous FIG. 5. Instantaneous l-PIV measurement in the x-y plane over a section of the rectangular microchannel centered on x40 mm ðx=w18:2Þ. The flow direction is from top to bottom, the volume flow rate is Q¼0:0555 lls 1(_ cwall 2s 1) and flow has been imposed for 50 min (Multimedia view) [URL: http://dx.doi.org/10.1063/1.4869476.1]. FIG. 6. l-PIV measurements in the y-z plane of a rectangular microchannel (w¼1.50 mm, d¼250 lm) at the flow rates indicated on the plots. The corresponding nominal wall shear rates are: (a) _ cwall 4:4s 1; and (b) _ cwall 89 s1. The flow direction is into the page. The superimposed dashed lines indicate the approximate measurement position used by Ober et al. 22 124101-4 Haward et al. Appl. Phys. Lett. 104, 124101 (2014) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 193.136.33.132 On: Wed, 21 May 2014 01:53:10 behavior, suggests a probable coupling between shearthinning, normal stress, and memory effects. Future attempts to shed light on the problem will initially focus on examining alternative fluids such as shear-thinning but relatively inelastic xanthan gum solutions, and constant viscosity but elastic dilute polymer solutions. This will help to clarify if any particular aspect of the non-Newtonian rheology is at the root. Also interesting will be to examine different wormlike micellar solutions such as the well-studied shear-banding cetylpyridinium chloride CPyCl/NaSal system. 19,20,22,24,27 Ober et al. 22 already presented three-dimensional velocimetry data for such a fluid in their rectangular microchannel. There were no abnormalities apparent in their data, but only one particular flow rate was tested, so the results are not necessarily conclusive. It may also be valuable to probe the fluid microstructure within the high velocity and low velocity regions of the microchannel. Flow-induced birefringence measurements may provide an effective qualitative method for doing this. Notwithstanding the mechanism or reason for the anomalous velocity distributions observed here with the CTAB solution, they clearly have some important consequences. It is normally taken for granted that the flow in such high aspect ratio channels as those used here is essentially twodimensional, i.e., invariant across the large dimension (w,as it is defined here). It is now clear that this assumption cannot be relied upon for all complex fluids. This has practical implications for the use of high aspect ratio microchannels as microfluidic rheometers. 25,28 It is essential the flow field be fully characterized in such devices over the range of shear rates to be measured if the results obtained are to be properly valid. These phenomena also need to be understood and accounted for in order to optimize industrial processes involving transport and mixing of complex fluids in micro to milli-scale geometries such as those encountered in inkjets, porous media, and lab-on-a-chip devices, for example. S.J.H. was supported by a Marie Curie International Incoming Fellowship within the 7th European Community Framework Programme. F.J.G.R. acknowledges financial support from FCT and FEDER through scholarship SFRH/BPD/69663/2010. P.B. and M.A.A. acknowledge financial support from the European Research Council (ERC) under the European Commission “Ideas” specific programme of the Seventh Framework Programme (Grant No. 307499). 1 R. G. Larson, The Structure and Rheology of Complex Fluids (Oxford University Press, New York, 1999). 2 J. F. Berret, Molecular Gels: Materials With Self-Assembled Fibrillar Networks (Springer, Dordrecht, 2006), pp. 667–720. 3 S. Ezrahi, E. Tuval, and A. Aserin, Adv. Colloid Interface Sci. 128–130, 77 (2006). 4 J. Yang, Curr. Opin. Colloid Interface Sci. 7, 276 (2002). 5 J. P. Rothstein, Rheology Reviews (The British Society of Rheology, Aberystwyth, 2008), pp. 1–46. 6 M. E. Cates, J. Phys. Chem. 94, 371 (1990). 7 N. A. Spenley, M. E. Cates, and T. C. B. McLeish, Phys. Rev. Lett. 71, 939 (1993). 8 R. Makhloufi, J. P. Decruppe, A. A€ ıt-Ali, and R. Cressely, Europhys. Lett. 32, 253 (1995). 9 R. Mair and P. T. Callaghan, Europhys. Lett. 36, 719 (1996). 10 S. Manneville, Rheol. Acta 47, 301 (2008). 11 P. D. Olmsted, Rheol. Acta 47, 283 (2008). 12 S. Lerouge and J. F. Berret, Adv. Polym. Sci. 230, 1 (2010). 13 M. Vasudevan, E. Buse, D. L. Lu, H. Krishna, R. Kalyanaraman, A. Q. Shen, B. Khomami, and R. Sureshkumar, Nature Mater. 9, 436 (2010). 14 N. Dubash, J. Cardiel, P. Cheung, and A. Q. Shen, Soft Matter 7, 876 (2011). 15 J. A. Pathak and S. D. Hudson, Macromolecules 39, 8782 (2006). 16 M. A. Fardin, D. Lopez, J. Croso, G. Gr egoire, O. Cardoso, G. H. McKinley, and S. Lerouge, Phys. Rev. Lett. 104, 178303 (2010). 17 P. Nghe, S. M. Fielding, P. Tabeling, and A. Ajdari, Phys. Rev. Lett. 104, 248303 (2010). 18 M. A. Fardin, T. J. Ober, C. Gay, G. Gr egoire, G. H. McKinley, and S. Lerouge, Europhys. Lett. 96, 44004 (2011). 19 S. J. Haward, T. J. Ober, M. S. N. Oliveira, M. A. Alves, and G. H. McKinley, Soft Matter 8, 536 (2012). 20 S. J. Haward and G. H. McKinley, Phys. Rev. E 85, 031502 (2012). 21 M. A. Fardin, T. J. Ober, C. Gay, G. Gr egoire, G. H. McKinley, and S. Lerouge, Soft Matter 8, 910 (2012). 22 T. J. Ober, J. Soulages, and G. H. McKinley, J. Rheol. 55, 1127 (2011). 23 T. Takahashi, M. Shirakashi, K. Miyamoto, and G. G. Fuller, Rheol. Acta 41, 448 (2002). 24 C. Masselon, J.-B. Salmon, and A. Colin, Phys. Rev. Lett. 100, 038301 (2008). 25 P. Nghe, G. Degr e, P. Tabeling, and A. Ajdari, Appl. Phys. Lett. 93, 204102 (2008). 26 C. D. Meinhart, S. T. Wereley, and M. H. B. Gray, Meas. Sci. Technol. 11, 809 (2000). 27 H. Rehage and H. Hoffmann, Mol. Phys. 74, 933 (1991). 28 C. J. Pipe, T. S. Majmudar, and G. H. McKinley, Rheol. Acta 47, 621 (2008). 124101-5 Haward et al. Appl. Phys. Lett. 104, 124101 (2014) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 193.136.33.132 On: Wed, 21 May 2014 01:53:10 CEFT Centro de Estudos de Fenómenos de Transporte Transport Phenomena Research Center