A PLIC-VOF-Based Simulation of Water-Organic Slug Flow Characteristics in a T-Shaped Microchannel

Abstract

A water-organic slug flow in a T-shaped microchannel was numerically studied due to its importance in the microreactor system. Various factors affecting the flow mode were studied, for example, channel width, fluid viscosity, interfacial tension, and inlet velocity. 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 results show that strong vortexes appear in both phases at the meeting sites of main and lateral channels where an organic slug is producing. Inlet velocity influences the slug length and flow mode greatly. The ratio between the slug lengths of two phases in the main channel is almost equal to the ratio between their inlet velocities. If the slug is produced, the interfacial tension and organic viscosity have less effect on the slug length for 200 μm microchannel. The slug producing rate is much higher in a narrow channel than that in a wide channel.

1. Introduction

Nowadays microreaction 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. Some of these devices even are available commercially and being tested by a growing community of researchers [1]. The spectrum of applications includes gas and liquid flow as well as gas-liquid and liquid-liquid multiphase flow. Recently, immiscible liquid-liquid slug flows in microchannels absorb extensive attentions due to the segmentation of the two fluids, in which the chemical reaction, mixing, and diffusion of the solutes can be enhanced greatly. On the other hand, microdroplets can be produced in such a liquid-liquid system in the microchannel under some certain conditions, for example, small slug, large interfacial tension and small contact angle between fluid and solid wall. The efficiency and properties of the slug flow in microreactors are mainly influenced by the channel geometry, viscosity ratio between two fluids, interfacial tension, and the flux of two liquids [2]. The various affecting parameters mean that optimizing the slug-flow system in micro-eactors needs extensive experimental work. Therefore, numerical study on such a system is indispensable to provide a reasonable and economical designing process.

As we know, the variety and complexity of flow phenomena clearly pose major challenges to the modeling approach. Accurate simulation of a multiphase flow problems with a moving interface lies in the obtaining of a sharp interface, which is quite difficult and a great challenge for the CFD researchers. Several techniques were developed in the last 25 years and can be classified into three main categories according to physical and mathematical approaches: capturing (Lagrangian), tracking (Eulerian), and combined methods. Reference [3] presents the detailed review of these methods. Among the various approaches of multiphase flow simulation, volume-of-fluid (VOF) [4, 5], level set method [6], front-tracking method [7], and lattice Boltzmann method [8] are known as useful tools. Further new developments on capturing interface include particle-mesh method [9] and CIP method [10].

The VOF method is widely adopted by in-house codes and built-in commercial codes. It is a popular interface tracking algorithm which has proved to be a useful and robust tool since its development decades ago, in which a color function C_k 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 < C_k < 1. However, a sharp interface is hard to obtain when solving the color function equation due to the numerical diffusion, which is the key point for this method. Therefore, the reconstruction of the interface is required. Several interface reconstruction technique were developed. for example, donor acceptor, flux line-segment model for advection and interface reconstruction, simple line interface calculation, and piecewise linear interface calculation. The historical perspective and comparison of these interface reconstruction techniques were introduced in detail in [11].

In the multiphase flow, interfacial tension force plays a crucial role. A popular model for the interfacial tension force is continuum surface force (CSF) model, which was developed by Brackbill in 1992 [12]. In 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 CSF model, the wall adhesion/wetting effect can be easily simulated by a simple boundary condition on the interfacial tension force, which is an important phenomenon in a microsystem [13].

In the present study, using PLIC-VOF method combined with CSF interface tension model, a water-organic slug flow in a T-shaped microchannel was numerically studied including the effects of channel size, viscosity of organic phase, interfacial tension and inlet velocities of two fluids.

2. Model System

Figure 1 shows the present model system. Water and organic phases were simultaneously injected into the microchannel from the inlets of main channel and lateral channel, respectively. The inlet velocities were indicated as u₁ for water and u₂ for organic phase. Both phases were assumed to be incompressible, Newtonian and immiscible at 1 atm, 25°C. The flow was treated as laminar one, which proved correct for the microflow [13]. The gravitational force was not taken into consideration due to its less importance in the microflow. The width of the channel was indicated as h, the length of the lateral channel was l, and it is equal to 5h. The length of the main channel was 4l.

Figure 1:

Computational model system.

In the present work, two-dimensional computation was considered. The computational domain (xl × yl) is shown in Figure 1. Here, xl = 4l and yl = l + h.300 × 120 rectangular grids were adopted. Since the T-shaped flow region occupied only part of the computational domain (xl × yl), the velocities outside the T-shaped region were set to zero manually in the computation.

