# VOF Modeling and Analysis of the Segmented Flow in Y-Shaped Microchannels for Microreactor Systems

Authors: Xian Wang; Hiroyuki Hirano; Gongnan Xie; Ding Xu
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

Microscaled devices receive great attention in microreactor systems for producing high renewable energy due to higher surface-to-volume, higher transport rates (heat or/and mass transfer rates), and other advantages over conventional-size reactors. In this paper, the two-phase liquid-liquid flow in a microchannel with various Y-shaped junctions has been studied numerically. Two kinds of immiscible liquids were injected into a microchannel from the Y-shaped junctions to generate the segment flow mode. The segment length was studied. The volume of fluid (VOF) method was used to track the liquid-liquid interface and the piecewise-liner interface construction (PLIC) technique was adopted to get a sharp interface. The interfacial tension was simulated with continuum surface force (CSF) model and the wall adhesion boundary condition was taken into consideration. The simulated flow pattern presents consistence with our experimental one. The numerical results show that a segmented flow mode appears in the main channel. Under the same inlet velocities of two liquids, the segment lengths of the two liquids are the same and depend on the inclined angles of two lateral channels. The effect of inlet velocity is studied in a typical T-shaped microchannel. It is found that the ratio between the lengths of two liquids is almost equal to the ratio between their inlet velocities.


---

## Body

## 1. Introduction

To save energy or decrease energy consumption, the development of new advanced technologies is desired. Among them, the microreaction technology is considered an attractive solution way to replace the traditional energy resources like fossil energy resources. Nowadays, micro-reaction technology is widely used due to the excellent mass and heat transfer properties as well as uniform flow patterns and residence time distributions of the small microfabricated reactors [1]. The microreaction technology has other merits of perfect mixing, effective heat exchange, and small reaction volume. The representations in better performance of heat exchange are higher surface-to-volume ratio and better heat transfer magnitudes compared to normal equipment. The better mixing is attributed to the laminar flow where fast diffusion and thereby perfect mixing within seconds are achieved. Extensive reviews of advances in microreaction technology can be found in [2–6].The spectrum of applications includes gas and liquid flow as well as gas-liquid and liquid-liquid multiphase flows. Recently, immiscible liquid-liquid slug flows in Y-shaped microchannels received extensive attention due to the segmentation of the two fluids, in which the chemical reaction, mixing, and diffusion of the solutes can be enhanced greatly [7–9]. On the other hand, microdroplets can be produced in such a liquid-liquid system in Y-shaped microchannels under some certain conditions. Such can be widely used in the areas of pharmacy, cosmetics, food, and various powder technologies. The efficiency of Y-shaped microchannels for various purposes depends on the segment lengths of both phases. For example, when the segment length is small, the enhanced reactions and small-sized particles can be obtained. Besides, there are various parameters affecting the segment length, such as the channel geometry, viscosity/density ratio between two fluids, interfacial tension, and the flux of two liquids. The various affecting parameters mean that optimizing the segment flow system in microchannels needs extensive experimental works, which would cause a large amount of human and financial resources. In this sense, numerical investigation (including simulation and prediction) on such flow system is indispensable to provide a reasonable and economical designing process for microreactor systems.The variety and complexity of multiphase flow phenomena clearly pose major challenges to the modeling computational approaches. Tang et al. [10] presented the detailed review of the computational methods for multiphase flow simulation. Among the various approaches of multiphase flow simulation, volume of fluid (VOF) [11, 12], level set method [13], front-tracking method [14], and lattice boltzmann method [15] are known as useful tools for modeling multiphase flow phenomena. Further new developments on capturing interface include particle-mesh method [16] and constrained interpolation profile (CIP) method [17].The VOF method is widely adopted by in-house codes and built-in in commercial codes. It is a popular interface tracking algorithm and has proved to be a useful and robust tool. The idea in VOF method is that a color functionCk is set as 1 for the kth fluid and 0 for the void, so that the interface locates in the computational cells in which 0<Ck<1. However, a sharp interface is hard to obtain when solving the color function equation due to the numerical diffusion. This is the key point for this method. Therefore, the reconstruction of the interface is required. Fortunately, several interface reconstruction techniques were developed successfully to fix the requirement, for example, Donor Acceptor, FLAIR (flux line-segment model for advection and interface reconstruction), SLIC (simple line interface calculation) and PLIC (piecewise linear interface calculation). The historical perspective and comparison of these interface reconstruction techniques were introduced in [18].In the multiphase flow, interfacial tension force plays a crucial role. A popular model for the interfacial tension force is the continuum surface force (CSF) model, which was first developed by Brackbill et al. in 1992 [19]. In the CSF model, the interfacial tension was interpreted as a continuous, three-dimensional effect across an interface, rather than as a pressure boundary condition on the interface. It is widely used in the multiphase simulation due to its easy implementation and high accuracy. On the other hand, using the CSF model, the wall adhesion/wetting effect can be conveniently simulated by a simple boundary condition on the interfacial tension force, which is an important phenomenon in a microsystem [20].In this study, the influence of two factors on the segment length of a liquid-liquid flow system in Y-shaped microchannels is numerically studied. These are the inclined angles of Y-shaped junctions and the inlet flow rate of fluids. In the present simulation, the volume of fluid (VOF) method was used to track the liquid-liquid interface and the piecewise-liner interface construction (PLIC) technique was adopted to get a sharp interface. The interfacial tension was simulated with continuum surface force (CSF) model.


