Two-phase gas-liquid flow regimes for smooth and ribbed rectangular ducts

In this research, adiabatic two-phase air-water flow was investigated, and results for smooth and ribbed rectangular ducts are presented here. The test fluids were air and water at approximately atmospheric conditions. Three ribs of different heights were used; the rib width and pitch were held constant. The ribs were positioned in the duct at three different locations to establish three different conditions: on the bottom wall (water side), on the top wall (air side) and on both the top and bottom walls. The flow regimes in the smooth and ribbed ducts, which were recorded with a video camera, were classified as plug, stratified, slug and wavy flow. The location of the ribs in the duct did not alter the shape of the flow regimes, but the regime boundaries were considerably changed (repositioned). The effects of using ribs of different heights on regime boundaries are represented with flow map diagrams and discussed in detail. Compared to the smooth duct, the ribbed duct had different regime boundary positions. Increasing the rib height initiated hydrodynamical instability at lower fluid velocities. These findings are relevant for the operation and design of pipes, boilers and heat exchangers.     

Thermal and friction characteristics of laminar flow through rectangular and square ducts with transverse ribs and wire coil inserts

The heat transfer and the pressure drop characteristics of laminar flow of viscous oil (195 < Pr < 525) through rectangular and square ducts with internal transverse rib turbulators on two opposite surfaces of the ducts and with wire coil inserts have been studied experimentally. Circular duct has also been used. The transverse ribs in combination with wire coil inserts have been found to perform better than either ribs or wire coil inserts acting alone. The heat transfer and the pressure drop measurements have been taken in separate test sections. Heat transfer tests were carried out in electrically heated stainless steel ducts incorporating uniform wall heat flux boundary conditions. Pressure drop tests were carried out in acrylic ducts. The flow friction and thermal characteristics are governed by duct aspect ratio, coil helix angle and wire diameter of the coil, rib height and rib spacing, Reynolds number and Prandtl number. Correlations developed for friction factor and Nusselt number have predicted the experimental data satisfactorily. The performance of the geometry under investigation has been evaluated. It has been found that on the basis of constant pumping power, up to fifty per cent heat duty increase occurs for the combined ribs and wire coil inserts case compared to the individual ribs and wire coil inserts cases in the measured experimental parameters space. On the constant heat duty basis, the pumping power has been reduced up to forty per cent for the combined enhancement geometry than the individual enhancement geometries.     

The Effects of duct height on heat transfer enhancement of a co-rotating type rectangular finned surface in duct

Experiments were performed to investigate the effect of duct height on heat transfer enhancement of a surface affixed with arrays (7 × 7) of short rectangular plate fins of a co-rotating type pattern in the duct. An infrared imaging system is used to measure detailed distributions of the heat transfer at the endwall along with the fin base. An infrared camera of TVS 8000 with 160 × 120 point In–Sb sensor was used to measure the temperature distributions in order to calculate the local heat transfer coefficients of the representative fin regions. Pressure drop and heat transfer experiments were performed for a co-rotating fin pattern varying the duct height from 20?50 mm. The friction factor calculated from the pressure drop shows that comparatively larger friction occurs for the smaller duct cases and the friction factor slowly decreases with increasing Reynolds number. The effect of duct height on the area-averaged heat transfer results show that heat transfer initially increases with duct height and then finally decreases with increasing the duct height. Detailed heat transfer analysis and iso-heat transfer coefficient contour gives a clear picture of heat transfer characteristics of the overall surface. The relative performance graph indicates that a 25 mm duct is the optimum duct height for the highest thermal performance. In addition, a significant thermal enhancement, 2.8?3.8 times the smooth surface, can be achieved at lower Reynolds number with a co-rotating fin pattern in the duct.     

Effects of rotation on flow in an asymmetric rib-roughened duct: LES study

We report on large-eddy simulations (LES) of fully-developed asymmetric flow in a duct of a rectangular cross-section in which square-sectioned, equally-spaced ribs oriented perpendicular to the flow direction, were mounted on one of the walls. The configuration mimics a passage of internal cooling of a gas-turbine blade. The duct flow at a Reynolds number Re = 15,000 (based on hydraulic diameter D_h and bulk flow velocity U₀) was subjected to clock-wise (stabilising) and anti-clock-wise (destabilising) orthogonal rotation at a moderate rotational number Ro = ΩD_h/U₀ = 0.3, where Ω is the angular velocity. The LES results reproduced well the available experimental results of Coletti et al. (2011) (in the mid-plane adjacent to the ribbed wall) and provided insight into the whole duct complementing the reference PIV measurement. We analyzed the effects of stabilising and destabilising rotation on the flow, vortical structures and turbulence statistics by comparison with the non-rotating case. The analysis includes the identification of depth of penetration of the rib-effects into the bulk flow, influence of flow three-dimensionality and the role of secondary motions, all shown to be strongly affected by the rotation and its direction.     

