Advances in High Energy Physics

Volume 2013 (2013), Article ID 732682, 6 pages

http://dx.doi.org/10.1155/2013/732682

Research Article

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

Xian Wang,1,2 Hiroyuki Hirano,² Gongnan Xie,³ and Ding Xu¹

¹State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace, Xi'an Jiaotong University, Xi'an 710049, China

²Department of Applied Chemistry, Okayama University of Science, Okayama 700-0005, Japan

³Engineering Simulation and Aerospace Computing (ESAC), School of Mechanical Engineering, Northwestern Polytechnical University, Xi'an 710072, China

Received 28 May 2013; Accepted 7 July 2013

Academic Editor: Carlo Cattani

Copyright © 2013 Xian Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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 function is set as 1 for the th fluid and 0 for the void, so that the interface locates in the computational cells in which . 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

Figure 1 shows the present model system. Two immiscible fluids were simultaneously injected into a Y-shaped microchannel from two lateral channels with the inlet velocities and , 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 are also shown in Figure 1. The widths of the lateral and main channels are the same and equal to 200 μm. The lengths of the lateral channels and the main channel are  h and  h, respectively. The inclined angles of two lateral channels (Y-shaped junction) are indicated as and , 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

Continuity equation

Momentum equations

where and . Subscripts 1 and 2 denote fluid 1 and 2, respectively. 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 is taken at (0, 1). In this study, the physical properties of the two fluids are assumed to be  kg/m³, Pa·s and  kg/m³, Pa·s. The interfacial tension is set to be  N/m. The magnitudes of inlet velocities are  m/s, and their directions are parallel to the lateral channels (Figure 1). Accordingly, the corresponding typical dimensionless groups for microchannel are and .

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 (or pressure correction ) 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 until the continuity equation is satisfied.

Since 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

, , and can be computed according to . Based on the values of , , 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 in its neighboring grids. Our previous work [23, 24] shows the detailed implementation and performance of our PLIC-reconstruction method.

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]:

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:

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

where is the unit normal to the interface and . The normal vector . 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:

where is the contact angle between fluid and solid wall. It is set to be 90° in the present computation. is the unit wall normal directed into the solid wall and 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. Figure 2 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  m/s, and . 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 at  μm,  m/s and . Red: cyclohexane; green: water.

Figure 3: Numerical simulation result for the segment flow in a Y-shaped microchannel at  μm,  m/s, and . Blue: cyclohexane; green: water.

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

Figure 4 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, is the channel width (Figure 4) and . and 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 as the result in Table 1 for the sake of a slight difference between and caused by some unavoidable computational errors.

Table 1: Relative segment length at various inclined angles of Y-shaped junction at  kg/m³,  Pa·s,  kg/m³,  Pa·s,  N/m, and  m/s.

w/h

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.62

(a)

(b)

(c)

(d)

Figure 4: Simulated segment flow mode in a Y-shaped microchannel at various inclined angles and .

It can be seen from Table 1, under current conditions, when is equal to 15°, the segment length is largest for both cases (and ). When is equal to 90°, although the segment length is not the smallest, it keeps to a small value for both cases ( and ). For the cases of (Table 1 column ), the value of affects the segment length greatly. When , the segment length decreases with the increase in . When , the segment length reaches its smallest value. When , the segment length increases with . However, when (Table 1 column ), the value of the inclined angle does not affect the slug length so much. Except the cases of and , 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 () or get the small slugs and droplets (), 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 ( and ) on segment lengths ( and ) are studied with a typical Y-shaped case of and (hereafter called the T-shaped microchannel). Figure 5 shows the simulated segment flow modes at and 0.02 m/s. is not changed and kept to 0.005 m/s. It can be seen that the segment length of phase 1, , increases with its inlet velocity , while decreases. The segment length can be controlled easily by the inlet flux of the liquids. Figure 6 shows the plot of against for the T-shaped microchannel. Here, is the ratio of the segment lengths of two phases, , and is the ratio of their inlet velocities , which is equal to their inlet volume flow rates. Here, different values of are obtained by different . From Figure 6, it can be seen that is almost linearly proportional to . 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  m/s, the segment flow mode cannot be obtained. Instead, a parallel flow mode appears in the T-shaped microchannel ( and ).

(a)  m/s

(b)  m/s

Figure 5: Simulated segment flow mode in a T-shaped microchannel at  m/s, and (a)  m/s; (b)  m/s.

Figure 6: The plot of against at  m/s for the T-shaped microchannel of and .

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 , and is kept to a small value at for both cases ( and ). 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 ( and ).

Nomenclature

: | Color function
: | Total external force per volume on the fluid (N/m³)
: | The width of microchannel (μm)
: | The length of the lateral channel (μm)
: | The length of the slug (μm)
: | The ratio between organic slug length and water slug length
: | Normal vector to the interface
: | Unit normal vector to the interface
: | Time (s)
: | Time interval (s)
: | Velocity vector (m/s)
: | Velocity component in direction (m/s)
: |
: | Velocity component in direction (m/s)
: | The length of solution domain (m)
: | The width of solution domain (m).

Greeks

: | A constant term in the line segment equation
| Contact angle (°)
: | The curvature of the interface (m⁻¹)
: | Viscosity of fluid (Pa·s)
: | Density of fluid (kg/m³)
: | Interfacial tension coefficient between water and organic fluid (N/m)
: | The angle between lateral channel and the extension central line of main channel in its opposite direction (°).

Subscripts

: | Indicators of cell-centered value
: | Orangic phase
: | Wall tangent component
: | Water phase, or wall nomal component
: | Components in and directions
: | Indicators for fluids 1 (water) and 2 (organic fluid).

