Wavelet decomposition method decoupled boiling/evaporation oscillation mechanisms over two to three timescales: A study for a microchannel with pin fin structure

Boiling/evaporation heat transfer in a microchannel with pin fin structure was performed with water as the working fluid. Simultaneous measurements of various parameters were performed. The chip wall temperatures were measured by a high spatial-time resolution IR image system, having a sensitivity of 0.02 °C. The flow pattern variations synchronously changed wall temperatures due to ultra-small Bi number. The wavelet decomposition method successfully identified the noise signal and decoupled various temperature oscillations with different amplitudes and frequencies. Three types of temperature oscillations were identified according to heat flux q and mass flux G. The first type of oscillation occurred at q/G < 0.62 kJ/kg. The approximation coefficient of wavelet decomposition decided the dominant cycle period which was ∼3 times of the fluid residence time in the microchannel, behaving the density wave oscillation characteristic. The detail coefficients of wavelet decomposition decided the dominant cycle period, which matched the flow pattern transition determined value well, representing the flow pattern transition induced oscillation. For the second type of oscillation, the wavelet decomposition decoupled the three oscillation mechanisms. The pressure drop oscillation caused the temperature oscillation amplitudes of 5–10 °C and cycle periods of 10–15 s. The density wave oscillation and flow pattern transition induced oscillation are embedded with both the pressure rise and decrease stages of the pressure drop oscillation. The third type of oscillation happened at q/G > 1.13 kJ/kg, having the density wave oscillation coupled with the varied liquid film evaporation induced oscillation. The liquid island, retention bubble induced nucleation sites and cone-shape two-phase developing region are unique features of microchannel boiling with pin fin structure. This study illustrated that pressure drop oscillation and density wave oscillation, usually happened in large size channels, also take place in microchannels. The flow pattern transition and varied liquid film evaporation induced oscillations are specific to microchannel boiling/evaporation flow.     

New contributions on laminar flow of inelastic non-Newtonian fluid in the two-dimensional symmetric expansion: Creeping and slowly moving flow conditions

The steady flow of generalized Newtonian fluid in a two-dimensional 1:3 sudden expansion was studied numerically. Finite volume method was applied to solve the momentum equations along with the continuity equation and the Power law rheological model within the laminar flow regime for a range of Reynolds number and Power law index values. The values of generalized Reynolds number, based on physical and rheological properties, upstream channel height and bulk velocity, were varied between 0.0001 ≤ Re_gen ≤ 10, while the Power law index values mapped the 0.60 ≤ n ≤ 1.40 range, allowing for the investigation of both shear-thinning and shear-thickening effects at creeping as well as slowly moving fluid flow conditions. We report accurate results of a systematic study with a focus on most important characteristics of recirculating fluid flow in the downstream section of sudden expansion geometry. It is shown that for the creeping flow regime there exist finite sized redevelopment length, extra pressure drop (Couette correction) and recirculation zones (also called as Moffatt vortices) that are influenced by the non-Newtonian viscous behaviour.     

Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal 1.030 mm circular channel

This research focuses on acquiring accurate flow boiling heat transfer data and flow pattern visualization for three refrigerants, R134a, R236fa and R245fa in a 1.030 mm channel. We investigate trends in the data, and their possible mechanisms, for mass fluxes from 200 to 1600 kg/m²s, heat fluxes from 2.3 kW/m² to 250 kW/m² at T_sat = 31 °C and ΔT_sub from 2 to 9 K. The local saturated flow boiling heat transfer coefficients display a heat flux and a mass flux dependency but no residual subcooling influence. The changes in heat transfer trends correspond well with flow regime transitions. These were segregated into the isolated bubble (IB) regime, the coalescing bubble (CB) regime, and the annular (A) regime for the three fluids. The importance of nucleate boiling and forced convection in these small channels is still relatively unclear and requires further research.     

An investigation into flow mode transition and pressure fluctuations for fluidized dense-phase pneumatic conveying of fine powders