## 2. Physical Model

Figure1 shows the present model system. Two immiscible fluids were simultaneously injected into a Y-shaped microchannel from two lateral channels with the inlet velocities u→01 and u→02, respectively. Both phases were assumed to be incompressible and Newtonian at 1 atm and 25°C. The flow was treated as laminar one, which proves correct for the microflow under current situations [1]. The gravitational force was not taken into consideration due to its unimportance in the microflow.Figure 1
Computational model system.The present computations are two-dimensional and steady, and the geometrical parameters of the computational domain(xl×yl) are also shown in Figure 1. The widths of the lateral and main channels h are the same and equal to 200 μm. The lengths of the lateral channels and the main channel are l=5 h and 8l=40 h, respectively. The inclined angles of two lateral channels (Y-shaped junction) are indicated as φ1 and φ2, respectively. An independent grid system of 800 × 220 cells is adopted for all simulations. The motivation of this study is to study the segmented flow at various inclined angles or inlet velocities of a typical Y-shaped microchannel with fixed inclined angles.


## 3. Governing Equations

The governing equations for the fluids are as follows.Color function
(1)
∂
C
∂
t
+
∂
(
u
C
)
∂
x
+
∂
(
v
C
)
∂
y
=
0
.Continuity equation
(2)
∂
u
∂
x
+
∂
v
∂
y
=
0
.Momentum equations
(3)
∂
u
∂
t
+
u
∂
u
∂
x
+
v
∂
u
∂
y
=
-
1
ρ
∂
p
∂
x
+
μ
ρ
(
∂
2
u
∂
2
x
+
∂
2
u
∂
2
y
)
+
1
ρ
f
sv
x
,
∂
v
∂
t
+
u
∂
v
∂
x
+
v
∂
v
∂
y
=
-
1
ρ
∂
p
∂
y
+
μ
ρ
(
∂
2
v
∂
2
x
+
∂
2
v
∂
2
y
)
+
1
ρ
f
sv
y
,
where ρ=ρ2+C(ρ1-ρ2) and μ=μ2+C(μ1-μ2). Subscripts 1 and 2 denote fluid 1 and 2, respectively. C is the color function, which is set to 1 for fluid 1 and 0 for fluid 2, respectively. On the interface of two fluids, the value of C is taken at (0, 1). In this study, the physical properties of the two fluids are assumed to be ρ1=1000 kg/m3, μ1=0.001Pa·s and ρ2=800 kg/m3, μ2=0.02Pa·s. The interfacial tension is set to be σ=0.05 N/m. The magnitudes of inlet velocities are |u→01|=|u→02|=0.005 m/s, and their directions are parallel to the lateral channels (Figure 1). Accordingly, the corresponding typical dimensionless groups for microchannel are Re=ρ1u01h/μ1=1.0 and Ca=μ12/(σρ1h)=10-4.The flow field was solved with finite difference method on a staggered mesh system, and the HSMAC (highly simplified marker and cell) algorithm [21] is employed to solve the pressure field and improve the velocity field. For HSMAC, only the local pressure value pi,jn+1 (or pressure correction pi,j′) is regarded as unknown in the Poisson equation of pressure and is solved by Newton-Raphson method. The velocities in the neighbor grids are corrected by pi,j′ until the continuity equation is satisfied.SinceC is not a continuous variable, (1) cannot be solved as usual. To get a sharp interface, the PLIC technique [22] is adopted, in which a cell-centered line segment is used to simulate the interface and the flux of fluid is computed geometrically according to the interface location and slope. The line segment can be described as
(4)nxx+nyy+η=0.nx, ny, and η can be computed according to C. Based on the values of nx, ny, and η, the flux of fluid 1, whose color function is equal to 1, can be calculated analytically. In the present computation, the outflow flux in a grid is computed, and the in-flow flux is corrected by the out-flow fluxes of the four neighboring grids. This is because the in-flow flux of a grid should be computed according to the values of C in its neighboring grids. Our previous work [23, 24] shows the detailed implementation and performance of our PLIC-reconstruction method.f
→
sv in (3) is the surface tension force per volume. Its direction is pointing into the concave side of the interface. Assuming the center of the interface curve is in fluid 1, for the CSF method, the surface tension is modeled according to [15]:
(5)f→sv=σκ∇ρρ1-ρ2=σκ∇C,
where σ is the interfacial tension coefficient and κ is the local curvature of the interface curve. To avoid the acceleration due to surface tension depending on the density, (5) is improved as follow:
(6)f→sv=ρ(ρ1+ρ2)/2σκ∇C
Thus, when (6) is applied on (3), the last term (i.e., acceleration due to surface tension) is independent of the density. The curvature κ is calculated from
(7)κ=-∇·n^=1|n|[(n|n|·∇)|n|-(∇·n)],
where n^ is the unit normal to the interface and n^=n/|n|. The normal vector n=∇C. The detailed computational equations of every term in (7) are shown in [20].The effect of wall adhesion at fluid interfaces in contact with rigid boundaries, which is very important in a microflow, can be estimated with CSF model by applying a boundary condition in (7) as follows:
(8)n^=(cosθca)n^w+(sinθca)n^t,
where θca is the contact angle between fluid and solid wall. It is set to be 90° in the present computation. n^w is the unit wall normal directed into the solid wall and n^t is the unit tangent normal pointing into fluid 1.


