Comparison of Co-Current and Counter-Current Flow Fields on Extraction Performance in Micro-Channels

Several applications such as liquid-liquid extraction in micro-fluidic devices are concerned with the flow of two immiscible liquid phases. The commonly observed flow regimes in these systems are slug-flow and stratified flow. The latter regime in micro-channels has the inherent advantage that separation of the two liquids at the exit is efficient. Recently extraction in a stratified counter-current flow has been studied experimentally and it has been shown to be more efficient than co-current flow. An analytical as well as a numerical method to determine the steady-state solution of the corresponding convection-diffusion equation for the two flow-fields is presented. It is shown that the counter-current process is superior to the co-current process for the same set of parameters and operating conditions. A simplified model is proposed to analyse the process when diffusion in the transverse direction is not rate limiting. Different approaches to determining mass transfer coefficient are compared. The concept of log mean temperature difference used in design of heat exchangers is extended to describe mass transfer in the system.


Introduction
The physical effects which govern the behavior of fluids flowing at the micro-scale and the macro-scale are different.Surface tension, viscous effects, energy dissipation and capillary action begin to dominate system behavior at the micro-scale.Micro-fluidics studies the behavior of the fluids at the micro-scale induced by these effects.It helps exploit the behavior at these scales for new applications by improving the efficiency of current processes.
Micro-fluidics results in Process Intensification through miniaturization.These systems are characterized by a dominance of viscous forces as compared to inertial forces; hence, typically low Reynolds numbers are encountered.Consequently, the flow regimes observed in these systems is laminar.Mixing in these micro-channels occurs primarily by molecular diffusion.The time needed for mixing by molecular diffusion is proportional to the square of the length of the diff usion path.The marked shortening of the diffusion path in a micro-channel results in relatively good mixing.
Multiphase liquid-liquid flows arise when two or more partially miscible or completely immiscible fluids are brought in contact and subjected to a pressure gradient.The resulting systems display different kinds of flow behavior, e.g.droplet, slug or stratified flow.These regimes depend on the relative flow rates of the fluid phases involved, the resulting interaction between interfacial and viscous forces and the wetting behavior of the channel walls.Different liquid-liquid two-phase flow patterns in micro-channels have been experimentally analysed, see Dessimoz et al. [1].
In the context of mass transfer across membranes Guo and Ho [2] have analysed an analytical solution based on separation of variables for co-current and counter-current flows in an annulus.The eigen-functions were obtained using a power series expansion.Here the velocity fields in the two chambers were decoupled as they were separated by a porous membrane.Electro-osmotic flows in micro channels have been frequently proposed as a way to overcome dispersion effects.This causes the velocity profile in the micro-channel to remain almost uniform except for a small region near the walls.The flow-field under these conditions can be approximated as an ideal plug-flow in the micro-channel.Liu et al. [3] and Gao et al. [4] have studied the electro-osmotic flow in a rectangular channel when one fluid is conducting and the other is non conducting as is usually the case in extraction.The flow-profiles were obtained analytically for the two fluids when the non conducting fluid was dragged by the viscous force of the conducting fluid.Wang et al. [5] have studied experimentally as well as theoretically how the interface level can be controlled in a two liquid system under the influence of a pressure drop and an electric field applied to one of the fluids.Their results indicate that the flow profiles in the two fluids can be approximated as plug-flows with a jump discontinuity at the interface for some operating conditions.
Recent studies have focused on liquid-liquid extraction in the stratified flow regime in the micro-channels.Here the two fluids flow side by side.This flow-pattern can be exploited to facilitate complete separation at the channel exit.The extraction of vanillin dissolved in water using toluene in micro-structured devices made of Poly Di Methyl Siloxane (PDMS) was studied experimentally by Fries et al. [6].Here the performance of segmented and stratified flow regimes were compared.LIF and micro-PIV measurements showed a laminar profile for stratified flow, whereas vortices in the slug were detected for segmented flow.The influence of channel width and therefore, the surface-to-volume ratio was investigated for stratified flow.There was a significant enhancement of mass transfer with decrease in the channel cross-sectional area for the stratified flow regime.
Three different fluid-flow patterns in a Y-shaped micro-channel, contact or stratified flow, segmented flow and emulsification were investigated in Okubo et al. [7].
Here a one-dimensional model for extraction assuming the interface to be at the centre of the channel was used to compare the model predictions with the experimental behavior.A two-dimensional flow-field taking into account the effect of the interface not being at the centre was analysed numerically by Žnidaršič-Plazl and Igor Plazl [8].They compare their model predictions with experimental results on steroid extraction.Most of the research in stratified flow has been when the fluids flow co-currently i.e. in the same direction.TeGrotenhuis et al. [9] has studied the counter-current mass transfer in a micro-channel when the two fluids are separated by a membrane.Here the diff usional resistance through the membrane was incorporated in the analysis.Recently it has been shown that by suitable modification of the channel surface it is possible to have counter-current flow in the micro-channels over a wide range of operating conditions, Aota et al. [10].They found that a maximum possible theoretical plate number of 4.6 is achievable in counter-current flow as opposed to co-current flow.The pressure drop characteristics in a counter-flow micro-channel have been investigated by Hibara et al. [11].The velocity profiles in a 100 micron channel with butylacetate and an aqueous phase flowing in a countercurrent manner was measured using micro-PIV, Aota et al. [12].
The counter-current flow is known to be more efficient in the context of heat exchanger networks.In this work the performance of co-current and counter-current flows in extraction is being studied with focus on micro-channels.Recently it has been experimentally shown that counter-current operation is possible in micro-channels.The primary objective of this work is to establish conditions under which the counter-current operation is superior to the co-current operation.To the best of our knowledge a theoretical analysis of this system has not been carried out.The main motivation is to show that improvements in the extraction performance are possible when the flow is counter-current as opposed to co-current.The convective diffusion equation is solved analytically for the co-current operation.This is a one dimensional model with diffusion being considered only in the direction transverse to the flow direction.Here axial dispersion effects are neglected.The counter-current system is solved numerically.The algorithm proposed exploits the features of the system.A lumped model is analysed where the concentration dependency on the flow direction alone is considered.It is shown that the counter-current flow performs better than the co-current flow.Different methods to compute the mass transfer coefficient as proposed in the literature are compared.

