Phase split of nitrogen/non-Newtonian fluid two-phase flow at a micro-T-junction

An experimental investigation has been undertaken to understand the phase split of nitrogen gas/non-Newtonian liquid two-phase flow passing through a 0.5 mm T-junction that oriented horizontally. Four different liquids, including water and aqueous solutions of carboxymethyl cellulose (CMC) with different mass concentrations of 0.1, 0.2 and 0.3 wt%, were employed. Rheology experiments showed that different from water, CMC solutions in this study are pseudoplastic non-Newtonian fluid whose viscosity decreases with increasing the shear rate. The inlet flow patterns were observed to be slug flow, slug–annular flow and annular flow. The fraction of liquid taken off at the side arm for nitrogen gas/non-Newtonian liquid systems is found to be higher than that for nitrogen gas/Newtonian liquid systems in all inlet flow patterns. In addition, with increasing the pseudoplasticity of the liquid phase, the side arm liquid taken off increases, but the increasing degree varies with each flow pattern. For annular flow, the increasing degree is much greater than those for slug and slug–annular flows.     

Characterization of diabatic two-phase flows in microchannels: Flow parameter results for R-134a in a 0.5 mm channel

An optical measurement method for two-phase flow pattern characterization in microtubes has been utilized to determine the frequency of bubbles generated in a microevaporator, the coalescence rates of these bubbles and their length distribution as well as their mean velocity. The tests were run in a 0.5 mm glass channel using saturated R-134a at 30 °C (7.7 bar). The optical technique uses two laser diodes and photodiodes to measure these parameters and to also identify the flow regimes and their transitions. Four flow patterns (bubbly flow, slug flow, semi-annular flow and annular flow) with their transitions were detected and observed also by high speed video. It was also possible to characterize bubble coalescence rates, which were observed here to be an important phenomena controlling the flow pattern transition in microchannels. Two types of coalescence occurred depending on the presence of small bubbles or not. The two-phase flow pattern transitions observed did not compare well to a leading macroscale flow map for refrigerants nor to a microscale map for air–water flows. Time averaged cross-sectional void fractions were also calculated indirectly from the mean two-phase vapor velocities and compared reasonably well to homogeneous values.     

Macro-to-microchannel transition in two-phase flow: Part 2 - Flow boiling heat transfer and critical heat flux

This part of the paper presents the current experimental flow boiling heat transfer and CHF data acquired for R134a, R236fa and R245fa in single, horizontal channels of 1.03, 2.20 and 3.04 mm diameters over a range of experimental conditions. The aim of this study is to investigate the effects of channel confinement, heat flux, flow pattern, saturation temperature, subcooling and working fluid properties on the two-phase heat transfer and CHF. Experimentally, it was observed that the flow boiling heat transfer coefficients are a significant function of the type of two-phase flow pattern. Furthermore, the monotonically increasing heat transfer coefficients at higher vapor qualities, corresponding to annular flow, signifies convective boiling as the dominant heat transfer mechanism in these small scale channels. The decreasing heat transfer trend at low vapor qualities in the slug flow (coalescing bubble dominated regime) was indicative of thin film evaporation with intermittent dry patch formation and rewetting at these conditions. The coalescing bubble flow heat transfer data were well predicted by the three-zone model when setting the dryout thickness to the measured surface roughness, indicating for the first time a roughness effect on the flow boiling heat transfer coefficient in this regime. The CHF data acquired during the experimental campaign indicated the influence of saturation temperature, mass velocity, channel confinement and fluid properties on CHF but no influence of inlet subcooling for the conditions tested. When globally comparing the CHF values for R134a in the 0.51-3.04 mm diameter channels, a peak in CHF peak was observed lying in between the 0.79 (Co ≈ 0.99) and 1.03 (Co ≈ 0.78) mm channels. A new CHF correlation has been proposed involving the confinement number, Co that is able to predict CHF for R134a, R236fa and R245fa in single-circular channels, rectangular multichannels and split flow rectangular multichannels. In summary, the present flow boiling and CHF trends point to a macro-to-microscale transition as indicated by the results presented in Ong and Thome (2011) [1].     