Thermohydraulics of turbulent flow through rectangular and square ducts with axial corrugation roughness and twisted-tapes with and without oblique teeth

The heat transfer and the pressure drop characteristics of turbulent flow of air (10,000 < Re < 100,000) through rectangular and square ducts with combined internal axial corrugations on all the surfaces of the ducts and with twisted-tape inserts with and without oblique teeth have been studied experimentally. The axial corrugations in combination with twisted-tapes of all types with oblique teeth have been found to perform better than those without oblique teeth in combination with axial corrugations. The heat transfer and the pressure drop measurements have been taken in separate test sections. Heat transfer tests were carried out in electrically heated stainless steel ducts incorporating uniform wall heat flux boundary conditions. Pressure drop tests were carried out in acrylic ducts. The flow friction and thermal characteristics are governed by duct aspect ratio, corrugation angle, corrugation pitch, twist ratio, space ratio, length, tooth horizontal length and tooth angle of the twisted-tape, Reynolds number and Prandtl number. Correlations developed for friction factor and Nusselt number have predicted the experimental data satisfactorily. The performance of the geometry under investigation has been evaluated. It has been found that on the basis of constant pumping power, up to 55% heat duty increase occurs for the combined axial corrugation and regularly spaced twisted-tape elements inserts with oblique teeth case compared to without oblique teeth twisted-tape inserts cases in the measured experimental parameters space. On the constant heat duty basis, the pumping power has been reduced up to 47% for the combined enhancement geometry than the individual enhancement geometries. However, full-length and short-length twisted-tapes with oblique teeth in combination with axial corrugations show only marginal improvements over the twisted-tapes without oblique teeth.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

Two-phase pressure drop and flow visualization of FC-72 in a silicon microchannel heat sink

The rapid development of two-phase microfluidic devices has triggered the demand for a detailed understanding of the flow characteristics inside microchannel heat sinks to advance the cooling process of micro-electronics. The present study focuses on the experimental investigation of pressure drop characteristics and flow visualization of a two-phase flow in a silicon microchannel heat sink. The microchannel heat sink consists of a rectangular silicon chip in which 45 rectangular microchannels were chemically etched with a depth of 276 μm, width of 225 μm, and a length of 16 mm. Experiments are carried out for mass fluxes ranging from 341 to 531 kg/m² s and heat fluxes from 60.4 to 130.6 kW/m² using FC-72 as the working fluid. Bubble growth and flow regimes are observed using high speed visualization. Three major flow regimes are identified: bubbly, slug, and annular. The frictional two-phase pressure drop increases with exit quality for a constant mass flux. An assessment of various pressure drop correlations reported in the literature is conducted for validation. A new general correlation is developed to predict the two-phase pressure drop in microchannel heat sinks for five different refrigerants. The experimental pressure drops for laminar-liquid laminar-vapor and laminar-liquid turbulent-vapor flow conditions are predicted by the new correlation with mean absolute errors of 10.4% and 14.5%, respectively.     

Gas–liquid pressure drop in vertical internally wavy 90° bend

Experiments of air water two-phase flow pressure drop in vertical internally wavy 90° bend have been carried out. The tested bends are flexible and made of stainless steel with inner diameter of 50 mm and various curvature radiuses of 200, 300, 400 and 500 mm. The experiments were performed under the following conditions of two-phase parameters; mass flux from 350 to 750 kg/m² s. Gas quality from 1% to 50% and system pressure from 4 to 7.5 bar. The results demonstrate that the effect of the above-mentioned parameters is very significant at high ranges of mass flow quality. Due to the increasing of two-phase flow resistance, energy dissipations, friction losses and interaction of the two-phases in the vertical internally wavy 90° bend the total pressure drops are perceptible about 2–5 times grater than that in smooth bends. Based on the mass and energy balance as well as the presented experimental results, new empirical correlation has been developed to calculate the two-phase pressure drop and hence the two-phase friction factor of the tested bends. The correlation includes the relevant primary parameter, fit the data well, and is sufficiency accurate for engineering purposes.     

