Scott S. H.
Tsai†
ab,
Jason S.
Wexler
b,
Jiandi
Wan‡
b and
Howard A.
Stone
*b
aSchool of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
bDepartment of Mechanical and Aerospace Engineering, Princeton University, NJ 08544, USA. E-mail: hastone@princeton.edu
First published on 5th October 2012
We describe a technique that measures ultralow interfacial tensions using paramagnetic spheres in a co-flow microfluidic device designed with a magnetic section. Our method involves tuning the distance between the co-flowing interface and the magnet's center, and observing the behavior of the spheres as they approach the liquid–liquid interface–the particles either pass through or are trapped by the interface. Using threshold values of the magnet-to-interface distance, we make estimates of the two-fluid interfacial tension. We demonstrate the effectiveness of this technique for measuring very low interfacial tensions, O(10−6–10−5) N m−1, by testing solutions of different surfactant concentrations, and we show that our results are comparable with measurements made using a spinning drop tensiometer.
Despite the utility of these fluid systems, quantifying ultralow interfacial tensions practically and accurately is challenging. Traditional force tensiometry methods such as the Wilhelmy Plate and Du Noüy Ring are not suitable because of difficulties resolving the small forces experienced by the probes at the interfaces.8 Measuring very low interfacial tensions using the spinning drop method is also not ideal since the droplets equilibrate very slowly, making it challenging to extract an equilibrium value.9 The variety of techniques being proposed, including the breaking thread method,9 the laser interface manipulation method,8 the method of measuring static profiles near single walls,10 the method of visual observation of thermal capillary waves,11 and the method of deforming emulsion drops with optical tweezers12 suggest that a simple way to measure ultralow interfacial tensions is desirable.
Here we describe a method for measuring ultralow interfacial tensions by applying external magnetic fields to paramagnetic spheres in a microfluidic co-flow. Our approach has the distinct advantage of reducing the amount of reagents required, as well as the possibility of integration into lab-on-a-chip devices. Although microfluidic tensiometers have been demonstrated before,13–16 all such devices have been developed for measuring typical O(10−3) N m−1 interfacial tensions. Our technique is designed for measuring ultralow interfacial tensions. Specifically, we obtain estimates of the interfacial tension by determining whether and how paramagnetic spheres break through the interface between co-flowing fluids.
Fig. 1 (a) A schematic diagram of the microfluidic tensiometer with a co-flow geometry, where the aqueous solution flows through the buffer inlet to center the beads, which are supplied through another inlet. Dodecane with 10 wt% Span 80 flows through the oil inlet. The permanent magnet draws paramagnetic beads towards, and possibly through, the two-fluid interface. (b) A magnified view of the region of influence of the magnet shows the geometric length scales of the system as well as the magnetization direction of the magnet, M = Mẑ. |
The non-aqueous (oil) solution enters through a side inlet to form a co-flow with the aqueous fluid and the beads. The non-aqueous solution consists of dodecane with 10 wt% sorbitan monooleate surfactant (Span 80), while the aqueous solution is a mixture of deionized water with sodium dodecyl sulfate (SDS) concentrations varying from 4–20 mM. Shum et al.19 reported interfacial tensions γ ≈ (10−6–10−5) N m−1 for similar mixtures, by spinning drop tensiometry. In our experiments, 25 wt% glycerol is added to the aqueous solution to more closely match the viscosities of the aqueous (μa) and non-aqueous (μna) phases, which are μa = 1.8 and μna = 2.7 mPa s, respectively.
A single flat neodymium iron boron magnet (K. J. Magnetics, Jamison, PA, USA), with magnetization M = 106 A m−1 and dimensions m × wm × hm = 9.5 × 4.8 × 1.6 mm, is placed adjacent to the microchannel, with the magnetization direction oriented through the thickness (Fig. 1(b)). A slot is cut in the PDMS next to the channel prior to bonding to the glass slide to allow close placement of the permanent magnet.
The perpendicular distance, , between the magnet center and the fluid-fluid interface in the microchannel (see Fig. 1(b)) is controlled by carefully adjusting the magnet's position with an antimagnetic tweezer (Sigma Aldrich, St. Louis, MI, USA). This distance is further fine-tuned by controlling the relative flow rates of the two co-flowing phases. Constant flow-rate syringe pumps (Harvard Apparatus, Holliston, MA, USA) provide the flow to each inlet, and the total flow rate in the system is always maintained at 8 μL min−1.
The experiments are imaged using bright-field microscopy (Leica Microsystems, Wetzlar, Germany) with 10× and 20× objectives and a high-speed camera (Vision Research, Wayne, NJ, USA). We use ImageJ software to measure the distance and track the particles in the microchannel.
In separate control experiments, we measure interfacial tensions using a Krüss SITE100 spinning drop tensiometer, and follow standard operating procedures. As described later, these measurements give values for the interfacial tension that are comparable to the estimates obtained with our new particle and co-flow approach.
Fig. 2 (a–b) A time-series of images taken at intervals of 5 ms shows 2.5 μm radius paramagnetic spheres, in a solution of 8 mM SDS in water with 25 wt% glycerol, approaching a fluid-fluid interface. The spheres (a) deform and pass through the interface into a solution of 10 wt% Span 80 in dodecane, or (b) become trapped at the two-fluid interface when the magnetic force is not strong enough to pull the spheres through the interface. Plots of particle positions in the directions (c) of the flow and (d) perpendicular to the flow versus time. In the flow direction, the particles move at approximately constant speeds even as they cross the interface (at t = 0 ms), until they slow down as they approach the channel wall. In the perpendicular direction, when the distance between the fluid-fluid interface and the magnet center is above a critical value, the spheres stop passing through the interface and remain trapped in the aqueous phase. |
We observe that the speed of the particles in the flow direction (Fig. 2(c)) is approximately constant even as the particles are being pulled through the interface, which indicates that the spheres are following the bulk fluid flow, and that the speed of the two co-flowing solutions have equilibrated to an average speed, uo, upstream from the location of our measurements. The particles only slow down as they approach the channel wall next to the permanent magnet. Thus, from a Lagrangian perspective, we can approximate the particles as moving only perpendicular to the fluid-fluid interface.
In the perpendicular direction (Fig. 2(d)), we observe that in the case = 3.04 mm, the particle we track does not pass through the interface; rather it remains at the interface upon contact. For lower values of the magnet distance the other particles we follow eventually pass through the interface. We also note that smaller values of the length result in more rapid deformation of the interface and a more prompt passage.
The origin of our coordinate system is the middle of the microchannel, parallel to the front edge of the permanent magnet, as depicted in Fig. 1(b). In our calculations, we non-dimensionalize all length scales, time scales, and the magnetic field via
(1) |
We obtain the full analytical solution to the magnetic field generated by the magnet with dimensions m × wm × hm and magnetization M = Mẑ using the approach described previously.22 Since the final expression is lengthy, and shows that the magnetic field is a function of the magnet's geometry, magnetization, and position relative to the origin, we will simply write the dimensionless expression as
Ĥ = f(X, Y, Z, Wm, Hm, Dc, L), | (2) |
Next, to obtain the particle trajectories we balance the most significant force contributions acting on the particles, which based on arguments proposed in the magnetophoresis literature, are magnetic and hydrodynamic forces.22,23 The magnetic force on a particle is given by22–24
(3) |
Fd = 6πaηa(u − v), | (4) |
(5) |
The particle's trajectory in dimensionless form is
(6) |
(7) |
(8) |
(9) |
We can make scaling and physical arguments to further simplify eqn (6) and (8). Fig. 3(a), which plots against the spatial positions X and Y, depicts large gradients at the magnet corners (X = 0 and X = 1). The large gradients decay quickly away from the corners. We can exploit the decay of the magnetic field gradient and the small parameter = O(10−3) in our system and write (see also Fig. 3(b)). This approximation agrees with experimental observations of constant particle speeds in the flow direction (Fig. 2(c)). We also make the assumption in our model that the beads remain in the mid-plane of the depth of the channel and that the vertical component of the magnetic force on the beads is negligible,
Fig. 3 (a) A three-dimensional plot of the flow-directed gradient of the squared magnetic field, , is shown for the region of the microfluidic channel that is influenced by the applied magnetic field. The magnitude of the gradient decays away from the magnet corners, such that when (b) the speed is calculated at Y = 0, the speed can be reduced to as becomes asymptotically small. (c) A three-dimensional plot of the gradient of the squared magnetic field perpendicular to the flow, , shows that the field gradient is strongest in the region adjacent to the magnet and its strength decreases with the distance away from the magnet. |
To take into account the presence of a fluid-fluid interface, we make use of the classical result of spheres settling towards an undeformable viscous fluid-fluid interface, as described by Brenner25 and Bart.26 As a sphere approaches a viscous two-fluid interface, the Stokes drag acting on the sphere is modified by a correction factor λ,
(10) |
(11) |
(12) |
With a fitting parameter χ, we numerically solve eqn (7) and (10) for L = 0.305 and L = 0.315 and compare with our experimental results in Fig. 4(a) and 4(b). The models agree with the experimental data until the particles reach the interface. Very close to the interface, the particles' trajectories follow the results of the modified model that accounts for the interface (eqn (10)). As expected, the modified model fails to capture the dynamics of the trajectory at the fluid-fluid interface since it assumes an undeformed interface, while in reality the interface is deformed by the spheres. We find the best fit of the magnetic susceptibility χ ≈ 8 × 10−3.
Fig. 4 The approach of paramagnetic spheres towards the fluid-fluid interface for cases where (a) L = 0.305 and (b) L = 0.315. Circles represent experimental data. The dotted line represents numerical results of particle trajectory without considering the fluid-fluid interface. The solid lines are the trajectory model taking into account the two-fluid interface. The experimental values agree well with the theory. Deviations are observed very close to the fluid-fluid interface as expected since the model assumes an undeformed interface, while in reality the spheres deform and pass through the interface. The trajectories are used to calibrate for χ, the magnetic susceptibility of the 2.5 μm radius beads, and we find χ ≈ 8 × 10−3. |
Fig. 5 The interface-to-magnet distance L plotted against the SDS concentration in the aqueous phase. Solid circles (•) represent experiments where the beads experience lower magnetic force and cannot overcome the interfacial tension between the two liquid phases. Empty circles (○) represent beads that break through the interface. Scale bar 20 μm. |
For each experiment of different SDS concentration, we take the pair of data points closest to the cross-over from trapping to passage, and numerically calculate the magnetic forces (eqn (3)) in the direction perpendicular to the interface, corresponding to the average values of L. We then take a force balance on a single sphere at the two-phase interface (see Fig. 6(a)), 2πaγsinθ ≈ Fm, where θ denotes the angle of the interface, measured from the horizontal line that intersects the center of the particle.
Fig. 6 (a) A schematic diagram that illustrates the force balance on a particle at the two fluid interface. (b) The average dimensionless distance L of a pair of data points closest to the cross-over threshold is used to calculate the magnetic field gradient and estimate the interfacial tension γ viaeqn (13). γ is plotted as SDS concentration is increased from 4–20 mM. Results from microfluidic (⊕) and spinning drop (⊞) experiments are shown in the figure. The inset shows images from spinning drop experiments for different SDS concentrations in the aqueous phase. Scale bar 5 mm. |
We assume that during the interfacial deformation, the maximum value of sinθ ≈ 1, and that the interaction (if any) between the surface of the particle and the surfactants is minimal and contained within the force balance. Thus we find
(13) |
The magnetic field gradient is evaluated at Xo = 0 and 0.5 for the pair of threshold dimensionless distance L to obtain the range of interfacial tensions for each set of experimental solutions. Median interfacial tension values (⊕) are plotted in Fig. 6(b) and the error bars indicate the range. The results from spinning drop tensiometry experiments (⊞), where non-aqueous phase droplets are immersed in the aqueous phase, are also plotted in Fig. 6(b).
We find good agreement between the two techniques in the range of 4–12 mM SDS in the aqueous phase. However, the spinning drop measurements become very difficult at above 8 mM SDS in the aqueous phase as microemulsions appear in the mixture and obstruct the images of the drops (see inset of Fig. 6(b)). We are also unable to study the time dependence of the spinning drop measurements because of the tendency of the drops to shift in position over time. These challenges further convince us of the value of our magnetic microfluidic tensiometry technique.
We find that all of the solutions we test using the microfluidic tensiometer result in interfacial tensions of O(10−6–10−5) N m−1 and that the interfacial tension γ decreases as the SDS concentration in water increases. Below SDS concentration of 4 mM in water, we are unable to pull any particles through the fluid-fluid interface, which indicates that the upper limit of our interfacial tension measurements is set by how much magnetic force can be applied to the particles.
Using a mathematical model, we calculate the magnetic force that acts on paramagnetic spheres corresponding to the threshold dimensionless distance L, and estimate the interfacial tension between the co-flowing fluids by balancing magnetic and interfacial forces on the spheres. We also vary the surfactant concentration in the aqueous fluid to demonstrate that the microfluidic tensiometer can measure a range of interfacial tensions, from O(10−6–10−5) N m−1. We validate the effectiveness of our interfacial tension measurement technique by comparing our results with measurements from a spinning drop tensiometer and obtain good agreement.
The ability to measure very low interfacial tensions is challenging, but also relevant to different applications. We present a microfluidic ultralow interfacial tensiometer that simplifies and miniaturizes the measurement of interfacial tensions so that only small fluid samples are necessary. The approach we demonstrate should also allow for the interfacial tension measurement of other two-fluid systems such as phase-separated colloid-polymer solutions and all-aqueous polymer solutions.
Footnotes |
† Present address: Department of Mechanical and Industrial Engineering, Ryerson University, Toronto, ON, M5B 2K3, Canada. |
‡ Present address: Microsystems Engineering, Rochester Institute of Technology, Rochester, NY, 14623, USA. |
This journal is © The Royal Society of Chemistry 2013 |