Pressure drop and heat transfer characteristics of boiling water in sub-hundred micron channel

The current work focuses on the pressure drop, heat transfer and stability in two phase flow in microchannels with hydraulic diameter of less than one hundred microns. Experiments were conducted in smooth microchannels of hydraulic diameter of 45, 65 μm, and a rough microchannel of hydraulic diameter of 70 μm, with deionised water as the working fluid. The local saturation pressure and temperature vary substantially over the length of the channel. In order to correctly predict the local saturation temperature and subsequently the heat transfer characteristics, numerical techniques have been used in conjunction with the conventional two phase pressure drop models. The Lockhart–Martinelli (liquid–laminar, vapour–laminar) model is found to predict the two phase pressure drop data within 20%. The instability in two phase flow is quantified; it is found that microchannels of smaller hydraulic diameter have lesser instabilities as compared to their larger counterparts. The experiments also suggest that surface characteristics strongly affect flow stability in the two phase flow regime. The effect of hydraulic diameter and surface characteristics on the flow characteristics and stability in two phase flow is seldom reported, and is of considerable practical relevance.     

Flow of mono-dispersed particles through horizontal bend

Pipeline slurry flow of mono-dispersed particles through horizontal bend is numerically simulated by implementing Eulerian two-phase model in FLUENT software. A hexagonal shape and Cooper type non-uniform three-dimensional grid is chosen to discretize the entire computational domain, and a control volume finite difference method is used to solve the governing equations. The modeling results are compared with the experimental data collected in 53.0 mm diameter horizontal bend with radius of 148.4 mm for concentration profiles and pressure drops. Experiments are performed on narrow-sized silica sand with mean diameter of 450 μm and for flow velocity up to 3.56 m/s (namely, 1.78, 2.67 and 3.56 m/s) and four efflux concentrations up to 16.28% (namely, 0%, 3.94%, 8.82% and 16.28%) by volume for each velocity. Eulerian model gives fairly accurate predictions for both the pressure drop and concentration profiles at all efflux concentrations and flow velocities.     

Flow pattern based correlations of two-phase pressure drop in rectangular microchannels

Numerous pressure drop correlations for microchannels have been proposed; most of them can be classified as either a homogeneous flow model (HFM) or a separated flow model (SFM). However, the predictions of these correlations have not been compared directly because they were developed in experiments conducted under a range of conditions, including channel shape, the number of channels, channel material and the working fluid. In this study, single rectangular microchannels with different aspect ratios and hydraulic diameters were fabricated in a photosensitive glass. Adiabatic water-liquid and Nitrogen-gas two-phase flow experiments were conducted using liquid superficial velocities of 0.06–1.0 m/s, gas superficial velocities of 0.06–72 m/s and hydraulic diameters of 141, 143, 304, 322 and 490 μm. A pressure drop in microchannels was directly measured through embedded ports. The flow pattern was visualized using a high-speed camera and a long-distance microscope. A two-phase pressure drop in the microchannel was highly related to the flow pattern. Data were used to assess seven different HFM viscosity models and ten SFM correlations, and new correlations based on flow patterns were proposed for both HFMs and SFMs.     

Numerical study of oscillating boundary layer flow over a flat plate using k–kL–ω turbulence model

Oscillating boundary layer flow over an infinite flat plate at rest was simulated using the k–k_L–ω turbulence model for a Reynolds number range of 32 ⩽ Re_δ ⩽ 10,000 ranging from fully laminar flow to fully turbulent flow. The k–k_L–ω model was validated by comparing the predictions with LES results and experimental results for intermittently turbulent and fully turbulent flow regimes. The good agreement obtained between the k–k_L–ω model prediction with the experimental and LES results indicate that the k–k_L–ω model is able to accurately simulate transient intermittently turbulent flow and as well as accurately predict the onset of turbulence for such oscillatory flows.     

Void fraction and pressure drop during external upward two-phase crossflow in tube bundles – part I: Experimental investigation

The present paper is the part I of a broad study concerning void fraction and pressure drop for air-water upward external flow across tube bundles. Experimental results were obtained for liquid and gas superficial velocities ranging from 0.02 to 1.50 m/s and 0.20 to 10.00 m/s, respectively. Void fraction measurements were performed for bubbly flow using a capacitive probe. The test section consisted of a triangular tube bundle counting with 19 mm OD tube and transverse pitch of 24 mm. Initially, the paper describes the test facility and the data regression and experimental procedures. Then, the pressure drop and void fraction measurements are validated based on tests for single-phase flow and quiescent liquid conditions, respectively. Finally, the experimental data are presented and analyzed. In the second part of this study (Part II), a literature review on predictive methods for void fraction and pressure drop is presented. Additionally, these methods are compared with the database presented in Part I and new predictive methods for void fraction and frictional pressure drop are proposed.     