## 4. Results and Discussion

### 4.1. Comparison with Experiment

The validity of our computational code was proved by comparing with our experimental result. Figure2 is the experimental photograph. The red liquid is cyclohexane (with red color dye) and the green one is water (with blue color dye). The channel width is 200 μm, the inlet velocities are u→01=u→02=0.025 m/s, and φ1=φ2=30∘. The interfacial tension is 0.025 N/m. The segment lengths of cyclohexane and water are about 0.67 and 0.7 mm, respectively. It should be noted that the exact segment length is very hard to be measured because the real interfaces shown in Figure 2 are not as clear as the computational one, which is due to the irregular contact angle caused by the complicated solid wall condition. Figure 3 shows the corresponding numerical result of our code under the same conditions with experiment. The blue and green liquids in Figure 3 are organic and aquatic phases, respectively. The computational segment lengths of cyclohexane and water are 0.68 and 0.71 mm. The images of slug generation near the Y-shaped junction, which is indicated in the circle in Figure 2, are quite similar for experimental and numerical results. Therefore, from both qualitative and quantitative points of view, the numerical and experimental results agree well.Figure 2
Experimental photograph for the segment flow in a Y-shaped microchannel ath=200 μm, |u→01|=|u→02|=0.025 m/s and φ1=φ2=30∘. Red: cyclohexane; green: water.Figure 3
Numerical simulation result for the segment flow in a Y-shaped microchannel ath=200 μm, |u→01|=|u→02|=0.025 m/s, and φ1=φ2=30∘. Blue: cyclohexane; green: water.


