Two-Phase Flow Regime Transitions in Microchannels: A Comparative Experimental Study

Even though numerical simulation of two-phase flows in microchannels has been attempted by many investigators, most efforts seem to have failed in correctly capturing the flow physics, especially the slug-flow regime characteristics. The presence of a thin liquid film in the order of 10 μm around the bubble (sometimes called gas pocket or gas slug) may be a contributing factor to the above difficulty. Typically, liquid films have a significant effect on the flow field. Thus, there is a strong motivation to employ numerical simulation methods to avoid some of the experimental difficulties. In this paper, the characteristics of two-phase slug flows in microchannels are calculated with the help of the volume-of-fluid (VOF) method. Formation of the slugs for different superficial velocities, capillary numbers, and gas volume fractions are investigated. The minimum mesh resolution required to capture the liquid film surrounding the gas bubble is reported employing a dynamic mesh adaptation methodology with interface tracking. Results are shown to be in good agreement with experimental data and empirical correlations.

We present a numerical simulation of two-phase flow in a three-dimensional cross-junction microchannel by using the lattice Boltzmann method (LBM). At first, we validated our LBM code with the velocity profile in a 3-dimensional rectangular channel. Then, we developed a lattice Boltzmann code based on the free energy model to simulate the immiscible binary fluid flow. The parallelization of the developed code is implemented on a PC cluster using the MPI program. The numerical results of two-phase flow in the microchannel reveal droplet formation process, which compares well with corresponding experimental results. The size of droplet decreases with increase of the flow-rate ratio and the capillary number. The movement of a droplet through the microchannel induces three-dimensional circulating flow inside the droplet. This complex flow is thought to enhance the mixing and reaction of reagents.

Nitrogen/non-Newtonian fluid two-phase upward flow in non-circular microchannels

The effects of gas-liquid inlet geometry and mixing method on adiabatic gas-liquid two-phase flow in a microchannel of 100 μm diameter have been investigated using a T-junction inlet with the same internal diameter as the microchannel. Two-phase flow patterns, void fraction, and friction pressure drop data obtained with the T-junction inlet were found to be significantly different from those obtained previously with a reducer inlet. For the T-junction inlet, the two-phase flow patterns in the microchannel were predominantly intermittent flows with short gas and liquid plugs/slugs flowing with nearly equal velocities. The void fraction data then conformed nearly to that of a homogeneous flow model, and the two-phase friction multiplier data could be described by the Lockhart–Martinelli correlation applicable to larger channels. However, when a reducer inlet was used previously and the diameter of the inlet section was much larger than that of the microchannel, an intermittent flow of long gas slugs separated by long liquid slugs became prevalent and the void fraction decreased to values far below the homogeneous void fraction. The differences in the two-phase flow characteristics between a T-junction inlet and reducer inlet were attributed to the differences in the gas bubble/slug generation mechanisms.

. This study investigates the geometric effects of 90-degree vertical elbows and flow configurations in two-phase flow. The study shows that the elbows make a significant effect on the transport characteristics of two-phase flow, which includes the changes in interfacial structures, bubble interaction mechanisms and flow regime transition. The effect of the elbows is characterized for global and local two-phase flow parameters. The global two-phase flow parameters include two-phase pressure, interfacial structures and flow regime transition. In order to characterize the frictional pressure drop and minor loss across the vertical elbows, pressure measurements are obtained across the test section over a wide range of flow conditions in both single-phase and two-phase flow conditions. A two-phase pressure drop correlation analogous to Lockhart-Martinelli correlation is proposed to predict the minor loss across the elbows. A high speed camera is employed to perform extensive flow visualization studies across the elbows in vertical upward, horizontal and vertical downward sections and modified flow regime maps are proposed. It is found that modified flow regime maps immediately downstream of the vertical upward elbow deviate significantly from the conventional flow regime map. A qualitative assessment of the counter-current flow limitation characteristics specific to the current experimental facility is performed. A multi-sensor conductivity probe is used to measure local two-phase flow parameters such as: void fraction, bubble velocity, interfacial area concentration and bubble frequency. The local measurements are obtained for six different flow conditions at ten measurement locations along axial direction of the test section. Both the vertical-upward and vertical-downward elbows have a significant impact on bubble distribution, resulting in, a bimodal distribution along the horizontal radius of the tube cross-section and migration of bubbles towards the inside of the elbow curvatures immediately downstream of the vertical-upward and vertical-downward elbows, respectively. The elbow effect decays further downstream of the elbow and bubbles migrate to more conventional distribution patterns. The axial transport of void fraction and interfacial area concentration shows that the elbows promote bubble disintegration. Preliminary comparisons between the interfacial area transport model and the experimental data for verticalupward and vertical downward section are also presented.