3. Computational Methods and Governing Equations

In the present computation, the immiscible interface of the two phases was tracked by volume-of fluid (VOF) method [4, 5] on Eulerian grids. The color function C was set as 1 for water and 0 for organic phase. The interface was reconstructed by piecewise liner (PLIC) technique which was initially proposed by Youngs [14]. The surface tension was computed by continuum surface force (CSF) model [12].

3.1. Dimensional Governing Equations

The dimensional governing equations are as follows.

Color function:

\frac{\partial C}{\partial t} + u \frac{\partial C}{\partial x} + v \frac{\partial C}{\partial y} = 0 .

(1)

Continuity equation:

\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0 .

(2)

Momentum equation:

\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = - \frac{1}{ρ} \frac{\partial p}{\partial x} + \frac{μ}{ρ} (\frac{\partial^{2} u}{\partial x^{2}} + \frac{\partial^{2} u}{\partial y^{2}}) + \frac{1}{ρ} f_{sv x}, \frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} = - \frac{1}{ρ} \frac{\partial p}{\partial y} + \frac{μ}{ρ} (\frac{\partial^{2} v}{\partial x^{2}} + \frac{\partial^{2} v}{\partial y^{2}}) + \frac{1}{ρ} f_{sv y},

(3)

where ρ = ρ₂ + C(ρ₁ – ρ₂) and μ = μ₂ + C(μ₁ – μ₂). Subscripts 1 and 2 stand for water (fluid 1) and organic phase (fluid 2), respectively. In the present computation, the physical properties ρ₁ = 1000 [kg/m³], μ₁ = 1.0 × 10⁻³ [Pa·s] were adopted for the water phase and ρ₂ = 800 [kg/m³], μ₂ = (1~20)μ₁ for the organic phase. f_sv in (3) is the total external forces per volume on the fluids. In this work, only interfacial tension force was considered as an external force.

The flow field was solved with finite difference method on a staggered mesh system. The highly simplified marker and cell (HSMAC) algorithm [15] was adopted to solve the pressure field and to improve the velocity field.

3.2. Interfacial Tension-CSF Model

f_sv in (3) is the surface tension force per volume, and the direction of f_sv is pointing into the concave side of the interface. Assuming the center of the interface curve is in fluid 1 (water), for the CSF method, the surface tension is modeled as follows [12]:

f_{sv} = σ κ \frac{\nabla ρ}{ρ_{1} - ρ_{2}} = σ κ \nabla C,

(4)

where σ is the interfacial tension coefficient [N/m] and κ is the local curvature of the interface curve [1/m]. To avoid the acceleration due to surface tension depending on the density, (4) is improved as

f_{sv} = \frac{ρ}{(ρ_{1} + ρ_{2}) / 2} σ κ \nabla C .

(5)

Thus, when (5) is applied on (3), the last term, that is, acceleration due to surface tension, is independent on the density. The curvature κ is calculated from

κ = - \nabla \cdot \hat{n} = \frac{1}{| n |} [(\frac{n}{| n |} \cdot \nabla) | n | - (\nabla \cdot n)],

(6)

where, $\hat{n}$ is the unit normal to the interface and $\hat{n} = n / | \hat{n} |$ . The normal vector n = ∇ C, and the cell-centered n_{i, j} is computed as follows:

\begin{matrix} n_{x i, j} = \frac{1}{4} (n_{x i, j} + n_{x i + 1 / 2, j + 1 / 2} + n_{x i - 1 / 2, j + 1 / 2} + n_{x i + 1 / 2, j - 1 / 2} + n_{x i - 1 / 2, j - 1 / 2}), \\ n_{y i, j} = \frac{1}{4} (n_{y i, j} + n_{y i + 1 / 2, j + 1 / 2} + n_{y i - 1 / 2, j + 1 / 2} + n_{y i + 1 / 2, j - 1 / 2} + n_{y i - 1 / 2, j - 1 / 2}) . \end{matrix}

(7)

For example;

\begin{matrix} n_{x i + 1 / 2, j + 1 / 2} = \frac{1}{2} (\frac{C_{i + 1, j + 1} - C_{i, j + 1}}{Δ x} + \frac{C_{i + 1, j} - C_{i, j}}{Δ x}), \\ n_{y i + 1 / 2, j + 1 / 2} = \frac{1}{2} (\frac{C_{i + 1, j + 1} - C_{i + 1, j}}{Δ y} + \frac{C_{i, j + 1} - C_{i, j}}{Δ y}) . \end{matrix}

