Pinch-off of microfluidic droplets with oscillatory velocity of inner phase flow.

When one liquid is introduced into another immiscible one, it ultimately fragments due to hydrodynamic instability. In contrast to neck pinch-off without external actuation, the viscous two-fluid system subjected to an oscillatory flow demonstrates higher efficiency in breaking fluid threads. However, the underlying dynamics of this process is less well understood. Here we show that the neck-thinning rate is accelerated by the amplitude of oscillation. By simply evaluating the momentum transfer from external actuation, we derive a dimensionless pre-factor to quantify the accelerated pinch-off. Our data ascribes the acceleration to the non-negligible inner fluid inertia, which neutralizes the inner phase viscous stress that retards the pinch-off. Moreover, we characterize an equivalent neck-thinning behavior between an actuated system and its unactuated counterpart with decreased viscosity ratio. Finally, we demonstrate that oscillation is capable of modulating satellite droplet formation by shifting the pinch-off location. Our study would be useful for manipulating fluids at microscale by external forcing.

In many natural systems, coordination between dynamic and geometrical parameters plays a crucial role in retaining the reliability and efficiency for system performance. For example, a beating heart expands and contracts periodically to pump blood through vessels1, thereby maintaining one’s life; most male frogs distend vocal sacs to intensify the calls, thereby attracting females and reinforcing the success in reproduction2. In principle, the reliability and efficiency usually arise from the periodic deformation in system configuration coordinated with dynamic parameters such as pressure gradient. Technologically, human beings benefit from this coordination in developing new techniques, for instance the invention of steam engine3 by which high-pressure steam expands to perform mechanical work. Recently, scientists also use the coordination strategy to control droplet generation by acoustically actuating45 or mechanically vibrating678 fluid systems. Subjected to an oscillatory flow, the bulbous drops pulsate with the rarefaction and compression of fluid pressure, followed by a sudden breakup due to a fast drainage in the fluid neck8. The resultant micro-droplets have various potential applications, ranging from material synthesis9, chemical reactions and mixing10 to biological assays11, drug encapsulation and delivery12 and cell cultures13. In these applications, highly-controllable droplet generation with high throughput is usually required. Note that the formation of satellite droplets is unfavourable in the interest of controllable generation of droplets, and the production frequency is often limited by neck thinning velocity. It is thus essential to have deep understanding toward the dynamics of droplet breakup to meet the demand of uniformity and high throughput.

The dynamics of Newtonian threads rupture without external forcing have been well established for both liquid-in-air and two-fluid systems14. In the vicinity of pinch-off, the neck profile is self-similar and the minimum neck radius R obeys scaling laws of , with exponent α determined by the asymptotic force balances under various circumstances14; t is the critical time at neck pinch-off and t is the time. Unlike pinch-off in air, viscous effect of the surrounding liquid is non-negligible in viscous two-fluid systems15. Stokes flow dominates two-fluid pinch-off when 15, where η, ρ and γ represent respectively the dynamic viscosity, volumetric density and interfacial tension, and subscript i and o stand for inner and outer phase, respectively. In Stokes regime, capillary effect is counteracted by inner and outer viscous dissipations, while fluid inertia is negligible, which gives a linear scaling of R with α = 1161718. The neck profile is featured by the asymmetric double-cone shape. The cone slopes only depend on the viscosity ratio λ (λ = η/η) of the two-fluid system16171819, independent of other parameters, such as nozzle diameter, interfacial tension, and density difference18. If the thinning neck radius can decrease down to the molecular scales, a transition to the thermal-fluctuation dominated pinch-off occurs when R < L2021. Here is the thermal length scale, with k and T being Boltzmann’s constant and temperature, respectively. The neck profile is symmetric and the liquid-liquid interface is rough due to the effects of thermal fluctuations. Both experimental2021 and theoretical22 results suggest the power α ≈ 0.42 for thermal-fluctuation regime. However, L is usually nanometric for simple liquids, where the effects of thermal fluctuations are hardly to be observed at laboratory-scale. As such, Stokes flow is the final asymptotic regime in viscous two-fluid pinch-off, provided the molecular scales are not reached.

