Motion of singles bubbles moving under a slightly inclined surface through stationary liquids

The influence of the liquid properties on the dynamical bubble shape and on the bubble motion has been investigated for bubbles moving under a downward facing inclined surface. The Morton number Mo varied from 2.59 × 10⁻¹¹ to 2.52 × 10⁺⁰¹. The Bond number Bo covered the range from 10 to 150 and the surface inclination angle θ was varied from 2° to 6°. To cover the wide range of Mo, several liquids such as glycerine, propanediol, water and isopropanol were used. The results have shown that the relation Fr = Fr(Bo, Mo, θ) is not adequate to describe the bubble motion, where Fr is the terminal Froude number. The choice of the terminal Reynolds number Re as the dependent parameter, allowed the clarification of the role of the Morton number on the bubble motion. At a given Bond number, the bubble Reynolds number decreases monotonously with the Morton number. Furthermore, an empirical correlation Re = Re(Bo, Mo, θ) is given that can be readily used in the mathematical modelling of bubble laden flows under solids.     

Experimental investigation of the onsets of gas and liquid entrainment from a small branch mounted on an inclined wall

The phenomena of the onsets of liquid entrainment and gas entrainment were investigated experimentally for the case of a flat plane with a circular outlet branch of diameter d (=6.35 mm) at the plane centre. This flat plane was situated in a large tank containing a stratified mixture of air and water under pressure (317 kPa for most experiments and 520 kPa for a few experiments) and at room temperature. The plane was inclined through various angles (θ) in increments of 30°, from the outlet branch orientation being vertically upward through the horizontal to vertically downward. For both onsets the vertical distance between the centre of the outlet branch and the undisturbed gas–liquid interface (h) was measured for various angles of inclination and Froude numbers. Both onsets were observed visually through a large viewing part of the test section. It was found that for both onsets there is a range of inclination angles where the onset h depended on θ and a range where the onset h essentially did not depend on θ. The data were correlated in terms of onset h/d, Froude number, and θ where there was dependence of onset h on the angle of inclination.     

Terminal velocities of clean and fully-contaminated drops in vertical pipes

Terminal velocities and shapes of drops rising through vertical pipes in clean and fully-contaminated systems are measured by using a high-speed video camera and an image processing method. Silicon oils and glycerol water solutions are used for the dispersed and continuous phases, respectively. Triton X-100 is used for surfactant. Clean and contaminated drops take either spherical, spheroidal or deformed spheroidal shapes when the diameter ratio λ is less than a critical value, λ_C, whereas they take bullet shapes for λ > λ_C (Taylor drops). The applicability of available drag and Froude number correlations is examined through comparisons with the measured data. Effects of surfactant on the shape and terminal velocity of a Taylor drop are also discussed based on the experimental data and interface tracking simulations. The conclusions obtained are as follows: (1) drag and Froude number correlations proposed so far give reasonable estimations of the terminal velocities of clean drops at any λ, (2) the terminal velocities of contaminated drops are well evaluated by making the viscosity ratio μ^* infinity in the drag correlation for clean drops in the viscous force dominant regime, (3) the effects of surfactant on the shape and terminal velocity of a Taylor drop become significant as the Eötvös number, Eo_D, decreases and μ^* increases, and (4) the reduction in surface tension due to the addition of surfactant would be the cause of the increase in the terminal velocity and elongation of a contaminated Taylor drop.     

Evaporation heat transfer of R-134a inside a microfin tube with different tube inclinations

An experimental investigation has been carried out to study the heat transfer characteristics during evaporation of R-134a inside a single helical microfin tube. The microfin tube has been provided with different tube inclination angles of the direction of fluid flow from horizontal, α. The experiments were performed for seven different tube inclinations, α, in a range of −90° to +90° and four mass velocities of 53, 80, 107 and 136 kg/m² s for each tube inclination angle during evaporation of R-134a. The results demonstrate that the tube inclination angle, α, affects the boiling heat-transfer coefficient in a significant manner. For all refrigerant mass velocities, the best performing tube is that having inclination angle of α = +90°. The effect of tube inclination angle, α, on heat-transfer coefficient, h, is more prominent at low vapor quality and mass velocity. An empirical correlation has also been developed to predict the heat-transfer coefficient during flow boiling inside a microfin tube with different tube inclinations.     

Upward and downward inclination oil–water flows