This paper presents the results of an ongoing investigation into the fluctuations of pressure signals due to solids–gas flows for dense-phase pneumatic conveying of fine powders. Pressure signals were obtained from pressure transducers installed along different locations of a pipeline for the fluidized dense-phase pneumatic conveying of fly ash (median particle diameter 30 μm; particle density 2300 kg/m³; loose-poured bulk density 700 kg/m³) and white powder (median particle diameter 55 μm; particle density 1600 kg/m³; loose-poured bulk density 620 kg/m³) from dilute to fluidized dense-phase. Standard deviation and Shannon entropy were employed to investigate the pressure signal fluctuations. It was found that there is an increase in the values of Shannon entropy and standard deviation for both of the products along the flow direction through the straight pipe sections. However, both the Shannon entropy and standard deviation values tend to decrease after the flow through bend(s). This result could be attributed to the deceleration of particles while flowing through the bends, resulting in dampened particle fluctuation and turbulence. Lower values of Shannon entropy in the early parts of the pipeline could be due to the non-suspension nature of flow (dense-phase), i.e., there is a higher probability that the particles are concentrated toward the bottom of pipe, compared with dilute-phase or suspension flow (high velocity), where the particles could be expected to be distributed homogenously throughout the pipe bore (as the flow is in suspension). Changes in straight-pipe pneumatic conveying characteristics along the flow direction also indicate a change in the flow regime along the flow.     

Origin of regime transition to turbulent flow in bubble column: Orifice- and column-induced transitions

The possible events during bubble formation on an orifice were investigated using a rectangular bubble column (30 cm × 30 cm × 100 cm). The gas flow rate through a single orifice was adjusted from 0.1 dm³/min to 5.0 dm³/min covering a high flow rate regime. At the high gas flow rate, the bubble formation process was complicated by diverse events, such as wake effect, channeling, and orifice-induced turbulent flow. The detachment period could be used to discern the bubble formation steps because it was strongly affected by the above events. The bubble size distribution around the orifice was also analyzed to gain a clearer understanding of the bubble formation process. Above the rate of 3.0 dm³/min through a single orifice, the detachment period converged to a value of 25 ms irrespective of the orifice diameter. The bubble size distribution also showed little difference in this range of gas flow rate. This could be explained by the development of turbulent flow around the orifice. A 0.15 m in-diameter bubble column was tested to investigate the effect of orifice-induced turbulent flow on the regime transition in which the homogeneous flow regime is converted into the heterogeneous flow regime in the column. Obvious distinction between the orifice- and column-induced transitions was observed.     

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.     

Influence of sudden contractions on in situ volume fractions for oil–water flows in horizontal pipes

Oil–water two-phase flow experiments were conducted in horizontal ducts made of Plexiglas® to determine the in situ oil fraction (holdup) by means of the closing valves technique, using mineral oil (viscosity: 0.838 Pa s at 20 °C; density: 890 kg m⁻³) and tap water. The ducts present sudden contractions from 50 mm to 40 mm i.d. and from 50 mm to 30 mm i.d., with contraction ratios of 0.64 and 0.36, respectively. About 200–320 tests were performed by varying the flow rates of the phases. Flow patterns were investigated for both the up- and downstream pipe in order to assess whether relevant variations of the flow patterns across the sudden contraction take place. Data were then compared with predictions of a specific correlation for oil–water flow and some correlations for gas–water flow. A drift-flux model was also applied to determine the distribution parameter.     

Technique for quantitative mapping of three-dimensional liquid–gas phase boundaries in microchannel flows

A diagnostic technique capable of characterizing interfaces between transparent, immiscible fluids is developed and demonstrated by investigating the morphology of liquid–gas interfaces in an adiabatic two-phase flow through a microchannel of 500 μm × 500 μm square cross section. Water seeded with 0.5 μm-diameter fluorescent polystyrene particles is pumped through the channel, and the desired adiabatic two-phase flow regime is achieved through controlled air injection. The diagnostic technique relies on obtaining particle position data through epifluorescent imaging of the flow at excitation and emission wavelengths of 532 nm and 620 nm, respectively. The particle position data are then used to resolve interface locations to within ±1 μm in the focal plane. By mapping the interface within individual focal planes at various depths within the channel, it is possible to determine the complete liquid–gas interface geometry across the channel cross section in a dynamic flow environment. Utilizing this approach, the liquid–gas phase boundaries of annular flows within a microchannel have been successfully characterized.     

Two-phase flow operational maps for multi-microchannel evaporators

The current paper presents new operational maps for several different multi-microchannel evaporators, with and without any inlet restrictions (micro-orifices), for the two-phase flow of refrigerants R245fa, R236fa, and R1234ze(E). The test fluids flowed in 67 parallel channels, each having a cross-sectional area of 100 × 100 μm². In order to emulate the power dissipated by active components in a 3D CMOS CPU chip, two aluminium microheaters were sputtered onto the back-side of the test section providing a 0.5 cm² each. Without any inlet restrictions in the micro-evaporator, significant parallel channel flow instabilities, vapor back flow, and flow maldistribution led to high-amplitude and high-frequency temperature and pressure oscillations. Such undesired phenomena were then prevented by placing restrictions at the inlet of each channel. High-speed flow visualization distinguished eight different operating regimes of the two-phase flow depending on the tested operating conditions. Therefore, the preferred operating regimes can be easily traced. In particular, flashing two-phase flow without back flow appeared to be the best operating regime without any flow and temperature instabilities.     