Except for the pinch-off of complex liquids in air2324 and zero-viscosity fluids (for example bubbles) inside another highly viscous liquids2526, the neck thinning dynamics of Newtonian fluids usually shows a universal behavior, which means that the neck profile and scaling law depends on neither initial nor boundary conditions. For instance, in liquid bridge experiments27, varying the stretching velocity can change the breakup locations28 and alter satellite droplet formation29, but the local dynamics of pinch-off remains unchanged and is identical to any pinch-off event involving that fluid282930. However, in microfluidics, the channel confinement and dynamic process (such as flow rates) complicate the pinch-off scenario and would alter the neck scaling31323334 due to the alteration in local force balance. Utilization of external forcing introduces more variables to be considered. Previous experiments show that perturbing the inner fluid pressure displays an enhancement in droplet breakup68. Assuming the minimum neck radius R still obeys the scaling law (C is the proportionality factor), the neck thinning velocity gives . Influencing either the proportionality factor C or the exponent α would induce variations in the dynamics of droplet pinch-off, and thus the neck thinning velocity V. Nevertheless, knowledge of pinch-off in microfluidics subjected to external forcing is limited, leaving several questions unaddressed like whether the external forcing would affect C or α or both and to what extent, and how it varies satellite droplet formation.

Here, we explore the dynamics of viscous two-fluid pinch-off in a co-flow microfluidic capillary device where the inner fluid velocity is oscillatory, actuated by mechanical perturbation (Fig. 1a). Our study lies into the Stokes flow regime, where the actuated neck radius still scales linearly with time prior to pinch-off (α = 1). However, the neck thinning velocity V is accelerated by the oscillation amplitude compared to that without oscillation (Fig. 1b–d). By estimating the pulsatile velocity of inner fluid at the ejection nozzle, we derive a dimensionless pre-factor to quantify the neck thinning velocity V, and rescale the perturbed pinch-off. The pulsatile velocity reinforces the transient inner fluid inertia, which, in turn, counteracts inner viscous stress, thus enhancing the pinch-off. Moreover, we find that shifts in pinch-off location, raised from the oscillation amplitude, impact the formation of satellite droplets. Our work has significance in microscale hydrodynamic systems with external actuation, such as surface-acoustic-wave5 and mechanically-perturbed678 microfluidics.

Figure 1

Neck pinch-off in microfluidics.

(a) Experimental setup. (b) Image of a drop discharged from ejection nozzle. D, D and Z are the nozzle diameter, minimum neck diameter and location of the minimum neck, respectively. Scale bar, 100 μm. (c,d) Snapshots of (c) unperturbed and (d) perturbed neck thinning versus time τ remaining to pinch-off (see supplementary Movies S1,S2,S3,S4,S5). In time period 0 < τ < 6.16 ms before pinch-off, the difference in neck diameter for perturbed case (ε= 2 mm) is larger than that for unperturbed case, showing the enhancement of pinch-off by mechanical perturbation. A suction of inner fluid towards the nozzle in (d) characterizes the last stage of neck pinch-off modulated by perturbation. The viscosity ratio of fluid system is λ = 0.157. Inner and outer fluid flow rates are Q = 0.3 mL h−1 and Q = 1 mL h−1, respectively. Scale bars, 100 μm. (e) Plot of minimum neck diameter D versus time τ for unperturbed two-fluid pinch-off with various viscosity ratios. Neck thins faster for fluid system with smaller viscosity ratio. (f) Scaling law for unperturbed neck pinch-off with different viscosity ratios. In Stokes regime, linear scaling gives . Inset: minimum neck diameter D versus time τ plotted in a log-log plane. The error comes from the limiting resolution in the captured image.

Materials and Methods

Experimental procedure