The effect of upward (+5°, +10°) and downward (−5°) pipe inclinations on the flow patterns, hold up and pressure gradient during two-liquid phase flows was investigated experimentally for mixture velocities between 0.7 and 2.5 m/s and phase fractions between 10% and 90%. The investigations were performed in a 38 mm ID stainless steel test pipe with water and oil as test fluids. High-speed video recording and local impedance and conductivity probes were used to precisely identify the different flow patterns. In both positive and negative inclinations the dispersed oil-in-water regime extended to lower mixture velocities and higher oil fractions compared to horizontal flow. A new flow pattern, oil plug flow, appeared at both +5° and +10° inclination while the stratified wavy pattern disappeared at −5° inclination. The oil to water velocity ratio was higher for the upward than for the downward flows but in the majority of cases and all inclinations oil was flowing faster than water. At low mixture velocities the velocity ratio increased with oil fraction while it decreased at high velocities. The increase became more significant as the degree of inclination increased. The frictional pressure gradient in both upward and downward flows was in general lower than in horizontal flows while a minimum occurred at all inclinations at high mixture velocities during the transition from dispersed water-in-oil to dual continuous flow.     

Natural convection in an inclined quadrantal cavity heated and cooled on adjacent walls

Natural convective flow and heat transfer in an inclined quadrantal cavity is studied experimentally and numerically. The particle tracing method is used to visualize the fluid motion in the enclosure. Numerical solutions are obtained via a commercial CFD package, Fluent. The working fluid is distilled water. The effects of the inclination angle, ? and the Rayleigh number, Ra on fluid flow and heat transfer are investigated for the range of angle of inclination between 0° ? ? ? 360°, and Ra from 10⁵ to 10⁷. It is disclosed that heat transfer changes dramatically according to the inclination angle which affects convection currents inside, i.e. flow physics inside. A fairly good agreement is observed between the experimental and numerical results.     

Experimental investigation of two-phase discharge from a stratified region through a small branch mounted on an inclined wall

Experimental data are presented for the mass flow rate and quality of two-phase discharge through a small branch of diameter d (=6.35 mm) attached normally to an inclined flat plane. The flat plane was situated in a large tank containing a stratified mixture of air and water under pressure (316 kPa) and at room temperature. The plane was inclined through various angles (θ) in increments of 30°, from the outlet-branch orientation being vertically upward through the horizontal to vertically downward. The bulk of the data correspond to seven inclination angles and two test-section-to-separator pressure differences (ΔP) of 11.0 and 115.5 kPa, and for each combination of θ and ΔP, the mass flow rate and quality were measured at different values of the interface level (h) between the onsets of gas and liquid entrainment. Four additional data sets were generated for other values of ΔP in order to confirm certain trends. Influences of these independent variables on the mass flow rate and quality are discussed and normalized plots are presented showing that the data can be collapsed for a wide range of conditions. Comparisons are made between the present data and previous correlations/models and new empirical correlations are formulated and shown to be capable of predicting the present data with good accuracy.     

Reynolds number, thickness and camber effects on flapping airfoil propulsion

The effect of varying airfoil thickness and camber on plunging and combined pitching and plunging airfoil propulsion at Reynolds number Re=200, 2000, 20 000 and 2×10⁶ was studied by numerical simulations for fully laminar and fully turbulent flow regimes. The thickness study was performed on 2-D NACA symmetric airfoils with 6-50% thick sections undergoing pure plunging motion at reduced frequency k=2 and amplitudes h=0.25 and 0.5, and for combined pitching and plunging motion at k=2, h=0.5, phase ?=90°, pitch angle θ_o=15° and 30° and the pitch axis was located at 1/3 of chord from leading edge. At Re=200 for motions where positive thrust is generated, thin airfoils outperform thick airfoils. At higher Re significant gains could be achieved both in thrust generation and propulsive efficiency by using a thicker airfoil section for plunging and combined motion with low pitch amplitude. The camber study was performed on 2-D NACA airfoils with varying camber locations undergoing pure plunging motion at k=2, h=0.5 and Re=20 000. Little variation in thrust performance was found with camber. The underlying physics behind the alteration in propulsive performance between low and high Reynolds numbers has been explored by comparing viscous Navier-Stokes and inviscid panel method results. The role of leading edge vortices was found to be key to the observed performance variation.     

Terminal velocity of a Taylor drop in a vertical pipe