Co-Current and Counter-Current Flow
We consider three different flow regimes for the analysis in this work: 1) co-current laminar flow, 2) co-current plug flow, and 3) counter-current plug flow.The stratified flow (fluids flow side by side as shown in Figure 1) of two liquid phases is analysed between two infinite horizontal plates extending to infinity in the z-direction.This assumption on the geometry helps us focus on the physics of the problem keeping the mathematics tractable.The distance between the plates (along the x-direction) is taken as H and the liquid-liquid interface is at distance h l (subscript l for laminar) from the lower plate.The flow is assumed to be in the y-direction.Figure 1 shows a schematic of the system being analyzed.

Laminar Flow
In the case of laminar or Poiseuille flow, the velocity profile is obtained assuming the flow to be steady, fully developed and the liquids to be incompressible.The velocity profiles of the system are governed by the equations for 0 , for .
These equations are subject to the conditions of no slip at the walls and continuity of elocity and shear stress at v

Co-Current Plug Flow
the interface which is located at h l .So In this case the velocities v 1 and v 2 have the same sign, and are constant within their phase.It is well established that plug flow behavior can be achieved in a micro channel using electro-osmosis [4].This eliminates axial dispersion effects which arise from laminar flow in micro-channels.The electric field can also be used to control the interface height for a fixed combination of flowrates.

H
(2) In the above the subscript 1 and 2 are used to denote the fluid in the first and second region respectively.The solutions to the above equation yield the velocity profiles of the two liquids as Equation ( 3).
In the case of plug flow, with equal velocity in both phases, we denote interface as h = h n .Here Here the imposed pressure gradient is denoted as P  .Both the fluids are subject to the same pressure drop.The flow behavior is hence similar to the Hagen-Poiseuille flow (parabolic in shape).The velocity field is continuous at the interface h l but its derivatives are discontinuous.A schematic of the velocity profile in the channel is shown in Figure 2. , When plug flow is assumed viscosity does not play any role in determining the interface position.If on the other hand we assume the flow to be laminar, then the interface position is determined by the viscosity of the two fluids.The experimentalist operates the system at fixed flowrates The are given by Equation ( 4).y These equations can be used to determine the pressure drop ∇P and the height h of the interface for a given combination of flow-rates and fluids.Alternatively, if the pressure drop and height of interface are specified, the velocity profiles in each liquid layer can be found and from this the flow rates can be determined.