Acknowledgments

The would like to acknowledge financial support for this work provided by the Advanced Research in Evolution Area by Okayama Prefecture Industrial Promotion Foundation and QOL Innovative Research (2012–2016) from the Ministry of Education, Culture, Sports, Science and Technology of Japan and the National Natural Science Foundation of China (no. 11242010).

References

W. Ehrfeld, V. Hessel, and H. Lowe, Microreactors, New Technology for Modern Chemistry, vol. D-69469, Wiley-VCH GmbH, Weinheim, Germany, 2000.

P. Watts and C. Wiles, “Recent advances in synthetic micro reaction technology,” Chemical Communications, no. 5, pp. 443–467, 2007.

C. Wiles and P. Watts, “Recent advances in micro reaction technology,” Chemical Communications, vol. 47, no. 23, pp. 6512–6535, 2011.

P. Watts and C. Wiles, “Micro reactors, flow reactors and continuous flow synthesis,” Journal of Chemical Research, vol. 36, no. 4, pp. 181–193, 2012.

J. M. Köhler, T. Henkel, A. Grodrian et al., “Digital reaction technology by micro segmented flow—components, concepts and applications,” Chemical Engineering Journal, vol. 101, no. 1–3, pp. 201–216, 2004.

T. Nisisako, T. Torii, and T. Higuchi, “Novel microreactors for functional polymer beads,” Chemical Engineering Journal, vol. 101, no. 1–3, pp. 23–29, 2004.

L. Wu, M. Tsutahara, L. S. Kim, and M. Ha, “Three-dimensional lattice Boltzmann simulations of droplet formation in a cross-junction microchannel,” International Journal of Multiphase Flow, vol. 34, no. 9, pp. 852–864, 2008.

A.-L. Dessimoz, L. Cavin, A. Renken, and L. Kiwi-Minsker, “Liquid-liquid two-phase flow patterns and mass transfer characteristics in rectangular glass microreactors,” Chemical Engineering Science, vol. 63, no. 16, pp. 4035–4044, 2008.

Y. Okubo, T. Maki, N. Aoki, T. Hong Khoo, Y. Ohmukai, and K. Mae, “Liquid-liquid extraction for efficient synthesis and separation by utilizing micro spaces,” Chemical Engineering Science, vol. 63, no. 16, pp. 4070–4077, 2008.

H. Tang, L. C. Wrobel, and Z. Fan, “Tracking of immiscible interfaces in multiple-material mixing processes,” Computational Materials Science, vol. 29, no. 1, pp. 103–118, 2004.

C. W. Hirt and B. D. Nichols, “Volume of fluid (VOF) method for the dynamics of free boundaries,” Journal of Computational Physics, vol. 39, no. 1, pp. 201–225, 1981.

B. D. Nichols, C. W. Hirt, and R. S. Hotchkiss, “SOLA-VOF: a solution algorithm for transient fluid flow with multiple free boundaries,” Los Alamos Scientific Laboratory Report LA-8355, 1980.

M. Sussman, P. Smereka, and S. Osher, “A level set approach for computing solutions to incompressible two-phase flow,” Journal of Computational Physics, vol. 114, no. 1, pp. 146–159, 1994.

G. Tryggvason, B. Bunner, A. Esmaeeli et al., “A Front-Tracking Method for the Computations of Multiphase Flow,” Journal of Computational Physics, vol. 169, no. 2, pp. 708–759, 2001.

M. R. Swift, W. R. Osborn, and J. M. Yeomans, “Lattice Boltzmann simulation of nonideal fluids,” Physical Review Letters, vol. 75, no. 5, pp. 830–833, 1995.

J. Liu, S. Koshizuka, and Y. Oka, “A hybrid particle-mesh method for viscous, incompressible, multiphase flows,” Journal of Computational Physics, vol. 202, no. 1, pp. 65–93, 2005.

T. Yabe, F. Xiao, and T. Utsumi, “The Constrained Interpolation Profile Method for Multiphase Analysis,” Journal of Computational Physics, vol. 169, no. 2, pp. 556–593, 2001.

W. J. Rider and D. B. Kothe, “Reconstructing Volume Tracking,” Journal of Computational Physics, vol. 141, no. 2, pp. 112–152, 1998.

J. U. Brackbill, D. B. Kothe, and C. Zemach, “A continuum method for modeling surface tension,” Journal of Computational Physics, vol. 100, no. 2, pp. 335–354, 1992.

V. Hessel, S. Hardt, and H. Lowe, Chemical Micro Process Engineering, Fundamentals, Modeling and Reactions, WILEY-VCH Verlag GmbH & Co.KGaA, Weinheim, Germany, 2004.

C. W. Hirt, B. D. Nichols, and N. C. Romero, “SOLA: a numerical solution algorithm for transient fluid flow -Addendum-,” Los Alamos Scientific Laboratory Report LA-5852, 1975.

D. L. Youngs, “Time-dependent multi-material flow with large fluid distortion,” in Numerical Methods For Fluid Dynamics, Edited By, K. W. Morton and M. J. Baines, Eds., pp. 273–285, Academic Press, New York, NY, USA, 1982.

X. Wang, H. Hirano, and G. N. Xie, “A PLIC-VOF based simulation of water-organic slug flow characteristics in a T-shaped microchannel,” Advances in Mechanical Engineering, vol. 2013, Article ID 987428, 8 pages, 2013.

X. Wang, H. Hirano, and N. Okamoto, “Numerical investigation on the liquid-liquid two phase flow in a Y-shaped microchannel,” The Australian and New Zealand Industrial and Applied Mathematics Journal, vol. 48, no. C, pp. C963–C976, 2008.