A scaling analysis based on the field equations for two phases and the jump conditions at the interface is carried out to deduce a balance of forces acting on a Taylor drop rising through stagnant liquid in a vertical pipe. The force balance is utilized to deduce a functional form of an empirical correlation of terminal velocity of the Taylor drop. Undetermined coefficients in the correlation are evaluated by making use of available correlations for two limiting cases, i.e. extremely high and low Reynolds number Taylor bubbles in large pipes. Terminal velocity data obtained by interface tracking simulations are also used to determine the coefficients. The proposed correlation expresses the Froude number Fr as a function of the drop Reynolds number Re_D, the Eötvös number Eo_D and the viscosity ratio μ^*. Comparisons between the correlation, simulations and experimental data confirm that the proposed correlation is applicable to Taylor drops under various conditions, i.e., 0.002 < Re_D < 4960, 4.8 < Eo_D < 228, 0 ? μ^* ? 70, 1 < N < 14700, −12 < log M < 4, and d/D < 1.6, where N is the inverse viscosity number, M the Morton number, d the sphere-volume equivalent drop diameter and D the pipe diameter.     

The effect of the VOF–CSF static contact angle boundary condition on the dynamics of sliding and bouncing ellipsoidal bubbles

The static contact angle is the only empiricism introduced in a Volume of Fluid–Continuum Surface Force (VOF–CSF) model of bubbly flow. Although it has previously been shown to have a relatively limited effect on the accuracy of velocity and shape predictions in the case of large gas bubbles sliding under inclined walls (e.g. Cook and Behnia, 2001), it may have a more determining influence on the numerical prediction of the dynamics of smaller ellipsoidal bubbles which were shown by Tsao and Koch (1997) to bounce repeatedly when sliding under inclined walls at certain wall inclinations. The present paper reports on the influence of surface tension and the static contact angle on the dynamics of an ellipsoidal air bubble of equivalent diameter D_e = 3.4 mm. The bubble Eötvös and Morton numbers are Eo = 1.56 and Mo = 2 × 10⁻¹¹ respectively. The computational results are achieved with a Piecewise Linear Construction (PLIC) of the interface and are reviewed with reference to experimental measurements of bubble velocity and interface shape oscillations recorded using a high speed digital camera. Tests are performed at plate inclination angles θ ∈ {10°, 20°, 30°, 45°} to the horizontal and computational models consider three static contact angles θ_c ∈ {10°, 20°, 30°}. The static contact angle has been found to have a significant effect on the bubble dynamics but to varying degree depending on the plate inclination. It is shown to promote lift off and bouncing when the plate inclination angle reaches 30° in a way that does not necessarily reflect experimental observations.     

Single-beam gamma densitometry measurements of oil–water flow in horizontal and slightly inclined pipes

Gamma densitometry is a frequently used non-intrusive method for measuring component volume fractions in multiphase flow systems. The application of a single-beam gamma densitometer to investigate oil–water flow in horizontal and slightly inclined pipes is presented. The experiments are performed in a 15 m long, 56 mm diameter, inclinable stainless steel pipe using Exxsol D60 oil (viscosity 1.64 mPa s, density 790 kg/m³) and water (viscosity 1.0 mPa s, density 996 kg/m³) as test fluids. The test pipe inclination is changed in the range from 5° upward to 5° downward. Experimental measurements are reported at three different mixture velocities, 0.25, 0.50 and 1.00 m/s, and the inlet water cut is varied from 0 to 1. The gamma densitometer is composed of radioactive isotope of Am-241 with the emission energy of 59.5 keV, scintillation detector [NaI(Tl)] and signal processing system. The time averaged cross-sectional distributions of oil and water phases are measured by traversing the gamma densitometer along the vertical pipe diameter. Based on water volume fraction measurements, water hold-up and slip ratio are estimated. The total pressure drop over the test section is measured and frictional pressure drop is estimated based on water hold-up measurements. The measurement uncertainties associated with gamma densitometry are also discussed. The measured water hold-up and slip ratio profiles are strongly dependent on pipe inclination. In general, higher water hold-up values are observed in upwardly inclined pipes compared to the horizontal and downwardly inclined pipes. At low mixture velocities, the slip ratio decreases as the water cut increases. The decrease is more significant as the degree of inclination increases. The frictional pressure drop for upward flow is slightly higher than the horizontal flow. In general, there is a marginal difference in frictional pressure drop values for horizontal and downwardly inclined flows.     

Bubble migration in a turbulent boundary layer