Hydrodynamics and energy consumption studies in a three-phase liquid circulating three-phase fluid bed contactor

The hydrodynamics and energy consumption have been studied in a cold flow, bubbling and turbulent, pressurized gas–liquid–solid three-phase fluidized bed (0.15 m ID × 1 m height) with concurrent gas–liquid up flow is proposed with the intention of increasing the gas hold up. The hydrodynamic behaviour is described and characterised by some specific gas and liquid velocities. Particles are easily fluidized and can be uniformly distributed over the whole height of the column. The effect of parameters like liquid flow rate, gas flow rate, particle loading, particle size, and solid density on gas hold up and effect of gas flow rate, solid density and particle size on solid hold up, energy consumption and minimum fluidization velocity has been studied. At the elevated pressures a superior method for better prediction of minimum fluidization velocity and terminal settling velocities has been adopted. The results have been interpreted with Bernoulli’s theorem and Richardson–Zaki equation. Based on the assumption of the gas and liquid as a pretend fluid, a simplification has been made to predict the particle terminal settling velocities. The Richardson–Zaki parameter n′ was compared with Renzo’s results. A correlation has been proposed with the experimental results for the three-phase fluidization.     

Multi-scale particle dynamics of low air velocity in a horizontal self-excited gas–solid two-phase pipe flow

The particle fluctuation velocities of a horizontal self-excited gas–solid two-phase pipe flow with soft fins near MPD (minimum pressure drop) air velocity are first measured by high-speed PIV in the acceleration and fully-developed regimes. Then orthogonal wavelet multi-resolution analysis and power spectrum are used to reveal multi-scale characteristics of particle fluctuation velocity. It is observed that the pronounced peaks of the spectra of axial and vertical fluctuation velocities appear in the range of low frequency near the bottom of pipe. These peaks of spectra become larger and their frequencies decrease by using fins. In the range of low frequencies (3–25 Hz), the wavelet components of the fluctuating energy of axial particle velocity make the main contribution accounting for 87% and 93% respectively for non-fin and using fins near the bottom of pipe. In the range of relatively high frequency (50–400 Hz), however, the wavelet components of using fins, accounting for about 49%, become smaller than that of non-fin, accounting for about 72%, in the suspension flow regime near the top of pipe. The skewness factor of axial particle fluctuation velocity indicates that the wavelet components follow the Gaussian probability distribution as the central frequency decreases.     

Saturated critical heat flux in a multi-microchannel heat sink fed by a split flow system

An extensive experimental campaign has been carried out for the measurement of saturated critical heat flux in a multi-microchannel copper heat sink. The heat sink was formed by 29 parallel channels that were 199 μm wide and 756 μm deep. In order to increase the critical heat flux and reduce the two-phase pressure drop, a split flow system was implemented with one central inlet at the middle of the channels and two outlets at either end. The base critical heat flux was measured using three HFC Refrigerants (R134a, R236fa and R245fa) for mass fluxes ranging from 250 to 1500 kg/m² s, inlet subcoolings from ?25 to ?5 K and saturation temperatures from 20 to 50 °C. The parametric effects of mass velocity, saturation temperature and inlet subcooling were investigated. The analysis showed that significantly higher CHF was obtainable with the split flow system (one inlet–two outlets) compared to the single inlet–single outlet system, providing also a much lower pressure drop. Notably several existing predictive methods matched the experimental data quite well and quantitatively predicted the benefit of higher CHF of the split flow.     

Experimental study of the effect of friction phenomena on actual and calculated inventory in a small-scale CFB riser