We conducted experiments in a co-flow capillary microfluidic device that was fabricated by aligning one round glass capillary with a taped nozzle (outer diameter of the nozzle D = 151.443 μm, Fig. 1a,b) inside another intact capillary (inner diameter 580 μm). Both inner and outer fluids were injected by syringe pumps (Longer Pump) into the microcapillary device, Fig. 1a. The outer fluid flow rate was kept as Q = 1 mL h−1, while inner flow rate Q was varied (0.02 mL h−1 ≤ Q ≤ 0.45 mL h−1). In the operation window of Q and Q, only dripping was observed, no jetting occurred. A mechanical vibrator (Pasco Scientific, Model SF-9324) perturbed the inner fluid microtubing (inner-wall radius R = 0.43 mm) in the gravitational direction with a displacement of . The perturbed microtubing was put nearly straight and meanwhile enabled to be displaced freely without deformation by the vibrator. The two ends of the microtubing were held still during experiments, one connected to the syringe needle and the other to the device (Fig. 1a). The distance L (Fig. 1a) between the perturbing location and nozzle was constant to be L = 40 cm. We studied the oscillation amplitude ε between 0 and 5 mm, with the upper bond limited by the vibrator, and examined frequency f ≤ 100 Hz. When frequency is higher than 100 Hz, enhancement of pinch-off decays significantly3536, probably due to two reasons: first, amplitude ε dissipates when frequency is too high; second, according to linear stability analysis14, when frequency is larger than a critical value, the cut-off wavenumber kR > 1 (k is the wavenumber and R is the unperturbed jet radius) and thus prevents the external perturbation from growing (see the Supplementary Information for validation). The flow was visualized, monitored and recorded (images and videos) by an inverted microscope (Nikon Eclipse TS100, Inverted Microscope) equipped with a high-speed camera (MotionPro® X4, IDT, Taiwan). Captured images and videos were analyzed by ImageJ.

Water-in-oil two-phase flow was examined in experiments, where outer fluid was silicone oil (viscosity: η = 492.875 mPa s), and inner fluids were various glycerol-water mixtures (viscosity: 0.992 mPa s ≤ η≤ 831.07 mPa s). This choice of two-fluid systems enables a wide range of viscosity ratios (10−3 < λ < 10), but narrow variations in interfacial tension (24.2 mN m−1 ≤ γ ≤ 32 mN m−1) and fluid density (0.998 g cm−3 ≤ ρ ≤ 1.261 g cm−3), so as to isolate the effect of viscosity ratio on pinch-off dynamics. Viscosity was measured by a viscometer (microVISCTM, RheoSense, Inc.), interfacial tension by a ring tensiometer (Surface Tensiometer 20, Cole-Parmer), and density by quantifying the volume of a known-mass fluid.

Estimate of the inner fluid oscillatory velocity V

With perturbation, the inner phase pressure is disturbed periodically, so that the velocity V of the inner fluid ejected from the nozzle pulsates with time. We have made a simple analysis to determine the pulsatile velocity V. Before accounting for the effect of mechanical perturbation, we first consider the unperturbed case, in which the inner phase flow rate Q is estimated by Poiseuille’s law37,

where R is inner-wall radius of the microtubing (R = 0.43 mm), η is the viscosity of inner fluid, P and P are pressures respectively at the perturbing location and the nozzle, and L is the length between the two locations. Based on Q, we evaluate the steady ejection velocity of the inner fluid as,

For perturbed case, we simply assume that the gravitational displacement of modifies the steady pressure P into the time-dependent P(t),

with inner fluid density being ρ and the gravitational acceleration being g. Since the linear pinching time-scale t (t ~ 1 ms, which is shown in the next section, Fig. 2) is much smaller than the perturbation period t (t≥ 10 ms for f ≤ 100 Hz in experiments), we assume that the Poiseuille’s law is still valid to estimate the fluctuating inner flow rate Q(t) in the linear pinch-off regime. When P is replaced by P(t) in Eq. (1), we obtain Q(t) as a function of time t,

Figure 2

Acceleration of pinch-off with oscillatory velocity.

(a,b) Effects of oscillation amplitude (a) and frequency (b) on pinch-off. Compared with unperturbed case, neck thins faster in the perturbed situation, but still displays linear thinning dynamics within the last ~1 ms (indicated by the dashed line) before pinch-off. The acceleration of pinch-off increases with oscillation amplitude in (a). By contrast, for different frequencies in (b), the enhancement of pinch-off shows weak dependence on the perturbation frequency. (c,d) Dynamics of perturbed neck thinning versus rescaled time for various oscillation amplitudes (c) and frequencies (d). Time τ is rescaled by multiplying a dimensionless pre-factor . In (a,c), Q = 0.45 mL h−1, Q = 1 mL h−1; in (b,d), Q = 0.3 mL h−1, Q = 1 mL h−1. Viscosity ratio, λ = 0.0387.

Now the oscillatory velocity V discharged from the nozzle is approximated to be , and hence in the form of

where is the dimensionless amplitude of the fluctuation term.