One-group interfacial area transport equation and its sink and source terms in narrow rectangular channel

The characteristics of two-phase flow in a narrow rectangular channel are expected to be different from those in other channel geometries, because of the significant restriction of the bubble shape which, consequently, may affect the heat removal by boiling under various operating conditions. The objective of this study is to develop an interfacial area transport equation with the sink and source terms being properly modeled for the gas–liquid two-phase flow in a narrow rectangular channel. By taking into account the crushed characteristics of the bubbles a new one-group interfacial area transport equation was derived for the two-phase flow in a narrow rectangular channel. The random collisions between bubbles and the impacts of turbulent eddies with bubbles were modeled for the bubble coalescence and breakup respectively in the two-phase flow in a narrow rectangular channel. The newly-developed one-group interfacial area transport equation with the derived sink and source terms was evaluated by using the area-averaged flow parameters of vertical upwardly-moving adiabatic air–water two-phase flows measured in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm. The flow conditions of the data set covered spherical bubbly, crushed pancake bubbly, crushed cap-bubbly and crushed slug flow regimes and their superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. Good agreement with the average relative deviation of 9.98% was obtained between the predicted and measured interfacial area concentrations in this study.     

Impact of aspect ratio on flow boiling of water in rectangular microchannels

In this paper we focus on the impact of varying the aspect ratio of rectangular microchannels, on the overall pressure drop involving water boiling. An integrated system comprising micro-heaters, sensors and microchannels has been realized on (1 1 0) silicon wafers, following CMOS compatible process steps. Rectangular microchannels were fabricated with varying aspect ratios (width [W] to depth [H]) but constant hydraulic diameter of 142 ± 2 μm and length of 20 mm. The invariant nature of the hydraulic diameter is confirmed through two independent means: physical measurements using profilometer and by measuring the pressure drop in single-phase fluid flow. The experimental results show that the pressure drop for two-phase flow in rectangular microchannels experiences minima at an aspect ratio of about 1.6. The minimum is possibly due to opposing trends of frictional and acceleration pressure drops, with respect to aspect ratio. In a certain heat flux and mass flux range, it is observed that the two-phase pressure drop is lower than the corresponding single-phase value. This is the first study to investigate the effect of aspect ratio in two-phase flow in microchannels, to the best of our knowledge. The results are in qualitative agreement with annular flow model predictions. These results improve the possibility of designing effective heat-sinks based on two-phase fluid flow in microchannels.     

Equal quality distribution of gas–liquid two-phase flow by partial separation method

This paper proposes a new method for equal quality distribution of gas–liquid two-phase flow by partial separate-phase distribution with a dual-header distributor. The upper and liquid (lower) headers are interconnected with five vertical downward arms. A gas–liquid two-phase mixture enters the distributor from the upper header where most of the liquid of the mixture is removed through the downward arms into the liquid header. Hence, firstly, the remaining gas-rich fluid can be uniformly distributed into the outlet branches, and then secondly, the liquid collected in the liquid header can be uniformly re-distributed into the individual outlet branches. Because both distribution processes are conducted in the condition of single or near single-phase flow, mal-distribution of the two-phase flow is essentially eliminated, and a satisfactory equal quality distribution of gas–liquid two-phase flow is reached. Experiments were conducted in an air–water two-phase flow test loop. The inner diameter of the inlet pipe was 60 mm, the superficial velocity ranges of gas and liquid were 3–32 m/s and 0.02–0.17 m/s respectively, and the quality ranged from 0.02 to 0.44. The flow pattern in the inlet pipe included stratified flow, wavy stratified, slug flow, and annular flow. The experimental results showed that this new method could significantly improve the distribution performance of the two-phase flow. The maximum quality deviation between each outlet branch and the inlet pipe is less than ±1% under the conditions of stratified, wavy stratified and slug flows in the upper header, and less than ±5% in annular flow.     

Experimental study on turbulent mixing of spray droplets in crossflow