(8)

Then, we get

n_{x} = \frac{1}{8 Δ x} (C_{i + 1, j + 1} + 2 C_{i + 1, j} + C_{i + 1, j - 1} - C_{i - 1, j + 1} - 2 C_{i - 1, j} - C_{i - 1, j - 1}),

(9)

n_{y} = \frac{1}{8 Δ y} (C_{i + 1, j + 1} + 2 C_{i, j + 1} + C_{i - 1, j + 1} - C_{i + 1, j - 1} - 2 C_{i, j - 1} - C_{i - 1, j - 1}),

(10)

where Δx and Δy are the dimensions of a grid in x and y directions.

Cell-centered tension force f_sv can be rewritten as

f_{sv} = \frac{\nabla ρ}{(ρ_{1} + ρ_{2}) / 2} σ n (- \nabla \cdot \hat{n}) .

(11)

The required face-centered values f_svx, f_svy in (3) are obtained by interpolating from the two nearest cell-centered values.

The effect of wall adhesion at fluid interfaces in contact with rigid boundaries, which is very important in a microflow, can be estimated easily in CSF model by applying a boundary condition in (11) as follow:

\hat{n} = (\cos θ_{ca}) {\hat{n}}_{w} + (\sin θ_{ca}) {\hat{n}}_{t},

(12)

where θ_ca is contact angle between fluids and solid wall. Contact angle is influenced by various factors, for example, the species of the fluids, the smoothness, and geometry of the solid wall. The static contact angle can be measured experimentally when the fluids are at rest. ${\hat{n}}_{w}$ is the unit wall normal directly into the solid wall, and ${\hat{n}}_{t}$ is the unit tangent normal pointing into fluid 1. For the points on the solid walls or near the solid walls (the case when the grid boundary does not overlap with the wall), the tangent component n_t can be computed using color function C reflected at the wall, and the wall component n_w can be computed in terms of n_t, that is, n_w = n_tcot⁻¹θ_ca. In the present study, the contact angle θ_ca was assumed to be 90° for the easiness of computation. Our work [16] shows the effect of contact angle on the slug shape.

3.3. PLIC Reconstruction of the Interface

To get a sharp interface, PLIC technique [14] was used to make reconstruction of the interface, 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 is described as

n_{x} x + n_{y} y + η = 0,

(13)

where n_x and n_y were computed according to (9) and (10), respectively.

Figure 2 shows four typical illustrations of the moving interface reconstruction in a computational grid. The dark area is full of fluid 1 and the light one fluid 2. Based on the sign of n_x, n_y, 16 reconstruction images can be obtained. Assuming β is the acute angle between the line segment and x-coordinate, then cot β = |n_x/n_y|. Assuming $cot α = cot β (Δ x / Δ y)$ , the four basic types shown in Figure 2, for which n_x > 0, n_y < 0, can be decided by the color function C and α as follows:

type (a): C ≤ 0.5cot α or C ≤ 0.5cot⁻¹ α,

type (b): C ≤ 1 – 0.5cot α,

type (c): C ≤ 1 – 0.5cot⁻¹ α,

type (d): C ≥ 1 – 0.5cot α or C ≥ 1 – 0.5cot⁻¹ α.

Figure 2:

Sample illustrations of the moving interface reconstruction in a computational grid.

Other 12 types can also be decided by the same way. Besides, there are two special cases: n_x = 0 and n_y = 0, that is, the interface is horizontal and vertical. They should be considered separately from the above 16 cases.

After deciding the type of the line segment, the flux of fluid 1 whose color function is equal to 1 can be calculated analytically. In the present computation, the out-flow flux in a grid was computed, and the in-flow fluxes was corrected by the out-flow fluxes in the four neighbor grids. This is because the in-flow flux should be computed according to the value of C in the neighbor grids. Figure 3 shows the example of flux calculation. Assuming the velocities on the four sides are all positive, that is, u_l, u_b, u_r, u_t > 0, the out-flow fluxes on the left and bottom sides of the grid f_l, f_b are zero. The ones on the right and top sides f_r, f_t are the shaded areas S_r and S_t, respectively, which can be computed geometrically by β and C. The volume of fluid 1 in the grid (dark area) is CΔxΔy for two dimensional computation. Then, the color function C was computed as follows:

C (i, j + 1) = C^{*} (i, j + 1) + \frac{f_{t}}{(Δ x Δ y)}, C (i, j - 1) = C^{*} (i, j - 1) + \frac{f_{b}}{(Δ x Δ y)}, C (i + 1, j) = C^{*} (i + 1, j) + \frac{f_{r}}{(Δ x Δ y)}, C (i - 1, j) = C^{*} (i - 1, j) + \frac{f_{l}}{(Δ x Δ y)}, C (i, j) = C^{*} (i, j) - \frac{(f_{t} + f_{b} + f_{r} + f_{l})}{(Δ x Δ y)},

(14)

where superscript * stands for the value before the current correction. It should be noted that, in the boundary girds, both the in-flow and out-flow fluxes at the boundary side should contribute to the flux change at the same time.

Figure 3:

The sketch of flux calculation in a computational grid.

Figure 4 shows the performances of our PLIC-reconstruction method on Zalesak's slotted cylinder test problem [17]. The solution region is 1 × 1 and grids are 200 × 200. Time step Δτ = 10⁻⁴. It shows clearly that the present PLIC reconstruction tracks the sharp boundary efficiently.

Figure 4:

Performance of VOF-PLIC reconstruction on Zalesak's slotted cylinder test problem using 200 × 200 grids.

4. Result and Discussion

4.1. Transient Flow Mode and Velocity Field

Figure 5 shows the transient flow modes at h = 200 μm, |u₁| = |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m, and μ₂ = 20.0 × 10⁻³ Pa·s. The light color fluid is water, and the dark one is the organic fluid. The initial condition is shown in the picture t = 0.0 s. For the multiphase flow computation, a reasonable initial condition for the color function C is very important. In the present work, the same initial condition were set for all the computational cases. At t = 0.08–0.16 s, the organic fluid flows into the main channel and the first organic slug (dark one) starts to be produced near the crossing of the lateral and main channels. At t = 0.24–0.32 s, the second organic slug is produced and the same situation occurs as the first one at t = 0.08–0.16 s. Such a process repeats with time proceeding. A water-organic alternating slug flow pattern in the main channel is obtained. In this paper, symbols L_w and L_o are used to indicate the lengths of water and organic slugs shown in Figure 5 (t = 0.40 s).

Figure 5:

The transient flow modes at h = 200 μm, |u₁| = |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa·s.

Figure 6 shows the local fluid streak lines and velocity vectors at t = 0.2–0.34 s for h = 200 μm, |u₁| = |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m, and μ₂ = 20.0 × 10⁻³ Pa · s. At t = 0.2 s, an old organic slug was just produced. The fluid streak lines are wavy and no vortex is founded. A complete image is refereed to picture t = 0.162 s in Figure 5. At t = 0.204–0.232 s, a new organic slug is growing. In this period, strong vortexes appear on the interface of two phases and inside the organic phase. As soon as the organic phase touches the bottom of the solid channel wall, the water-organic interface on the right-hand side becomes perpendicular to the wall at once due to θ_ca = 90°, which happens at t = 0.236–0.280 s. In this period, the vortex appears mainly in the water phase. At t = 0.320 s, vortexes disappear and the streak lines become wavy again. At t = 0.340 s, the next organic slug starts to be produced and the old one flows downstream. The process from t = 0.200 s to 0.320 s shown in Figure 6 repeats in the main channel. An organic slug was produced in about 0.12 s under current conditions.

Figure 6:

Local fluid streak lines and velocity vectors at t = 0.2–0.34 s for h = 200 μm, |u₁| = |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa·s.

4.2. The Effects of Channel Width h, Interfacial Tension, σ and Organic Viscosity μ₂

The channel width h affects the length of slug to a large extent. The cases of h = 500, 200, 100 and 50 μm at |u₁| = |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa · s are studied here. The results show that for the microchannels with h = 500, 200, 100, and 50 μm, the corresponding lengths of organic slugs L_o (Figure 5) are 1970, 730, 353, and 177 μm, respectively. The corresponding producing rate of organic slug N [s⁻¹], that is, the number of organic slug produced in a second, is about 3, 7, 21 and 39, respectively. It is concluded that with the decrease in the channel width h, the organic slug length decreases, and the producing rate of organic slug increases greatly.