Two-phase flow characteristics in gas–liquid microreactors

Multiphase chemical microreactors require a detailed knowledge of the flow conditions inside the reaction system. This paper reports flow visualization measurements of the two-phase gas–liquid flow pattern and the liquid velocity distribution inside liquid plugs of an intermittent flow. Rectangular cross-section silicon microchannels with hydraulic diameters between 187.5 and 218 μm are fabricated. Laser Induced Fluorescence (LIF) is used to determine the flow pattern. To analyze the influence of the liquid properties and the channel diameter on the two-phase flow pattern, we present flow regime maps using different channel geometries and fluids. A universal flow pattern map based on dimensional analysis is presented. In contrast to microchannel flows, a great number of correlations for flow characteristics for multiphase flow in (round) pipes with diameters >1 mm exist. We compare our experimental results from optical flow visualizations in microreactors with common flow correlations and regime maps for macro- and microchannels. The recirculation motion in the liquid segments of an intermittent gas–liquid flow is analyzed using micron-resolution particle image velocimetry (μPIV). The velocity distribution influences the mixing and the mass transport towards the reactive phase interface dealing with two-phase chemical reactions. For straight microchannels hardly any mass transport over the center line is quantified. For enhanced mixing geometrical adaptations are suggested.     

Flow regime development analysis in adiabatic upward two-phase flow in a vertical annulus

In this work radial and axial flow regime development in adiabatic upward air-water two-phase flow in a vertical annulus has been investigated. Local flow regimes have been identified using conductivity probes and neural networks techniques. The inner and outer diameters of the annulus are 19.1 mm and 38.1 mm, respectively. The equivalent hydraulic diameter of the flow channel, D_H, is 19.0 mm and the total length is 4.37 m. The flow regime map includes 1080 local flow regimes identifications in 72 flow conditions within a range of 0.01 m/s < 〈j_g〉 < 30 m/s and 0.2 m/s < 〈j_f〉 < 3.5 m/s where 〈j_g〉 and 〈j_f〉 are, respectively, superficial gas and liquid velocities. The local flow regime has been classified into four categories: bubbly, cap-slug, churn-turbulent and annular flows. In order to study the radial and axial development of flow regime the measurements have been performed at five radial locations. The three axial positions correspond to z/D_H = 52, 149 and 230, where z represents the axial position. The flow regime indicator has been chosen as some statistical parameters of local bubble chord length distributions and self-organized neural networks have been used as mapping system. This information has been also used to compare the results given by the existing flow regime transition models. The local flow regime is characterized basically by the void fraction and bubble chord length. The radial development of flow regime shows partial and complete local flow regime combinations. The radial development is controlled by axial location and superficial liquid velocity. The radial flow regime transition is always initiated in the center of the flow channel and it is propagated towards the channel boundaries. The axial development of flow regime is observed in all the flow maps and it is governed by superficial liquid velocity and radial location. The prediction results of the models are compared for each flow regime transition.     

Visualization of flow boiling of liquid nitrogen in a vertical mini-tube

The flow boiling patterns of liquid nitrogen in a vertical mini-tube with an inner diameter of 1.931 mm are visualized with a high-speed digital camera. The superficial gas and liquid velocities are in the ranges of 0.01–26.5 m/s and 0.01–1.2 m/s, respectively. Four typical flow patterns, namely, bubbly, slug, churn and annular flow are observed. Some interesting scenes about the entrainment and liquid droplet deposition in the churn and annular flow, and the flow reversal with the indication of negative pressure drop, are also presented. Based on the visualization, the two-phase flow regime maps are obtained. Compared with the flow regime maps for gas–water flow in tubes with similar hydraulic diameters, the region of slug flow in the present study reduces significantly. Correspondingly, the transition boundary from the bubbly flow to slug flow shifts to higher superficial gas velocity, and that of churn to annular flow moves to lower superficial gas velocity. Moreover, time-averaged void fraction is calculated by quantitative image-digitizing technique and compared with various prediction models. Finally, three kinds of oscillations with long-period and large-amplitude are found, possible explanation for the oscillations is given by comparing the instantaneous flow images with the data of pressure, mass flux and temperature recorded synchronously.     