Although we ignore several effects that might influence the transient velocity V, such as shaking and deformation of the microtubing, and fluid acceleration in the microtubing, the simplified form of Q(t) in Eq. (4) (as well as V in Eq. (5)) does characterize essentially how oscillation amplitude ε, in association with material properties (ρ and η) and geometrical parameters (R and L), modulates the ejection velocity. We verify Eq. (4) by comparing the experimental measurements to theoretical predictions of the inner fluid discharged volume (see supplementary Figs S1 and S2). Moreover, we will show in the next section that pinch-off is enhanced by the vibration in the gravitational direction, while no enhancement is identified when the vibration is parallel to the horizontal plane (see supplementary Fig. S3). This difference suggests the dominant effect of gravity in impacting neck pinch-off. As a result, the above analysis paves a way to evaluate the average momentum transferred into the fluid system, to quantify the acceleration of pinch-off with perturbation, to rescale the perturbed neck pinch-off, and to account for the non-negligible inner fluid inertia in perturbed case, which are shown in the next section.

Results and Discussion

Scaling of classical viscous two-fluid pinch-off

We start from the case without oscillation but with the influence of viscosity ratio λ. Keeping silicone oil as the outer fluid, and choosing various glycerol-water mixtures as the inner fluid, we can tune viscosity ratio λ over three magnitudes (10−3 < λ < 10). Fluid system with smaller λ means less viscous inner fluid in the context, and vice versa. Figure 1e compares unperturbed minimum neck diameter D (Fig. 1b, D = 2R) versus time τ (τ = t − t, time remaining to neck pinch-off) for different viscosity ratios λ. Fluid system with smaller viscosity ratio pinches off faster (Fig. 1e) due to the less inner viscous resistance. In the vicinity of pinch-off, a linear scaling of the neck diameter is presented (inset of Fig. 1f). In this linear pinch-off scenario, two-fluid Stokes flow dominates, where inner axial viscous stress (with the kinematic estimate being , z the length scale of the pinching neck and u the axial velocity of fluid flowing out of the neck), outer radial viscous stress , and capillary pressure compete with each other15. This viscous-capillary force balance yields the scaling of neck diameter as (see the Supplementary Information for derivations), where is viscous-capillary velocity, and F( is a dimensionless proportionality factor that can be evaluated from linear stability analysis16171838. Normalizing D by the nozzle diameter D, we obtain a linear scaling of the minimum neck diameter in the form of

with viscous time scale . Based on the analysis by Tomotika38, we evaluate F( for different viscosity ratios (see the Supplementary Information for the estimate of F(λ)), and replot the data in Fig. 1f, which is in good agreement with previous study of viscous two-fluid drop pinch-off15161718.

Perturbation accelerated pinch-off

Now we consider the fluid system with fixed viscosity ratio (λ = 0.0387) under periodic perturbation. When undergoing pinch-off subjected to an oscillatory flow, the droplet initially grows in its size until the neck appears. Afterwards, a sudden suction of the inner fluid, caused by the pulsatile velocity V (Eq. (5)), triggers and subsequently accelerates the neck pinch-off (exampled in Fig. 1d, and see supplementary Movies S1,S2,S3,S4,S5). The above process is identical to any individual pinch-off with oscillation (see supplementary Movie S6)). As a result, actuated neck thinning velocity V is faster than the case without actuation, while the linear neck thinning regime features both cases in the vicinity of pinch-off (Fig. 2a,b). The neck thinning velocity V increases apparently with the oscillation amplitude ε (Fig. 2a). However, the pinching velocity weakly depends on the oscillation frequency for f ≤ 100 Hz (Fig. 2b), probably due to the mismatch between the two time scales, the linear pinching time scale t (t ~ 1 ms; see dashed line in Fig. 2a,b) and the perturbation period t (t≥ 10 ms for f ≤ 100 Hz).

We now quantify the acceleration of perturbed pinch-off by estimating the extra momentum contribution from mechanical vibration. After simply calculating the root mean square of the fluctuation term in Eq. (5), , we extrapolate that the average momentum is increased by a ratio of due to mechanical vibration, and hypothesize that the neck thinning velocity is accelerated by a ratio of as an estimate. To verify this hypothesis, we replot Fig. 2a,b by rescaling D versus the product of time τ and the dimensionless factor . After rescaling, all data essentially collapse onto a single master curve in the vicinity of pinch-off (Fig. 2c,d). The enhancement of pinch-off comes thus from the additional momentum transportation induced by mechanical vibration, and is proportional to the dimensionless amplitude ε.

Rescaling perturbed neck thinning with various viscosity ratios

In subsequent tests, we focus on the influence of viscosity ratio on fluid systems with oscillation. In comparison to unactuated case, the oscillation-enhanced pinch-off has a faster neck thinning rate for all fluid systems examined, as shown in Fig. 3a. Note that Fig. 2b only indicates a dependence of neck thinning velocity on oscillation amplitude ε, but Fig. 3a further clarifies the dependence on the dimensionless amplitude ε, in that larger enhancement of pinch-off is observed for fluid system with the same ε but smaller viscosity ratio λ (smaller η in experiments). Therefore, Eq. (6) fails in scaling the perturbed neck thinning with arbitrary oscillation amplitudes and viscosity ratios. From previous discussions, neck thinning velocity is accelerated by a ratio of . We then rescale the neck shrinkage behavior via multiplying the right hand side of Eq. (6) by the dimensionless pre-factor :

Figure 3

Influence of oscillatory flow on two-fluid pinch-off with various viscosity ratios.

(a) Minimum neck diameter versus time for perturbed and unperturbed pinch-offs with various viscosity ratios. (b) Rescaling the perturbed neck pinch-off in linear pinching regime, . Inset: representative data in the linear scaling regime. (c) Comparison of inner fluid radial Reynolds numbers for perturbed and unperturbed cases. Perturbed Reynolds number Re (circle) increases sharply during the last stage of pinch-off, and thus surpasses unperturbed Re (square), indicating an enhancement of inner fluid inertia by perturbation. (d) An example of equivalent neck thinning behavior. Perturbed system (λ = 0.157, filled triangle) displays almost the same pinch-off dynamics as an unperturbed system with less viscous inner fluid (λ = 0.1305, circle) in the linear pinch-off regime. (e) Equivalent viscosity ratio λ as a function of dimensionless amplitude ε and viscosity ratio λ of the perturbed system. For the calculated domain, 0 ≤ ε ≤ 1, 10−2 ≤ λ ≤ 102, and the resultant equivalent viscosity 10−4 < λ < 102.

as shown in Fig. 3b (data rescaled from Fig. 3a). To highlight this rescaling, we replot representative data close to pinch-off in the inset of Fig. 3b, where all data falls onto a single master straight line with unit-slope, in excellent agreement with Eq. (7).

Non-negligible inner fluid inertia

Having quantified the accelerated pinch-off, we now explore the origin of the acceleration in pinch-off with oscillatory flow. For conventional pinch-off scenario, two-fluid Stokes flow dominates near the singularity. However, oscillation modifies the force balance via disturbing the flow velocity, therefore affecting the whole pinch-off dynamics. A recent study proposed that, when pinching in air, fluid thread passes through a plethora of transient regimes before transferring into the final inertial-viscous regime39. Likewise, various regime transitions also exist in two-fluid pinch-off scenario before the final Stokes regime arrives, as indicated by the deviation of slope from unity when τ > 1 ms (Fig. 3a), especially for cases with high viscosity ratios (for example λ = 1.6862, open square in Fig. 3a). Oscillatory flow, however, smooths the original slope-varying curve into a more uniform unit-slope form (Fig. 3a, filled square for λ = 1.6862 with oscillation) by impacting the transient local force balance.

We hypothesize that oscillation amplifies the effects of axial inner fluid inertia. During pinch-off in Stokes regime, inner fluid thread of radius R and length z thins radially with velocity V, while fluid flows axially out of the thread with velocity u. Local volume-flux conservation requires . Provided 15, u is thus estimated as . Since neck thinning velocity V is accelerated by the oscillation amplitude, u increases with ε, so does the inner axial fluid inertia . To verify this hypothesis, we compare the local axial Reynolds number Re in the neck for perturbed and unperturbed cases, respectively. Having and (), we assume that the axial is of the same magnitude of the radial . Figure 3c displays the results of radial Re versus time τ during neck pinching (λ = 0.0387). Close to pinch-off, a sharp increase of Re occurs for the perturbed neck, while undisturbed Re varies more gently. The increased Re indicates an undoubted increase of inner fluid inertia due to perturbation. It is now probably invalid to ignore the effects of oscillation-enhanced fluid inertia in determining the neck thinning dynamics. Actually, the increased fluid inertia neutralizes the inner fluid viscous dissipation.

We now elucidate the counterbalance between the enhanced inner fluid inertia and inner viscous stress. Attributed to oscillation-enhanced inner fluid inertia, the effective inner dissipation is weakened roughly to be (compared to unperturbed case of ). Because of the counterbalance, the perturbed pinch-off displays an equivalent neck thinning dynamics to the unperturbed case with smaller inner viscosity η (smaller viscosity ratio λ too). Preliminary, an observation (Fig. 3a) supports this speculation; for example, fluid system with viscosity ratio λ = 0.6643 disturbed by ε0 = 3 mm (filled circle, Fig. 3a) has the similar neck thinning dynamics to the unperturbed “less viscous” system with λ = 0.157 (open triangle, Fig. 3a).

To quantify this “less viscous” performance due to perturbation, we explore the equivalent viscosity ratio λ corresponding to a perturbed system with viscosity ratio λ (λ < λ). This can be done by calculating from Eq. (7). Take system with λ = 0.157 perturbed by ε0 = 3 mm for example, we find out λ1 = 0.1305. Experiments show that, close to pinch-off, the perturbed system with λ = 0.157 pinches equivalently to the undisturbed one with λ = 0.1305 (Fig. 3d). Furthermore, an equivalence between λ and λ with arbitrary ε is developed (Fig. 3e) based on the formula: . According to above analysis, we can even modulate a glycerol-water mixture undergoing “less viscous” pinch-off dynamics than distilled water by exerting appropriate oscillation amplitude (Fig. 4). The equivalence of pinch-off dynamics indicates potential in breaking highly viscous liquid in an easier manner by applying an oscillatory flow. Our findings show that dynamics of neck thinning in the linear regime is controlled by both fluid properties and oscillation amplitude. By characterizing the combined effects of viscosity ratio and dimensionless amplitude (Eq. (7) and Fig. 3e), we can predict the dynamics of linear neck pinch-off.

Figure 4

Dynamics of distilled water (crossed diamond) and 47 wt.% glycerol + 53 wt.% water (star) pinch-off in silicone oil.

Although glycerol-water mixture pinches slower (open star) than pure water (crossed diamond), perturbing the mixture by appropriate amplitude (filled star) have the mixture pinch faster than water.

Shifts of pinch-off location with perturbation

Besides neck thinning dynamics, pinch-off location also changes with applied oscillation. When breaking up, the liquid thread that connects the main drop and the residual fluid can pinch off either at its front or rear side4041 (Fig. 5a, and see arrows in Fig. 5b–f). Compared with unactuated case (Fig. 5b), applied oscillation shifts pinch-off location upstream towards the nozzle (Fig. 5c–f, and see supplementary Movies S1,S2,S3,S4,S5) due to the pulsatile velocity V. More remarkably, while front pinch-off occurs before rear pinch-off for unperturbed case (Fig. 5b), large driving amplitude ε can introduce the reversed case (Fig. 5d–f), between which a symmetric neck forms with front and rear pinch-offs occurring simultaneously (Fig. 5c). Meanwhile, oscillation can also induce a transition of pinch-off location from outside the nozzle (Fig. 5b–d) to inside the nozzle (Fig. 5f).

Figure 5

Shifts of pinch-off location with oscillation.

(a) Schematic showing two possible pinch-off locations: front pinch-off and rear pinch-off. (b–f) Snapshots of pinch-off location transition and satellite droplet formation. In (b–f), the first and second rows display the first and second neck pinch-off (arrows), respectively; the third row shows an instant for satellite droplet formation; the last row presents the instant when all satellites are formed. The less visible satellites are indicated by circles and the numbers below. Time value on snapshots indicates the time after the first pinch-off. (b) ε0 = 0, front pinch-off occurring before rear pinch-off. (c) ε0 = 1 mm, symmetric neck formation. (d–f) Rear pinch-off occurring before front pinch-off. (d) ε0 = 2 mm, pinch-off outside the nozzle. (e) ε0 = 3 mm, pinch-off at the nozzle. (f) ε0 = 5 mm, pinch-off inside the nozzle. (g) Evolution of the location of minimum neck Z versus time (t − t) during the pinch-off. t is the critical time for the first pinch-off. Data are measured from cases in (b–f). (h) Transition of pinch-off location with dimensionless amplitude ε. Regions (b,d,f) are separated by boundaries (c,e). ε = 0.110 for the symmetric neck formation (c); εc = 0.225 for pinch-off at the nozzle (e).

We map the pinch-off locations in Fig. 5g by tracing the location of minimum neck Z (Fig. 1b) versus time (t − t) before pinch-off. Two boundaries are emphasized in Fig. 5g: symmetric neck formation (Z bifurcates, pentagon) and pinch-off at the nozzle (Z = 0 when t = t, circle). The transition of pinch-off location is the result of momentum transport. Because axial velocity u of inner fluid scales as , it is reasonable to assume that the two boundaries correspond to certain critical values of ε. Experimentally, we find that, in the tested range of viscosity ratios, ε= 0.110 and ε = 0.225 (Fig. 5h) are the two critical values for the symmetric neck formation and pinch-off at the nozzle, respectively. These results provide a method to precisely manipulate pinch-off location by using oscillatory velocity.

Being susceptible to the details of breakup4041, satellite droplet formation varies with the shift of pinch-off location. The multi-breakup of liquid thread generates multiple satellite and subsatellite droplets (circles in Fig. 5b–f) owing to the self-repeated neck formation4243. With increased oscillation amplitude, the number of thread fragmentation rises, inducing more satellites and subsatellites (Fig. 5b–f); for instance, as many as ten satellites are generated with ε= 5 mm (Fig. 5f). In addition, the size distribution of satellite drops relies on the symmetry of neck formation. Symmetric pinch-off induces symmetric size distribution (Fig. 5c); asymmetrical neck formation otherwise renders the size distribution asymmetric (Fig. 5b,d,e). Astonishingly, rear pinch-off inside the nozzle produces satellite drops with descending size distribution (Fig. 5f), which are similar to those generated by the tip-multi-breaking, a recently reported droplet breakup mode4445. Despite their non-uniformity, the satellites depend on oscillation, which may provide an additional handle to harness the formation of droplets with various volumes, a required feature in some applications. For example, in multi-volume droplet digital polymerase chain reaction (MV-dPCR)46, droplets with various volumes enable simultaneous measurements of a sample at different copies per droplet. Compared to single-volume digital PCR, MV-dPCR achieves higher detection reproducibility, wider dynamic range and better resolution while reducing the total number of droplets/wells required for the measurements464748.

Concluding Remarks

In conclusion, we have experimentally examined the dynamics of pinch-off in viscous two-fluid systems with oscillatory velocity actuated by mechanical perturbation. Attributed to the oscillatory flow, an enhanced suction of inner fluid towards the nozzle occurs in the last stage of pinch-off. In this scenario, the inner fluid thread thins radially faster prior to breakup compared to the unactuated case, still in a linear way though. The enhancement of pinch-off by external actuation depends on the oscillation amplitude. We rescale the actuated neck radius by a dimensionless pre-factor , where ε is the dimensionless amplitude. Actuating the fluid system modulates the local force balance via increasing the effects of inner fluid inertia. Therefore, the enhanced fluid inertia is non-negligible and responsible for the accelerated neck thinning. Meanwhile, the actuated pinch-off displays a “less viscous” behavior due to the counterbalance between the enhanced inertia and the inner viscous resistance. Such a “less viscous” performance would enable robust control over breaking viscous liquid filament more easily. Moreover, the enhanced inner fluid inertia shifts the pinch-off location upstream towards the nozzle, which afterwards affects satellite droplet formation. By quantifying the relationship between the oscillation amplitude and the pinch-off location, novel control over satellites formation by oscillation could probably be developed.

Beyond technical benefit in modulating droplet formation resulted from our study, this work raises several questions of hydrodynamic interest but remaining to be addressed. Among them, two issues are of the most importance. One is that if the oscillation effect could dominate capillary effects, a new regime might occur, which may be similar to the thermal fluctuation dominant regime2021. Another is that the existence of transient regimes during two-fluid pinch-off from the initial to final regime remains unexplored. Deeper understanding towards these issues calls for experimental studies in association with numerical simulations and theoretical explanations.

Additional Information

How to cite this article: Zhu, P. et al. Pinch-off of microfluidic droplets with oscillatory velocity of inner phase flow. Sci. Rep. 6, 31436; doi: 10.1038/srep31436 (2016).