Hydrodynamics of Liquid-Liquid Flows in Micro Channels and Its Influence on Transport Properties: A Review

Hydrodynamics plays a major role in transport of heat and mass transfer in microchannels. This includes flow patterns and flow regimes in which the micro-channels are operated. The flow patterns have a major impact the transport properties. Another important aspect is the pressure drop in micro-channels. In the present review, the experimental and Computational Fluid Dynamics (CFD) studies covering all the above aspects have been covered. The effect of geometrical parameters like shape of channel, channel size, material of construction of channels; operating parameters like flow velocity, flow ratio and fluid properties have been presented and analyzed. Experimental and analytical work of different pressure drop models has also been presented. All the literature related to influence of flow patterns on transport properties like volumetric mass transfer coefficients (VMTC) and heat transfer coefficients (HTC) have been presented and analyzed. It is found that most works in Liquid-Liquid Extraction (LLE) systems have been carried out in slug flow and T-junctions. Models for coupled systems of flow and mass transfer have been presented and works carried out for different coupled systems have been listed. CFD simulations match experimental results within 20% deviations in quantitative and qualitative predictions of flow phenomena for most research articles referred in this review. There is a disparity in prediction of a generalized regime map and a generalized regime map for prediction of flow patterns for various systems would need the help of Artificial Intelligence.