Counter-Current Plug Flow
In this case the velocities v 1 and v 2 have an opposite sign.We take 1 and 2 0 v  0 v  .Counter-current flow can be theoretically simulated using a combination of Poiseulle flow and a Couette flow.To generate a clear separation of the two phases the interface must be located at The electric field can be manipulated to increase or decrease the velocity for a fixed flow rate.This can be used to control the interface position "h".Experimentally counter-current flow has been achieved by surface modifications of the micro-channels [10].
For the sake of mathematical simplicity the velocity in each of the phases is assumed to be uniform across the transverse direction in counter-current flow.

Mass Transfer in Extraction
The mass transfer behavior in stratified flow of a liquid-liquid extraction system in a micro-channel is now analysed.Here we consider the flow of a solute in the first fluid which is being extracted by the second fluid.The concentration in fluid 1, respectively 2 is represented by C 1 , respectively C 2 .Considering steady-state operation with convection in the y-direction and diffusion in the x-direction we obtain the equations which govern the behavior of the system as for 0 , for .
Here the expressions of v 1 , v 2 take on distinct values for laminar, co-current and counter-current flows.At the interface we have, for , for , at the walls we have 1 2 0 for 0, 0 0 for , 0, and at the inlet y = 0 we have The film interface conditions result in a discontinuous concentration profile, while keeping the mass flux continuous.If K > 1, C 2 remains below the value of C 1 at the interface.When K < 1 the reverse is true and the second fluid extracts the solute out of the first strongly.Here the concentration C 1 is depleted at the interface and we obtain a larger C 2 concentration.In our computations we use C in = 1 mol/m 3 .
In the co-current flow (superscript co) the concentrations of the outlet streams both tend to an equlibrium and this limits the extraction.In the counter-current flow (superscript cc) this limitation does not exist and hence the performance is much better.
For co-current (laminar or plug flow), the conservation of mass states that under steady-state conditions.This is valid for a long channel when the two exiting streams are in equilibrium.Hence For counter-current flow the overall mass balance gives This is used to check the numerical solution.The mixed cup average concentration at a particular "y" is given by

Analytical Solution for Co-Current Plug Flow
The convection diffusion equation can be solved analytically and elegantly under the assumptions of 1) the co-current Plug Flow Regime (PFR) when the velocity in the two fluids is uniform (v 1 and v 2 are constant), and 2) a constant transverse diffusion coefficient (D 1 and D 2 ).We start with non dimensionalizing the equations with respect to their characteristic lengths and initial concentrations, which gives the dimensionless form as 2 1 where , 1, 2.
For simplicity, we drop the superscript * from now on.We seek the solution C i (x, y) in the form g i (y) f i (x).Substituting this in (8) gives This results in an eigen value problem in the x direction whose solution is The boundary conditions at x = 0, 1 yield a = c = 0.At x = h, the boundary condition We seek b and d to be non-zero.This yields the characteristic equation which determines the eigen values λ as the solution to where . The eigen functions corresponding to the nth eigen value is It has been shown in [13] that this system is self-adjoint in the inner product when the velocities in the two fluids are equal i.e., v The eigen functions are normalized with respect to this inner product and the constants b n , d n are obtained as The solution for the y dependency is For co-current extraction the initial condition is C 1 = 1 for 0 < x < h and C 2 = 0 for h < x < 1.The coefficient k n can be obtained from the initial condition as Since the boundary conditions are homogeneous Neumann in the x-direction, λ = 0 is also an eigen value, which corresponds to n = 1.The eigen-function corresponding to this is the equilibrium solution and is given by The complete solution to the convection diffusion equation is hence It was found that it was sufficient to take the first fifty terms in the summation in the above solution to obtain convergence.This implies that the eigen-value problem (14) is solved for the first 50 roots.Care must be taken to ensure that no roots are missed and no roots are calculated more than once.This analytical solution is used to validate the numerical code based on the method of lines with a second order finite difference scheme in the transverse direction (x).The numerical method was used to determine the concentration profiles in the laminar flow regime.
In Figure 3 we show how the analytical solution for the co-current plug flow based on the separation of variables (solid line) validates the predictions of the numerical code (points) based on finite differences.It is seen that both the average concentration as well as the concentration at a fixed y using the two approaches agree quantitatively.
The numerical code was then used to simulate the behavior for the laminar flow profile in co-current mode.
Here the velocity profile obtained in Equation ( 3) is used to simulate the laminar behavior.The comparison of the cup-averaged concentration profiles obtained using the laminar flow and the plug flow behavior in a microchannel is shown in Figure 4.For the plug flow simulation the average velocity of the laminar flow is used.It is seen that the extraction performance under laminar flow conditions is superior to that of the plug flow conditions.