Macro-to-microchannel transition in two-phase flow: Part 1 - Two-phase flow patterns and film thickness measurements

The classification of macroscale, mesoscale and microscale channels with respect to two-phase processes is still an open question. The main objective of this study focuses on investigating the macro-to-microscale transition during flow boiling in small scale channels of three different sizes with three different refrigerants over a range of saturation conditions to investigate the effects of channel confinement on two-phase flow patterns and liquid film stratification in a single circular horizontal channel (Part 2 covers the flow boiling heat transfer and critical heat flux). This paper presents the experimental two-phase flow pattern transition data together with a top/bottom liquid film thickness comparison for refrigerants R134a, R236fa and R245fa during flow boiling in small channels of 1.03, 2.20 and 3.04 mm diameter. Based on this work, an improved flow pattern map has been proposed by determining the flow patterns transitions existing under different conditions including the transition to macroscale slug/plug flow at a confinement number of Co ≈ 0.3-0.4. From the top/bottom liquid film thickness comparison results, it was observed that the gravity forces are fully suppressed and overcome by the surface tension and shear forces when the confinement number approaches 1, Co ≈ 1. Thus, as a new approximate rule, the lower threshold of macroscale flow is Co = 0.3-0.4 while the upper threshold of symmetric microscale flow is Co ≈ 1 with a transition (or mesoscale) region in-between.     

Experimental investigation of vapor bubble growth during flow boiling in a microchannel

Experiments were conducted to analyze flow boiling characteristics of water in a single brass microchannel of 25 mm length, 201 μm width, and 266 μm depth. Different heat flux conditions were tested for each of two different mass flow rates over three different values of inlet fluid temperature. Temporal and spatial surface temperature profiles were analyzed to show the relative effect of axial heat conduction on temperature rise along the channel length and the effect of flow regime transition on local surface temperature oscillation. Vapor bubble growth rate increased with increasing wall superheat. The slower a bubble grew, the further it was carried downstream by the moving liquid. Bubble growth was suppressed for increased mass flux while the vapor bubble was less than the channel diameter. The pressure spike of an elongating vapor bubble was shown to suppress the growth of a neighboring bubble by more than 50% of its volume. An upstream progression of the Onset of Bubble Elongation (OBE) was observed that began at the channel exit and progressed upstream. The effects of conjugate heat transfer were observed when different flow regime transitions produced different rates of progression for the elongation sequence. Instability was observed at lower heat fluxes for this single channel experiment than for similar studies with multiple channels.     

The characteristics of gas/liquid flow in large risers at high pressures

Although most of the work reported on two-phase flows are limited to small pipe diameters, two-phase flow in large risers are increasingly being encountered in the petroleum and nuclear industries. In the present work, a wire mesh sensor was employed to obtain void fraction and bubble size distribution data and visualizations of steam/water flow in a large vertical pipe (194 mm in diameter) at 46 bar. For comparison purposes, measurements were made at similar phase velocities and physical properties to a dataset for nitrogen/naphtha flow in a similar-sized riser. There exist significant differences between both sets of data. Churn-turbulent flow is observed in the present work instead of slug flow, and this differs from the intermittent and semi-annular flow patterns reported for nitrogen/naphtha data. The mean void fraction of the nitrogen/naphtha data is higher than that of the present steam/water data due to the differences in purity in the liquid phases. Furthermore, core peak distributions are observed for the present work in contrast to the flatter profiles deduced for the nitrogen/naphtha using a power law relationship.     