Table 1 shows the length of organic slug L_o [μm] under various interfacial tensions σ and organic viscosities μ₂ at h = 200 μm and |u₁| = |u₂| = 5 mm/s. It can be seen for two cases σ = 1.0 × 10⁻² N/m, μ₂ = 10.0 × 10⁻³ and 20.0 × 10⁻³ Pa · s that the parallel flow is founded; that is, there is no slug produced in the mail channel and a slightly wavy interface appears at the horizontal centerline of the main channel. Except for the above two cases, the slug can be produced and the slug lengths have no large difference (maximum difference 130 μm). Therefore, as long as slugs can be produced in the water-organic system, the organic slug length is not affected by σ and μ₂ so much for the channel of h = 200 μm. It should be mentioned that the position where the organic slug is produced is different. For large σ and small μ₂, the organic slug is produced just at the crossing of the lateral and main channels. For small σ and large μ₂, this position lies in the middle or downstream of the main channel and away from the crossing; that is, there exists a parallel flow before an organic slug is produced.

Table 1:

The length of organic slug L_o (μm) under various interfacial tensions σ (N/m) and organic viscosities μ₂ (Pa·s). Other conditions are: h = 200 μm and $| u_{1} | = | u_{2} | =$ 5 mm/s, θ_ca = 90°, ρ₁ = 1000 kg/m³, μ₁ = 1.0 × 10⁻³ Pa·s and ρ₂ = 800 kg/m³.

L_o (μm)	μ₂ (Pa·s)
L_o (μm)	1.0 × 10⁻³	2.0 × 10⁻³	5.0 × 10⁻³	10.0 × 10⁻³	20.0 × 10⁻³

σ (N/m)
0.05	624	648	612	608	730
0.03	636	658	612	730	670
0.01	600	730	636	Parallel flow	Parallel flow

4.3. The Effect of Inlet Velocity of Water u₁

In this section, the effects of the inlet velocity of water phase on the organic slug length (L_o) and water slug length (L_w) are studied. Although it is reasonable to study the effect of inlet volume flow rate [m³/s], due to the two-dimensional computation, the cross section area cannot be evaluated correctly. Since the lateral and main channel have the same-sized cross section, the inlet velocity ratio of two phases is equal to the ratio of their inlet volume flow rates, and we use the inlet velocity to study the effect of inlet volume flow rate. To show this effect clearly, the main channel length is increased to 8l in this section.

Figure 7 shows the flow modes under various inlet velocities of water |u₁| at t = 0.8 s, h = 200 μm, |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa · s. It can be seen that the inlet velocity of water affects flow mode and slug lengths greatly. With |u₁| increasing, the organic slug length decreases and the spacing between neighboring organic slugs increases. Particular, at |u₁| = 40 mm/s, there is no slug produced; instead of it, a parallel flow mode appears and the water phase occupies most of the channel. Figure 8 shows the slug lengths L_o and L_w [μm] against |u₁|[mm/s]. Figure 9 shows the ratio of slug lengths of two phases ( $\bar{L} = L_{o} / L_{w}$ ) against the ratio of inlet velocities of two phases ( $\bar{u} = | u_{2} | / | u_{1} |$ ). Other conditions in Figures 8 and 9 are same with these in Figure 7. It can be seen that organic slug length L_o decreases with the increase in |u₁| while water slug length L_w increases (Figure 8). From Figure 9, we can see $\bar{L}$ is almost linearly proportional to $\bar{u}$ . That means thatthe ratio between the slug volumes of two phases in the main channel is almost equal to the ratio of their inlet volume flow rates.

Figure 7:

Flow modes under various inlet velocities of water |u₁| at t = 0.8 s, h = 200 μm, |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa·s.

Figure 8:

Slug lengths L_o and L_w [μm] against |u₁| [mm/s]. Other conditions are: h = 200 μm, |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa·s.

Figure 9:

$\overset{—}{L}$ against $\overset{—}{u}$ . Other conditions are: h = 200 μm, |u₁| = |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa·s.

Figure 10 shows the producing rate of organic slug N [s⁻¹] against |u₁| [mm/s]. It can be seen that N increases rapidly with |u₁| increasing when |u₁| is in some range. For the current condition, |u₁| < 40 mm/s.

Figure 10:

Producing rate of organic slug N [s⁻¹] against |u₁| [mm/s]. Other conditions are: h = 200 μm, |u₂| = 5 mm/s, σ = 5.0 × 10⁻² N/m and μ₂ = 20.0 × 10⁻³ Pa·s.

5. Conclusion