Experimental investigation of R-134a and R-245fa two-phase flow in microchannels for different flow conditions

Numerical Investigation of Heat Transfer Enhancement due to Two-Phase Flow in Microchannel

Slug flow is the representative flow regime of two-phase flow in micro channels. It is well known that the thin liquid film formed around the confined vapor bubble plays an important role in micro channel heat transfer. In the present study, experiments are carried out to clarify the effects of parameters that affect the formation of the thin liquid film in micro channel slug flow. Laser focus displacement meter is used to measure the thickness of the thin liquid film. Air, ethanol, water and FC-40 are used as working fluids. Circular tubes with five different diameters, D = 0.3, 0.5, 0.7, 1.0, 1.3 mm, and square channels with two different sizes, 0.3 × 0.3 and 0.5 × 0.5 mm, are used. It is confirmed that the liquid film thickness is determined only by capillary number at small capillary numbers. However, the effect of inertial force and flow acceleration cannot be neglected as capillary number increases. The effect of cross sectional shape is also investigated. Experimental correlation for the adiabatic liquid film thickness in circular tubes based on capillary number, Reynolds number and Weber number is proposed. When viscous boundary layer is thin, liquid film thickness is limited by the viscous boundary layer thickness. Thus, in order to develop precise flow boiling models in micro tubes, it is important to consider the effects of inertial force and boundary layer thickness on the liquid film thickness.

Motivated by systematic CO2 evaporation experiments which recently became available (J. Pettersen, “Flow vaporization of CO2 in microchannel tubes,” Doctor technicae thesis, Norwegian University of Science and Technology, 2002), the present work constitutes an exploratory investigation of isothermal flow of CO2 near its liquid–vapor critical point through a long 5 μm diameter microchannel. A modified van der Waals constitutive model—with properties closely approximating those of “real” near-critical CO2—is incorporated in a two-dimensional lattice Boltzmann hydrodynamics model by embedding a dimensionless parameter X, with X→1 denoting the “real” fluid. The hydrodynamic phenomena resulting by imposing a constant pressure gradient along a periodic channel are investigated by considering two regimes in tandem: (1) transition from bubbly to annular flow with a liquid film formed at the channel walls and (2) destabilization of the liquid film by the Kelvin–Helmholtz instability. Due to numerical constraints, intrinsic modeling errors are introduced and are shown to be associated with discrepancies in the relative vapor–liquid interfacial thickness, which is expressed by X. The effects of these errors are investigated both theoretically and numerically in the physical limit X→1. Numerically determined flow patterns compare qualitatively well with direct visualization results obtained by Pettersen. Overall, the characteristics of isothermal near-critical two-phase flow in microchannels can be reproduced by the appropriate modification of the thermophysical properties of CO2.

We exhibit a single interval map (called a universal map) that admits all those orbit patterns which are available in the first Sharkovsky class. An interval map is said to be in the first Sharkovsky class if every periodic point of it is a fixed point. This provides a way to find universal maps in the class of contractions on intervals. We also characterize all such universal maps in the first Sharkovsky class.