Accurate information concerning riser inventory in a fluidized bed is required in some applications such as the calcium looping process, because it is related to the CO₂ capture efficiency of the system. In a circulating fluidized bed (CFB), the riser inventory is normally calculated from the riser pressure drop; however, the friction and the acceleration phenomena may have a significant influence on the total riser pressure drop. Therefore, deviation may occur in the calculation from the actual mass. For this reason the magnitude of the friction and the acceleration pressure drop in the entire riser is studied in small-scale risers. Two series of studies were performed: the first one in a scaled cold model riser of the 10 kW_th facility, and the second one in the 10 kW_th fluidized bed riser under process conditions. The velocities were chosen to comply with the fluidization regimes suitable for the calcium looping process, namely, the turbulent and the fast. In cold-model experiments in a low-velocity turbulent fluidization regime, the actual weight (static pressure drop) of the particles is observed more than the weight calculated from a recorded pressure drop. This phenomenon is also repeated in pilot plant conditions. In the cold-model setup, the friction and acceleration pressure drop became apparent in the fast fluidization regime, and increased as the gas velocity rose. Within calcium looping conditions in the pilot plant operation, the static pressure drop was observed more than the recorded pressure drop. Therefore, as a conservative approach, the influence of friction pressure drop may be neglected while calculating the solid inventory of the riser. The concept of transit inventory is introduced as a fraction of total inventory, which lies in freefall zones of the CFB system. This fraction increases as gas velocity rises.     

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.     

Numerical evaluation of turbulence models for dense to dilute gas–solid flows in vertical conveyor

A two-fluid model (TFM) of multiphase flows based on the kinetic theory and small frictional limit boundary condition of granular flow was used to study the behavior of dense to dilute gas–solid flows in vertical pneumatic conveyor. An axisymmetric 2-dimensional, vertical pipe with 5.6 m length and 0.01 m internal diameter was chosen as the computation domain, same to that used for experimentation in the literature. The chosen particles are spherical, of diameter 1.91 mm and density 2500 kg/m³. Turbulence interaction between the gas and particle phases was investigated by Simonin's and Ahmadi's models and their numerical results were validated for dilute to dense conveying of particles. Flow regimes transition and pressure drop were predicted. Voidage and velocity profiles of each phase were calculated in radial direction at different lengths of the conveying pipe. It was found that the voidage has a minimum, and gas and solid velocities have maximum values along the center line of the conveying pipe and pressure drop has a minimum value in transition from dense slugging to dilute stable flow regime. Slug length and pressure fluctuation reduction were predicted with increasing gas velocity, too. It is shown that solid phase turbulence plays a significant role in numerical prediction of hydrodynamics of conveyor and the capability of particles turbulence models depends on tuning parameters of slip-wall boundary condition.     

An experimental investigation of two-phase air–water flow through a horizontal circular micro-channel

This paper is a continuation of the authors’ previous work. Two-phase air–water flow experiments are performed in a horizontal circular micro-channel. The test section is made of a fused silica tube with an inner diameter of 0.15 mm and a length of 104 mm. The flow phenomena, which are liquid/unstable annular alternating flow (LUAAF), liquid/annular alternating flow (LAAF), and annular flow, are observed and recorded by a high-speed camera mounted together with a stereozoom microscope. A flow pattern map is presented in terms of the phase superficial velocities and is compared with those of other researchers obtained from different working fluids. Image analysis is performed to determine the void fraction, which increases non-linearly with increasing volumetric quality. It is revealed that the two-phase frictional multiplier data show a dependence on flow pattern rather than mass flux. Based on the present data, a new pressure drop correlation is proposed for practical applications. According to the present study, in general the data for the two-phase air–water flow characteristics are found to comply with those of working fluids other than air–water mixtures.     

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.     

CFD–DEM simulation of the pneumatic conveying of fine particles through a horizontal slit

Fine particles play a significant role in many industrial processes. To study the dynamic behavior of fine particle and their deposition in rock fractures, the pneumatic conveying of fine particles (approximately 100 μm in diameter) through a small-scale horizontal slit (0.41 m × 0.025 m) was studied, which is useful for the sealing technology of underground gas drainage in coal mining production. The CFD–DEM method was adopted to model the gas-particle two-phase flow; the gas phase was treated as a continuum and modeled using computational fluid dynamics (CFD), particle motion and collisions were simulated using the DEM code. Then, the bulk movement of fine particles through a small-scale horizontal slit was explored numerically, and the flow patterns were further investigated by visual inspection. The simulation results indicated that stratified flow or dune flow can be observed at low gas velocities. For intermediate gas velocities, the flow patterns showed pulsation phenomena, and dune flow reappeared in the tail section. Moreover, periodic flow regimes with alternating thick and sparse stream structures were observed at a high gas velocity. The simulation results of the bulk movement of fine particles were in good agreement with the experimental findings, which were obtained by video-imaging experiments. Furthermore, the calculated pressure drop versus gas velocity profile was investigated and compared with relative experimental findings, and the results showed good agreement. Furthermore, the particle velocity vectors and voidage distribution were numerically simulated. Selected stimulation results are presented and provide a reference for the further study of fine particles.     

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).     

