The effects of hydrophilic surfactant concentration and flow ratio on dynamic wetting in a T-junction microfluidic device
Chris P. Tostado,Jianhong Xu,Guangsheng Luo
DOI: https://doi.org/10.1016/j.cej.2011.05.043
IF: 15.1
2011-01-01
Chemical Engineering Journal
Abstract:Highlights ► The contact angle of an oil drop on a PMMA surface in water ranges depending on the CTAB surfactant concentration. ► The PMMA surface transitions from being oleophilic (<<90°) to oleophobic near a CTAB concentration of 1.0 × 10 −2 % w/w. ► Microchannel flows showed different wetting patterns at different surfactant concentrations due to CTAB adsorption. ► CTAB alters the surface energy of PMMA, changing the wetting state and allowing for the production of O/W and W/O flows. ► CTAB alone (no binding) cannot lower the water surface tension and interfacial tension enough to induce a phase inversion. Abstract An initial study was conducted on the influence of the hydrophilic surfactant concentration and two-phase flow-rate ratio on the aqueous-phase wetting ability (dynamic contact angle) in a T-junction PMMA microfluidic device. Cetyltrimethyl ammonium bromide (CTAB) was used as the hydrophilic surfactant for the deionized water aqueous phase while the organic phase was corn oil. Static contact angles and interfacial tensions were measured for systems with surfactant concentrations ranging from 0% to 5.0 × 10 −1 % (w/w). Static contacts angles for an oil drop on PMMA in water show a transition from <90° to >90° at CTAB concentrations near 1.0 × 10 −2 % (w/w), suggesting a change in the surface energy and a transition from a slightly hydrophobic to hydrophilic PMMA surface for the given systems. Microchannel flow experiments were conducted to study dynamic wetting character, and special attention was paid to the CTAB concentration range between 3.0 × 10 −3 % (w/w) and 3.0 × 10 −2 % (w/w) to focus on the region in which a wetting transition was observed in the static case. Dynamic contact angles and wetting responded more to changes in the aqueous phase flow rate than to changes in the organic (continuous) phase flow rate, and receding contact angles (plug tails) changed from >90° to <90° at a W:O flow ratio near unity. Different dynamic wetting patterns were produced at different CTAB surfactant concentrations because of the dynamic adsorption of hydrophilic surfactant at the channel wall interface. At concentrations below the critical micelle concentration, CTAB can be used to dynamically alter the surface energy of a PMMA microchannel and change its wetting state allowing for the production of both O/W and W/O type flows. The wetting state (i.e. the extent to which each phase wets the channel wall) is crucial for numerous lab on a chip applications. Keywords T-junction microchannel CTAB Phase inversion Dynamic wetting Contact angle 1 Introduction The production of two-phase liquid–liquid drop and plug flow in microchannels has become an important hallmark of microfluidic research because of its robust application to different industries including the biotechnology, pharmaceutical and chemical industries [1–5] . Liquid–liquid two-phase flows also provide unique opportunities for studying novel fluid dynamic and interfacial phenomena [6–8] . Numerous studies have been conducted in the past on different aspects of two-phase flows of immiscible liquids in various microfluidic channel geometries including the T-shaped, cross-flow, flow focusing and head-on junctions [9–14] . However, an in-depth understanding of not only the methods for controlling the two-phase flow patterns produced in microchannels, but also the ways in which these flows and their wetting characteristics can be manipulated, is crucial for the advancement of microfluidic technologies which rely on transfer phenomena and chemical reactions [15] . Although the study of two-phase flow instabilities and phenomena is not new [16] , one particular topic – surfactant-controlled drop and plug flow – has been the subject of numerous studies and has also now been applied to numerous microfluidic systems. In these systems, the Bond number and Reynolds number are both usually ≪1. The values of the Bo and Re signify the ratios between gravitational/interfacial and inertial/viscous forces respectively. Thus gravitational and inertial effects are insignificant in comparison to interfacial and viscous forces in microfluidic systems. Furthermore, because of the relative importance of interfacial phenomena, a recurrent theme within numerous studies involving two-phase flow is the central role which fluid wetting plays in determining the flow regime [10,17] . In an attempt to better understand and explain the role that wetting plays, many studies including those mentioned above have used interfacial tension and static contact angle data as a starting point for determining the wetting characteristics observed for a particular system. Both parameters are related to the surface energy by the Young equation: (1) γ L ( cos θ ) = γ SV − γ SL where γ L is the interfacial tension of the liquid in the ambient fluid (e.g. air), γ SV is the surface tension of the solid in equilibrium with the liquid vapor in the ambient air, γ SL is the interfacial tension of the solid–liquid interface and θ is the contact angle [18] . The interfacial tension, γ is measured in units of force per unit distance (N/m) or units of energy per unit area and at constant temperature, volume, and mass, is thermodynamically defined by: (2) γ = ∂ F ∂ A T , V , n where F is the free energy of the system, and A is the interfacial area [19] . Surfactant-controlled drop and plug formations and deformations have been studied in the references. Dreyfus et al. [6] reported on the wetting characteristics and drop formation of deionized water/tetradecane-Span 80 systems in both T-junction and cross-flow microfluidic devices. At various oil to water flow ratios, they observed different flow regimes similar to those obtained in previous studies [7] . They concluded that the wetting properties of the fluids are mainly influenced by the presence and phase location of the surfactant, the entry conditions of the fluids, and the dimensions of the channels, that determine the flow regime and channel wetting characteristics obtained within a microfluidic device. These results were confirmed by Xu et al. [9,10] with different experimental setups, emphasizing the importance of the surfactant in determining the flow and channel wetting (e.g. O/W or W/O) within the same microfluidic device. In most of these studies, surfactants have been added in high amounts (∼2.0% w/w) to stabilize the microfluidic flows, preventing or reducing the coalescence of drops and plugs. Directly measured interfacial tension data has often been combined with static contact angle data to determine wetting regimes in which the walls of the microchannel exhibit hydrophobic or hydrophilic character or even a combination of the two [20] . Fig. 1 shows the change in liquid–liquid interfacial tension and static contact angle over a wide CTAB surfactant concentration range. The system – corn oil and deionized water with different concentrations of the hydrophilic surfactant CTAB – shows the interfacial tension decreases and the contact angle (which an oil droplet makes on a poly(methyl methacrylate) surface in water) increases as the concentration of CTAB is raised. These two inverse trends, consistent with previous studies, are particularly interesting in the surfactant concentration range between 5.0 × 10 −3 % (w/w) and 3.0 × 10 −2 % (w/w) because a transition takes place in which the wetting character between the organic phase and solid changes from lipophilic (contact angle <90°) to lipophobic (contact angle >90°) ( Fig. 1 inset). A potential dynamic state analog to this transition is the phase inversion, a phenomenon commonly observed in microchannels in which the dispersed and continuous phases of flow are reversed, e.g. from oil-in-water to water-in-oil flow. [10] In liquid–liquid flows, oil-in-water to water-in-oil (and vice versa) phase inversions have recently become the focus of attention because of their potential capability of enhancing and controlling mixing processes as well as other industrial applications [15,21] . However, much is still unknown regarding the energetics and conditions surrounding phase inversions and several views exist regarding the extent to which liquid–liquid interfacial tensions and liquid–solid interfacial tensions play a role in causing them. In this initial study, we report on the plug formations and wetting character of the two-phase flows produced within a microchannel for a particular liquid–liquid system, which has not been previously studied, with special emphasis on the surfactant concentration range in which a static wetting transition was observed. One of our objectives in conducting dynamic flow experiments within a microchannel was to determine if any dynamic transition occurs analogously to the one observed in the static case. Drop and plug formations were closely observed at different concentrations within the contact angle transition region and special attention was paid to the receding and advancing contact angles formed by the aqueous phase plugs. Finally, in an attempt to explain the role of surfactant-controlled liquid–liquid interfacial tension on two-phase flow in microchannels, an equation of state approach was used to estimate the interfacial energies of various plug formations. 2 Experimental 2.1 Microfluidic device A simple T-junction microfluidic device was used for all flow experiments and the channel geometry and dimensions are shown in Fig. 2 . Two channels were constructed on a PMMA chip using a computer numerically controlled (CNC) end-mill device. The main channel dimensions were 200 μm × 300 μm for the depth and width, and the channel's total length was 45 mm. The organic phase passed through the main channel while the aqueous phase flowed through the side channel, which intersected with the main channel at 10 mm and possessed a depth, height and length of 50 μm and 200 μm and 12 mm, respectively. Both phases were pumped into the microchannels using two separate micro-syringe pumps. 2.2 Materials Corn oil (with viscosity of ∼67.0 mPa s at 20 °C) and deionized water (with viscosity of ∼1.0 mPa s at 20 °C) were used for the organic and aqueous phases. Different concentrations of hydrophilic surfactant cetyltrimethyl ammonium bromide (CTAB) were dissolved into the aqueous phase. 2.3 Experimental design and analysis Viscosities of the organic and aqueous phases were measured with a viscometer at 20 °C. Viscosities of the aqueous phase solutions did not vary significantly (<1.00%) with the addition of CTAB. The liquid–liquid interfacial tensions and static contact angles were measured using a DataPhysics OCA20 device with a measurement range of 0.01–2000 mN/m and an error of ±2%. The measured interfacial tension between the oil and water phase without any surfactant was 24.0 mN/m. This value is consistent for ordinary corn oil (∼98% triglycerides), which has not been purified further [22] . The pure corn oil and deionized water surface tensions at 20 °C were 33.28 and 71.85 mN/m, respectively, also consistent with the literature [23,24] . An inverted sessile drop method was used to show the contact angle a small drop (<0.1 μl) of corn oil makes on a poly(methyl methacrylate) (PMMA) surface in the presence of water with different concentrations of hydrophilic CTAB ranging from 0% to 5.0 × 10 −1 % (w/w). Flow experiments within the T-junction microfluidic device were performed for the same concentration range with special emphasis on the region between 3 × 10 −3 % (w/w) and 3 × 10 −2 % (w/w) to focus on the transition of the wettability of PMMA. W:O flow rate ratios ranged from 0.1 to 10. Microchannel flow experiments were executed and monitored in a manner similar to previous studies performed by Xu et al. [9,10] . The photos of two-phase flow patterns and plug formations were taken at 95 frames/s with a 0.04 ms exposure time using a high-speed CCD camera with a 1280 × 1024-pixel resolution attached to an Olympus U-TV0-5xC-3 microscope at 10× magnification. Additionally, the advancing and receding dynamic contact angles of aqueous plugs at high CTAB concentrations (>1 × 10 −2 ) were measured using DataPhysics SCA20 contact angle measurement software Ver. 3.12.11. 3 Results and discussion 3.1 Surfactant concentration, capillary number and plug length As shown in Fig. 1 , the surfactant concentration range of interest lies between 5.0 × 10 −3 % and 3.0 × 10 −2 % (w/w). At higher (>3.0 × 10 −2 % w/w) concentrations of CTAB, the interfacial tension approaches a minimum value of approximately 0.6 mN/m and does not change at higher concentrations past the critical micelle concentration (CMC) of ∼4.0 × 10 −2 % (w/w). The contact angle the oil drop makes on the PMMA surface increases with surfactant concentration, reaching a maximum value of approximately 163° at concentrations above 3.0 × 10 −2 % (w/w). CTAB concentrations of interest were categorized as lower or higher with a transition concentration at 1.0 × 10 −2 % (w/w). (1) At lower CTAB concentrations below 1.0 × 10 −2 % (w/w), capillary numbers, Ca = μν / γ , where μ and ν are the viscosity and flow rate velocity of the continuous phase, respectively and γ is the liquid–liquid interfacial tension, ranging roughly from 2 × 10 −3 to 2 × 10 −2 , are lower because of higher liquid–liquid interfacial tension values and two-phase flow and wetting character is largely influenced by interfacial forces. (2) At transition concentration of 1.0 × 10 −2 % (w/w), Ca values are ranging from 3 × 10 −3 to 3 × 10 −2 depending on the oil phase flow rate, which ranged from 1 to 10 μl/min. (3) At higher concentrations of CTAB above 1.0 × 10 −2 % (w/w), Ca values are above 1 × 10 −2 . In this region, liquid–liquid interfacial tension decreases (capillary numbers grow larger) and viscous forces influence the plugs flow and the dynamic wetting character much more noticeably. Fig. 3 shows the effects of surfactant concentration and capillary number on the aqueous plug length. The aqueous phase flow rates are fixed at 1, 5 and 10 μl/min, with oil phase flow rates varying from 1 to 50, 1 to 200, and 1 to 400 μl/min, respectively. For all three aqueous phase flow rates, plug length varies dramatically at Capillary numbers below 2.0 × 10 −2 . Above a Capillary number of 2 × 10 −2 , plug lengths become much more uniform and are not greatly affected by the continuous flow rate or the surfactant concentration. The Ca number at which this transition from largely variable to more uniform plug lengths agrees with the model proposed by De Menech et al. [25] . A Ca value around 2.0 × 10 −2 is thought to be the point at which a transition occurs from interfacial dominated flow in which flow rates primarily dictate plug sizes, to shear-viscous dominated flow in which the viscosity ratio between phases dominates. The effects of surfactant concentration and capillary number on two-phase flow patterns and dynamic wetting properties were discussed at different surfactant concentration regions respectively in the following sections. 3.2 Lower CTAB concentrations—partial wetting flow At CTAB concentrations between 0% and 5.0 × 10 −3 % (w/w), well below the CMC, water plugs of varying sizes are formed. Tables 1–3 show the flow types formed at different CTAB concentrations, with low W:O flow ratios ranging from 0.1 to 1.0. Plug sizes decrease with the increase of capillary number. At concentration as low as 3.0 × 10 −3 % (w/w), a partial wetting plug flow is observed in which the aqueous phase plugs cling sporadically to the PMMA walls of the microchannel ( Fig. 4 ). The phenomenon is observed even more as the concentration of surfactant increases. Although the water plugs could cling to different walls of the microchannel, one side of the microchannel tended to be favored over the others. This is most likely due to an uneven distribution of surfactant molecules on the PMMA surface. As the CTAB concentration increases to 5.0 × 10 −3 % (w/w), the water plug profiles become less stable and water tails begin to form, which also wet the PMMA walls. Flow instability reaches a climax at a CTAB concentration of 5.0 × 10 −3 % (w/w) and a water flow rate of 10 μl/min in our experiments. When Ca ≥ 0.02, no partial wetting by the aqueous phase was observed and small plugs or water droplets were produced. Changes to the W:O flow ratio did not lead to significant changes in wetting character at the same surfactant concentration. For lower CTAB concentrations, the channel wetting character and water-in-oil type flows produced resemble the corresponding static contact angle data. The existence of partial wetting can be attributed to the fact that PMMA wall is neither a strongly hydrophobic or hydrophilic surface. The increase in the tendency of the plugs to temporarily stick to the PMMA walls as surfactant concentration increases. This tendency implies that a change takes place at a solid–liquid or the liquid–liquid interface (or some combination) sufficient enough to alter the interfacial energy and affect wetting character. 3.3 Transition CTAB concentration—unstable flow patterns Plug formations are noticeably less stable at the transition CTAB concentration of 1.0 × 10 −2 % (w/w). Aqueous plugs, in addition to exhibiting the partial wetting type character observed at lower CTAB concentrations, also possess prominent trailing tails, which significantly obfuscate the back end profile. Typical flow patterns and flow rates are shown in Table 4 . Instability increases with the increase in aqueous phase flow rate and reaches an apparent climax at an aqueous flow rate of 10 μl/min when the W:O flow ratio was 10:2 μl/min and the capillary number was 6.4 × 10 −3 . Here a spontaneous phase inversion occurred. Although a phase inversion was not observed for any other experiments with different flow conditions, experiments with conditions closely resembling those at which the inversion occurred exhibited flow patterns, which indicated an impending inversion. Fig. 5 shows images of different flow patterns with a higher W:O flow ratio of 10:1, and lower Ca of 3.2. × 10 −3 ( Fig. 5 A) as well as a lower W:O flow ratio of 10:3 and higher Ca of 9.6 × 10 −3 , ( Fig. 5 B). Instability preceded the stable phase inversion with long aqueous plugs coalescing and blocking (building up pressure) the flow of oil at the T-junction ( Fig. 5 A), while immediately after the phase inversion at a slightly lower W:O flow ratio and higher capillary number, the flow pattern returns to an unstable state and neither the oil phase nor the aqueous phase completely wets the channel walls. Videos of these flow patterns have been provided Similarly to the low CTAB concentration regime, the transition CTAB concentration regime showed that different dynamic wetting patterns were produced at different CTAB surfactant concentrations due to surfactant binding at the channel wall. At concentrations below the critical micelle concentration, CTAB can be used to dynamically alter the surface energy of a PMMA microchannel and change its wetting state allowing for the production of both O/W and W/O type flows. The wetting state (i.e. the extent to which each phase wets the channel wall) is crucial for numerous lab-on-a-chip applications including wall-coated catalysis [26] and plug-based microreactions [27] . 3.4 Higher CTAB concentrations—dynamic wetting and shear driven instabilities Beyond the transition CTAB concentration of 1.0 × 10 −2 % (w/w), flow patterns begin to stabilize and reflect the increasing contribution of viscous forces to plugs/drops formation. At these concentrations, aqueous phase plugs exhibit dynamic advancing and receding contact angles which were measured with OCA 20 contact angle software to infer wetting states ( Fig. 6 ). The different flow types encountered at higher CTAB concentrations can be divided into four wetting regimes: stratified wetting, dynamic wetting plug-flow, non-wetting plug-flow, and shear deformed plug flow, as shown in Fig. 7 . In all these flow regimes, the aqueous phase plug head shape and the advancing contact angle changed very little (advancing contact angle ∼170°). The main difference is the change in plug tail shape and the receding contact angle. The first three wetting categories were encountered at CTAB concentrations slightly above the transition concentration at 2.0 × 10 −2 % (w/w) and 3.0 × 10 −2 % (w/w). Table 5 shows the images of typical flow patterns for different flow ratios and Ca values observed at CTAB concentration 2.0 × 10 −2 % (w/w). Stratified wetting closely resembled the instability encountered at high flow ratios (>5) of the transition CTAB concentration in that neither phase completely wets the channel walls. The flow cycled through different states of small aqueous plugs, which grow into longer wetting aqueous plugs until the plugs eventually completely wet the PMMA channel walls ( Fig. 7 A). This cycle was most likely the result of an artificial pressure difference caused by the syringe pump as observed by Cabral and Hudson [28] . For ratios between 5 and 1, stable dynamic plug flow ensued similar to liquid–gas flows observed by Cubaud et al. [29] . Plugs varied in length from ∼1800 μm at flow ratio 5 to slightly higher than 900 μm at a flow ratio of 1. The wetting character is also affected by the aqueous phase flow rate. Fig. 8 shows the different flow types produced at 2.0 × 10 −2 % (w/w) CTAB concentration for aqueous flow rates of 1.0 and 5.0 μl/min, respectively. At each aqueous phase flow rate, plug lengths decreased in size according to the organic continuous flow rate (Ca value), but wetting trends tended to remain the same when the flow ratio is lower than 1. While the wetting character changed from no wetting to dynamic wetting by changing the aqueous phase flow rate at a fixed organic continuous flow rate (results are similar for 3.0 × 10 −2 %, not shown). The receding contact angles of at least 10 different plugs for each experimental condition are shown in Table 5 , along with those for CTAB concentration of 3.0 × 10 −2 % (w/w) are shown in Fig. 9 . The receding contact angles for cases in which the aqueous phase exhibited stratified wetting or when drop diameter was smaller than the channel width were assumed to be 0° and 180°, respectively. Also, in cases of stratified wetting, the elongated plugs observed before complete wetting ensued were used to determine the advancing contact angle. It can be seen that the receding contact angle at the tail end of the plug decreasing significantly with flow ratio. And there a critical flow ration of 1 to distinguish the different flow regimes. At flow ratios lower than 1, non-wetting stable plug flow (receding contact angles >90 °C) was observed. While flow ratios were larger than 1, stable dynamic wetting plug flow (receding contact angles <90 °C) was observed. These results are mainly caused by the dynamic adsorption of surfactant molecules on solid–liquid and liquid–liquid interface in microchannels. The effects of surfactant concentration and two-phase flow ratio on the dynamic adsorption behavior of surfactant molecules should be investigated in the future work. Shear-deformed plug flow was observed at CTAB concentrations above 1.0 × 10 −1 % (w/w). In this regime, distinct bullet-shaped aqueous plugs and drops of various sizes were formed, similar to those described in the literature, most likely due to contact angle hysteresis and the Bretherton effect [15,30,31] . Despite having lengths longer than the channel width, the plugs were surrounded by a layer of oil ( Table 6 ). Kolb and Cerro [32,33] , for example, have done extensive research on the motion of two-phase gas–liquid flow in square microchannels and found experimentally that a transition in radial symmetry occurs at Ca ∼0.1 ( Fig. 10 ). In our experiments (although the channels dimensions are not exactly square), a transition in cross-sectional area may occur at Ca value slightly lower than 0.1 at 0.093 ( Table 6 ) due to the difference in liquid phase viscosities (W:O viscosity ratio was 1:67). 3.5 Plug shape analysis The diversity of the plug shapes encountered throughout the flow experiments merited a more in-depth examination of the energies of formation for the different flow types. The drop/plug shape is determined by many factors, such as liquid–solid and liquid–liquid interfacial tensions, the contact angle, and hydrodynamic–viscous effects. In this section, we used the static energy estimates method [34] to discuss the influence of surfactant-controlled liquid–liquid interfacial tension on two-phase flow at higher surfactant concentrations within the microchannel. Fig. 11 shows the interfacial energy estimates of flows at varying flow ratios with different combinations of surface and interfacial tensions at a constant PMMA surface tension assuming no surfactant adsorption. The estimates were obtained by using a method similar to the one used by Kwok et al. [34] for calculating the solid–liquid interfacial tensions. For the system used here, at the CMC of CTAB the lowest combination of water surface tension and liquid–liquid interfacial tension measured was 37.0 and 0.6 mN/m, respectively. The solid surface tension is 38.5 mN/m without surfactant adsorption. The interfacial areas of the plug formations encountered at CTAB concentration 3.0 × 10 −2 % (w/w) were calculated using a MATLAB program. Four different assumptions were made regarding the wetting of the flows: complete wetting by oil (case A), partial wetting by both phases (case B), or complete wetting by water (case C and D), as shown in Fig. 12 . Fig. 12 A shows that the partial wetting state by the aqueous plugs is the lowest in energy when water surface tension and liquid–liquid interfacial tension are near their lower limit values ( Fig. 13 ) of 37.0 and 0.6 mN/m, respectively. Thus, the CTAB surfactant's influence over the surface tension of water and liquid–liquid interfacial tension could change the two-phase flow from complete wetting by oil to partial wetting by both phases, which is consistent with the experimental results ( Fig. 7 B). Substituting lower water surface tension and liquid–liquid interfacial tension values −35.0 and 0.1 mN/m, respectively – into the interfacial energy model, however, shows that complete wetting by water (i.e. phase inversion) becomes the lowest energy state ( Fig. 12 B). Above the CMC value of 4.0 × 10 −2 , the surface tension of water and the water-oil interfacial tension reach a lower limit, which does not change upon increasing the surfactant concentration. Because the lowest experimentally obtained values for the surface tension of water and liquid–liquid interfacial tension were above 35.0 and 0.1 mN/m, CTAB alone cannot lower the water surface tension and liquid–liquid interfacial tension sufficiently enough to induce a phase inversion. In liquid–liquid multi-phase flow systems, however, surfactant concentration and its ability to alter both the liquid–liquid and liquid–solid interfacial tension should still be considered. This method does not take into account changes in the solid surface's interfacial tension due to the presence of surfactant. Furthermore, as with all static models, the viscosity ratio between phases and viscous forces are not taken into account. So there are several weaknesses in the model's applicability. Experimental verification of the interfacial energy estimates requires an in-depth study on the dynamic adsorption of surfactant on the liquid–liquid and solid–liquid interfaces in microchannels and its influence on the dynamics of two-phase flow, which should be investigated in the future work. 4 Conclusions The interfacial tensions and static contact angles of different corn oil-water + CTAB liquid–liquid flows on PMMA were measured across CTAB concentrations ranging from 0% to 5.0 × 10 −1 % (w/w) and a static transition in wetting character (from a contact angle of <90° to >90° for an oil drop on PMMA in water) was observed near a CTAB concentration of 1.0 × 10 −2 % (w/w). Dynamic flow experiments were conducted in a PMMA T-junction microfluidic device to determine if an analogous dynamic wetting transition occurs. At lower CTAB concentrations below 1.0 × 10 −2 % (w/w), stable W/O flow patterns formed but eventually turned into a partial wetting type flow with increasing CTAB concentration. Instability in plug formation climaxed at the transition concentration and a spontaneous phase inversion occurred due to imbalances between interfacial and viscous forces rather than energetic stability. At concentrations higher than the transition concentration, stable plug formations ensued but varied in type. Four different distinguishable flow regimes were observed: stratified wetting, dynamic wetting plug-flow, non-wetting plug-flow, and shear deformed plug flow. The first three types of flow were all observed at concentrations just slightly higher than the transition concentration while shear-deformed plug flow was observed at concentrations above 1.0 × 10 −1 % (w/w). The dynamic receding contact angles of plugs formed at concentrations 2.0 × 10 −2 % and 3.0 × 10 −2 % (w/w) for various flow ratios and capillary numbers were recorded and show that the receding contact angle (at the plug tail) decreases with the increase W:O flow ratio. And there a critical flow ration of 1 to distinguish the different flow regimes. Interfacial energy estimates verify the importance of the solid surface energy and liquid–liquid interfacial tension in determining the flow type. However, the dynamics of the drop/plug formation for systems with surfactant are affected by many factors, such as liquid–solid and liquid–liquid interfacial tensions, the channel surface character, and hydrodynamic–viscous effects and an in-depth study on the dynamic adsorption of surfactant on the liquid–liquid and solid–liquid interfaces in microchannels and its influence on the dynamics of two-phase flow will be investigated in future work. Acknowledgements We gratefully acknowledge the support of the National Natural Science Foundation of China ( 20525622 , 20806042 , 20876084 ), A Foundation for the Author of National Excellent Doctoral Dissertation of PR China (FANEDD: 201053 ) and SRFDP ( 20090002110070 ) on this work. References [1] T. Nakashima M. Shimizu M. Kukizaki Particle control of emulsion by membrane emulsification and its applications Advanced Drug Delivery Review 2000 47 56 [2] H. Song R.F. Ismagilov Millisecond kinetics on a microfluidic chip using nanoliters of reagents Journal of the American Chemical Society 2003 14613 14619 [3] X. Chen H. Wu C. Mao G.M. Whitesides A prototype two-dimensional capillary electrophoresis system fabricated in poly(dimethylsiloxane) Analytical Chemistry 2002 1772 1778 [4] S.J. Harrison V. Skelton Chemical and biochemical microreactors Trends in Analytical Chemistry 2000 389 395 [5] O. Ozer Different emulsion systems for drug and cosmetic delivery European Journal of Pharmaceutical Sciences 1 2007 S11 14 [6] R. Dreyfus P. Tabeling H. Williame Ordered and disordered patterns in two-phase flows in microchannels Physical Review Letters 2003 144505 144509 [7] T. Thorsen R. Roberts F. Arnold S. Quake Dynamic pattern formation in a vesicle-generating microfluidic device Physical Review Letters 2001 4163 4166 [8] L. Shui J. Eijkel A. van de Berg Multiphase flow in microfluidic systems—control and applications of droplets and interfaces Advanced Colloid Interface Science 2007 35 49 [9] J.H. Xu G.S. Luo S.W. Li G.G. Chen Shear force induced monodisperse droplet formation in a microfluidic device by controlling wetting properties Lab on a Chip 2006 131 136 [10] J.H. Xu S.W. Li J. Tan Y.J. Wang G.S. Luo Controllable preparation of monodisperse O/W and W/O emulsions in the same microfluidic device Langmuir 2006 7943 7946 [11] J.H. Xu S.W. Li J. Tan Y.J. Wang G.S. Luo Preparation of highly monodisperse droplet in a T-junction microfluidic device AIChE Journal 9 2006 3005 3010 [12] J. Tan J.H. Xu S.W. Li G.S. Luo Drop dispenser in a cross-junction microfluidic device: scaling and mechanism of break-up Chemical Engineering Journal 2–3 2008 306 311 [13] S.L. Anna N. Bontoux H.A. Stone Formation of dispersions using “flow focusing” in microchannels Applied Physics Letters 2003 364 366 [14] C.P. Ody C.N. Baroud E. de Langre Transport of wetting liquid plugs in bifurcating microfluidic channels Journal of Colloid and Interface Science 2007 231 238 [15] S.A.K. Oskooei D. Sinton Partial wetting gas–liquid segmented flow microreactor Lab on a Chip 2010 1732 1734 [16] G.I. Taylor The formation of emulsions in definable fields of flow Proceedings of the Royal Society of London Series A 858 1934 501 523 [17] W.H. Tan S. Takeuchi Monodisperse alginate hydrogel microbeads for cell encapsulation Advanced Materials 2007 2696 2701 [18] M.E. Schrader Young–Dupre revisited Langmuir 1995 3585 3589 [19] P.G. De Gennes F. Brochard-Wyart D. Quéré Capillary and Wetting Phenomena: Drops, Bubbles Pearls and Waves 2004 Springer New York [20] W. Li Z.H. Nie H. Zhang C. Paquet M. Seo P. Garstecki E. Kumacheva Screening the effect of surface energy on microchannels on microfluidic emulsification Langmuir 15 2007 8010 8014 [21] L. Shui A. van de Berg J. Eijkel Interfacial tension controlled W/O and O/W 2-phase flows in microchannel Lab on a Chip 2009 795 801 [22] G.A. Gaonkar Interfacial tensions of vegetable oil/water systems: effect of oil purification Journal of the American Oil Chemists’ Society 8 1989 1090 1092 [23] A. Halpern The surface tension of oils Journal of Physical Chemistry B 6 1949 895 897 [24] K. Holmberg Handbook of Applied Surface and Colloid Chemistry vol. 1 2002 John Wiley and Sons Hoboken [25] M. De Menech P. Garstecki F. Jousse H.A. Stone Transition from squeezing to dripping in a microfluidic T-shaped junction The Journal of Fluid Mechanics 2008 141 161 [26] E.V. Rebrov A. Benrenguer-Murcia H.E. Skelton B.F.G. Johnson A.E.H. Wheatley J.C. Schouten Capillary microreactors wall-coated with mesoporous titania thin film catalyst supports Lab on a Chip 9 2009 503 506 [27] D.L. Chen L. Li S. Reyes D.N. Adamson R.F. Ismagilov Using Three-phase flow of immiscible liquids to prevent coalescence of droplets in microfluidic channels: criteria to identify the third liquid and validation with protein crystallization Langmuir 23 2007 2255 2260 [28] J.T. Cabral S.D. Hudson Microfluidic approach for rapid multicomponent interfacial tensiometry Lab on a Chip 2006 427 436 [29] T. Cubaud U. Ulmanella C.M. Ho Two-phase flow in microchannels with surface modifications Fluid Dynamics Research 2006 772 786 [30] F.P. Bretherton The motion of long bubbles in tubes The Journal of Fluid Mechanics 1961 166 188 [31] P. Aussilous D. Quere Quick deposition of a fluid on the wall of a tube Physics of Fluids 10 2000 2367 2371 [32] W.B. Kolb R.L. Cerro Coating the inside of a capillary of square cross section Chemical Engineering Science 9 1991 2181 2195 [33] W.B. Kolb R.L. Cerro The motion of long bubbles in tubes of square cross section Physics of Fluids A 7 1993 1549 1557 [34] D.Y. Kwok A. Leung C. Lam A. Li R. Wu A.W. Neumann Low-rate dynamic contact angles on poly(methyl methacrylate) and the determination of solid surface tensions Journal of Colloid and Interfacial Science 1998 44 51