Bubbly-to-cap bubbly flow transition in a long-26 m vertical large diameter pipe at low liquid flow rate

The concurrent upward two-phase flow of air and water in a long vertical large diameter pipe with an inner diameter (D) of 200 mm and a height (z) of 26 m (z/D = 130) was investigated experimentally at low superficial liquid velocities from 0.05009 to 0.3121 m/s and the superficial gas velocities from 0.01779 to 0.5069 m/s. The resultant void fractions range from 0.03579 to 0.4059. According to the observations using a high speed video camera, the flow regimes of bubbly, developing cap bubbly and fully-developed cap bubbly flows prevailed in the flows. The developing cap bubbly flow appeared as a flow regime transition from bubbly to fully-developed cap bubble flow in the vertical large diameter pipe. The developing cap bubbly flow changes gradually and lasts for a long time period and a wide axial region in the flow direction, in contrast to a sudden transition from bubbly to slug flows in a small diameter pipe. The analysis in this study showed that the flow regime transition depends not only on the void fraction but also on the axial distance in the flow and the pipe diameter. The axial flow development brings about the transition to happen in a lower void fraction flow and the increase of pipe diameter causes the transition to happen in a higher void fraction flow. The measured void fraction showed an N-shaped axial changing manner that the void fraction increases monotonously with axial position in the bubbly flow, decreases non-monotonously with axial position in the developing cap bubbly flow, and increases monotonously again with axial position in the fully-developed cap bubbly flow. The temporary void fraction decrease phenomenon in the transition region from bubbly to cap bubbly flow can be attributed to the formation of medium to large cap bubbles and their gradual growth into the maximum size of cap bubble and/or cluster of large cap bubbles in the developing cap bubbly flow. In order to predict the N-shaped axial void fraction changing behaviors in the flow regime transition from bubbly to cap bubbly flow, the existing 12 drift flux correlation sets for large diameter pipes are reviewed and their predictabilities are studied against the present experimental data. Although some drift flux correlation sets, such as those of Clark and Flemmer (1986) and Hibiki and Ishii (2003), can predict the present experimental data with reasonable average relative deviations, no drift flux correlation set for distribution parameter and drift velocity can give a reliable prediction for the observed N-shaped axial void fraction changing behaviors in the region from bubbly to cap bubbly flow in a vertical large diameter pipe.     

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.     

Experimental study on the direct contact condensation of the steam jet in subcooled water flow in a rectangular channel: Flow patterns and flow field

Direct contact condensation (DCC) of steam jet in subcooled water flow in a channel was experimentally studied. The main inlet parameters, including steam mass flux, water mass flux and water temperature were tested in the ranges of 200–600 kg/(m² s), 7–18 × 10³ kg/(m² s), 288–333 K, respectively. Two unstable flow patterns and two stable flow patterns were observed via visualization window by a high speed camera. The flow patterns were determined by steam mass flux, water mass flux and water temperature, and the relationship between flow patterns and flow field parameters was discussed. The results indicated that whether pressure or temperature distributions on the bottom wall of channel could represent different flow patterns. And the position of pressure peak on the bottom wall could almost represent the condensation length. The upper wall pressure distributions were mainly dependent on steam and water mass flux; and the upper wall temperature distributions were affected by the three main inlet parameters. Moreover, the bottom wall pressure and temperature distributions of different unstable flow patterns had similar characteristics while those of stable flow patterns were affected by shock and expansion waves. The underlying cause of transition between different flow patterns under different inlet parameters was reflected and discussed based on pressure distributions.     

The effects of channel diameter on flow pattern,void fraction and pressure drop of two-phase air–water flow in circular micro-channels

Two-phase air–water flow characteristics are experimentally investigated in horizontal circular micro-channels. Test sections are made of fused silica. The experiments are conducted based on three different inner diameters of 0.53, 0.22 and 0.15 mm with the corresponding lengths of 320, 120 and 104 mm, respectively. The test runs are done at superficial velocities of gas and liquid ranging between 0.37–42.36 and 0.005–3.04 m/s, respectively. The flow visualisation is facilitated by systems mainly including stereozoom microscope and high-speed camera. The flow regime maps developed from the observed flow patterns are presented. The void fractions are determined based on image analysis. New correlation for two-phase frictional multiplier is also proposed for practical applications.     