With increasing computer capabilities, numerical modeling of two-phase flows has developed significantly over the last few years. Although there are two main categories, namely 'one' fluid and 'two' fluid methods, the 'one' fluid methods are more commonly used for tracking or capturing the interface between two fluids. Level set (LS), volume-of-fluid (VOF), front tracking, marker-and-cell (MAC) and lattice Boltzmann (LB) methods are all 'one' fluid methods. It is clear that there is no perfect method; each method has advantages and disadvantages which make it more appropriate for one kind of problem than for others. For instance, a LS method will accurately compute the curvature and the normal to the interface, but tends to loss mass which is physically incorrect. On the other hand, a VOF method will conserve mass up to machine precision, but the computation of the curvature and normal to the interface is not as accurate. In order to minimize the disadvantages of these methods, several authors have used two or more methods together to model two-phase flows. This is the case for the CLSVOF (Couple Level Set Volume Of Fluid) method, where LS and VOF are coupled together in order to better capture the interface. In CLSVOF, the level set function is used to compute the interface curvature and normal to the interface, while the volume of fluid function is used to capture the interface. For two-phase flows in microchannels, surface tension forces play an important role in determining the dynamics of bubbles whereas gravitational forces are generally negligible. Also it is very important to consider the interaction between the boundaries and the fluids by prescribing or computing the correct contact angle between them. The commercial CFD code FLUENT allows the use of static constant angles, or the use of User Defined Functions (UDF) to compute the dynamic contact angles. It is inappropriate to use a static contact angle to model cases involving moving contact lines. For such cases a dynamic contact angle scheme should be implemented. In this study, FLUENT was used to model adiabatic and diabatic, time dependent two-phase flows. Since FLUENT already contains a VOF method, a LS method was implemented and coupled with VOF into FLUENT via UDFs. Furthermore, since the LS function, used to compute the surface tension force, ceases to be a signed distance to the interface even after one time step, a re-initialization equation was solved after each time step. This involved using a fifth order WENO (Weighted Essentially Non Oscillatory) scheme to discretize the space derivatives (otherwise oscillations of the interface occurred), and a first order Euler method for the time integration. In another part of the study, a 3D dynamic contact angle model based on volume fraction, interface reconstruction, and experimentally available advancing and receding static contact angles was also developed and implemented into FLUENT via UDFs. Several validations for the developed CLSVOF method and dynamic contact angle model are presented in this thesis, these includes a static bubble, a bubble rising in a stagnant liquid for Morton numbers ranging from 102 to 10-11, droplet deformation due to a vortex flow field, droplets spreading over a wall under the gravity effect and droplets sliding over a wall due to gravity. These validations demonstrated the high accuracy and the stability of our methods for modeling these phenomena. A heat and mass transfer model was also implemented into the commercial CFD code FLUENT for simulating of boiling (and condensation) heat transfer. Several simulations were presented with water and R134a as working fluids. The influence of the contact angle and the wall superheat was also studied.

We have experimentally studied a two-phase flow in a microchannel with a height of 150 μm and a width of 20 mm. Different liquids have been used, namely, a purified Milli-Q water, an 50% aqueous-ethanol solution, and FC-72. Before and after the experiment, the height of the microchannel was controlled, as well as the wettability of its walls and surface tension of liquids. Using the schlieren method, the main characteristics of two-phase flow in wide ranges of gas- and liquid-flow rates have been revealed. The flow regime-formation mechanism has been found to depend on the properties of the liquid used. The flow regime has been registered when the droplets moving along the microchannel are vertical liquid bridges. It has been shown that, when using FC-72 liquid, a film of liquid is formed on the upper channel wall in the whole range of gas- and liquid-flow rates.

The effect of channel diameter on adiabatic two-phase flow characteristics in microchannels

A numerical investigation is carried out to study hydrodynamic behavior and thermal performance of slug flow in serpentine microchannel with different curvature ratio and wall properties. By using coupled Level Set and Volume of Fluid (CLSVOF) method implemented in Fluent, gas-liquid interface shows good consistency with experimental results especially for slug flow. Distribution of sectional velocity profile in curved microchannel with various curvature ratios was investigated. On the basis of fluid flow in microchannel with non-slip boundary condition, by means of user defined function (UDF), influence of slip velocity on hydrodynamic characteristics was performed numerically. For a better prediction on thermal performance, further analysis of combined effect of wall properties and curvature ratio on fluid flow characteristics was also conducted. As a consequence, the increasing slip velocity on channel wall can reduce hydraulic resistance. Both the contact angle and the curvature ratio have an appreciable influence on thermal performance of gas-liquid two-phase flow in serpentine microchannel.

Analysis of two-phase flow regimes and parameter distribution in 5 x 5 rod bundle using image processing techniques