The wall void peaking distribution observed in an upward turbulent bubbly boundary layer along a flat plate is generated by bubbles that move towards the plate, come into contact with the wall and then slide along it. This transverse ‘migration’ has been studied using flow visualization, high speed video and particle tracking techniques to measure the trajectories of mono-disperse air bubbles at very low void fractions. Investigations have been performed at four Reynolds numbers in the range 280 < Re_θ < 3000, covering both the laminar and turbulent regimes, with mono-disperse bubbles of mean equivalent diameter between 2 mm and 6 mm. Lagrangian statistics calculated from hundreds of trajectories show that the migration only occurs in the turbulent regime and for bubble diameters below some critical value: 3.5 mm < d_eqcrit < 4 mm. Above this size (We > 3), the interface deformation is such that bubbles do not remain at the wall, even when they are released at the surface. Also, bubble migration is shown to be non-systematic, to have a non-deterministic character in the sense that trajectories differ significantly, to increase with Reynolds number and to take place on a short time scale. A series of experiments with isolated bubbles demonstrates that these results are not influenced by bubble–bubble interactions and confirm that two-way coupling in the flow is limited. Flow visualizations show that the migration originates with the capture of bubbles inside the large turbulent structures of the boundary layer (‘bulges’). The bubbles begin to move towards the wall as they cross these structures, and the point at which they reach the wall is strongly correlated with the position of the deep ‘valleys’ which separate the turbulent ‘bulges’. The analysis of the mean Lagrangian trajectories of migrating bubbles confirms these observations. Firstly, the average time of migration calculated from these trajectories coincides with the mean transit time of the bubbles across the structures. Secondly, once the trajectories have been scaled by this transit time and the boundary layer thickness δ, they all have the same form in the region y/δ < 0.4, independent of the Reynolds number.     

Viscous oil–water flows in a microchannel initially saturated with oil: Flow patterns and pressure drop characteristics

Immiscible viscous liquid–liquid two-phase flow patterns and pressure drop characteristics in a circular microchannel have been investigated. Water and silicone oil with a dynamic viscosity of 863 mPa s were injected into a fused silica microchannel with an inner diameter of 250 μm. As the microchannel was initially filled with the silicone oil, an oil film was found to always form and remain on the microchannel wall. Different flow patterns were observed and classified over a wide range of water and oil flow rates. A flow pattern map is presented in terms of Re, Ca, and We numbers. Two-phase pressure drop data have also been collected and analyzed to develop a simple correlation for slug, annular and annular-droplet flow patterns in terms of superficial water and oil velocities.     

Exact solutions for the flow of a generalized Oldroyd-B fluid induced by a constantly accelerating plate between two side walls perpendicular to the plate

The unsteady flow of an incompressible generalized Oldroyd-B fluid induced by a constantly accelerating plate between two side walls perpendicular to the plate has been studied using Fourier sine and Laplace transforms. The obtained solutions for the velocity field and shear stresses, written in terms of the generalized G and R functions, are presented as sum of the similar Newtonian solutions and the corresponding non-Newtonian contributions. For α = β = 1 and λ_r → λ these solutions are going to the corresponding Newtonian solutions. Furthermore, the solutions for generalized Maxwell fluids as well as those for ordinary Oldroyd-B and Maxwell fluids, performing the same motion, are also obtained as limiting cases of our general solutions. In the absence of the side walls, namely when the distance between the two walls tends to infinity, the solutions corresponding to the motion over an infinite constantly accelerating plate are recovered. For λ_r → 0 and β → 1, these solutions reduce to the known solutions from the literature. Finally, the effect of the material parameters on the velocity profile is spotlighted by means of the graphical illustrations.     

Hydrodynamics and boiling phenomena of water droplets impinging on hot solid

The collision of single water droplets with a hot Inconel 625 alloy surface was investigated by a two-directional flash photography technique using two digital still cameras and three flash units. The experiments were conducted under the following conditions: the pre-impact diameters of the droplets ranged from 0.53 to 0.60 mm, the impact velocities ranged from 1.7 m/s to 4.1 m/s, and the solid surface temperatures ranged from 170 °C to 500 °C. When a droplet impacted onto the solid at a temperature of 170 °C, weak boiling was observed at the liquid/solid interface. At temperatures of 200 or 300 °C, numerous vapor bubbles were formed. Numerous secondary droplets then jetted upward from the deforming droplet due to the blowout of the vapor bubbles into the atmosphere. No secondary droplets were observed for a surface temperature of 500 °C at the low-impact Weber numbers (∼30) associated with the impact inertia of the droplets. Experiments using 2.5-mm-diameter droplets were also conducted. The dimensionless collision behaviors of large and small droplets were compared under the same Weber number conditions. At temperatures of less than or equal to 300 °C, the blowout of vapor bubbles occurred at early stages for a large droplet. At a surface temperature of 500 °C, the two dimensionless deformation behaviors of the droplets were very similar to each other.     

Low-Reynolds-number motion of a deformable drop between two parallel plane walls