Effects of channel dimension,heat flux,and mass flux on flow boiling regimes in microchannels

Experiments are conducted with a perfluorinated dielectric fluid, Fluorinert FC-77, to investigate the effects of channel size and mass flux (225–1420 kg/m²s) on microchannel flow boiling regimes by means of high-speed photography. Seven different silicon test pieces with parallel microchannels of widths ranging from 100 to 5850 μm, all with a depth of 400 μm, are considered. Flow visualizations are performed with a high-speed digital video camera while local measurements of the heat transfer coefficient are simultaneously obtained. The visualizations and the heat transfer data show that flow regimes in the microchannels of width 400 μm and larger are similar, with nucleate boiling being dominant in these channels over a wide range of heat flux. In contrast, flow regimes in the smaller microchannels are different and bubble nucleation at the walls is suppressed at a relatively low heat flux for these sizes. Two types of flow regime maps are developed and the effects of channel width on the flow regime transitions are discussed.     

Interaction between Taylor bubbles rising in stagnant non-Newtonian fluids

The interaction between Taylor bubbles rising in stagnant non-Newtonian solutions was studied. Aqueous solutions of carboxymethylcellulose (CMC) and polyacrylamide (PAA) polymers were used to study the effect of different rheological properties: shear viscosity and viscoelasticity. The solutions studied covered a range of Reynolds numbers between 10 and 714, and Deborah numbers up to 14. The study was performed with pairs of Taylor bubbles rising in a vertical column (0.032 m internal diameter) filled with stagnant liquid. The velocities of the leading and trailing bubbles were measured by sets of laser diodes/photocells placed along the column. The velocity of the trailing bubble was analysed together with the liquid velocity profile in the wake of a single rising bubble (Particle Image Velocimetry data obtained from the literature). For the less concentrated CMC solutions, with moderate shear viscosity and low viscoelasticity, the interaction between Taylor bubbles was similar to that found in Newtonian fluids. For the most concentrated CMC solution, which has high shear viscosity and moderate viscoelasticity, a negative wake forms behind the Taylor bubbles, inhibiting coalescence since the bubbles maintain a minimum distance of about 1D between them. For the PAA solutions, with moderate shear viscosity but higher viscoelasticity than the CMC solutions, longer wake lengths are seen, which are responsible for trailing bubble acceleration at greater distances from the leading bubble. Also in the PAA solutions, the long time needed for the fluid to recover its initial shear viscosity after the passage of the first bubble makes the fluid less resistant to the trailing bubble flow. Hence, the trailing bubble can travel at a higher velocity than the leading bubble, even at distances above 90D.     

Two-phase flow patterns of a top heat mode closed loop oscillating heat pipe with check valves (THMCLOHP/CV)

This research is aimed at studying the two-phase flow pattern of a top heat mode closed loop oscillating heat pipe with check valves. The working fluids used are ethanol and R141b and R11 coolants with a filling ratio of 50% of the total volume. It is found that the maximum heat flux occurs for the R11 coolant used as the working fluid in the case with the inner diameter of 1.8 mm, inclination angle of ?90?, evaporator temperature of 125?C, and evaporator length of 50 mm. The internal flow patterns are found to be slug flow/disperse bubble flow/annular flow, slug flow/disperse bubble flow/churn flow, slug flow/bubble flow/annular flow, slug flow/disperse bubble flow, bubble flow/annular flow, and slug flow/annular flow.     

Flow pattern characterization for R-245fa in minichannels: Optical measurement technique and experimental results

An optical measurement method using image processing for two-phase flow pattern characterization in minichannel is developed. The bubble frequency, the percentage of small bubbles as well as their velocity are measured. A high-speed high-definition video camera is used to measure these parameters and to identify the flow regimes and their transitions. The tests are performed in a 3.0 mm glass channel using saturated R-245fa at 60 °C (4.6 bar). The mass velocity is ranging from 100 to 1500 kg/m² s, the heat flux is varying from 10 to 90 kW/m² and the inlet vapor quality from 0 to 1. Four flow patterns (bubbly flow, bubbly–slug flow, slug flow and annular flow) are recognized. The comparison between the present experimental intermittent/annular transition lines and five transition lines from macroscale and microscale flow pattern maps available in the literature is presented. Finally, the influence of the flow pattern on the heat transfer coefficient is highlighted.     