The fluid-fluid flow plays an important role in micro-reactor for the enhancement of heat, mass, and momentum transfer. In the present work, we numerically studied water-organic flow in a T-shaped microchannel which is often met in the microreactor. Several important factors affecting the flow mode were studied, for example, the channel width, inlet velocity, fluid viscosity, and interfacial tension. These factors influence the flow mode, slug length, slug producing rate in the microchannel, especially the inlet velocity. It is also found that when a slug is producing, strong vortexes appear in both fluids. After a slug is produced, it moves along the channel, and there is no vortexes founded in the downstream of the main channel. The numerical results are expected to be a guide for various experimental works.

Footnotes

Nomenclature

Acknowledgments

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

References

Ehrfeld

Hessel

, and Lowe

, Microreactors, New Technology for Modern Chemistry, WILEY-VCH, Weinheim, Germany, 2000.

Harries

Burns

J. R.

Barrow

D. A.

, and Ramshaw

, “A numerical model for segmented flow in a microreactor,” International Journal of Heat and Mass Transfer, vol. 46, no. 17, pp. 3313–3322, 2003.

Tang

Wrobel

L. C.

, and Fan

, “Tracking of immiscible interfaces in multiple-material mixing processes,” Computational Materials Science, vol. 29, no. 1, pp. 103–118, 2004.

Hirt

C. W.

and Nichols

B. D.

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

Nichols

B. D.

Hirt

C. W.

, and Hotchkiss

R. S.

, “SOLA-VOF: a solution algorithm for transient fluid flow with multiple free boundaries,” Los Alamos Scientific Laboratory Report LA-8355, 1980.

Sussman

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

Tryggvason

Bunner

Esmaeeli

, “A front–tracking method for the computations of multiphase flow,” Journal of Computational Physics, vol. 169, no. 2, pp. 708–759, 2001.

Swift

M. R.

Osborn

W. R.

, and Yeomans

J. M.

, “Lattice Boltzmann simulation of nonideal fluids,” Physical Review Letters, vol. 75, no. 5, pp. 830–833, 1995.

Liu

Koshizuka

, and Oka

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

10.

Yabe

Xiao

, and Utsumi

, “The constrained interpolation profile method for multiphase analysis,” Journal of Computational Physics, vol. 169, no. 2, pp. 556–593, 2001.

11.

Rider

W. J.

and Kothe

D. B.

, “Reconstructing volume tracking,” Journal of Computational Physics, vol. 141, no. 2, pp. 112–152, 1998.

12.

Brackbill

J. U.

, “A continuum method for modeling surface tension,” Journal of Computational Physics, vol. 100, no. 2, pp. 335–354, 1992.

13.

Hessel

Hardt

, and Lowe

, Chemical Micro Process Engineering, Fundamentals, Modeling and Reactions, WILEY-VCH, Weinheim, Germany, 2004.

14.

Youngs

D. L.

, “Time-dependent multi-material flow with large fluid distortion,” in Numerical Methods for Fluid Dynamics, Morton

K. W.

and Baines

M. J.

, Eds., pp. 273–285, Academic Press, New York, NY, USA, 1982.

15.

Hirt

C. W.

Nichols

B. D.

, and Romero

N. C.

, “SOLA: a numerical solution algorithm for transient fluid flow -Addendum-,” Los Alamos Scientific Laboratory Report LA-5852, 1975.

16.

Wang

Hirano

, and 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, pp. C963–C976, 2008.

17.

Zalesak

S. T.

, “Fully multidimensional flux-corrected transport algorithms for fluids,” Journal of Computational Physics, vol. 31, no. 3, pp. 335–362, 1979.

A PLIC-VOF-Based Simulation of Water-Organic Slug Flow Characteristics in a T-Shaped Microchannel

Abstract

1. Introduction

2. Model System

3. Computational Methods and Governing Equations

3.1. Dimensional Governing Equations

3.2. Interfacial Tension-CSF Model

3.3. PLIC Reconstruction of the Interface

4. Result and Discussion

4.1. Transient Flow Mode and Velocity Field

4.2. The Effects of Channel Width h, Interfacial Tension, σ and Organic Viscosity μ2

4.3. The Effect of Inlet Velocity of Water u1

5. Conclusion

Footnotes

Nomenclature

Acknowledgments

References

4.2. The Effects of Channel Width h, Interfacial Tension, σ and Organic Viscosity μ₂

4.3. The Effect of Inlet Velocity of Water u₁