### 4.2. The Effect of Inclined Angles of Y-Shaped Junction

Figure4 shows three typical images of simulated segment flow modes in the Y-shaped microchannels. Table 1 shows the relative segment length at various inclined angles of Y-shaped junction. Here, h is the channel width (Figure 4) and w=(w1+w2)/2. w1 and w2 are the segment lengths of liquid 1 and 2, respectively (Figure 4). The two values should be same theoretically since the inlet fluxes are the same. We used average value w as the result in Table 1 for the sake of a slight difference between w1 and w2 caused by some unavoidable computational errors.Table 1
Relative segment lengthw/h at various inclined angles of Y-shaped junction at ρ1=1000 kg/m3, μ1=0.001 Pa·s, ρ2=800 kg/m3, μ2=0.02 Pa·s, σ=0.05 N/m, and |u→01|=|u→02|=0.005 m/s.
φ
2
w/h
φ
1
=
0
∘
φ
1
=
φ
2
15°
6.70
4.32
30°
4.13
3.46
45°
3.64
3.44
60°
3.30
3.90
75°
3.50
3.45
90°
3.55
3.45
150°
3.60
3.62Simulated segment flow mode in a Y-shaped microchannel at various inclined anglesφ1 and φ2.
(a)
φ
1
=
φ
2
=
45
°
(b)
φ
1
=
φ
2
=
150
°
(c)
φ
1
=
0
°
,
φ
2
=
45
°
(d)
φ
1
=
0
°
,
φ
2
=
90
°It can be seen from Table1, under current conditions, when φ2 is equal to 15°, the segment length is largest for both cases (φ1=0∘and φ1=φ2). When φ2 is equal to 90°, although the segment length is not the smallest, it keeps to a small value for both cases (φ1=0∘ and φ1=φ2). For the cases of φ1=0∘ (Table 1 column φ1=0∘), the value of φ2 affects the segment length greatly. When φ2<60∘, the segment length decreases with the increase in φ2. When φ2=60∘, the segment length reaches its smallest value. When φ2>60∘, the segment length increases with φ2. However, when φ1=φ2 (Table 1 column φ1=φ2), the value of the inclined angle does not affect the slug length so much. Except the cases of φ1=φ2=15∘ and 60∘, the segment lengths for other cases are almost the same.Thus, the segment length can be controlled by adjusting the inclined angles of Y-shaped junction. T-shaped junction tends to generate more interfaces of two liquids (θca=90∘) or get the small slugs and droplets (θca≠90∘), which is very important for the mixing, chemical reaction, and making micro/nanosized particles using microreactor.


### 4.3. The Effect of Inlet Flow Rate