Numerical Solution Counter-Current Flow
For the counter-current plug flow the convective-diffusion equations are solved numerically.Two challenges arise in this and need to be addressed.These are, 1) the jump discontinuity in concentrations at the interface and 2) the inlet of the two fluid streams being at the two end points.The latter renders the system a boundary value problem.The numerical algorithm we use for solving the steady-state convection diffusion equation for extraction under counter-current flow is now described: 1) The channel length is divided into Ny grids in the "y" direction.The values of the solute concentrations at the interface on the fluid1 side are assumed.
2) The values of the concentrations at the interface on the fluid 2 side are obtained using the equilibrium condition.
3) Now the convection diffusion equation in each fluid is solved using the method of lines.This is possible as we have a Dirichlet boundary condition at one end (the interface) and a Neuman condition at the other end (wall) with known inlet conditions.Here a second order scheme is used to discretise the equations in the transverse direction and the equations are integrated along the axial direction.
4) After the solutions are obtained the fluxes at the Ny grid points are calculated in each fluid.The difference in the fluxes at the interface has to be zero.This condition is used to iterate on the concentrations at the interface on fluid 1 till convergence is achieved using a Newton-Raphson technique.
The above algorithm is implemented in Matlab.The cup-mixed average concentration profile along the axis obtained using the above algorithm is shown in Figure 5.

Simplified Model for Co-Current and Counter-Current Flow
To obtain a quick physical insight into the behavior obtained in the two flow-regimes of co-current and countercurrent flow, a simplified model is proposed in this Section.It is valid under the assumptions of a very small height H of the channel (as prevailing in micro-channels), and large diffusion coefficients D i .Under these conditions the concentration variation in the direction transverse to the flow can be neglected and the evolution of the average concentration along the axial direction is governed by ordinary differential equations.For simplicity we assume the velocity profile to follow plug flow.

Co-Current Flow
The simplified equations of mass balance are now given by Here l represents an overall mass transfer coefficient (between the two phases).We now define The solution to the above two equations is given by

Counter-Current Flow
The simplified equations are now given by Introducing and as before, the solution is given by We now describe how the mass transfer coefficient k l can be estimated for a system experimentally.
For Co-current flow, the simplified Equations (18a) and (18b) can be rearranged to yield The rate at which mass is transferred when the concentration drops to c 1 or c 2 in the system is Using these equations, we obtain Rearranging and eliminating the terms containing the interface position "h" using (22b) we obtain At y = L, the exit The logarithmic mean concentration difference is defined as In a similar manner, k l can be calculated for countercurrent flow, using the simplified Equations (20a) and (20b).Following the procedure for co-current flow it can be seen that The mass transfer coefficient k l can also be defined using the driving force for extraction to be the deviation from equilibrium value, see Dessimoz [1].This gives the mass transfer coefficient as

Results and Discussion
In order to compare the results of our simulations and to be consistent with the literature, and evaluate the performance of a specific micro-channel set-up, we introduce some characteristic quantities.These are now defined.