Introduction
Microchannels have been an area of interest since the last two decades in various applications including petroleum industries, emulsification, fine and speciality chemicals (like nitration reactions), extraction processes, biochemical applications etc. Flow patterns play a vital role in the transport properties (heat and mass transfer) in non-reactive and reactive systems. Conventional equipments like agitated vessels, rotating discs, impinging jets (for mass transfer) [1] and jet impinge heat exchange mechanisms used for cooling gas turbine components, heat treatments, [2] preheating in glass manufacturing [3] and many other conventional heat transfer equipments [4] have order of magnitude lower transport properties than those offered by microchannels. Flow patterns are formed due to interplay between inertial and surface tension forces. The gravitational forces are of less importance in the microchannel based devices and systems. The viscosities differences of the continuous and dispersed phases give rise to several different flow patterns. The dimensionless numbers used in investigations are: Reynolds number: Re = ρvD µ ; Capillary number: Ca = µv σ ; Weber number: We = Re*Ca; Ohnesorge number: Oh = µ √ Dσρ . Bond number: Bo = ∆ρgd 2 /σ. The same dimensionless numbers have been used for the inlet junction flow patterns. Excellent reviews have been published on different flow patterns as well as its effect on transport properties both in gas-liquid (GL) and liquid-liquid (LL) systems. [5,6]. However, reviews on flow patterns of LL systems focus on major flow patterns like slug-flow, parallel flow and transition flow. It is important to note that the viscosity of the continuous and dispersed phase gives rise to various different flow patterns [7,8]. Also, the flow patterns depend on the bubble formation and breakup during the tubing, squeezing and dripping regimes and transitions between them. Further, reviews [5] have focused mostly on the effect of geometrical and operating parameters, material of construction of microchannels but the effect of surfactants, effect of hydrophobicity and hydrophilicity on walls, effects of wall films on flow patterns have not been thoroughly discussed. The present work focuses on reviewing all the prominent literature works (including experimental, numerical works) on flow patterns and regime maps. First, state of the art of experimental and CFD studies describing flow patterns in liquid-liquid flows have been analyzed. This includes flow patterns in different geometries and geometrical parameters like shape, size, material of construction, operating parameters like flow rates, flow ratios, fluids properties of the fluids under consideration, measurement techniques used by researchers to understand internal circulations in slugs of both phases are also presented. Secondly, most important flow regime maps signifying the maps based on forces and based on dimensionless groups (individual and combination of different dimensional numbers have been presented. Finally, pressure drop in microchannels (experimental data) of various authors has been presented for different two phase systems and geometries. The corresponding models and their limitations have been presented. Finally, the review lists out the various liquid-liquid reactive and non-reactive systems carried out by researchers. Pressure drop studies has been covered well in recent review by Al-Azzawi et al. [5]. In the present work, the most prominent works have been reported. The prominent works showing effect of the flow patterns on transport properties (specifically Volumetric Mass Transfer Ceofficients (VMTC)) have been included. Figures 1 and 2 show the entire structure of the manuscript. Figure 1 focuses on the experimental studies while Figure 2 focuses on the CFD studies.

Hydrodynamics in Microchannels
Micro-channels provide advantages like high surface to volume ratio that ensures better control on size of the droplets. A typical experimental set-up for carrying out experiments to study flow patterns, mass transfer, pressure drop etc. is shown in Figure 3. The two immiscible liquids (aqueous and organic, from reservoirs (1 and 2) are introduced by continuously operating high-precision piston pumps (3 and 4) to a symmetric 120 • Y-junction or T-junction made of materials like Teflon, glass, steel etc. (5). The pumps are controlled with the help of a computer to adjust the flow rate. The capillary contactor (6) is made of poly-tetra-fluoro-ethylene (PTFE) and was attached directly downstream of the Y-or T junction. The length of the capillary micro-channel can be adjusted easily. Thus, the setup was very flexible with respect to the residence time used. At the end of the contacting stage, a flow splitter (7), was used to separate the phases.  (8,9) sampling bottles. (Reproduced with permission from [9]) the two outlets of the splitter were collected in the sampling bottles (8,9) for further analysis. With permission from Kashid et al. [1].
Typical configurations for micro-channels are shown for in Figure 4. The micro-channel geometries are different in cross sectional shape, main channel shape and mixing junction configuration. Regarding the cross sectional shape of the micro-channel, micro-channels can be classified into two main groups: circular and non-circular micro-channels, the latter can be square, rectangular, trapezoidal, etc. Different types of junction configurations are used at the micro-channel entrance. The most popular mixing junctions are T-junction, Y-junction, cross-shaped junction, cross flow T-junction and concentric junction. The cross sectional dimension of the inlet channels can be different from those of the main channel. Different junction configurations can produce different stable flow regimes. For example, a micro-channel with a Y-shaped junction is more appropriate for producing a stable liquid-liquid parallel flow as compared to that with a T-shaped junction [9]. Droplet flow in a micro-channel can be achieved by using a concentric inlet or cross-shaped junction. Further, details into dispersion in micro-channels with different junction structures can be found in the literature [10][11][12]. In the forthcoming sections, we will discuss about the flow patterns, pressure drop in micro-channels using T and Y junctions and review the current experimental and CFD modeling works.

Experimental Studies: Flow Patterns, Flow Regimes, Pressure Drop
In this section we present the various experimental works by the researchers for studying the different flow patterns occurring in micro-channels and how they have been represented using regime maps. Measured pressure drop for different two phase systems at different operating conditions and effect of material characteristics and operating parameters are also presented. Research works involving CFD simulations depend on mathematical modelling to track the interface and reproducing the experimentally observed flow patterns for different geometric and operating conditions have been presented.
Experimental research works mainly in published literature have mainly focused on flow patterns, internal circulations inside slugs [13] (in slug flows), flow regimes and pressure drop. Some analytical models for pressure drops have also been presented. These models have been validated with experimental works in literature.

Flow Patterns
Flow patterns are divided into two parts: the first and most significant is flow patterns in the main channel and second is the flow patterns in the junction. Flow patterns in microchannels have been studied explicitly by many researchers [1,7,[14][15][16]. In micro-channels distinctly visible flow patterns are annular flow, parallel flow, slug flow, droplet flow and deformed interface flow provided by the authors [6,7,15,16]. Other flows which may form are slug-drop, slug-dispersed etc. depending on the operating parameters like flow rates of the two fluids. The viscosity of the fluids play a major role in determining some flow like rivulet and serpentine flow (Yagodnitsyna et al. [7]) instead of annular flows. Flow patterns at the junctions are carried out by a few researchers [11,[17][18][19][20][21].
It should be noted that the flow patterns in main channels depend on which regime exists in the junctions. Different patterns are observed at cross-junctions and cross Tjunctions. The patterns of cross-junctions and cross-T junctions are Tubing/threading, dripping and jetting. Three different systems with different viscosities and surface tensions namely, water-kerosene, water-paraffin-oil and paraffin-oil-castor-oil were considered by researchers [7]. They observed rivulet flow and serpentine flow in addition to the flow observed mentioned previously.
The different types of flows in the main channel which are essentially due to the interplay between the shear, inertial and surface tension forces are as follows (Refer Figure 5): 1.
Droplet flow: In this type of flow the dispersed phase flows in the form of droplets. The droplets have diameter quite lesser than the hydraulic diameter and have continuous renewal of surface due to the continuous phase flow (Figure 5a).

2.
Slug flow: This is the most observed flow in micro-channels. When two miscible fluids are introduced into the junction, one fluid passes downstream while other penetrated over opposite side of the junction. This mutual displacement causes the slug flow structure (Kashid et al. [14]). In this type of flow, the dispersed phase flows in the continuous phase in the form of slugs. These slugs have a diameter and a slug length little lesser/same than the hydraulic diameter of the channel and if lesser due to the film thickness. The slugs are separated from each other with a specific distance from each other. There are internal circulations within the slugs which aid in decreasing the boundary layers near the interface of the dispersed and continuous phases thus promoting mass transfer. In addition to this, there are circulations in between the two slugs which promote diffusion ( Figure 5b).

3.
Slug Plug flow: This type of flow is described differently by different authors. Plugs as liquid segments in symmetrical plug which does not necessarily wet the channel wall have been defined by Yagodnitsyna et al. [7]. The plug length is greater than the slug length. 8.
Rivulet flow: This new flow pattern was observed by Yagodnitsyna et al. [7]. The authors defined the rivulet flow pattern to be a parallel flow when one of the liquids forms a rivulet which wets channel wall and its shape along the channel length can be warped (Figure 5g). 9.
Serpentine flow: This new flow pattern was observed by Yagodnitsyna et al. [7]. In this type of flow, a periodical wavy interface (like a twisted tape) is observed. The periodical wavy interface is similar to the overlapping longitudinal and traverse mode vibrations in Melde's experiment and is stationary from beginning to end (Figure 5h). 10. Dispersed Flow: This type of flow might occur only if there is sufficient amount of mixing due to obstructions etc. We will see this flow in greater detail in the forthcoming sections. The types of flows in the junctions are given below Wu et al. [11]: (explain in detail regarding the flow regimes) 1. Tubing/threading: In this type of flow, organic liquid flows from the junction in the form of an annular core in the main channel with minor interface deformations. Threading has a cylindrical organic core while tubing has approximate square cross section with round corners (Figure 5i) 2. Dripping: In this type of flow the organic core tends to largely obstruct the intersection which constricts the aqueous (water in the case considered by the authors) flow. Pressure is exerted by water coming in from both sides (i) squeezing pressure and (ii) viscous shear. Both these forces overcome the interfacial shear to send the organic droplet into main channel in the form of slug. (Figure 5j,k) Cubaud and Mason [17]. Inertial forces have no significant effect on the droplet formation due to proximity of Re<1 for both phases. Droplet size variation depends on Ca number. For example, for Ca = 0.019 droplet The types of flows in the junctions are given below Wu et al. [11]: (explain in detail regarding the flow regimes)

1.
Tubing/threading: In this type of flow, organic liquid flows from the junction in the form of an annular core in the main channel with minor interface deformations. Threading has a cylindrical organic core while tubing has approximate square cross section with round corners (Figure 5i Dripping: In this type of flow the organic core tends to largely obstruct the intersection which constricts the aqueous (water in the case considered by the authors) flow. Pressure is exerted by water coming in from both sides (i) squeezing pressure and (ii) viscous shear. Both these forces overcome the interfacial shear to send the organic droplet into main channel in the form of slug. (Figure 5j,k) Cubaud and Mason [17]. Inertial forces have no significant effect on the droplet formation due to proximity of Re < 1 for both phases.
Droplet size variation depends on Ca number. For example, for Ca = 0.019 droplet size varies upto 15%. While the variation is lower for Ca = 0.011 is higher than 15% [22].

3.
Transition regime: This regime is transition from dripping to sqeezing regime. In transition regime the Ca numbers are in range (0.009 < Ca < 0.013) and Re < 1. Hence inertial forces have no significant effect on flow patterns 4.
Sqeezing regime: In this regime the inertial forces play a big role.

5.
Jetting: In this type of flow the organic core thread falls in certain distance into the outlet main micro-channel before it breaks up into slugs/droplets. At jet breakup, interfacial shear is overcome by viscous shear. (Figure 5l).
In the upcoming subsections, the effect of different parameters like geometrical parameters like shape, size of microchannels, material of construction of channels, type of junction for flow mixing and inlet, fluid properties of the dispersed and continuous phases like the viscosity ratio, interfacial tension and the effect of other factors like addition of agents or application of vacuum to stabilize flow patterns have been presented. Table 1 gives the literature works on the flow patterns while Figures 6 and 7 gives the chart showing the flow patterns due to the sensitivity of several parameters involved.   Role of Shape, Size, Geometry and Material Properties of Microchannels Svenssen et al. [25] investigated the effect of square and rectangular geometries on the flow patterns for a water-nhexane (where n-hexane was organic phase and water was aqueous phase) system in glass microchannels (T-junctions). Slug flow in dripping regime was observed. The flow patterns remain same with increase in size for same flow velocities but effect the transport properties which will be discussed in Section 3. Kashid and Kiwi-Minsker [24] studied the effect of flow velocities with T-junctions of different shapes like square and trapezoidal for variation of flow velocities from 20 to 150 mm/s. For velocities up to 250 mm/s slug flow was observed while for velocities up to 500 mm/s slug-drop flow was observed. Further, from 500 to 1100 mm/s deformed-interface flow was observed while for 1200 mm/s to 2000 mm/s parallel flow was observed. While the flow patterns were same for both shapes the effect on transport properties were better in square than trapezoidal shapes. Further, the authors performed studies in caterpillar microchannels which gave flow patterns like slug flow, slug-dispersed flow and dispersed flow and gave very good results as far as transport properties are concerned.
Kashid and Kiwi-Minsker [24] carried out experimental investigations in 4 different types of geometries namely T-square, Trapezoidal, Y-rectangular, Concentric microchannels. The authors observed that hydraulic diameter of the microchannels and volume fraction of the dispersed phase plays an important role in addition to dimentionless numbers.
Tsaoulidis et al. [31], carried out experimental investigations on 0.2 mm and 0.27 mm micro-channels with two different materials glass and Teflon. Further, the authors used Ionic Liquids (IL) as the dispersed phase and found 3 distinct flow patterns namely annular, plug and drop flow.

The Role of Fluid Properties
Sugiura et al. [32] have considered tri-olein water system. They observed that the droplet diameter was independent of surface tension (σ = 3.6 to 17.5 mN/m) for lower flow velocity ranges. Yagodnitsyana et al. [29] carried out experimental studies of Ionic Liquid (IL) (1-Butyl-3Methyl-imidazolium bis (trifluoromethylsulfonyl)amide)-Water. The authors found a relation between plug velocity (plug as defined by Kovalev et al. [4]) and bulk velocity. U plug = 2.34U 1.087 bulk . The authors investigated the influence of bulk velocity on plug lengths and found that plug length increased with flowrate ratio increase while decreased with increase in bulk flow.
Yagodnitsyna et al. [30] carried out investigations in a T-shaped microchannel using 3 sets of LL systems: Kerosene-water; paraffin-oil-water; castor-oil-water. The authors obseved different flow patterns like rivulet flow instead of annular flow due to the wettability at the walls by the liquids of systems investigated. Serpentile flow was observed for higher Weber numbers. The authors concluded that flow pattern maps were different for different fluidic systems. Further, absence of abundant information on influence of parameters like density, viscosity, surface tension, contact angle on flow patterns were highlighted by the authors. Most research works used systems where one of the liquids did not wet the fluid walls. The authors concluded that rivulets had less velocities than plugs based on a relation between velocity of plugs and bulk velocities. For the systems considered the relation was U plug = 1.21U bulk . Further, the circulations grow linearly with plug velocity.
Viscosity of both phases play an important role in flow patterns both in the continuous and dispersed phases (velocities inside droplets, slugs/plugs). Su et al. [23] found that when dispersed phase viscosities are higher than continuous phase viscosities, faster flow transition took place from slug flow to parallel flow and resulted in smaller slugs. Kovalev et al. [8] carried out experimental investigations in T-junctions. Effects of continuous phase viscosity on water slug lengths have been studied for very low capillary number (2.32 × 10 −3 -2.32 × 10 −2 ) by Garstecki et al. [33]. They observed marginal change in the slug length with an increase in silicon oil viscosity from 10 mPas to 100 mPas.
Kovalev et al. [8] found different flow patterns for low viscosity ratios upto 0.001. The six distinct flow patterns were Plug flow, plug flow with suspended particles, droplet flow, throat annular flow, parallel flow and slug flow ( Figure 8). The definition of slug flow (as defined by the authors) is different than the definition of literature works. Here slug flow is periodical flow of dispersed phase segments which have assymetrical shapes and wet at least one channel wall (as shown in Figure 8). According to the authors the plug flow is the one which have diameter and length little less than the hydraulic diameter of column. Further, the authors defined dropet flow as flow of almost spherical droplets. Interesting features of plug flows have been reported by the authors (a) The plugs took different shapes bullet like, dumble like shape with possibility of microdroplet breakoff. (b) The velocity distribution inside the plugs were found to be three dimensional in nature in contrast to two dimensional nature generally observed. (c) Flow structures evolved downstream in contrast to that observed in literature. (d) Velocity distribution in plug-plug interface was closely associated with each other. Dimensionless parameters (product of flow ratio and Capillary number) were defined to define the plug flow.  Figure 9 shows how plugs and slugs are distinguished according to their lengths as observed by Foroughi and Kawaji, [34]. If average length of the water segments was equal to or less than 5 channel diameters the flow was classified as slug flow. The author have kept the dispersed phase flow rate constant while increasing the flow rate of continuous phase. With increase in continuous phase flow the flow patterns changed from drop to plug to slug. With further increase in water flow sausage shaped annular patterns were noticed. With highest flow rates, water droplets surrounding the oil core were observed.
The authors developed a correlation for slug length as a function of flow ratio and capillary number and also related plug velocity as a function of bulk velocity. Vansteene et al. [26], carried out experimental investigations in in T and Y junctions in glass microchannels in the dripping regime. The authors developed an empirical correlation for normalized droplet diameter (NDD) for high viscosity ratios. NDD was a function of flow ratio and capillary numbers. The empirical model was found to control flow parameters as a means for optimization of droplet production frequency and specific interfacial area. The authors developed a correlation for slug length as a function of flow ratio and capillary number and also related plug velocity as a function of bulk velocity. Vansteene et al. [26], carried out experimental investigations in in T and Y junctions in glass microchannels in the dripping regime. The authors developed an empirical correlation for normalized droplet diameter (NDD) for high viscosity ratios. NDD was a function of flow ratio and capillary numbers. The empirical model was found to control flow parameters as a means for optimization of droplet production frequency and specific interfacial area.

Internal Circulations Inside Slugs/Plugs, Droplets Using Velocimetry Techniques
Internal circulations play an important role in enhancement of transport properties and reaction rates [30].
Kovalev et al. [8] used Particle Tracking Velocimetry (PTV) to understand the hydrodynamics inside plugs. The authors found important findings that the particle distribution inside plugs were different than the results obtained from flow visualization specifically the dumble like interface shapes at high plug velocities ( Figure 9).
The flow inside the plugs were essentially found to be three dimensional than the two dimensional flows often reported in literature. Further, internal circulations inside the slugs (dispersed phase) as well as outside the slugs (continuous phase in between two slugs) have been studied using μ-PIV (as seen in Figure 10).

Internal Circulations Inside Slugs/Plugs, Droplets Using Velocimetry Techniques
Internal circulations play an important role in enhancement of transport properties and reaction rates [30].
Kovalev et al. [8] used Particle Tracking Velocimetry (PTV) to understand the hydrodynamics inside plugs. The authors found important findings that the particle distribution inside plugs were different than the results obtained from flow visualization specifically the dumble like interface shapes at high plug velocities ( Figure 9).
The flow inside the plugs were essentially found to be three dimensional than the two dimensional flows often reported in literature. Further, internal circulations inside the slugs (dispersed phase) as well as outside the slugs (continuous phase in between two slugs) have been studied using µ-PIV (as seen in Figure 10). Li et al. [28] have studied internal flow in droplets for varying Bond numbers (0.003 < Bo < 0.007) and Capillary numbers (4.05 x10 −4 < Ca < 1.62 × 10 −4 ) to see the effect of viscosity ratio and flow on droplet size and droplet spacing and in turn on internal circulations. The authors used µ-PIV in their investigations. The investigations were based on dimensionless droplet size (ratio of droplet length to inlet channel width) and droplet spacing (ratio of distance between two droplets to inlet channel width).

Effect of Droplet Spacing on Internal Circulations
At low inlet flowrates of the continuous and dispersed phases, the authors observed very low velocities in the region between the droplets and having reverse circulations. As the flow rates increased the circulation velocities dissapeared with an increase in droplet spacing. The influence of the wall on droplet movement was found to decrease with an increase in droplet spacing.

Effect of Viscosity Ratio on Internal Circulations
Li et al. [28] observed the effect of viscosity ratio on the internal circulations inside droplets. For lowest flowrate of their investigation, Ca = 2.43 × 10 −5 ; a single clockwise eddy was observed in the upper part of the droplet while for highest flowrate, Ca = 9.71 × 10 −4 ; two clockwise eddies were observed at the upper part of the droplet. The authors clarified this due to the increase in the shear stress between the droplets and dispersed phase when continuous phase velocity decreases. This makes movement of droplets in continuous phase difficult and the droplets are dragged by the continuous phase resulting in higher internal circulations.

Effect of Droplet Size
Li et al. [28] observed that when the droplet size is same as the channel width two counter rotating eddies are observed on either side of the droplet while when it is 1.5 times the channel width, the four eddies are observed with zero velocities at the core of the eddies.

Effect of Curved Channels on Internal Circulations
Liu et al. [35] studied the internal flow fields in droplets, eddy dynamics behind droplets in curved channels to study the effect of Ca, viscosity ratio, droplet size, channel curvature and interfacial tension (by addition of surfactants). The authors observed that the above parameters except Capillary number reduced oil film gap resulting in dissapearance of eddies behind droplets. The internal circulations reduce from two circulations to single circulation with decrease in droplet size.

Effect of Flow Ratio
Garstecki et al. [33] also investigated the effect of flow ratio on slug flow formation mechanism and slug length. It was found that for flow ratio (q < 1) there was a marginal effect on slug length. However, for higher values of flow ratios (q > 1) slug lengths were found to increase. The flow ratio, was found to be an important parameter in transition of slug detachment mechanism. Flow ratio variation is an important parameter for impacting the transport properties like mass transfer coefficients. Due to increase in slug length by increase in flow ratio the transport properties like volumetric mass transfer coefficients are affected due to increase in interfacial area. This will be explained in greater detail in Section 3.

Role of Surfactants
Svensson et al. [25] found the effect of surfactants to two phase flows in square glass microchannels. They reported that the surfactant molecules adsorb on the glass wall hindering the flow between the organic phase slugs and the channel wall. This increases the resistance and effects the transport properties (specifically mass transfer). However, some authors have reported improvement in transport properties which will be discussed in Section 3. It is reported that the presence of surfactants enhance transport properties by 1.44 times.

Role of Vacuum at Microchannel Outlet
Parallel flow has received lesser importance due to its transition to slug flow. Recently, Li et al. [27], performed experiments to investigate stable parallel flows for Re < 0.01. The authors used Y microchannels and applied vacuum at the outlet to have stable parallel flow. The authors showed that the stability of parallel flow was due to the reduction in fluid-solid interaction at the walls due to application of vacuum. Most importantly, no surface modifications, additional channel structures or surfactants were required to get stable parallel flows in contrast to earlier works. The authors also observed transition of flow patterns by changing the inlet pressures of aqueous (water) and organic (oil) phases. The authors reported that when inlet pressure of water phase (P in,w ) was much higher than that of oil phase (P in,o ). When P in,o was higher than P in,w droplet flow was observed. When P in,o was further increased, parallel flow was observed while for much higher inlet pressures single phase oil flow was observed.

Flow Regimes and Regime Maps
Flow regimes are maps which determine the type of flow by taking into account the various parameters like geometric, dynamic and fluid properties. Researchers, [11,15] have used plots of continuous v/s dispersed phase velocities or dimensionless numbers (product of Weber and Ohnesorge numbers or Weber vs Capillary numbers) for continuous and dispersed phases plotted against each other to define regime maps ( Figure 11A,B). The regime maps based on the domination of different forces have also been shown ( Figure 11C,D). Finally, Figure 11E shows a generalized regime map which compares the data by various authors. The final regime map considers LLE systems which have a transfer component unlike other regime maps which are just two phase LL systems. Interestingly, Yagodnitsyna et al. [7] observed different flow regimes than the regular flow regimes reported by other authors.
Yagodnitsyna et al. [7] carried out dimensionless analysis to find that the regimes depended on two dimensionless number Weber number and Ohnesorge number with the relation We ε/2 * Oh β/2 . The data showed that ε = 2 and β = 1. By plotting the dimensionless number relation for continuous v/s dispersed phase Yagodnitsyna et al. [7] plotted a universal regime map ( Figure 11A) for main channels. The map consists of Weber number times Ohnesorge number for dispersed phase as its ordinate while Weber number times Ohnesorge number for continuous phase is its abscissa. The map covers all the flow regimes of the main channels namely parallel, plug, slug, droplet, rivulet, serpentine flows. The authors claim that this map can predict the flow regimes for any type of system and flow in micro-channels. The authors have presented their data as well as the data from Zhao et al. [16] for this regime map. Similarly, for micro-channel junctions (T and cross T junctions), Wu et al. [11] have predicted the flow patterns as shown in Figure 11B. The regime map is plotted with Weber number of continuous phase (We c ) as its ordinate while Capillary number for dispersed phase as its abscissa. The regime map ( Figure 11B) shows the transition lines of tubing-threading, dripping and dripping jetting with good predictive accuracy. However, physical properties have not been accurately captured in this regime map. Zhang et al. [12] carried out experimental investigations for flow patterns and mass transfer using Toluene-water; toluene-sulfuric acid and ethylacetate-water systems. The authors also carried out measurements to study the effect of temperature on the performance of micro-channels for the different systems. The fluid systems were chosen so that they could simulate reacting processes like nitration, sulfonation and acid catalysed reactions and have different fluid properties enabling parametric studies. The authors observed a lot of change in mixing behaviour for the operating range considered. Composite flow regime maps ( Figure 11E) were presented in the form of Ca c Re c 0.5 v/s Ca D 0.7 Re D 0.5 . The experimental data as well as data from literature was compared with the regime map and good prediction of flow patterns was found.

Studies in Pressure Drop in Microchannels
In this section, we present single and two phase pressure drops in the form of measurements and analytical models predicted by different authors [31,[36][37][38]. Table 2 explain the studies of pressure drop. Studies in pressure drops have been extensively covered in detail in review by Al-Azzawi et al. [5]. Hence, we try to cover this section briefly.
Experimental measurements of single phase pressure drops in micro-channels were carried out for 3 liquid mixtures Fu, [36]. The investigations were carried out for uniform, converging and diverging sections considering the T and Y junctions that are present in all micro-channel systems. The authors observed a linear relationship between pressure drop and volumetric flow rates for uniform, converging and diverging micro-channels. The maximum pressure drops were highest for diverging Y micro-channels (~100 kPa).
Kashid and Agar, [38] carried out two phase pressure drop measurements in Y-junction micro-channels for a set of different flow velocities in slug flow regime. The system chosen was Water-cyclohexane. The total pressure drop was given by the sum of hydrodynamic and capillary pressure drop. The hydrodynamic pressure drop was given by the Hagen Poiseuille equation while the capillary pressure drop was given by Young Laplace equation. The analytical model equation for total pressure drop of the main channel are given in Table 2. The pressure drop is a summation of individual pressure drops of water and cyclohexane multiplied by its volume fraction. As per our earlier discussion, it has been observed that various researchers have found film formation in these type of microchannels. Hence, the authors applied two approaches for pressure drop, one with film and one without film. The variation of pressure drop was carried out with different parameters like inlet flow rate, flow ratio and slug flow velocity.
Su et al. [23] carried out experimental investigations of pressure drop in microchannels and observed that hydrodynamics was characterized by Poiseulle number which is product of friction factor and Reynolds number (Re). Pressure drop was investigated with single phase fluid flow with different viscosities. The authors concluded that conventional theory could predict the pressure drop in microchannels with good accuracy.
Tsaoulidis et al. [31] carried out measurements for the two phase pressure drop for different volume fractions of IL and compared it with the single phase IL pressure drop. The authors found that the two phase pressure drop was always higher than the single phase pressure drop of IL with identical flowrates considered for comparison. Further, the two phase pressure drop increases with increase in volume fraction of IL. The trend of pressure drop however depended on material of construction of the channels. When glass micro-channels were used it was observed that at lower flow rates the two phase pressure drop was lower than the single phase pressure drop. Authors used the correlation by Jovanović et al. [39] for comparison with their measured values.
Ladosz and von Rohr [37] proposed two models to calculate pressure drop of liquidliquid slug flow in square micro-channels. Two different models for pressure drop were proposed (a) Moving film model and (b) no film model. The authors approximated a droplet surrounded by liquid film as annular flow for modelling the slug flow in the presence of film. The authors found good agreement of both models with toluene water system while for silicon-Oil-water system poor agreement was obtained. Pressure drop for two channel widths (0.2 mm and 0.4 mm) were measured and the analytical model developed ( Table 2) was validated. The authors found good agreement for two phase measurements with water-toluene system while with silicon-Oil-water system poor agreement was observed.
Assumptions: 1. Constant dynamic contact angle 2 Constant slug length 3 Equal number of slugs of water and cyclohexane under similar operating conditions. 4 Neglecting end effects 5 For uniform micro-channels, fully developed velocity profiles are obtained at x/L = 0.3. 6 In the flow rate calculations, water injection into the channel was assumed to have a negligible effect on the single-phase flow rate of oil through the needle 7 Hagen-Poiseuille equation was assumed to be valid 8 the dynamic viscosity of the mixture depends on the viscosity of each phase, their volume fractions, the temperature, and the magnitude of the dispersion of one phase in the other. Limitations: 1 CFD simulations for the geometries considered for experimental conditions have not been carried out 2 A single phase pressure drop instead of a two phase pressure drop was used for validation with experiments Findings: 1 An expression for two phase pressure drop across a Y-junction micro-channel for slug flow conditions with and without film formation 2 Power required per unit volume for most efficient conventional extractors was found to be around an order of magnitude higher than the micro-channels 3 An approximate single phase pressure drop model for liquid-liquid mixtures was derived and validated with experimental results 4 Analytical velocity and pressure profiles for converging and diverging sections were derived and compared with numerical (CFD) results. A good agreement was observed 5 Inlet configurations (T and Y junctions) did not affect the flow patterns or pressure drop 6 The flow patterns depend on the material of construction of the channels. In glass micro-channels plug flow or intermittent flow with IL was observed while in Teflon micro-channels annular flow was observed at low volume fractions while plug and drop flow were observed at higher volume fractions 7 Material also plays an important role in the determination of pressure drop. Two different Teflon micro-channels exhibit different pressure drops due to different plug length at same operating conditions 8 Modeling film thickness is very important for good prediction of overall pressure drop. An improved film thickness model gave good predictions for pressure drop when compared with experimental results 9 Proposed two models for evaluation of two phase pressure drop in square micro-channels. No film model and moving film model 10 The models showed a system specific behavior with good agreement with water-toluene system while high deviations with SiliconOil-water system. Hence, the models work well with fluids with similar viscosities 11 Pressure drop increases with an increase 12 Lockhart Martinelli and homogeneous were used to interpret pressure drop in microchannels and it was observed that homogeneous models do not provide good fit with experimental data. 13 A new correlation having two parameters that depend on microchannel type have been introduced 14 Pressure drops depend closely on type of micro channel material, the type of liquid that is introduced first and liquid flow rates 15 Prediction of flow regime map and flow patterns 16 Pressure drop could be predicted accurately using Lochart-Martinelli model for the system considered by the authors 17 Pressure drop is highly dependent on slug size 18 It is important to correctly measure interfacial pressure drop 19 Film velocity has negligible influence on pressure drop.
The authors attributed this to lower film thickness of the film than the one assumed for analytical model development [41]. For both systems considered there was an increase of pressure drop with flow rate for both channels. For example, for silicon-oil-water system, the pressure drop increased from 50 to 330 kPa for corresponding inlet velocities of 0.01 m/s and 0.14 m/s for 0.2 mm channel diameters, the increase was restricted from 10 to 60 kPa for corresponding inlet velocities of 0.01 m/s and 0.14 m/s for 0.4 mm channel diameters. Similarly values for toluene-water systems for pressure drop were 30 and 225 kPa for velocities of 0.02 and 0.27 m/s respectively for 0.2 mm channel diameters and 5 and 30 for 0.4 mm channel diameters.
Studies on two phase pressure drop carried by authors [37,38] show the following findings: 1.
Two phase pressure drops involve hydrodynamic and interfacial/capillary pressure drops in microchannels 2.
Both experimental investigations and analytical models show that two phase pressure drop is higher than single phase pressure drop 3.
Modelling two phase pressure drop involves careful consideration of two aspects (a) film near the wall and (b) the movement of the film 4.
Two phase pressure drop depends on (i) the fluids/system under consideration (ii) volume fraction of dispersed phase (iii) material of construction of capillaries.
Most of the points mentioned above have been considered by the researchers [37] but the model presented by the authors is valid only for fluids with similar viscosities and show large deviations for fluids with large viscosity differences. Figure 12A-C shows the experimental data of variation of single/two phase pressure drops w.r.t Re for different authors. Figure 12A shows the variation of pressure drop as a function of Re (5 ≤ Re ≤ 120) for water as working fluid. The channel diameter is 675 µm. A linear relationship is observed. The pressure drop values range from 5 kPa to 100 kPa respectively. Figure 12B shows the two pressure drop variation with respect to Reynolds numbers. Systems with similar viscosities have been considered. Channel diameters range from 400 to 500 µm respectively. It can be observed that the measured pressure drops for [37][38][39] are similar till Re = 40 while deviations upto 50% between measurements of [37,39] are observed after that. Deviations upto 50% can be observed. Probable reasons are the type of junction being used and the difference in channel diameter magnitudes. Figure 12C shows the pressure drop of Water-IL system by Tsaoulidis et al. [31]. The Re is below 0.5 due to high viscosities of IL. The pressure drops are however higher ranging from 20 to 180 kPa. The researchers have validated their data with respective models as shown in Table 2. Deviations of 10-20% have been shown by the authors for their respective models. The pressure drop mode1 Jovanović et al. [39] was compared with another pressure drop model of Kashid and Agar [38] and found that the latter over-predicted the experimental data of former. Tsaoulidis et al. [31] have used both T and Y junctions but at very low Reynolds numbers. They have also developed their own model which was seen to show good agreement with experimental data with deviations of around 10%.

CFD Studies
The CFD studies in LL systems have been focused on flow patterns, namely, slug flow, parallel flow and transition flow. The effect of operating parameters has also been shown by Raj et al. [42]. However, as per Raj et al. [42], some flow patterns which were observed in experiments could not be reproduced in their 3D simulations. Pressure drops using two phase CFD simulations have not been reported based on the simulation results in the literature.  [39] in Y junction micro-channel of diameter 498 µm Water-Toluene system by Ladosz and Rudolf von Rohr, [37] in Y junction micro-channel of diameter 400 µm (C) Variation of pressure drop for Re < 0.5 for Water-Ionic-liquid system for channel diameter of 270 µm by Tsaoulidis et al. [31].
In the upcoming subsections, different methods to handle deformable interfaces in two phase flows have been described. The studies of two systems mainly water-cyclohexane and water-toluene [14,[42][43][44] have been presented. Tables 3 and 4 gives the CFD details of prominent research works. Different strategies to handle moving interfaces are classified into moving mesh methods, fixed mesh methods and a combination of both. In the first method, a moving mesh is used to track the interface. As the topology of the interface changes, the mesh is adjusted in accordance with the change. Methods falling under this strategy are called front tracking methods. The second approach used a fixed (Eulerian) mesh, while the interface is tracked using different procedures, e.g., special markers or functions. Such methods are called front capturing methods; they include marker-and-cell (MAC), surface capturing and volume capturing approaches.
In front tracking methods, the interface is resolved directly. It always lies between two neighbouring mesh elements. As the interface moves, the neighbouring mesh elements deform in order to track the movement of the interface. For the description of the interface movement, an additional explicit equation is necessary. Usually the latter can be derived from an integral consideration of transport phenomena around the interface. It is worth mentioning that, in front tracking methods, each phase is considered separately, by using an individual set of equations describing the flow field. The front tracking methods are extremely accurate and robust, but rather complex to implement. Difficulties arise when two additional equations are needed to describe the mesh movement in three dimensions, because the second equation cannot be readily found in terms of the physical description of the system under study.
One of the oldest methods to describe interfacial movement is the marker and cell (MAC) method (Harlow and Welch, [45]). Here, massless particles are used to capture the motion of one fluid and thus the motion of the interface. Although its implementation is easy and it is capable of resolving the interface accurately, it fails to describe significant deformations of the interface. This is mainly because of the necessary redistribution of the markers, due to deformation or, in case where the interface expands, the required addition of new particles. Besides, this method can be very demanding in terms of computing power and time for 3D simulations.
One of the most popular surface capturing method is the level set (LS) method [46,47]. Here, a function is used to locate the interface, which takes positive and negative values on different sides of the interface and zero at the interface. The interface is therefore called zero LS. The surface tension is modelled using the method by Brackbill et al. [48]. An overview of the different LS methods is presented by Osher and Fedkiw, [49]. This method is conceptually simple and easy to implement. Its main drawback is the loss of mass (or volume), especially for significantly deformed interfaces van sind Annaland, [50].
The volume of fluid (VOF) method introduced by Hirt and Nichols, [51] represents a typical volume capturing method. The basic idea of the VOF method is the definition of a volume function which takes values between zero and one for the cells containing the interface.    A detailed review on the VOF methods is published by Rider and Kothe, [60]. This method is capable of handling problems with significant interface topology change and does not suffer from mass (or volume) losses. The extension of the method to 3D simulations is straightforward and no specific algorithm is needed for the case of merging or break up of the interface. However, the interface is smeared out and merges "numerically". The main drawback of this method is the inherent numerical smearing.
An alternative method is suggested by thegroup of researchers [61,62]. They used a combination of both front tracking and front capturing methods, whereby a fixed grid is used to describe the motion of the fluid flow and another moving grid with a low dimension is used to track the interface. Since for both phases a fixed grid is used, they are treated together, by solving a single set of governing equations for the whole flow field. The complexity of this method is due to the necessary dynamic remeshing of the moving sub-grid and mapping of the data coming from the moving grid to the fixed one. Additionally, for the merging and breakage of the interface, a special sub-grid algorithm is required. Most researchers have used the volume capturing (VOF) or surface capturing (LS) methods with a combination of both. The equations and boundary conditions used for VOF and Levelset equations along with continuity and momentum equations have been given below.
Continuity Equation: Momentum Equation: where σ = surface tension, κ n = surface curvature at the fluid interface computed from the divergence of the unit surface normal. F in the above equation is defined as: n asn = n |n| , n = ∇∅ = 0 VOF Equation: The boundary conditions used for most of the research articles considered in the present review are as follows: 1.
Outlet-Constant Pressure The mathematical equations used are as follows: Inlet: For no slip and zero flux: Outlet: The parametric sensitivity carried out by various authors in the subsections below:

Effect of Channel Geometry, Fluid Properties and Flowrates
Han and Chen [58] carried out 3D numerical simulations of merged droplet formation in squeezing regime in a double T-junction. Authors investigated the effects of flowrate of continuous phase, continuous phase viscosity and interfacial tension. The authors found that the droplet formation involved three steps: filling of the dispersed phase flow in the T-junction, blocking of the T-junction and eventually breakup of the droplet. During this period there was a periodic pressure cycle resulting in pressure drop. The pressure cycle period was equal to the droplet breakup period. The authors found that when Ca number increased the slug lengths decreased while when the viscosity increased the droplet breakup time decreased and so did the pressure cycle time. However reverse trend was observed when surface tension was increased.
Yang et al. [56] carried out 2D CFD simulations in OpenFoam Opensource CFD code in a system consisting of a main channel and system of capillaries for seperating dispersed phase. Test system of Water-Toluene was used for this purpose. The seperation was made in such a way that when dispersed and continuous phases travel together the dispersed phase is drawn into the side capillaries and complete seperation occurs by applying high pressure at the continuous phase outlet (or the main channel). Thus the dispersed phase passes through the capillary. The authors defined the pressure range in which the system can be operated by calculating lower limit of pressure drop by Hagen-Poiseulle equation and the upper limit by the Young-Laplace equation.
Sripadaraja et al. [59] carried out 2D numerical investigations for squeezing regime (0.0028 < Ca < 0.008) and squeezing to dripping regime. Three different geometries namely, cross junction, coflow and T-junction were considered for investigation. The authors found that the effect of surfactants on droplet generation is nominal. With increase in flow velocity, Ca increases. For Ca = 0.002, droplet generation increased by 16%. Similar to Han and Chen, [58], the authors defined a 3 step approach for droplet formation which the authors named Filling, necking and lag stage. Lan et al. [55], carried out 3D numerical simulations for two different geometries Co-tube microchannel and Flow focusing shape microchannel. A modified Level-Set method was used to resolve the interface movement due to the mass loss problem in original Level-Set method. The authors observed decrease in droplet sizes with increase in flow. Further, strong internal circulations were observed during the initial stages of droplet formation resulting in increase in transport property rates. However, the strength of the circulations decrease after droplet formation.

Effect of Flow Regime
Azarmaneesh and Farhadi, [22] studied droplet formation for 1 < Re < 25 for sqeezing regimes and stable and unstable jetting regimes. New vortices were created after droplet detachment when Re < 1. Droplet formation steps during squeezing regime were shown for Re = 0.8 and Ca = 0.003. At dripping regime the droplet formation was seen to be affected only by balance of shear stress and surface tension while inertia forces also influence the motion, changing droplet size and some new vortices in droplet. Kostinont et al. [54] carried out numerical investigations using the properties of glass. The channel structure was modifed by having guidelines at the centerline of the channels so that the interface is always at the center. The authors tried to get transition of parallel to slug transition by increase in flow velocities thus affecting Re numbers. DeMenech et al. [57] numerically investigated the mechanism of droplet breakup in squeezing, dripping and jetting regimes using the PhaseField method. The breakup process was found to be driven chiefly by buildup of pressure upstream of an emerging droplet. The capillary number was found to have a negligible effect on droplet size and dynamics of breakup. The dripping regime was influenced by the contrained geometry and modified the scaling law derived for size of droplets from interfacial and viscous stress balance. The jetting regime was found to set at high flow rates and higher values of Ca number.
The research works for CFD studies have considered flows in the slug flow regime. A numerical model was developed to study the slug flow patterns under different flow conditions and fluid properties Kashid et al. [14]. The authors validated their model with the experimental results. The authors found that CFD can be a very efficient tool for design of micro-extractors. Further, internal circulations inside the slugs (dispersed phase) as well as outside the slugs (continuous phase in between two slugs) were studied.
VOF based CFD methodology to simulate slug flow in micro-channels was used by some researchers Kashid et al. [44]. The authors considered two cases (i) when there is no film between the wall and the slugs and (ii) when the slug is seperated from the wall by a liquid film situations of slug flow were considered and validated. The authors found that VOF method described slug flow generation and its characteristics very well. The slug flow generation was found to depend on slug flow velocity and interfacial velocity. The authors also concluded that slug length can be preserved upto a channel diameter of 1 mm.

Effect of Channel Wettability
Yamasaki et al. [53] studied stable parallel flow for Water-Octane system using the VOF model coupled with Continuum-Surface Force (CSF) model. The Capillary number and Weber numbers chosen for these studies were far less than 1. The authors tried to simulate parallel flow which would aid the seperation of the two layers easily and reduce extra downstream equipment for seperation. 3D CFD studies of drop/slug formation in micro-channels were undertaken Raj et al. [42]. The drop-shape slug shape, slug lengths, slug formation was predicted with Volume of Fluid (VOF) interface tracking method along with Navier-Stokes governing equations for continuity and momentum and compared with measurements of Garstecki et al. [33]. The authors developed a CFD model to see the effect of wall adhesion on flow regimes and slug lengths. The authors claimed the model to be robust and that the model can be used to simulate mixing, transport processes and chemical reactions. However, the authors suggested that extensive and reliable measurements need to be carried out for simulating and predicting film thicknesses correctly. This is due to the effects of grid refinement near the walls, surface wettability and roughness of film thickness critically. Further, effect of surface tension at low velocity ranges have been found be independent on droplet diameter Sugiura et al. [32].
CFD studies were performed by Fu, [36] to understand the velocity and pressure distribution in the configurations of the micro-channels considered. The numerical velocity profiles for different horizontal locations were compared with analytical solution. For x/L = 0.6 and 0.4 for converging and diverging micro-channels results were in good agreement. Similarly, pressure distribution profiles were also compared with analytical solution. The converging pressure fell slowly near the entrance while decreased rapidly near the outlet. For diverging section, pressure fell rapidly near the entrance and slowly near the exit.

Influence of Flow Patterns on Transport Property
In this section, an effort has been made to summarize the findings of various authors specifically, the geometric parameters (like size of microchannel), operating parameters (flow velocity) on VMTC. While for most of the authors it was slug flow the major findings lie in the systems chosen for studies and parameteric studies by various authors.

Experimental Studies
Major amount of experimental work by researchers [9,10,16,[63][64][65][66][67][68][69][70][71][72] have been carried out in slug flow regime as seen from Table 5. The information on the geometry, operating parameters, experimental techniques, flow regimes of major LL systems, the findings, limitations and assumptions of the research works have been presented in Table 5. 3.1.1. Effect of Geometric, Operating Parameters, Wallfilms and Fluid Properties on VMTC Kashid et al. [1] carried out studies in mass transfer systems for three different systems 1. Water-iodine-kerosene; 2. Kerosene-acetic-acid-water and 3. Water-succinic-acid-nbutanol at constant flow velocities. The authors found that the VMTC's varied from (0.32 to 0.98 s −1 ) for channel diameters of 500 µm while it dropped to (0.13-0.32 s −1 ) for channel diameters of 1000 µm for water-iodine-kerosene system. The authors found that internal circulations increase with flow velocity. Hence, VMTC increases due to higher internal circulations. The authors, also, found that flow ratio (q) has a positive effect on VMTC. For water-succinic-acid-n-butanol system the extraction efficiencies decreased by around 10% with increase in slug flow velocity for constant values of internal diameter of the capillary. Zhao et al. [16] have conducted experimental investigations with water/n-butanol/succinic acid with micro-channels having internal diameter of 300-600 µm and length 45 mm. The authors observed an increase of VMTC values to about 4 times when the channels size was decreased from 600 µm to 300 µm. Further, this was observed for Re numbers in the range of 500-600. Other researchers, [6,9] have operated their respective systems for low Reynolds numbers (Re < 20) for toluene-trichloroacetic acid-water system for channel diameter of 400 µm and found that the VMTC's less than 1. Other researchers, Di Miceli Raimondi et al. [70] have carried out measurements of VMTC's in diameters of around 300 µm. The VMTC's obtained by the authors were around 8.44 s −1 for 210 µm while it decreased to around 2.65 s −1 for 300 µm micro-channels. The Re range in such systems were less than 50.
Dessimoz et al. [9] have used T-shaped and Y-shaped micro-channels. The authors used T-shaped micro-channel for generating slug flow while Y-shaped micro-channels were used for generating parallel flow. The authors have used hexane-water system for generation of slug flow while they added sodium hydroxide in the aqueous phase. Hexanetrichloroacetic acid-water-sodium hydroxide system was used due to higher interfacial forces responsible for the generation of slugs. Also addition of sodium hydroxide decreased the viscous forces and interfacial forces further dominated the system due to which slug flow could be formed. Xu et al. [69] carried out investigations on alkaline hydrolysis reactions for T-junction for slug flow (similar to Dessimoz et al. [9] and found increase in VMTC's with increase in NAOH concentrations Ghaini et al. [66] studied the effect of wall film on transport property. The authors used the fast instantaneous reaction. Matsuoka et al. [73] studied the influence of wall films during slug flow in micro-channels. The authors found the variation of VMTC with Reynolds number for a fixed velocity and different hydraulic diameters (0.6 mm to 2 mm). Maximum VMTC obtained was around 0.7 s −1 for a constant channel length. The VMTC was found to be a strong function of slug flow velocity and decreased due to resistance due to wall film.
Di Miceli Raimondi et al. [70] carried out measurements for dispersed phase VMTCs using the water toluene system. The maximum inlet velocities measured were up to 300 mm/s and high values of transport properties are obtained.
Arsenjuk et al. [71] carried out experimental investigations for a kerosene-glycerolwater system to see the effect of film thickness on transport property. The authors found that film thicknesses were higher for systems having a transfer component than systems not having transfer component. This influences the transport property. This study deserves special significance since the presence of transfer component affects two phase investigations and affects the flow properties.
Sattari-Najafabadi et al. [74] have also performed investigations to study the effect of shape of microchannels on the VMTC. Sattari-Najafabadi et al. [74] have found the impact of geometric configuration i.e., change in channel width (W) and type of junction. The authors used the T-junction, cross junction and cross T-junction to understand the effect on VMTC. Further, the channel width (W) to depth (H) ratio (W/H) was varied (1 ≤ W/H ≤ 2). The flow ratio (q) was varied from 0.5 to 2. The authors concluded that best VMTC's can be obtained for a cross-junction with flow ratio of 0.5.
Sattari-Najafabadi and Esfahany, [75] carried out experimental investigations to find the effect of surfactant on VMTC during two phase flow of a LLE system. The authors observed that VMTC obtained in presence of surfactant was 1.44 times better than VMTC without surfactant.

Comparison of VMTC's of Microchannels and Conventional Equipments
Singh et al. [76] carried out experimental investigations to compare the performance of different microchannels with conventional extraction equipments. The authors considered three different micro-channels namely T-junction micro-channels, serpentine microchannels and split and recombine micro-channels for experimentation. The authors observed that the process intensification possible by using micro-channels has the potential to reduce solvent inventory 10 times and correspondingly reduce health hazard and plant foot print. Process intensification can occur only for shorter contact times since for longer contact times results were similar for conventional extractors and microchannels.
Sahu et al. [77] carried out experimental investigations and compared the performance of continuous micro-channels v/s batch systems. Experiments were performed both in batch and micro-channels. They considered square Y and serpentine micro-channels for experiments. The authors observed superior performance of micro-channels than batch systems. This is attributed to strong internal circulations induced by shear between slugs and wall.

Evaluation of VMTC's for Systems Involving IL's, DES as Solvents
Bai et al. [10] considered water-Rhodamine6G-ionic liquid system to find the VMTC in microchannels of width 1 mm and depth 0.3 mm. Two types of geometric configurations mainly crossing junction and T-junction were used for investigation. Empirical correlations to predict VMTC were obtained from experimental data. The VMTC was found to be a strong function of flow ratio and Reynolds number. The measured and predicted VMTC's matched well with less deviations.
Tsaoulidis et al. [31] carried out studies of Dioxo-Uranium VI extraction from aqueous nitric acid using ILs. T-junction micro-channels were used. Slug flow was found to be the flow pattern during their study. The transport property/efficiency of the microchannel was influenced by the presence of low initial concentration of nitric acid.
Li and Angeli, [72] carried out experimental investigations for extraction of lanthanides. The authors used IL's to extract lanthanides from nitric acid. The authors observed that the parameters that affected the VMTC were flow ratio, length of microchannel and initial concentration of nitric acid.
He et al. [78] have carried out investigation on extraction of rare earths Pr/Nd from chloride solutions with Deep Eutectic Solvents (DES) using Y-junction serpentine microchannels. Jiang et al. [79] have carried out separation of lanthanum (La) and Cerium (Ce) from chlorides in serpentine micro-extractors and rectangular cross-sections. Sen et al. [80] used serpentile micro-channels for extraction of Uranium from Tributyl phosphate in dodecane (aqueous phase). Annular flow was obtained for phase ratio of 2 while after 2 the flows were in droplet regime. All the investigations involving rare earths had Re values in the range Re < 25 and VMTC values ranged from 0.05 to 0.17 s −1 .    Findings: 1. Liquid-liquid extraction is carried out in capillaries and VMTC's were found to be two orders of magnitude higher than conventional equipment's. 2. The improvement of mass transfer rates is due to two phenomena. A. the large interfacial areas which enhances the mass transfer rates and b. the internal circulations inside the slugs which promote surface renewal at the phase interface. 3. The study showed that capillary contactors can be a powerful tool for laboratory studies of liquid-liquid extraction 4. VMTC's have a very strong dependence on parameters like slug flow velocity, flowrates of organic and aqueous phases and capillary diameter. The VMTC's are constant after a certain slug flow velocity is reached. 5. VMTC's decreased with an increase in capillary diameter 6. Inlet locations play an important role in the mass transfer process 7. Decrease in height/depth of the channel causes a increase in VMTC 8. Correlations of VMTC's were developed 9. Found regime transition maps based on Capillary and Reynolds numbers for parallel and slug flow systems 10. The authors found that the main parameter for the determination of VMTC's for parallel flow was the initial concentration of water while for slug flow it was the total linear velocity of fluids (dependent on total flow rate through the channel) 11. The authors focused on the wetting behavior at the capillary walls due to which the true interfacial area is different than the one determined by the conventional snapshot method of determination. Hence, chemical methods need to be applied to such systems 12. Found the interfacial area for systems with and without wall films. Effect of slug flow velocity on interfacial showed that for systems with film 13. An inert carrier gas can be used to enhance mass transfer under certain circumstances and to tune the separation behavior of a capillary separator 14. Secondary flow patterns are developed due to the inert carrier gas in addition to the water slugs that exist in liquid-liquid flow which leads to an enhanced surface area between liquid phases 15. Extraction efficiency and VMTC decreased with increased volume flow rate ratio, total flow velocity, and the inner diameter of the main channel. 16 Initial concentration of aqueous phase has an effect on extraction efficiency. 17. VMTC decreases with channel length for an alkaline hydrolysis reaction for a micro-channel with slug flow 18. Overall extraction rate is a linear function of VMTC and reaction is mass transfer controlled 19. Slug velocity is an important factor for enhancing the VMTC due to increase in internal circulations within the slugs 20 Determination of VMTC's by experimental measurements and validating a correlation obtained from numerical simulations 21 The flow patterns inside the droplets confined to the walls causes enhancement in mass transfer compared to droplets that are not confined 22. Different flow patterns slug, annular, dispersed and inverted dispersed flow patterns were observed for experimental investigations of liquid-liquid extraction systems for three different flow orientations 23. VMTC decreases with increase in capillary diameter 24. Vertical downflow orientation gives higher VMTC's for both lower and higher diameter capillary than the other orientations 25. Suggested an innovative co-current tubular contactor that is trumpet shaped with short horizontal section with long vertical up and downflow sections 26. Compared the performance of micro-channels with conventional extractors and found that micro-channels have the potential to reduce solvent inventory and plant area due to miniaturization 27. The wall film in liquid-liquid slug flow is continuously renewed and thus actively participates in mass transfer. Hence this should not be neglected while designing micro-reactors with liquid-liquid slug flow28. Experimental investigations in micro-channels with two different inlet junctions (T-type and cross type) showed that droplet based microsystems provides higher VMTC's between water and IL for micro-channel configurations 29. VMTC at droplet formation stage is 3-4 times higher than at the latter stage due to strong disturbance inside droplet during formation 30. Found that under stable Taylor flow conditions in a 1 mm inner diameter glass channel, the organic phase (dodecane) formed liquid droplets and the aqueous phase(water) formed liquid slugs with a liquid film enclosing the organic phase droplets 31. Compared the performance of micro-channels with batch extractors and found that micro-channels can reduce the residence times to significant extent and make the processes continuous 32 have not been carried out to compare the flow patterns 5. The article does not consider the variation of VMTC due to change in channel diameter and type of inlet mixing (restricted to T-junction inlet). 6. The article is restricted to slug flow patterns in microchannels 7. The article is restricted to experimental investigation of effect of wall film on mass transfer between continuous and dispersed phases 8. The article is restricted to droplet flow regime and droplet formation only 9. The study does not consider the effect of wall film between the slugs which might affect the VMTC 10. The present article is restricted to prediction of flow patterns while measurement of VMTC's for these systems and operating parameters were not performed. 12. Pressure drop studies when the VMTC's were increased due to flow rates are not reported in the present study. Assumptions: 1 Outlet flow from steel outlet would be higher than the one from PTFE outlet due to higher affinity of aqueous phase for steel than affinity for organic phase of PTFE 2 The combination of a flow splitter and capillary has been assumed to act as a mixer-settler device as in conventional liquid-liquid systems 3 Properties of two immiscible liquids are constant for all the experiments 4 For cases where there is no formation of wall film, the slugs were considered to be of rectangular shape for calculation of interfacial area 5 Fluid is Newtonian and incompressible.

Comprehensive Regime Map for LL Systems Having Transferable Component
In Section 2.2, the regime maps presented by various researchers for two phase systems have been discussed. The flow regime maps are based on combination of dimensionless numbers which is the interplay between the inertial, viscous and surface tension forces. However, in Section 3.1. Arsenjuk et al. [71] have found that the presence of a transferable component affects the dynamics of two phase LL systems. The product of Weber number and Ohnesorge number for both continuous and dispersed phase seems to be a good way to represent the flow regimes. An effort has been made to collate the data of the LL systems having a transfer component and see if they adhere to the regime map provided by [7]. Figure 13 shows the plot of flow regimes. Yagodnitsyna et al. [7] (Figure 7) have shown a very comprehensive universal flow regime map which however deviates from the findings of the various authors in slug flow regime. During our analysis, we find a disparity in flow patterns shown by several authors. The present analysis shows that many of the literature works fall in the area of parallel flow in the regime map of Yagodnitsyna et al. [7] while in experimental works of researchers, they are in slug flow Kashid et al. [6].   Further, the flow pattern map of Zhang et al. [12] takes into consideration different powers of Reynolds and Capillary numbers and also claim to have a good fit of the literature.

CFD Studies of Influence of Flow Patterns on Transport Property
In this section, we investigate the mass transfer studies carried out by various researchers. In this section, the approaches available in the published literature for solving coupled equations have been presented. They are a. Moving boundary method and b. Front capturing method. Table 6 provides the geometric details and operating parameters of the research works while Table 7 provides details of CFD softwares, mesh, methods used for the research works considered.
A typical schematic of the moving boundary approach has been shown in Figure 14. The velocity of the moving wall is same as the slug velocity. Equal stress condition across the interface is considered for the momentum. Convection diffusion equations are solved for each phase. No flux boundary conditions is implemented across moving walls while the flux continuity and thermodynamic equilibrium boundary conditions are implemented with the concentration equations. To evaluate the VMTC, the simulation process decoupled the hydrodynamics of twophase flow from the mass transfer phenomena and was then divided into two steps. First, the steady-state flow field was established, considering the velocity to be time-invariant. Then, all the components in the solution were initialized in the organic and aqueous phases, respectively, and the flow time was reset to zero. The mass transfer equations were calculated with a time step of 0.001 s and 10 iterations at each time step by analysing data obtained with different time steps and iterations. For the initial conditions, the flow field at t = 0 s calculated in the first step of the computational procedure is the steady-state flow field. Since species 1 is only patched to a certain concentration value in phase 1, the appropriate initial and final boundary conditions of the concentration field are provided.
Many industrial separation processes occur in fluid-fluid systems including droplets/ bubbles in a continuous phase (e.g., solvent extraction, gas scrubbing, waste water treatment in bio reactors). For these processes interfacial mass transfer represents a crucial phenomenon. To describe it properly, a coupled problem including momentum and mass transport at and around moving interfaces must be solved. Such problems are very difficult, since the deformation of the interface influences interfacial mass transfer related to the handling of boundary conditions at the moving interface. Due to the assumption of the thermodynamic equilibrium, a concentration jump appears at the interface which is numerically demanding and difficult to implement into the front capturing methods. Besides, the component flux continuity should also be fulfilled at the interface. There are just few papers dedicated to this problem.
Bothe et al. [87] used the VOF method to model a gas-liquid system (air bubbles rising in water). They transformed the mass transfer equation to    The flow is two dimensional 5 Flow is fully developed and liquids are immiscible with each other 6 Liquid-liquid interface is assumed to be flat for the sake of simplicity 7 Marangoni stresses at the interface are neglected and flow field is assumed to be independent of concentration field and the interfacial flux continuity boundary condition to D l ∇ C l ·n = HD g ∇ C g ·n (13) where n is the unit vector normal to the interface, C is concentration, → u is velocity, D is molecular diffusivity and where indices l and g denote gas and liquid phases and H is the distribution coefficient (Henry's constant).
Yang and Mao, [92] carried out CFD simulations using the LS method for a rising droplet in a continuous liquid phase (water-succinic acid-butanol). The equation for mass transfer was modified in their work and is given below Here, H is the distribution for the liquid-liquid system considered and H ε (∅) is the regularized Heaviside function to avoid the numerical instabilities at the interface. It is given by Haroun et al. [93] suggested a VOF-based modelling approach for reactive absorption systems. They considered the film flow of a liquid flowing along a corrugated surface in counter-current with the laminar gas stream. The mass transport equation has been modified to: where Here, j indicates component number, f is volume fraction, W j is the production term related to chemical reaction. Haroun et al. [93] claimed that Equations (21) and (22) account for the jump conditions at the smeared-out interface, where the volume fraction gradients are usually very steep.
Kenig et al. [94] performed CFD simulations with coupled mass and momentum equations were solved with the incorporation of the boundary conditions in the concentration equation.
∂C d ∂t ∂C c ∂t where, subscript c denotes continuous phase while subscript d denotes dispersed phase, H D denotes the distribution coefficient, α 1 , α 2 are parameters, n is the unit vector normal to the interface. The values of α 1 , α 2 are adjusted in such a way that the boundary conditions are solved only at the interface. The authors have considered a droplet rising in an aqueous phase in a water-acetic-acid-toluene system and mass transfer from toluene to water. Two phase simulations with and without mass transfer have been carried out. CFD simulations using a mass transfer model that predicts mass transfer across the interface were carried out using COMSOL 3.5a. The model has been validated with liquid-liquid and gas-liquid systems with satisfactory qualitative results.
Following are the boundary conditions for moving boundary approach: 1.
Axial and radial velocities of the walls are equal to slug/plug velocity 2.
No slip boundary condition at channel wall 3.
Equal stresses at the interfaces between two phases 4.
Flux continuity boundary condition at the interface 5.
Zero flux set at the channel wall 6.
Periodic boundary condition to satisfy flux continuity at the interface.
In Section 2.2, CFD studies were discussed in which governing equations for mass and momentum were solved with interface tracking techniques like VOF. However, the most important parameter for transport property (mass transfer) is the VMTC and hence mass transfer from dispersed phase to continuous phase should be obtained by solving concentration governing equations. However, boundary conditions of thermodynamic equilibrium and flux continuity as given by Equations (2) and (3) should also need to be incorporated while the system of equations is being solved. The approaches mentioned in the previous subsection have been used by a few researchers and will be discussed in this section. Details of investigations are provided in Tables 6 and 7. Kashid et al. [65] performed simulations with the model provided by Yang and Mao [92] and used the VOF method for defining the interface. The authors observed internal circulations that keep on renewing at interface and increase the concentration gradient between two slugs which enhances the rate of mass transfer. The authors claimed that the model can be used for any liquid-liquid system with a fixed interface location. No free surface model equations were used for tracking the interface. To summarize, the authors have used the concept of moving boundaries for their numerical investigations.
Di Miceli Raimondi [89] have solved two dimensional simulations using their own inhouse code. The simulations were carried out in two steps. First, the hydrodynamic behaviour of the system is computed by means of the code. In the second step the concentration field of this solute present in the dispersed phase is computed by a code in Matlab 7.4. The numerical method used by the authors is an interface capturing technique without any interface reconstruction. The authors observed three flow structures (a) internal vortices in both continuous and dispersed phases and thin film that wets the wall, (b) several recirculation nodes inside droplets and (c) without any recirculation loops in continuous phase since droplet is not confined enough. The flow structures were found to be a function of Reynolds and capillary numbers. Several simulations were performed for different dispersed phase velocity, channel width, dispersed phase holdup. A correlation has been developed for the computation of VMTC based on the data. The authors have compared the correlation with other correlations [95,96] and the predictions were in good agreement.
Zhang et al. [86] performed CFD simulations to predict mass transfer characteristics with and without alkaline hydrolysis reactions during slug flow in micro-channels. The species concentration equation was used for computing the concentration. The moving boundaries approach was used. The predicted VMTC values were compared with experimental values and a close agreement was observed. The authors claimed that the model could be applied to various liquid-liquid systems with fixed interfaces. The authors also claimed that this model can be coupled with free surface models to capture complex hydrodynamics near the interface.
Cito et al. [90] studied wall mass transfer rates generated by capillary driven flow in a micro-channel. Simulations were performed for 2D rectangular and circular microchannels. The effect of Reynolds number, channel geometry and contact angle on average wall mass transfer rates were observed. A correlation for Sherwood number Sherwood and Wei [97] relating the above parameters was developed from the simulated data. The effect of wall curvature was considered negligible in the present work due to high value of Schmidt number. For Re < 500 the Sherwood number had a weak dependence on contact time. The dependence on time corresponds to unsteady mass diffusion problems.
Tsaoulidis and Angeli, [82] carried out experimental investigations for extraction of Uranium in straight microchannels. The internal diameter was varied from (0.5 ≤ D ≤ 2 mm). The authors found the VMTC in small channels was in the range 0.049 to 0.29 s −1 . A correlation for VMTC as a function of Re, Ca and (Internal dia (ID)/length of channel (D/L)) was developed and predictions were compared with the experimental measurements as well as other models available in the literature. The authors also investigated the effect of input parameters like internal diameter and mixture velocity on efficiency and found that output parameters like efficiency decreased with increase in ID. The authors have also carried out numerical simulations with commercial software COMSOL 4.3. The domain taken by the authors is a dispersed elongated structure which the authors mention as plug (instead of the general definition of slug by previous authors) while the surrounding continuous phase fluid is mentioned as the slug. Momentum and continuity equations were solved for this domain along with convection-diffusion equations. The interface is considered stationary and a moving wall boundary condition has been provided. The numerical results for the velocity fields were compared with velocity patterns from Particle Image Velocimetry (PIV). Images captured by High Speed camera were used to calculate output parameters like the film thickness, length of unit cell, plug length, plug velocity. The data generated from the images were used as boundary conditions of numerical simulations. Extraction efficiencies were compared with numerical results and a good agreement was obtained.
Gómez-Pastora et al. [88] performed experimental and numerical investigations of solute mass transfer in Y-Y micro-channels as a function of coupled influence of both flow patterns and mass transfer kinetics. The influence of ratio between residence and diffusion times and volumetric ratio between fluid phases was determined for 3 geometric variations Symmetric device: Y junction inlet and outlet with equal inlet and outlet diameters for phase 1 and phase 2. Configuration 1: Y junction inlet and outlet with inlet and outlet diameter of phase 2 four times higher than the ones of phase 1 3. Configuration 2: Y junction inlet and outlet with inlet and outlet diameter of phase 1 four times higher than the ones of phase 2. Figure 15 shows the schematic of the micro-channels used by the authors. Two factors (a) concentration factor (defined as the ratio of outlet concentration in phase 2 to outlet concentration of solute in phase 1) and separation factor (similar to conversion and defined as ratio of difference between inlet and outlet concentration of phase 1 to inlet concentration of phase 1) were used to quantify performance. The authors found that separation factor influenced the values of solute concentration in both phases which influences mass flow of solute and channel outlet. Two configuration of Y-Y channels were suggested for different applications. The authors recommended that systems with separation factor less than 1 can be used for applications where high solute flows are required (configuration 1 by authors can be used) while separation factors greater than 1 can be used for processes where high concentration factors are needed (configuration 2 can be used).

Influence of Flow Patterns on Transport Property: Heat Transfer
In the present section we briefly report the experimental and numerical studies on two phase liquid-liquid heat transfer. Liquid-liquid two phase flow is a preferred method over gas-liquid two phase flow due to the dispersed phase liquid having higher thermal conductivity and hence higher heat transfer coefficients. Further, in the slug flow regime, the continuous phase flow is disrupted by the presence of the dispersed phase which contributes to enhancement of heat transfer coefficient. Also, another reason for enhancement is the recirculations inside the dispersed phase. Major parameters influencing liquid-liquid two phase flow deals with influence of different parameters like fluid properties and geometrical dimensions.

Experimental Studies of Heat Transfer in Microchannels Involving Liquid-Liquid Two Phase Flow
Though plenty of studies in single phase systems in microchannels have been reported [4], a few experimental studies for two phase liquid-liquid systems have been carried out by researchers [98][99][100]. Systems used were water-light mineral [98]; siliconoilwater [99] and Water-hexadecane [100]. All the authors have observed slug flow in their investigations. Asthana et al. [98] have carried out experimental investigations in 100 × 100 µm six turn serpentine microchannel of square cross section. Increase in heat transfer due to two phase flow (water-mineral oil) of 4 times compared to pure water have been observed. The authors also observed a significant increase in pressure drop due to heat transfer. Giolla Eain et al. [99] studied the effect of decreasing slug lengths and increase in film thicknesses on heat transfer rates in a 1.5 mm circular channel under constant wall heat flux boundary conditions in an Re range 2.3 < Re < 92.1. The authors observed that heat transfer increased by decrease in slug lengths while decreased by increase in film thickness. Correlations for both heat transfer and film thickness as a function of dimension-less numbers (Ca and We) were proposed by the authors. Dai et al. [100] also carried out experimental investigations in a vertical 2 mm circular microchannel under constant wall heat flux boundary conditions in an Re range 15 < Re < 48. The authors found an empirical correlation for the same. Abdollahi et al. [101] carried out both experimental and numerical investigations in microchannels with hydraulic diameter of 2 mm, a oil-water system for constant heat flux boundary conditions. The authors determined correlations for friction factor and heat transfer coefficients and found an increase of 700% in heat transfer rates over single phase flow.

CFD Studies of Heat Transfer in Microchannels Involving Liquid-Liquid Two Phase Flow
A few numerical investigations have also been performed to study two phase flow. Vivekanand and Raju [102] carried out 2D axisymmetric numerical investigations in 100 µm channels using the VOF method. The investigation was carried out for both constant wall heat flux and constant wall boundary conditions. The authors observed increase in Nusselt (Nu) number upto 180% and 210% increase in two phase as compared to single phase. Urbant et al. [103] have found increase in heat transfer coefficients for water droplets in oil in circular microchannels of diameter 0.1 µm. The recirculations were seen in both phases resulting in the increase in HTC. Fischer et al. [104] have observed 400% increase in heat transfer rates for a Water-Silicon Oil system as compared to single phase fluids using 2D axisymmetric geometry of 0.1-1 mm diameter microcapillaries. Che et al. [105] did simulations with 2D geometry for microchannels and found that heat transfer increases with decrease in length of the slugs. Che et al. [106] carried out 3D large eddy simulations (LES) for 0.2-0.8 mm (characteristic length) microchannels for water-mineral oil system. An increase in Nu number was observed due to the two phase flow conditions due to recirculations present in both continuous and dispersed phase. Bandara et al. [107] carried out CFD simulations to investigate heat transfer characteristics (4.9 < Re < 21.9). The authors concluded that the heat transfer rates depended on the Ca numbers, slug sizes, film thickness and contact angles. The authors used correlations available in the literature (for GL flows) to find the Nu number and concluded that correlations needed to be developed specifically for LL flows. Presently two correlations for Nu number (one by numerical studies by Dai et al. [100] and other by experimental studies by Giolla Eain et al. [99] exist for LL flows. The equations are given in Table 8.

Other Aspects: Film Thickness, Peclet Number Effect, Systems Considered, Scope for Different Geometries
During heat transfer in two phase LL flows a thin film is formed between the slugs and the walls when flows are in slug flow regime. Techniques used in the literature consist of Micro-PIV [98] and optical macroscopy [99]. The heat transfer rate has an inverse relationship with film thickness and is decreased with an increase in film thikness. Most authors have used in the past the correlations of GL flows in LL flows to determine film thickness until Giolla Eain et al. [99] found a correlation for two phase flows based on their experimental data. The authors defined film thickness in terms of Ca number and We number δ In case of rectangular channels, heat transfer from the wall to center of liquid slug is further due to the space between the curvature of the liquid slug and the corners of the rectangular channel since heat transfer in this region occurs due to conduction only. Further, the heat transfer in LL flows will also depend on increase or decrease of flow which will alter the size of the slugs and hence increase heat transfer rates.
Another important parameter affecting heat transfer is the Peclet (Pe) number which is the ratio of heat transfer by convection to that by diffusion. Che et al. [105] in their numerical investigations found that increasing Pe number increases the Nu number but decreases the heat transfer process.

Authors Correlation
Dai et al. [100] 1 Nu FS Nu W = 5 + 1.5 Fluid properties and size and shape of the microchannels would play a vital role in heat transfer in LL flows. The literature review shows that mostly Oil-water systems have been used for both experimental and numerical studies. Shape of the microchannel would have an effect on the heat transfer process as seen in the description of rectangular channels. Further, 3D simulations would give more opportunities for studying the shape of the channels numerically.

Conclusions
The conclusions are summarized in Figure 16. The brief description is provided as follows: 1.
Flow patterns namely parallel, slug, annular, annular-slug, slug-dispersed have been reported by authors but some new patterns like rivulet and serpentine patterns have also been observed for liquid-liquid systems in main channels of microchannels 2.
Flow patterns in junctions (cross or T-junctions) are also summarized as tubing, dripping and jetting. Several works are available but a connection between the flow patterns in the junctions and the main channels needs to be made while creation of regime maps 3.
Regime maps have been presented by authors with dimensionless numbers like product of Reynolds numbers and Capillary numbers or product of Weber and Ohnesorge numbers. However, there is no consensus on the particular regime maps and no regime map is able to predict the works of the entire literature correctly 4.
Pressure drop studies are in the form of experimental investigations and analytical models. The pressure drop models are seen to predict well with the experimental data for similar viscosity systems 5.
Few CFD models focusing on flow patterns have been found for liquid-liquid systems. The models have shown good predictive capabilities with experimental results. However, Capillary number plays an important role in the predictions. For certain Capillary numbers the CFD models have been unable to predict the flow patterns as those of experimental observations suggesting more amount of research works to be carried out in this area.

6.
Mass transfer studies in experimentation have shown the effect of operating variables and properties and geometric dimensions. VMTC is found to increase with an increase in flow velocity, presence of inerts, presence of surfactants, increase in flow ratio. VMTC however decreases with increase in characteristic channel dimension (like channel width or diameter) for specific channel length. 7.
For cases of extraction of rare earths or uranium, the residence times are found to be less than 20 s which is an order of magnitude less than conventional batch or continuous processes 8.
The VMTC's for test systems involving conventional solvents or systems involving IL's or DES's are two orders of magnitude higher than the VMTC's of conventional equipment's. However, commercialization needs stacking up of a number of microchannels in series and hence needs substantial research in this direction 9.
For reactive systems, the measured MTC's are two to three times the ones predicted by film theory due to the enhancement in mass transfer during reaction 10. CFD studies have been carried out majorly with stationary interfaces with specified slug lengths obtained from experimental measurements. Boundary conditions for the interface are also specified from experimental measurements. In contrast, coupled system of equations involving interface capturing methods, pressure-velocity and concentration equations need to be solved. 11. All investigations in heat transfer studies have shown higher heat transfer rates from 180% to 700% as compared to single phase flow. It is important to note that lower sizes (~0.1 mm) give lower percentage increase in heat transfer rate (180%) while higher sizes (~2 mm) give higher percentage increase in heat transfer rate. 12. Most CFD studies have been carried out in 2D or 2D axisymmetric mode while CFD simulations need to carried out in 3 dimensions to get better insight of qualitative and quantitative results.

1.
Several works have been focused on understanding the flow patterns. Few researchers have focussed on flow patterns at the junctions while many of them have focused on the main channels only. In main channels, few researchers observed new flow patterns which are different than the obvious flow patterns of parallel, slug, annular and dispersed flow or a combination of these. These patterns need to be studied more extensively by research fraternity similar to the work of researchers Yagodnitsyna et al. [7], needs to be taken up by more researchers. Numerical simulations for these cases needs to be carried out which for two phase flows for various systems with different fluid properties and operating parameters and such models need to be validated with the rich experimental data available in the literature both in junctions and main channels 2.
Commercial applications need to be found where LLE can be carried out in microchannels. The regime in which such operations would be carried out will pave the way for efficient equipment design. Presently, from the reviewed literature, extraction of rare earth metals seems to be the most promising application that can be commercialized in micro-channels. 3.
There have been very extensive and comprehensive studies for predictions of regime maps by some researchers. A few researchers, claim universal regime maps. However, the individual investigations in literature (for reactive and non-reactive systems) for micro-channels does not quite comply with the flow regime maps except the regime map by Zhang et al. [12]. This might be due to the use of different material of constructions for experimental investigations. Further, several authors claim combination of different dimensional numbers. One of the recommendations would be to form a consortium of experts for development of comprehensive regime maps which can be used for any system which currently seems to be a challenge. Numerical simulations with the power of volume capturing methods would play a vital role in preparing such universal maps 4.
Surprisingly, most of the experimental works compare the VMTC's in micro-channels with conventional equipment's and the VMTCs are two orders of magnitude higher than ones in conventional extractors. But inspite of such various advantages of microchannels, commercial application of micro-channels are not reported in literature. Gaps in commercialization of micro-channel extraction units need to be found and research work needs to be taken up in that direction 5.
The advantages of CFD simulations in understanding flow patterns and mass transfer in micro-channels have not be harnessed to its full potential. More number of works of CFD with green solvents like ionic-liquids and Deep Eutetic Solvents (DES) need to be taken up and validated against the experimental results both in terms of hydrodynamics and mass transfer 6.
Mass transfer models need to be used for predictions for VMTC's with micro-channels as well for the various systems that have been discussed in the present article (both experimental and two phase CFD simulations). Surprisingly, mass transfer models have not been used extensively for predictions of interfacial area and VMTC's. Softwares like Ansys Fluent/Comsol/in-house-code have been used by authors who have developed CFD models. The surface is however considered constant in most of the cases. Movement of the surface is important to capture the Marangoni instabilities which enhances mass transfer. 7.
Very few works in terms of coupled equations of front tracking methods (Level-Set/VOF) with mass transfer models with Navier-Stokes equations exist in the literature for drops or bubbles moving in continuous liquids or liquid films on solid surfaces. However, in these models the distribution coefficient is used as a transformation parameter and thus assumed to be constant. However, in reality, the distribution coefficient depends on state variables and suffers significant changes. Hence, robust mass transfer models also needs to be developed before they can be applied to microchannels for prediction of VMTC's and interfacial areas. Ideally, a comprehensive mathematical model taking account of the interface deformation, interfacial concentration jump and variable distribution coefficient needs to be developed. 8.
More test systems need to be studied for understanding heat tranfer rates in liquidliquid systems with a wide range of sizes 9.
Higher range of Re numbers need to be considered (100 < Re < 1000) for understanding the effect of various parameters on heat transfer rates. Informed Consent Statement: Not applicable.

Data Availability Statement:
Since this is review article the data is taken from the literature.

Conflicts of Interest:
The authors declare no conflict of interest.  Figure 15, m H 2 Channel depth for microchannel as depicted in Figure 15 Table 7,m Heat capacity ratio in  Table 7,τ Contact time of aqueous and organic phases in micro-channel