Characterisation of gas-liquid two-phase flow in minichannels with co-flowing fluid injection inside the channel,part I: unified mapping of flow regimes

Present research highlights the potential of apparatuses with integrated minichannel packings to intensify gas-liquid-solid contacting. Especially an operation of these devices within the Taylor flow regime gained extraordinary attention due to its excellent heat and mass transfer and the segmented flow characteristics. However, criteria for flow regime transitions are mainly developed from water-similar fluids and are contradictory which hinders uniform flow regime prediction.This work presents a systematic analysis of adiabatic gas-liquid downflow in a square minichannel of 1.0 mm hydraulic diameter. In the mixing zone located within the flow channel, gas was injected into the co-flowing liquid by so-called capillary injectors with variable inner diameter (0.184, 0.317, 0.490 mm). Experiments were conducted using water, water-glycerol, and water-ethanol mixtures to cover a broad range of material properties. The gas and liquid superficial velocities were varied between 9.81·10^-4…2.72 m/s and 1.7·10⁻⁴…0.80 m/s, respectively. Taylor flow, Taylor-annular flow, annular flow, churn flow, and bubbly flow were observed. Using the Pi-theorem, 8 significant dimensionless groups dictating the flow transition were identified, namely u_{G, s}/u_{L, s}, Re_G, Re_L, We_G, We_L, Θ*, d_{In, CI}/d_h, and d_{Ou, CI}/d_h. Based on more than 1500 experimental data, criteria for the regime transitions of Taylor flow are provided. The derived flow regime map shows good agreement for all applied liquids and for the two larger injector geometries.     

Investigation of slip boundary conditions in the T-shaped microchannel

Flow regimes and mixing performance in a T-type micromixer at high Reynolds numbers were studied by numerical solution of the Navier–Stokes equations. The Reynolds number was varied from 1 to 1000. The cross section of the mixing channel was 100 μm × 200 μm, and its length was 1400 μm. The transverse inlet channels were symmetric about the mixing channel, and their cross-section was 100 μm × 100 μm, and the total length was 800 μm. Five different flow regimes were identified: (i) stationary vortex-free flow (Re < 5); (ii) stationary symmetric vortex flow with two horseshoe vortices (5 < Re < 150); (iii) stationary asymmetric vortex flow (150 < Re < 240); (iv) non-stationary periodic flow (240 < Re < 400); and (v) stochastic flow (Re > 400). Maximum mixing efficiency is obtained for nonstationary asymmetric vortex flow. In this case, an S-shaped vortex structure is formed in the flow field. The effect of the slip conditions on the flow pattern and mixing efficiency is studied. The slip length varied from 1 to 70 μm in the calculations. It was shown that the mixing can be controlled by hydrophobic coating.     

The effect of aspect ratio in counter-current gas-liquid bubble columns: Experimental results and gas holdup correlations

It is generally admitted that the gas holdup is independent of the column dimensions and gas sparger design if three criteria are satisfied: the diameter of the bubble column is larger than 0.15 m, gas sparger openings are larger than 1–2 mm and the aspect ratio is larger than 5. This paper contributes to the existing discussion; in particular, the effect of the aspect ratio (within the range 1–15) in a counter-current gas-liquid bubble column has been experimentally studied and a new gas holdup correlation to estimate the influence of aspect ratio, operation mode and working fluid on the gas holdup has been proposed. The bubble column, equipped with a spider gas sparger, is 5.3 m in height, has an inner diameter of 0.24 m; gas superficial velocities in the range of 0.004–0.23 m/s have been considered, and, for the runs with water moving counter-currently to the gas phase, the liquid has been recirculated at a superficial velocity of −0.0846 m/s. Filtered air has been used as the gaseous phase in all the experiments, while the liquid phase has included tap water and different aqueous solutions of sodium chloride as electrolyte. Gas holdup measurements have been used to investigate the flow regime transitions and the global bubble column hydrodynamics. The counter-current mode has turned out to increase the gas holdup and destabilize the homogeneous flow regime; the presence of electrolytes has resulted in increasing the gas holdup and stabilizing the homogeneous flow regime; the aspect ratio, up to a critical value, has turned out to decrease the gas holdup and destabilize the homogeneous flow regime. The critical value of the aspect ratio ranged between 5 and 10, depending on the bubble column operation (i.e., batch or counter-current modes) and liquid phase properties. Since no correlation has been found in the literature that can correctly predict the gas holdup under the investigated conditions, a new scheme of gas holdup correlation has been proposed. Starting from considerations concerning the flow regime transition, corrective parameters are included in the gas holdup correlation to account for the effect of the changes introduced by the aspect ratio, operation mode and working fluid. The proposed correlation has been found to predict fairly well the present experimental data as well as previously published gas holdup data.     