Centrifugal spray injected at various angles in gas crossflow has been studied experimentally using PIV visualization system and image-processing techniques. Experiments were carried out inside a rectangular duct (95 mm × 95 mm in cross-section) at ambient temperature and pressure, with different gas Reynolds numbers (vary from 12,900 to 45,000) and three injection angles (60°, 90° and 120°). The spray angle of the centrifugal nozzle is 80°, with D₃₂ of 80 μm. The instantaneous images of droplets distribution and the values of the line-averaged D₃₂ at different positions on the cross-sections along the flow field for each condition were obtained, and their flow field configurations were achieved. Quantitative assessments of mixing degree between two phases for different injection angles were determined using a spatial unmixedness parameter. It is found that the addition of droplets into the gas crossflow enhanced the turbulence intensity of the gas crossflow and caused different-scale vortices. The flow field structure, to a great extent, is dependent on the injection angle. The entrainment and centrifugal force of large vortex lead to uneven droplet distribution and moreover influence the mixing of droplets and gas crossflow. A better mixing result can be obtained with the injection angle of nozzles of 60°.     

Liquid jet breakup for non-circular orifices under low pressures

This study uses a high-speed visualization technique to investigate the breakup process and flow behavior of low pressure water jets issued from non-circular orifices including square, triangular, and rectangular shapes. These orifices have approximately the same sectional areas. Stability curve and Ohnesorge chart are employed to make a comparison with circular jets discharged from a circular orifice of the same sectional area. The analysis is carried out for gauge pressures varying from 0.1 psi to 70 psi with small pressure steps corresponding to a range from 0.7 kPa to 482.6 kPa in metric units. Axis-switching phenomenon is observed and analyzed through calculating the axis-switching wavelength and oscillation frequency for rectangular jets. It is found that results for circular jets agreed well with classic theory. Non-circular jets demonstrate enhanced instabilities as a whole compared to circular jets. The different behaviors of non-circular jets are reasonably explained by Rayleigh’s oscillation theory. Axis-switching and aspect-ratio effect in rectangular jets is found to slow down the increase of breakup-length in the Rayleigh breakup regime. Square and triangular jets are more susceptible to wind effects and they are more unstable especially at higher pressure conditions. This can be concluded from the shorter breakup-length and narrower transitional region from the Rayleigh regime to the wind-induced regime as compared to the circular and rectangular jets. Axis-switching wavelength of the rectangular jets is found to increase linearly with increasing jet velocity and oscillation frequency decreases correspondingly.     

An investigation of a model of the flow pattern transition mechanism in relation to the identification of annular flow of R134a in a vertical tube using various void fraction models and flow regime maps

In the present study, new experimental data are presented for literature on the prediction of film thickness and identification of flow regime during the co-current downward condensation in a vertical smooth copper tube having an inner diameter of 8.1 mm and a length of 500 mm. R134a and water are used as working fluids in the tube side and annular side of a double tube heat exchanger, respectively. Condensation experiments are done at mass fluxes of 300 and 515 kg m^?2 s^?1. The condensing temperatures are between 40 and 50 °C; heat fluxes are between 12.65 and 66.61 kW m^?2. The average experimental heat transfer coefficient of the refrigerant HFC-134a is calculated by applying an energy balance based on the energy transferred from the test section. A mathematical model by Barnea et al. based on the momentum balance of liquid and vapor phases is used to determine the condensation film thickness of R134a. The comparative film thickness values are determined indirectly using relevant measured data together with various void fraction models and correlations reported in the open literature. The effects of heat flux, mass flux, and condensation temperature on the film thickness and condensation heat transfer coefficient are also discussed for the laminar and turbulent flow conditions. There is a good agreement between the film thickness results obtained from the theoretical model and those obtained from six of 35 void fraction models in the high mass flux region of R134a. In spite of their different valid conditions, six well-known flow regime maps from the literature are found to be predictive for the annular flow conditions in the test tube in spite of their different operating conditions.     

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 side ratio on fluid flow and heat transfer from rectangular cylinders using the PANS method

A computational study of heat transfer from rectangular cylinders is carried out. Rectangular cylinders are distinguished based on the ratio of the length of streamwise face to the height of the cross-stream face (side ratio, R). The simulations were performed to understand the heat transfer in a flow field comprising separation, reattachment, vortex shedding and stagnation. The Partially-Averaged Navier–Stokes (PANS) modeling approach is used to solve the turbulent flow physics associated and the wall resolve approach is used for the near wall treatment because of the flow separation involved. The simulations were performed using a finite volume based opensource software, OpenFOAM, at Reynolds number (Re) = 22,000 for rectangular cylinder at constant temperature kept in an air stream. Two critical side ratios were obtained, R = 0.62 and 3.0. At R = 0.62, the maximum value of the drag coefficient (C_d) = 2.681 was observed which gradually reduced by 54% at R = 4.0. The base pressure coefficient and global Nusselt number also attained the maximum value at R = 0.62 and from R = 2.5 to 3.0 a sharp discontinuous increase by 140% in the Strouhal number was observed. At R = 0.62, it was observed that the separated flow reattaches at the trailing edge after rolling over the side face and therefore increases the overall Nusselt number. The phase averaging was also performed to analyze the unsteady behavior of heat transfer.     