In this section, the effects of the inlet velocities of two phases (u→01 and u→02) on segment lengths (w1 and w2) are studied with a typical Y-shaped case of φ1=0∘ and φ2=90∘ (hereafter called the T-shaped microchannel). Figure 5 shows the simulated segment flow modes at |u→01|=0.001 and 0.02 m/s. |u→02| is not changed and kept to 0.005 m/s. It can be seen that the segment length of phase 1, w1, increases with its inlet velocity |u→01|, while w2 decreases. The segment length can be controlled easily by the inlet flux of the liquids. Figure 6 shows the plot of w- against u- for the T-shaped microchannel. Here, w- is the ratio of the segment lengths of two phases, w1/w2, and u- is the ratio of their inlet velocities |u→01|/|u→02|, which is equal to their inlet volume flow rates. Here, different values of u- are obtained by different |u→01|. From Figure 6, it can be seen that w- is almost linearly proportional to u-. This means the ratio between the segment volumes of two phases is almost equal to the ratio of their inlet volume flow rates. It should be noted that under current conditions, when |u→01|≥0.04 m/s, the segment flow mode cannot be obtained. Instead, a parallel flow mode appears in the T-shaped microchannel (φ1=0∘ and φ2=90∘).Simulated segment flow mode in a T-shaped microchannel at|u→02|=0.005 m/s, and (a) |u→01|=0.001 m/s; (b) |u→01|=0.02 m/s.
(a)
|
u
→
01
|
=
0.001 m/s
(b)
|
u
→
01
|
=
0.02 m/sFigure 6
The plot ofw- against u- at |u→02|=0.005 m/s for the T-shaped microchannel of φ1=0∘ and φ2=90∘.


## 4.1. Comparison with Experiment

The validity of our computational code was proved by comparing with our experimental result. Figure2 is the experimental photograph. The red liquid is cyclohexane (with red color dye) and the green one is water (with blue color dye). The channel width is 200 μm, the inlet velocities are u→01=u→02=0.025 m/s, and φ1=φ2=30∘. The interfacial tension is 0.025 N/m. The segment lengths of cyclohexane and water are about 0.67 and 0.7 mm, respectively. It should be noted that the exact segment length is very hard to be measured because the real interfaces shown in Figure 2 are not as clear as the computational one, which is due to the irregular contact angle caused by the complicated solid wall condition. Figure 3 shows the corresponding numerical result of our code under the same conditions with experiment. The blue and green liquids in Figure 3 are organic and aquatic phases, respectively. The computational segment lengths of cyclohexane and water are 0.68 and 0.71 mm. The images of slug generation near the Y-shaped junction, which is indicated in the circle in Figure 2, are quite similar for experimental and numerical results. Therefore, from both qualitative and quantitative points of view, the numerical and experimental results agree well.Figure 2
Experimental photograph for the segment flow in a Y-shaped microchannel ath=200 μm, |u→01|=|u→02|=0.025 m/s and φ1=φ2=30∘. Red: cyclohexane; green: water.Figure 3
Numerical simulation result for the segment flow in a Y-shaped microchannel ath=200 μm, |u→01|=|u→02|=0.025 m/s, and φ1=φ2=30∘. Blue: cyclohexane; green: water.


## 4.2. The Effect of Inclined Angles of Y-Shaped Junction

Figure4 shows three typical images of simulated segment flow modes in the Y-shaped microchannels. Table 1 shows the relative segment length at various inclined angles of Y-shaped junction. Here, h is the channel width (Figure 4) and w=(w1+w2)/2. w1 and w2 are the segment lengths of liquid 1 and 2, respectively (Figure 4). The two values should be same theoretically since the inlet fluxes are the same. We used average value w as the result in Table 1 for the sake of a slight difference between w1 and w2 caused by some unavoidable computational errors.Table 1
Relative segment lengthw/h at various inclined angles of Y-shaped junction at ρ1=1000 kg/m3, μ1=0.001 Pa·s, ρ2=800 kg/m3, μ2=0.02 Pa·s, σ=0.05 N/m, and |u→01|=|u→02|=0.005 m/s.
φ
2
w/h
φ
1
=
0
∘
φ
1
=
φ
2
15°
6.70
4.32
30°
4.13
3.46
45°
3.64
3.44
60°
3.30
3.90
75°
3.50
3.45
90°
3.55
3.45
150°
3.60
3.62Simulated segment flow mode in a Y-shaped microchannel at various inclined anglesφ1 and φ2.
(a)
φ
1
=
φ
2
=
45
°
(b)
φ
1
=
φ
2
=
150
°
(c)
φ
1
=
0
°
,
φ
2
=
45
°
(d)
φ
1
=
0
°
,
φ
2
=
90
°It can be seen from Table1, under current conditions, when φ2 is equal to 15°, the segment length is largest for both cases (φ1=0∘and φ1=φ2). When φ2 is equal to 90°, although the segment length is not the smallest, it keeps to a small value for both cases (φ1=0∘ and φ1=φ2). For the cases of φ1=0∘ (Table 1 column φ1=0∘), the value of φ2 affects the segment length greatly. When φ2<60∘, the segment length decreases with the increase in φ2. When φ2=60∘, the segment length reaches its smallest value. When φ2>60∘, the segment length increases with φ2. However, when φ1=φ2 (Table 1 column φ1=φ2), the value of the inclined angle does not affect the slug length so much. Except the cases of φ1=φ2=15∘ and 60∘, the segment lengths for other cases are almost the same.Thus, the segment length can be controlled by adjusting the inclined angles of Y-shaped junction. T-shaped junction tends to generate more interfaces of two liquids (θca=90∘) or get the small slugs and droplets (θca≠90∘), which is very important for the mixing, chemical reaction, and making micro/nanosized particles using microreactor.