Effect of diameter on two-phase pressure drop in narrow tubes

The effect of tube diameter on two-phase frictional pressure drop was investigated in circular tubes with inner diameters of 0.6, 1.2, 1.7, 2.6 and 3.4 mm using air and water. The gas and liquid superficial velocity ranges were 0.01-50 m/s and 0.01-3 m/s, respectively. The gas and liquid flow rates were measured and the two-phase flow pattern images were recorded using high-speed CMOS camera. Unique flow patterns were observed for smaller tube diameters. Pressure drop was measured and compared with various existing models such as homogeneous model and Lockhart-Martinelli model. It appears that the dominant effect of surface tension shrinking the flow stratification in the annular regime is important. It was found that existing models are inadequate in predicting the pressure drop for all the flow regimes visualized. Based on the analysis of present experimental frictional pressure drop data a correlation is proposed for predicting Chisholm parameter “C” in slug annular flow pattern. For all other flow regimes Chisholm’s original correlation appears to be adequate except the bubbly flow regime where homogeneous model works well. The modification results in overall mean deviation of pressure drop within 25% for all tube diameters considered. This approach of flow regime based modification of liquid gas interaction parameter appears to be the key to pressure drop prediction in narrow tubes.     

Gas/shear-thinning liquid flows through pipes: Modeling and experiments

In chemical and oil industry gas/shear-thinning liquid two-phase flows are frequently encountered. In this work, we investigate experimentally the flow characteristics of air/shear-thinning liquid systems in horizontal and slightly inclined smooth pipes. The experiments are performed in a 9-m-long glass pipe using air and three different carboxymethyl cellulose (CMC) solutions as test fluids. Flow pattern maps are built by visual observation using a high-speed camera. The observed flow patterns are stratified, plug, and slug flow. The effects of the pipe inclination and the rheology of the shear-thinning fluid in terms of flow pattern maps are presented. The predicted existence region of the stratified flow regime is compared with the experimental observation showing a good agreement. A mechanistic model valid for air/power-law slug flow is proposed and model predictions are compared to the experimental data showing a good agreement. Slug flow characteristics are investigated by the analysis of the signals of a capacitance probe: slug velocity, slug frequency, and slug lengths are measured. A new correlation for the slug frequency is proposed and the results are promising.     

Flow regime transition criteria for two-phase flow at reduced gravity conditions

Flow regime transition criteria are of practical importance for two-phase flow analyses at reduced gravity conditions. Here, flow regime transition criteria which take the frictional pressure loss effect into account were studied in detail. Criteria at reduced gravity conditions were developed by extending an existing model from normal gravity to reduced gravity conditions. A comparison of the newly developed flow regime transition criteria model with various experimental datasets taken at microgravity conditions showed satisfactory agreement. Sample computations of the model were performed at various gravity conditions, such as 0.196, 1.62, 3.71 and 9.81 m/s² corresponding to micro-gravity and lunar, Martian and Earth surface gravity, respectively. It was found that the effect of gravity on bubbly–slug and slug–annular (churn) transitions in a two-phase flow system was more pronounced at low liquid flow conditions, whereas the gravity effect could be ignored at high mixture volumetric flux conditions. While for the annular flow transitions due to flow reversal and onset of droplet entrainment, higher superficial gas velocity was obtained at higher gravity level.     

Effects of the gas phase molecular weight and bubble size on effervescent atomization