The motion of a three-dimensional deformable drop between two parallel plane walls in a low-Reynolds-number Poiseuille flow is examined using a boundary-integral algorithm that employs the Green’s function for the domain between two infinite plane walls, which incorporates the wall effects without discretization of the walls. We have developed an economical calculation scheme that allows long-time dynamical simulations, so that both transient and steady-state shapes and velocities are obtained. Results are presented for neutrally buoyant drops having various viscosity, size, deformability, and channel position. For nearly spherical drops, the decrease in translational velocity relative to the undisturbed fluid velocity at the drop center increases with drop size, proximity of the drop to one or both walls, and drop-to-medium viscosity ratio. When deformable drops are initially placed off the centerline of flow, lateral migration towards the channel center is observed, where the drops obtain steady shapes and translational velocities for subcritical capillary numbers. With increasing capillary number, the drops become more deformed and have larger steady velocities due to larger drop-to-wall clearances. Non-monotonic behavior for the lateral migration velocities with increasing viscosity ratio is observed. Simulation results for large drops with non-deformed spherical diameters exceeding the channel height are also presented.     

Two-phase flow pressure change subject to sudden contraction in small rectangular channels

This study investigates the pressure change and flow pattern subject to the influence of sudden contractions. The air and water mixture flows from small rectangular channels (2 × 4, 2 × 6, 4 × 4 and 4 × 6 mm, respectively) into a 2 mm diameter tube. The total mass flux (G) ranges from 100 to 700 kg/m² s with gas quality (x) being varied from 0.001 to 0.8. In general the contraction pressure change increases with the rise of mass flux, and gas quality but an unique deflection of contraction pressure change pertaining to liquid vena contracta at a very low gas quality is encountered at a very low gas quality in the 4 ± 6 mm test section. For a low gas quality, elongated bubble prevails after the contraction, yet the size of elongated bubbles is reduced when the aspect ratio is increased. Comparisons amid the pressure change data of this study and available literature data with the predictions of existing model/correlations indicate that none of them can accurately predict the data. It is found that the influence of surface tension and outlet tube size, or equivalently the Bond number plays a major role for the departure of various models/correlations. Among the models/correlations being examined, the homogeneous model shows a little better than the others. Hence by taking account the influences of gas quality, Bond number, Weber number and area contraction ratio into the homogeneous model, a modified homogeneous correlation is proposed that considerably improves the predictive ability over existing correlations with a mean deviation of 30% to all the data.     

Fluid forces on a very low Reynolds number airfoil and their prediction

This paper presents the measurements of mean and fluctuating forces on an NACA0012 airfoil over a large range of angle (α) of attack (0-90°) and low to small chord Reynolds numbers (Re_c), 5.3 × 10³-5.1 × 10⁴, which is of both fundamental and practical importance. The forces, measured using a load cell, display good agreement with the estimate from the LDA-measured cross-flow distributions of velocities in the wake based on the momentum conservation. The dependence of the forces on both α and Re_c is determined and discussed in detail. It has been found that the stall of an airfoil, characterized by a drop in the lift force and a jump in the drag force, occurs at Re_c ? 1.05 × 10⁴ but is absent at Re_c = 5.3 × 10³. A theoretical analysis is developed to predict and explain the observed dependence of the mean lift and drag on α.     

Three-dimensional numerical simulation of sedimenting drops inside a vertical channel

Numerical simulation of sedimenting deformable drops inside a vertical channel has been performed at finite Reynolds numbers. The channel is confined by two vertical walls in x-direction and is periodic in y- and z-directions. Results are obtained using a finite difference/front-tracking method. The main dimensionless parameters are the Reynolds number, the Bond number and ratio of the length of channel to the diameter of drops. The effect of these parameters on lateral migration of drops is investigated. It is found that the wall repulsion is the main mechanism of the lateral migration of the drop, and drop migrates toward the channel axis. When the Reynolds number is relatively low, two different lateral migration regimes are observed: migration with monotonic approach and migration with damped oscillations. These regimes are affected by the dimensionless parameters. When the Bond number increases, the oscillations of drop around the centerline of channel are stronger and drop reaches the channel centerline in a larger period. Results of lateral migration of one drop are consistent with perturbation theory, and two-dimensional numerical simulations performed by Feng et al. (1994). The drag coefficient has also been calculated, and effect of various parameters has been discussed. Two drops interaction is similar to that observed by Feng et al. (1994) for two-dimensional circular cylinders. Results are consistent with experiments performed by Wu and Manasseh (1998). Simulations of four sedimenting drops show that depending on the relative size of drops, they either fall in two rows or they form a single horizontal layer and settle with a unique velocity.