## 4.3. The Effect of Inlet Flow Rate

In this section, the effects of the inlet velocities of two phases (u→01 and u→02) on segment lengths (w1 and w2) are studied with a typical Y-shaped case of φ1=0∘ and φ2=90∘ (hereafter called the T-shaped microchannel). Figure 5 shows the simulated segment flow modes at |u→01|=0.001 and 0.02 m/s. |u→02| is not changed and kept to 0.005 m/s. It can be seen that the segment length of phase 1, w1, increases with its inlet velocity |u→01|, while w2 decreases. The segment length can be controlled easily by the inlet flux of the liquids. Figure 6 shows the plot of w- against u- for the T-shaped microchannel. Here, w- is the ratio of the segment lengths of two phases, w1/w2, and u- is the ratio of their inlet velocities |u→01|/|u→02|, which is equal to their inlet volume flow rates. Here, different values of u- are obtained by different |u→01|. From Figure 6, it can be seen that w- is almost linearly proportional to u-. This means the ratio between the segment volumes of two phases is almost equal to the ratio of their inlet volume flow rates. It should be noted that under current conditions, when |u→01|≥0.04 m/s, the segment flow mode cannot be obtained. Instead, a parallel flow mode appears in the T-shaped microchannel (φ1=0∘ and φ2=90∘).Simulated segment flow mode in a T-shaped microchannel at|u→02|=0.005 m/s, and (a) |u→01|=0.001 m/s; (b) |u→01|=0.02 m/s.
(a)
|
u
→
01
|
=
0.001 m/s
(b)
|
u
→
01
|
=
0.02 m/sFigure 6
The plot ofw- against u- at |u→02|=0.005 m/s for the T-shaped microchannel of φ1=0∘ and φ2=90∘.


## 5. Conclusion

Microchannels are gradually used in microreactor systems with microreaction technology, where a typical application includes liquid-liquid multiphase flows. In the present work, a two-phase liquid-liquid segmented flow in a microchannel with various Y-shaped junctions has been studied numerically by means of volume of fluid (VOF) approach. Using liner interface construction (PLIC) technique, a sharp interface between two phases was simulated. The influence of the inclined angle of Y-shaped junction and the inlet velocity on the segmented flow length has been observed.A segment flow mode is found in the Y-shaped microchannel when the contact angle is set to 90°. Under current situations, the segment length is largest atφ2=15∘, and is kept to a small value at φ2=90∘ for both cases (φ1=0∘ and φ1=φ2). This means using T-shaped junction, small slug, and more interfaces of two fluids can be obtained. The ratio between the segment lengths of two liquids is almost equal to the ratio of their inlet velocities for the case of Y-shaped microchannel (φ1=0∘ and φ2=90∘).



---

*Source: 732682-2013-08-04.xml*