The effervescent atomization from an industrial Coker feed nozzle is compared for two different gas densities (air and mixed gas of 81.4 vol.% helium/18.6 vol.% nitrogen) at equivalent operating temperatures. The application is to observe the similarity of lab tests using air at 20 °C to the industrial process using steam at 300-400 °C. The effects of operating conditions, such as gas to liquid mass ratio, mixing pressure and void fraction on the flow regime, bubble size, and droplet size distribution were also examined in this study. The experiments were performed using mixtures of water with air or mixed gas, which resulted in gas to liquid mass ratios ranging from 1% to 4%.Stroboscopic back scattered imagery (SBSI) indicates that the average bubble size inside the nozzle conduit is similar when air and water are used as the process fluids, when compared to the case when mixed gas and water are used as the process fluids. Under similar conditions, the Phase Doppler Particle Anemometer (PDPA) data indicate that the droplet size in the spray is similar when using either mixed gas or air as the atomization gas.Experimental results obtained by high-speed video shadowgraphy (HSVS) indicate that the flow pattern inside the nozzle feeding conduit was slug flow with a tendency to attain annular flow with increased air to liquid mass ratios. Thus, from the experimental results it is evident that the smaller molecular weight of the mixed gas versus air (8.4 versus 29) does not significantly reduce the bubble (<±10% difference) and droplet size (<±1.5% difference), indicating a weak dependence of the gas phase density on two-phase atomization. This confirms that laboratory experiments on effervescent nozzles using air have reliable similarity to systems that use high temperature steam for the gas phase.     

Numerical modeling and experimental investigation of gas–liquid slug formation in a microchannel T-junction

Gas–liquid two-phase flow in a microfluidic T-junction with nearly square microchannels of 113 μm hydraulic diameter was investigated experimentally and numerically. Air and water superficial velocities were 0.018–0.791 m/s and 0.042–0.757 m/s, respectively. Three-dimensional modeling was performed with computational fluid dynamics (CFD) software FLUENT and the volume of fluid (VOF) model. Slug flow (snapping/breaking/jetting) and stratified flow were observed experimentally. Numerically predicted void fraction followed a linear relationship with the homogeneous void fraction, while experimental values depended on the superficial velocity ratio U_G/U_L. Higher experimental velocity slip caused by gas inlet pressure build-up and oscillation caused deviation from numerical predictions. Velocity slip was found to depend on the cross-sectional area coverage of the gas slug, the formation of a liquid film and the presence of liquid at the channel corners. Numerical modeling was found to require improvement to treat the contact angle and contact line slip, and could benefit from the use of a dynamic boundary condition to simulate the compressible gas phase inlet reservoir.     

Flow regime transition criteria for two-phase flow in a vertical annulus

In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug–bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been checked against air–water as well as boiling flow regime maps. In all the cases, the new developed model shows better predicting capabilities than the existing correlations most used in literature.     

Flow patterns of n-hexadecane–CO2 liquid–liquid two-phase flow in vertical pipes under high pressure

To promote a better understanding of liquid–liquid two-phase flow behavior, particularly under high pressure, flow patterns of n-hexadecane–CO₂ liquid–liquid two-phase upward flow in vertical stainless steel pipes were experimentally investigated. Observations were made in two 0.0015 m I.D. pipes of different lengths (0.068 m and 0.5 m) under high pressure varying from 10.3 to 29.6 MPa using a high pressure visualization system. The total flow rate was fixed at 2.0 × 10⁻⁶ m³/min, while the flow rate ratio (φ) varied from 0.05 to 19. Bubbly flow, plug flow, slug flow, annular flow, and near-one-phase flow regions were found in both pipes, while stratified flow was observed only in the 0.068 m pipe. Flow pattern maps were constructed in the flow rate ratio versus pressure graph, which demonstrates significant impacts of flow rate ratio, pipe length, and pressure on flow patterns. These impacts are discussed in detail. To the authors’ best knowledge, this work is the first attempt to observe complex liquid–liquid two-phase flow behavior with flow pattern transitions under high pressure, and contributes to a better understanding of liquid–liquid two-phase flow behavior.