Investigation of two-phase flow pattern,liquid holdup and pressure drop in viscous oil–gas flow

An experimental investigation was carried out on viscous oil–gas flow characteristics in a 69 mm internal diameter pipe. Two-phase flow patterns were determined from holdup time-traces and videos of the flow field in a transparent section of the pipe, in which synthetic commercial oils (32 and 100 cP) and sulfur hexafluoride gas (SF₆) were fed at oil superficial velocities from 0.04 to 3 m/s and gas superficial velocities from 0.0075 to 3 m/s.     

Numerical analysis of turbulent convection heat transfer from an array of perforated fins

Numerical investigation is made for three-dimensional fluid flow and convective heat transfer from an array of solid and perforated fins that are mounted on a flat plate. Incompressible air as working fluid is modeled using Navier–Stokes equations and RNG based k ? ? turbulent model is used to predict turbulent flow parameters. Temperature field inside the fins is obtained by solving Fourier’s conduction equation. The conjugate differential equations for both solid and gas phase are solved simultaneously by finite volume procedure using SIMPLE algorithm. Perforations such as small channels with square cross section are arranged streamwise along the fin’s length and their numbers varied from 1 to 3. Flow and heat transfer characteristics are presented for Reynolds numbers from 2 × 10⁴ to 4 × 10⁴ based on the fin length and Prandtl number is taken Pr = 0.71. Numerical computations are validated with experimental studies of the previous investigators and good agreements were observed. Results show that fins with longitudinal pores, have remarkable heat transfer enhancement in addition to the considerable reduction in weight by comparison with solid fins.     

Experimental investigation of the effects of blowing conditions and Mach number on the unsteady behavior of coolant ejection through a trailing edge cutback

The present paper shows the results of an experimental investigation into the unsteadiness of coolant ejection at the trailing edge of a highly loaded nozzle vane cascade. The trailing edge cooling scheme features a pressure side cutback with film cooling slots, stiffened by evenly spaced ribs in an inline configuration. Cooling air is also ejected through two rows of cylindrical holes placed upstream of the cutback. Tests were performed with a low inlet turbulence intensity level (Tu₁ = 1.6%), changing the cascade operating conditions from low speed (M_2is = 0.2) up to high subsonic regime (M_2is = 0.6), and with coolant to main stream mass flow ratio varied within the 0.5–2.0% range. Particle Image Velocimetry (PIV) and flow visualizations were used to investigate the unsteady mixing process taking place between coolant and main flow downstream of the cutback, up to the trailing edge. For all the tested conditions, the results show the presence of large coherent structures, which presence is still evident up to the trailing edge. Their shape and direction of rotation change with injection conditions, as a function of coolant to mainstream velocity ratio, strongly influencing the thermal protection capability of the injected coolant flow. The Mach number increase is only responsible for a positioning of such vortical structures closer to the wall, while the Strouhal number almost remains unchanged.     

Experimental investigation into transient pressure pulses during pneumatic conveying of fine powders using Shannon entropy

This paper presents the results of an ongoing investigation into transient pressure pulses using Shannon entropy. Pressure fluctuations (produced by gas–solid two-phase flow during fluidized dense-phase conveying) are recorded by pressure transducers installed at strategic locations along a pipeline. This work validates previous work on identifying the flow mode from pressure signals (Mittal, Mallick, & Wypych, 2014). Two different powders, namely fly ash (median particle diameter 45 μm, particle density 1950 kg/m³, loosely poured bulk density 950 kg/m³) and cement (median particle diameter 15 μm, particle density 3060 kg/m³, loosely poured bulk density 1070 kg/m³), are conveyed through different pipelines (51 mm I.D. × 70 m length and 63 mm I.D. × 24 m length). The transient nature of pressure fluctuations (instead of steady-state behavior) is considered in investigating flow characteristics. Shannon entropy is found to increase along straight pipe sections for both solids and both pipelines. However, Shannon entropy decreases after a bend. A comparison of Shannon entropy among different ranges of superficial air velocity reveals that high Shannon entropy corresponds to very low velocities (i.e. 3–5 m/s) and very high velocities (i.e. 11–14 m/s) while low Shannon entropy corresponds to mid-range velocities (i.e. 6–8 m/s).