Microsoft Word - 70-508-1-CE vetted COMPLETE cleaned_from the author IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 1 STUDY OF LIQUID-LIQUID SLUG BREAK UP MECHANISM IN A MICROCHANNEL T-JUNCTION AT VARIOUS MODIFIED WEBER NUMBER JITKAI CHIN Department of Chemical and Environmental Engineering, University of Notitngham Malaysia Campus, Jalan Broga, 43500 Semenyih, Selangor Darul Ehsan, Malaysia. Jit-Kai.Chin@nottingham.edu.my ABSTRACT: The formation of immiscible liquid slugs shows various potential applications in different fields because of absolute control of the species both spatially and temporally, as well as enhanced mixing efficiency. Although applications of droplet/slug based microfluidics have been widely demonstrated, the fundamental physics governing slug/ droplet break-up remains to be an area of active research. This study defines an effective Weber Number (Weeff) that characterizes the interplay of interfacial tension, shear stress and dynamic pressure in the micro-channel. Qualitative results show that, at relatively low Weeff, hydrodynamic force plays an important role in slug break- up at the micro-channel T-junction. At intermediate Weeff, both shear force and hydrodynamic force play an equally important role in slug break up. At relatively high Weeff, shear force is dominant in slug break up at the T-junction. ABSTRAK: Pembentukan slug cecair likat tidak bercampur mempunyai potensi di dalam pelbagai bidang kerana daya kawalan mutlak spesis tersebut secara spasial dan temporal, serta kecekapan pencampuran. Walaupun aplikasi titisan/ likat mikrofluida telah didemonstrasi secara meluas, kajian fizik asas yang menentukan pemecahan likat/ titisan masih menjadi bahan kajian yang aktif. Kajian ini mentakrifkan nombor Weber efektif (Weeff) yang mengaitkan ketegangan antara permukaan, tekanan ricih dan tekanan dinamik di dalam saluran mikro. Keputusan kualitatif menunjukkan bahawa, pada Weeff yang rendah, daya hidrodinamik memainkan peranan penting dalam pemecahan likat di persimpangan T saluran micro. Pada Weeff perantaraan, daya ricih dan daya hidrodinamik memainkan peranan yang penting dalam pemecahan slug. Pada Weeff yang relatifnya tinggi, daya ricih adalah dominan di dalam pemecahan likat di persimpangan T. KEYWORDS: liquid-liquid slug break-up, modified Weber Number, hydrodynamic force 1. INTRODUCTION Slug forms spontaneously when two immiscible liquids are brought into contact in a T-junction micro-channel at appropriate flow rates. The potential application of slug-based microfluidics has been demonstrated since years ago and attracted attention from researchers because of the absolute control on species both spatially and temporally, as IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 2 well as internal recirculation in the slug itself, induced by radial shearing along the channel enhances mass transfer that translates into a decrease of mixing time. Although the technique has been widely demonstrated in various applications [1-6], the discussion in the mechanism of slug break up in an immiscible liquid-liquid phase is seldom reported. The present study experimentally observes the dynamic evolution of droplet interfaces and flow patterns that precede a slug break up event in a T-junction microchannel. Interplay of interfacial tension, shear stress and pressure drop in driving the mechanism of break up at various flow rates were investigated. 2. THEORY Slug break up has been studied for decades by carrying out experiments in millimetre scaled capillary tubes. Although the scale discussed in the paper is larger than micrometer, some of the models and the ensuing understanding are relevant to the work. Deformation of a droplet is affected by the external flow pattern and physical properties of both the disperse phase and the continuous phase [7]. Force of inertia, shear stress, interfacial tension and viscosities ratio are the important parameters that act on the droplet which result in its deformation and break up. Hence the Weber Number (We) was proposed to describe the phenomenon. Investigation on droplet end-pinching and break up mechanism in large scale horizontal flow show that breakup of the two-phases flow is due to both shear stress and external pressure field applied by the continuous phase [8]. The imposed shear flow and interfacial tension act in opposite directions, in which imposed shear flow acts to deform the disperse phase while the interfacial tension tends to maintain its shape. The shear stress induces changes in the shape of the extruding disperse phase, causing the asymmetrical shape in disperse phase. The parameters laid the foundation for droplet and slug formations, which are simplified into two dimensionless numbers. The Weber Number, We and Capillary Number, Ca, respectively are defined as: force linterfacia force inertia NumberWeber = γ ρ Du We 2 = (1) tensionlinterfacia force viscous NumberCapillary = γ µu Ca = (2) Several models are proposed for the disintegration of droplet and slug break up in a microchannel configuration. The size of droplet produced under external shear force is a product of the balance between the shear forces with Laplace pressure [9], based on the derivation of viscosity, leading to equation 3. G r µ γ ~ (3) In which the shear rate is estimated from Hagen Poiseuille equation, IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 3 d v G gap 2 ~ (4) Further investigations show [2-3] that the slug formation at a microchannel T-junction is independent of Ca but highly dependent of the flow ratio of the disperse phase to the continuous phase. Their results also show that capillarity is no longer valid when Ca >> 1 x 10 -3 , in which slug break up occurs further downstream in the main channel but not at the T-junction. Further research reveals that slug break up is affected by viscosity of the continuous phase. These work shows that different shear forces applied at the interface might affect the range of each regime [3]. Extensive work also suggests that hydrodynamic force is also one of the important parameters in producing slugs in microchannels, rather than shear stress acting at the liquid-liquid interface [10], which draws conclusion that liquid-liquid slug break up is due to higher magnitudes of hydrodynamic force. However, the experiment only focuses on moderate operating range in which steady slug formation occurs, neglecting slug formation in the disperse phase laminar flow region in the main channel downstream. Notwithstanding, the result is consistent with the findings of Ismagilov et al. [2-3], in which Ca is not an appropriate parameter in characterisation of slug formation in microchannel since the flow dynamic is the most important factor in driving slug break up. In which case, the utilisation of We as a characterisation parameter is recommended. In addition, an effective viscosity is required in a two-phase system as both phases contribute to the shear stress. Hence, a modified We is used for the characterisation, involving the effective viscosity, ρ; continuous phase velocity, vcont ; as well as the hydrodynamic diameter, dhyd; which is shown in Equation 5. γ ρ hydconteff eff dv We 2 = (5) 2.1 Interfacial Tension As disperse phase slug grows in the main channel, the surface area of the disperse phase rear interface in the T-junction varies with time because of the different degree of tangential stress on it by the upstream continuous phase flow. The undefined contact line and unpredictable shape of the rear interface result in spatial and temporal interfacial force changes. Hence, it is calculated based on the surface area estimated from images taken during experiment. Assuming the shape of the slug is limited by the dimension of the channel and the cross sectional area of the microchannel equivalent to a bowl shape, the interfacial forces are equivalent to: ( ) 2 2 1 ~ 2 h hhwA w P π γ +−=       −∆ (6) (7) IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 4 ( ) ( )       +−−       +−      +− 2 2*~ 2 2* 22 ~ 2 2 h hhw w F h hhw hw F πγ π γ γ γ (8) (9) 2.2 Shear Force Similar to calculation of interfacial forces, the contact area between the immiscible phases is important in calculating the shear force. Calculation of the shear force is difficult because it depends on the contact between both phases, based on images taken during the experiment. Fτ= µGA (10) 2.3 Hydrodynamic Force The hyrodynamic force including both the pressure drop along the main channel well as capillary force across interface, is shown in Equation 11. Cdishconth PPPP ∆+∆+∆=∆ ,, (11) Since it is not possible to measure pressure force along the flow direction, Hagen Poiseuille equation is used in estimating pressure drop along the channel in the flow direction, including both disperse phase slug and continuous phase slug, as shown in Equation 12 and Equation 13. From the pressure drop, forces acting at the rear interface are calculated by multiplying the pressure by the cross sectional area, as shown in Fig. 1. contcont hyd cont disdis hyd dis toth conthdishtoth L D Q L D Q P PPP µ π µ π 44, ,,, 128128 +=∆ ∆+∆=∆ (12) (13) Capillary pressure drop is calculated based on Lucas-Washburn equation, including the differences between the receding contact angle and the advancing contact angle, located at the front interface and the rear interface of the slug respectively, which are included in Equation 14.               +−        +=∆ a h a w r h r w Cap rrrr P θθθθγ cos 2 cos 2 cos 2 cos 2 (14) IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 5 Fig. 1: Schematic description of liquid-liquid slug phase flow along the main channel. It is equivalent to a capillary flow in which the hydrodynamic pressure drop along the channel caused by each phase and the pressure drop across the interphase is calculated. The total hydrodynamic force is a combination of hydrodynamic pressure drop and capillary pressure drop over the cross sectional area of the channel, given in Equation 15. ( ) ( )[ ]                             +−        ++         + =∆ ∆+∆= a h a w r h r w contcont hyd cont disdis hyd dis Cap captotht rrrr L D Q L D Q AP PPAF θθθθγ µ π µ π cos 2 cos 2 cos 2 cos 2 128128 44 , (15) 3. EXPERIMENT Before the experiment, the properties of the chemicals are investigated. The experiment is carried out in a glass microchannel network fabricated with conventional wet etching method involving Fluoric Acid (HF) [11]. The channel has a bowl-shape dimension and consists of 2 simple T-junctions. However, only a junction was involved in this experiment. The top dimension of the channel is equivalent to 245 µ m while the height of the channel is equivalent to 50 µm. Tetradecane and acetonitrile (Sigma, 99% purity) are used as disperse phase and continuous phase respectively in this experiment. During experiment, the disperse phase is delivered from the side channel while the continuous phase is delivered into the main channel using standard syringe (BD Plastic), driven by a syringe pump (NE-1000, New Era Pump System Inc.). IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 6 3.1 Method The microchannel network was mounted onto a light inverted microscope platform before acetonitrile and tetradecane were injected into the device, driven by a syringe pump. The flowrate into the microchannel network is controlled from the digital panel on the syringe pump. Velocities for disperse phase and continuous phase range from 1.5 – 9.0 mm/s and 4.8 – 21 mm/s. Slug formed spontaneously when the liquid meet each other at the T-junction. Images of the slug break up were recorded using a high speed camera system (RedLake Inc) mounted onto the microscope, as shown in Figure 2. Images are taken at rates ranging from 250 – 5000 fps, depending on the frequency of slug formation. Image J, an image processing software, was used for analysis. Different flow ratios were controlled by using different syringe sizes. For example for a flow ratio of 0.27, a syringe size of 2.5 ml was used as disperse phase while syringe size of 5 ml was used for the continuous phase. Experiments were carried out with flow ratio of 0.35, 0.27 and 0.15 respectively, as shown in Table 1. The flowrates were then adjusted accordingly while the process of image recording was repeated. Fig. 2: Schematic representation of (a) The setup for this experiment; (b) The glass microfluidic device with the continuous phase flowing along the main channel while the disperse phase from the side arm. Both phases meet at the T-junction flowing downstream towards the outlet. (c) The dimension of the channel assumed to be a bowl shape, in which the cross sectional area is equivalent to a rectangular shape with a quarter circle at both ends. Disperse phase syringe, ml Continuous phase syringes, ml Flow ratio 2 5 0.35 2 10 0.27 1 5 0.14 Table 1: Combinations of different BD standard plastic syringe size giving different flow ratios IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 7 4. RESULTS AND DISCUSSION Surface tensions of acetonitrile and tetradecane are 28.46 mN/m and 25.97 mN/m respectively while the interfacial tension of acetotnitrile/tetradecane system is 15.84 mN/m. The measurements were carried out in a digital tensiometer (K10 ST, Kruss Instrument). Dynamic viscosities of tetradecane and acetonitrile are 2.27 and 0.38 mPaS respectively [12]. Advancing and receding contact angles of the tetradecane slug are measured from images taken during experiment. The ranges of angles are 125 o – 150 o and 40 o – 100 o respectively, depends on flowrate and flow ratio. Slug lengths were measured from images taken throughout the experiment. The discussion below excludes flow ratio at 0.42 because of limitation of the high speed camera in which the images are not sufficiently sharp for measurement and calculation. 4.1 At Relatively Low Weeff (Weeff < 2.0 x 10 -4 ) At relatively low Weeff and across different flow ratio, disperse phase was observed to have a clean snap off at the corner of the side and the main channel, as shown in Fig. 3. Symmetry of the disperse phase into the main channel is disturbed when it almost occupies the cross sectional area of the main channel (Fig. 3a). The disperse phase then fully occupies the cross sectional area of the main channel and blocks it while elongating in the main channel (Fig. 3b). Throughout the process, rear interface in the junction is forced to change from convex in to a concave shape (Fig. 3c), indicating pressure build up in upstream of the continuous phase. The pressure built up increases continuously while the rear interface is being displaced downstream. As the interface approaches the corner of the T-junction and just before breaking up, a thin but visible disperse phase connecting neck was generated between the disperse phase in the main channel and the disperse phase in the side channel (Fig. 3d). However, the plane view thickness and the length of this connecting neck are relatively small and a clean snap off at the junction is reasonably assumed. Due to the relatively short area available for shear and the explicit change in the shape of the rear disperse phase interface in the main channel, it is clearly shown that hydrodynamic plays an important role in slug break up at relatively low Weeff flow. Fig. 3: Images of slug break up captured by high speed camera at Weeff =8.8322 x 10 -5 and flow ratio of 0.35. The interval between images is 80 ms. Clean snap-off takes place during slug break-up. 3(a) 3(b) 3(c) 3(d) IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 8 4.2 At Intermediate Weeff (2.0 x 10 -4 < Weeff < 5.5 x 10 -4 ) As the total flow rate is further increased, the disperse phase slug growth in the main channel is approximating to that at relatively low Weeff (Fig. 4a, 4b, 4c) but the increase flow rate in the disperse phase causes a connecting neck of the disperse phase to evolve into a short trail of thin laminar flow before breaking up a short distance downstream of the T-junction (Fig. 4d). The thin laminar flow layer provides some surface area for shear to take place. Analyses from images show that both hydrodynamic force and shear stress play important role in slug break up. Fig. 4: Slug break up at intermediate Weeff , Weeff = 5.0 x 10 -4 and flow ratio = 0.35. Images were taken at interval of 20 ms. Small connecting neck exists before slug break- up. 4.3 At Relatively High Weeff (Weeff > 5.5 x 10 -4 ) As flow rates of both disperse phase and continuous phase are increased, the break up mechanism varies from previous cases. The continuous phase pushes the extruding disperse phase in the junction towards main channel downstream, before it can fully occupy the entire cross sectional area of the main channel. This causes increase in total contact area between both phases (Fig. 5a, 5b). Compared to low and intermediate Weeff, the change of the shape of the rear interface in the T-junction was not obvious and the plane view rear interface in the T-junction and downstream main channel did not approximate to the top width of the microchannel (Fig. 5c), at any time throughout the break up process. Shear interaction between the disperse phase and the continuous phase in the main channel causes elongation of the disperse phase laminar flow region in the main channel while the slug breaks off from the laminar flow region (Fig. 5d). Fig. 5: Slug break up at relatively high Weeff , Weeff = 1.18 x 10 -3 and flow ratio = 0.35. Images were taken at interval of 36 ms. A large disperse phase laminar flow region formed downstream in the main channel. Slug breaks up at the tip of the region in the main channel. 4(a) 4(b) 4(c) 4(d) 5(a) 5(b) 5(c) 5(d) IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 9 The experimental data and calculation is summarised in Table 2. It shows that hydrodynamic force is an important factor in slug break up at relatively low Weeff. At intermediate Weeff, both shear stress and hydrodynamic force are important for slug break up. At relatively high Weeff, the break up is dominated by shear force between both phases in the main channel. Flow ratio = 0.35 Weeff 8.83 x 10 -5 2.28 x 10 -4 4.34 x 10 -4 9.89 x 10 -4 Shear force, Fτ 3.50 x 10 -11 1.67 x 10 -10 2.42 x 10 -10 4.51 x 10 -10 Hydrodynamic force, Fh+c 1.89 x 10 -5 5.00 x 10 -5 7.46 x 10 -5 8.36 x 10 -5 Ratio of forces, Fτ/ Fh+c 1.86 x 10 -6 3.34 x 10 -6 3.24 x 10 -6 5.39 x 10 -6 Flow ratio = 0.27 Weeff 1.69 x 10 -4 2.22 x 10 -4 5.91 x 10 -4 8.88 x 10 -4 Shear force, Fτ 3.90 x 10 -11 1.04 x 10 -10 1.65 x 10 -10 2.40 x 10 -10 Hydrodynamic force, Fh+c 1.68 x 10 -5 2.44 x 10 -5 3.09 x 10 -5 3.21 x 10 -5 Ratio of forces, Fτ/ Fh+c 2.32 x 10 -6 4.26 x 10 -6 5.34 x 10 -6 7.48 x 10 -6 Flow ratio = 0.15 Weeff 7.25 x 10 -5 1.13 x 10 -4 1.72 x 10 -4 2.88 x 10 -4 Shear force, Fτ 7.32 x 10 -11 9.04 x 10 -10 1.15 x 10 -10 2.69 x 10 -10 Hydrodynamic force, Fh+c 4.44 x 10 -5 6.21 x 10 -5 8.05 x 10 -5 7.55 x 10 -5 Ratio of forces, Fτ/ Fh+c 1.65 x 10 -6 1.46 x 10 -6 1.43 x 10 -6 3.56 x 10 -6 Table 2: calculation of force ratio according to Weeff in different syringe combination giving different flow ratios In each flow ratio, Weeff is increased by increasing in overall flowrates in both syringes. In Flow ratio = 0.35, when Weeff < 2.5 x 10 -4 , tetradecane slugs have clean snap off taken place at the corner of T-junction. As the Weeff is increased, a disperse phase small trail starts to form at the corner. A disperse phase laminar flow region formed in the main channel as Weeff is further increased. Throughout the experiment, phenomena of clean snap off until failure of the slug formation at the T-junction is observed. In flow ratio = 0.27, the experiment starts with clean snap off of slug formation. However, a small trail is formed as increment in Weeff is carried out. Further increment in Weeff causes the tetradecane slugs to be produced from the laminar flow region formed in the main channel. In third experiment, in which 1ml syringe is used in disperse phase while 5ml syringe is used in continuous phase to give flow ratio of 0.14, tetradecane slug formation takes place across a wide range from 7.25 x 10 -5 until 1.72 x 10 -4 . A trail is formed at the corner for slug formation as Weeff is increased to 2.88 x 10 -4 . Further increase in Weeff results in the phenomenon failure at Weeff = 7.23 x 10 -4 , however meaningful results are limited due to capability of the high speed camera. IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 10 The preliminary study shows a general trend that the phenomena of slug formation in a T- junction, ie clean snap off, formation of a small trail at the corner and formation of a disperse phase laminar flow region depends on the Weeff but independent of flow ratio. Validation is carried out by plotting disperse phase slug length against Weeff, as shown in Fig. 6. It is found that slug length approaches infinity when Weeff is extremely low. As Weeff increases gradually, the slug length generated drops significantly until Weeff ~1 x 10 - 4 . Further increase in Weeff has negligible effect on the slug length. The plots are similar to the trend in Hagen Poisuelle equation, with t showing that the slug length as an inverse function of pressure drop in the main channel. Fig. 6 confirms the suitability of utilising Weeff as characterising parameter especially at relatively low flow rate. Fig. 6: Plotting of relative disperse phase slug length to Weeff. 5. CONCLUSION Understanding the fundamental physics of slug formation in a microchannel T- junction is important in exploring further potential application of the phenomenon into other fields. However the flow dynamic is always complicated and difficult to model with simulation any software available. Although the experiment was carried out only with a single immiscible liquid-liquid system, the qualitative result and simplified analysis method still give a few important points about slug break up phenomenon in a microchannel T-junction. From the experimental work, it can be concluded that at relatively low Weeff, (Weeff < 2.0 x10 -4 ), in which disperse phase slug clean snap off occurs, the ratio of shear force to the hydrodynamic force is always less than 3.0 x 10 -6 . hydrodynamic force plays an IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 11 important role in disperse phase slug break up at relatively low Weeff with a clean snap off at the T-junction. At intermediate Weeff, (Weeff = 2.5 - 5 x10 -4 ), the force ratio is approximately 3.0 – 5.0 x 10 -6 . A small trail of disperse phase extruding into the main channel, resulting in both hydrodynamic force and shear force playing important role in slug break up. In the case when Weeff > 5 x10 -4 , disperse phase slug breaks up from a laminar flow region in the main channel with the force ratio at approximately 4.5 – 6.0 x 10 -6 . large contact surface area of the disperse phase laminar flow region in the main channel provides large surface area for shear interaction between both phases. Further analysis justifies the utilization of Weeff as a characterizing parameter. 6. FUTURE WORK Experimentations with different immiscible liquid-liquid systems are being carried out in verifying the model and establishing the theory for slug break up across different flow ratio and viscosity ratios. In addition, the use of a high speed camera which gives higher resolution and higher frames per second, fps, are useful in analysis. ACKNOWLEDGEMENT The author would like to acknowledge Ministry of Higher Education Malaysia (MOHE) Fundamental Research Grant Scheme (FRGS/2/2010/TK/UNIM/03/3) for funding the research project. REFERENCES [1] Burns, J. R., and Ramshaw C., The intensification of rapid reactions in multiphase systems using slug flow in capillaries, Lab-on-a-chip, 2001, 1, pp. 10–15 [2] Song H., Tice J. D., and Ismagilov R. F., A microfluidic system for controlling reaction networks in time, Angw. Chem., 2003, pp. 768–772 [3] Tice J. D., Lyon A. D., and Ismagilov R. F., Effects of Viscosity on Droplet Formation and Mixing in Microfluidic Channels, Analytica Chimica Acta 507, 2004, pp. 73–77 [4] Niu, X.Z., Zhang, B., Marszalek, R.T., Ces, O., Edel J.B., Klug, D.R. and deMello, A.J. Droplet-based compartmentalization of chemically separated components in two-dimensional separations, Chem. Commun., 2009, 7(41), pp. 6159-6161 [5] Mazutis, L. And Griffiths, A.D., Preparation of Monodisperse Emulsions by Hydrodynamic Size Fractionation, Applied Physics Letters, 2009, 95, pp. 204103- 04105 [6] Srinivasan, V., Pamula, V.K. and Fair, R.B., An Integrated Digital Microfluidic Lab-on-a- chip for Clinical Diagnostics on Human Physicological Fluids, Lab-on-a-chip, 2004, 4, pp. 310-315 [7] Hinze J., Fundamentals of the hydrodynamic mehcanism of splitting in dispersion processes, AIChE, 1955, 1, pp. 289–295 [8] Zhang D. F., and Stone H. A., Drop formation in viscous flows at a vertical capillary tube, Phys. Fluids, 1997, 9, pp. 2234–2242 IIUM Engineering Journal, Vol. 12, No. 2, 2011 Chin 12 [9] Thorsen T., Roberts. R. W., Arnold F. H., and Quake S., Dynamic Pattern Formation in a Vesicle-Generating Microfluidic Device, Phy. Rev. Lett., 2001, pp. 4163–4166 [10] Garstecki P., Fuerstman M., Stone H. A., Whitesides G. M., Formation of droplets and bubbles in a microfluidic T-junction – Scaling and mechanism of break up, Lab-on-a-chip, 2006, 437 [11] Che-Hisn Lin, Gwo-Bin Lee, Yen-Heng Lin, and Guan-Liang Chang, A fast prototyping process for fabrication of microfluidic systems on soda-lime glass, J. Micromech. Microeng., 2001, 11, pp. 726–732 [12] Perry’s Chemical Engineers’ Handbook, 1997 NOMENCLATURE Re Reynolds number We Weber Number Ca Capillary Number Eff Effective ρ density, kg/m 3 L characteristic length, m u flow velocity, m/s µ viscosity, kg/m 2 s r radius, m h height of the channel, m w width of the channel, m Dhyd Hydraulic Diameter, m G shear rate, Cont Continuous Phase Disp Dispersed Phase γ Interfacial tension Q Volumetric Flowrate, kg/m 3