Characteristic Quantities
The first is the efficiency E, defined in terms of the mixed cup concentrations as where 2 is the concentration of the solute in the second region after equilibrium is attained, and typically eq C 2, 0 = in ."E" is a measure of how close the exiting stream is to equilibrium.The overall residence time t res for co-current flow is defined as For a given length L, a unique residence time and an extraction efficiency E(L) is obtained.E = 1 corresponds to the situation when the exiting streams are in equilibrium and no further separation can take place.
The second characteristic which can describe the system is the extraction ratio E r .It represents the fraction of the amount of solute that has been fed to the system which is removed by the second fluid.This is defined as Note that for co-current flow Figure 6 shows the dependency of concentration along  is estimated by minimizing the least squares error of the exit concentration from the two models.It can be seen that the prediction by the convective diffusion model and the simplified model agree well throughout the length of the curve.We would like to avoid the least squares estimation of the mass transfer coefficient since this is a statistical procedure and involves fitting a parameter.For this the concentration profiles are determined solving the convective diffusion equation.The exit concentration results are used as experimental results and the log mean concentration difference is used to determine the mass transfer coefficient as described earlier.The predictions of the concentration profiles of the convective diffusion equation and the simplified model are shown in Figures 8(a) and 8(b).It is seen that the concentration profiles are in better agreement for the counter-current flow as opposed to the co-current flow.The prediction of the simplified model using the log mean concentration difference is better when the diffusion coefficient is higher and this is seen when we compare Figure 8(c) with Figure 8(a).Figure 8(d) compares the mass transfer coefficient predicted by the log mean concentration difference and the deviation from equilibrium (Equation ( 28)).We see that the mass transfer coefficient at short lengths of the channels are high.This drops down sharply first and then slowly as the length is increased.This is to be expected since for short lengths the driving force is high and this results in a larger mass transfer coefficient.As the length increases the driving force decreases as we approach equilibrium and this lowers the mass transfer coefficient.
Figure 9(a) shows comparison of co-current concentration profiles obtained from the convection diffusion model and the simplified model.Here the mass transfer coefficient is calculated using the driving force as the deviation from the equilibrium concentration Equation (28).Here again when the diffusivity is high we see a much better match between the axial profiles (Figure 9(b)).At the same higher diffusivity, profiles using mass transfer coefficient calculated from the deviation from equilibrium matches better than those obtained using the log mean concentration difference.(compare Figure 9(b) and 8(c)).Figure 9(c) shows a comparison of k l values calculated using the above two approaches.It is observed that mass transfer coefficients from Equation (28) are much higher than that using log mean difference at all the lengths of the channel.This difference is very high at the shorter lengths and it decreases as the length increases as observed in Figure 8(d).

Conclusions
In this work we have compared the performance of the co-current and countercurrent modes of operation in an  extraction in a micro channel.This theoretical study examines the convective diffusion equation as well as a simplified model.An analytical solution is proposed for the co-current system based on the method of separation of variables.The counter-current system is solved numerically.For this a novel method which exploits the feature of the system is proposed.It is based on iterating on the interface values of the concentration.The study is numerical and shows how concepts from heat exchanger design can be extended to mass transfer systems.It is shown that the counter-current extraction is more efficient than the co-current operation.It is found that the different definitions of the mass transfer coefficient predict the behavior of the system accurately.

Figure 1 .
Figure 1.Micro-channel layout.1,2 describe regions occupied by the two fluids.These are separated at the interface h.The direction of velocity determines the flow type.(a) Co-current; (b) Counter-current.

Figure 2 .
Figure 2. Schematic picture of stratified flow system showing all variables.The velocity profile shown is for μ 1 > μ 2 , .0   P the point where the velocity is zero.Alternatively combination of electro osmotic flow with Poiseulle flow can give a counter-current flow when there are two immiscible liquids as the electric field affects the flow of only one of the two fluids.The electric field can be manipulated to increase or decrease the velocity for a fixed flow rate.This can be used to control the interface position "h".Experimentally counter-current flow has been achieved by surface modifications of the micro-channels[10].For the sake of mathematical simplicity the velocity in each of the phases is assumed to be uniform across the transverse direction in counter-current flow.

Figure 3 .
Figure 3.Comparison of the analytical solution and the numerical computation for co-current plug-flow.Here H = 4 × 10 −4 m, h = 2 × 10 −4 m, Q 1 = Q 2 = 14.2857 × 10 −6 m 2 /s, K = 4.14, D 1 = 7.4 × 10 −9 m 2 /s, D 2 = 5.64 × 10 −8 m 2 /s.Full line is the average concentration for x > h and the dashed line the concentration at x = 0.75 H from the analytical computation.The triangles are the corresponding computed values with the numerical computation, where a grid with 100 points in the transverse direction per phase is considered.

Figure 4 .
Figure 4. Comparison of the laminar flow versus PFR with fixed Q 1 and Q 2 .Full line is laminar flow with h = 0.42 H, K = 1/4.14,dashed line PFR with same interface position.