Heat transfer enhancement in a tube with equilateral triangle cross sectioned coiled wire inserts

The heat transfer and pressure drop were experimentally investigated in a coiled wire inserted tube in turbulent flow regime. The coiled wire has equilateral triangular cross section and was inserted separately from the tube wall. The experiments were carried out with three different pitch ratios (P/D = 1, 2 and 3) and two different ratio of equilateral triangle length side to tube diameter (a/D = 0.0714 and 0.0892) at a distance (s) of 1 mm from the tube wall in the range of Reynolds number from 3500 to 27,000. Uniform heat flux was applied to the external surface of the tube and air was selected as fluid. The experimental results obtained from a smooth tube were compared with those from the studies in literature for validation of experimental set-up. The use of coiled wire inserts leads to a considerable increase in heat transfer and pressure drop over the smooth tube. The Nusselt number rises with the increase of Reynolds number and wire thickness and the decrease of pitch ratio. The highest overall enhancement efficiency of 36.5% is achieved for the wire with a/D = 0.0892 and P/D = 1 at Reynolds number of 3858. Consequently, the experimental results reveal that the best operating regime of all coiled wire inserts is detected at low Reynolds number, leading to more compact heat exchanger.     

Experimental study on local characteristics of oil–water dispersed flow in a vertical pipe

The local flow characteristics of oil–water dispersed flow in a vertical upward pipe were studied experimentally. The inner diameter and length of the test section are 40 mm and 3800 mm, respectively. A double-sensor conductivity probe was used to measure the local interfacial parameters, including interfacial area concentration, oil phase fraction, interfacial velocity, and oil drops Sauter mean diameter. The water flow rates varied from 0.12 m/s to 0.89 m/s, while the oil flow rates ranged from 0.024 m/s to 0.198 m/s. Typical radial profiles of interfacial area concentration, oil phase fraction, interfacial velocity, and oil drops Sauter mean diameter are presented. An interesting phenomenon is that the local and cross-section-averaged interfacial area concentrations display concave change with water flow rate under constant oil flow rate. The physical mechanism of such a variation is discussed in details.     

Use of zero-net-mass-flow for separation control in diffusing S-duct

Flow control using zero-net-mass-flow jets in an S-shaped diffusing duct was investigated. Experiments were conducted in a channel flow facility at a Reynolds number, Re = 4.1 × 10⁴ with particle image velocimetry measurements in the symmetry plane of the duct. In the natural configuration, separation of the boundary layer occurs in a region of the duct with an high degree of curvature. A stability analysis of the wall normal base flow at the location of the applied control is presented and estimates the most effective frequency of the actuator. Time-averaged velocity fields show total reattachment of the boundary layer using active flow control.     

One-dimensional interfacial area transport of vertical upward bubbly flow in narrow rectangular channel

The design and safety analysis for miniature heat exchangers, the cooling system of high performance microelectronics, research nuclear reactors, fusion reactors and the cooling system of the spallation neutron source targets requires the knowledge of the gas–liquid two-phase flow in a narrow rectangular channel. In this study, flow measurements of vertical upward air–water flows in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm were performed at seven axial locations by using the imaging processing technique. The local frictional pressure loss gradients were also measured by a differential pressure cell. In the experiment, the superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. The developing two-phase flow was characterized by the significant axial changes of the local flow parameters due to the bubble coalescence and breakup in the tested flow conditions. The existing two-phase frictional multiplier correlations such as Chisholm, 1967, Mishima et al., 1993 and Lee and Lee (2001) were verified to give a good prediction for the measured two-phase frictional multiplier. The predictions of the drift-flux model with the rectangular channel distribution parameter correlation of Ishii (1977) and several existing drift velocity correlations of Ishii, 1977, Hibiki and Ishii, 2003 and Jones and Zuber (1979) agreed well with the measured void fractions and gas velocities. The interfacial area concentration (IAC) model of Hibiki and Ishii (2002) was modified by taking the channel width as the system length scale and the modified IAC model could predict the IAC and Sauter mean diameter acceptably.