Experimental study on the effect of drag reducing polymer on flow patterns and drag reduction in a horizontal oil–water flow

In this study, a HMW anionic co-polymer of 40:60 wt/wt NaAMPS/acrylamide was used as a drag reducing polymer (DRP) for oil–water flow in a horizontal 25.4 mm ID acrylic pipe. The effect of polymer concentration in the master solution and after injection in the main water stream, oil and water velocities, and pipe length on drag reduction (DR) was investigated. The injected polymer had a noticeable effect on flow patterns and their transitions. Stratified and dual continuous flows extended to higher superficial oil velocities while annular flow changed to dual continuous flow. The results showed that as low as 2 ppm polymer concentration was sufficient to create a significant drag reduction across the pipe. DR was found to increase with polymer concentration increased and reached maximum plateau value at around 10 ppm. The results showed that the drag reduction effect tends to increase as superficial water velocity increased and eventually reached a plateau at U_sw of around 1.3 m/s. At U_sw > 1.0 m/s, the drag reduction decreased as U_so increased while at lower water velocities, drag reduction is fluctuating with respect to U_so. A maximum DR of about 60% was achieved at U_so = 0.14 m/s while only 45% was obtained at U_so = 0.52 m/s. The effectiveness of the DRP was found to be independent of the polymer concentration in the master solution and to some extent pipe length. The friction factor correlation proposed by Al-Sarkhi et al. (2011) for horizontal flow of oil–water using DRPs was found to underpredict the present experimental pressure gradient data.     

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.     

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.     

Gas dispersion in horizontal pulp-fibre-suspension flow

Dispersion of gas into pulp-suspension horizontal flow was investigated downstream of 90° tees for ranges of fibre mass concentrations (0–3.0%), superficial liquid/pulp velocities (0.5–5.0 m/s) and superficial gas velocities (0.11–0.44 m/s) based on a gas mixing index, derived from the standard deviation of cross-sectional local gas holdup obtained from electrical resistance tomography. Mixing for dilute suspensions was similar to that for water, but differed significantly for higher suspension concentrations. Mixing worsened with increasing fibre mass concentration for the bubble flow regime, likely due to dense fibre networks in the core of the pipe causing bubbles to congregate near the wall. When buoyancy was significant, gas uniformity improved with increasing pulp concentration, since robust fibre networks caused liquid/pulp slugs to flow at the top of the pipe, whereas stratified flow was approached at lower concentrations. Mixing was less dependent on superficial liquid/pulp velocity at higher pulp concentrations, due to less variation in flow regimes.     

Influences of peripherally-cut twisted tape insert on heat transfer and thermal performance characteristics in laminar and turbulent tube flows

Effects of peripherally-cut twisted tape insert on heat transfer, friction loss and thermal performance factor characteristics in a round tube were investigated. Nine different peripherally-cut twisted tapes with constant twist ratio (y/W = 3.0) and different three tape depth ratios (DR = d/W = 0.11, 0.22 and 0.33), each with three different tape width ratios (WR = w/W = 0.11, 0.22 and 0.33) were tested. Besides, one typical twisted tape was also tested for comparison. The measurement of heat transfer rate was conducted under uniform heat flux condition while that of friction factor was performed under isothermal condition. Tests were performed with Reynolds number in a range from 1000 to 20,000, using water as a working fluid. The experimental results revealed that both heat transfer rate and friction factor in the tube equipped with the peripherally-cut twisted tapes were significantly higher than those in the tube fitted with the typical twisted tape and plain tube, especially in the laminar flow regime. The higher turbulence intensity of fluid in the vicinity of the tube wall generated by the peripherally-cut twisted tape compared to that induced by the typical twisted tape is referred as the main reason for achieved results. The obtained results also demonstrated that as the depth ratio increased and width ratio decreased, the heat transfer enhancement increased. Over the range investigated, the peripherally-cut twisted tape enhanced heat transfer rates in term of Nusselt numbers up to 2.6 times (turbulent regime) and 12.8 times (laminar regime) of that in the plain tube. These corresponded to the maximum performance factors of 1.29 (turbulent regime) and 4.88 (laminar regime).