Laminar bubbly flow in an open capillary channel in microgravity

The study of a bubbly laminar two-phase flow in an open capillary channel under microgravity conditions was conducted aboard the sounding rocket, Texus-45. The channel consists of two parallel plates of width b = 25 mm and distance a = 10 mm. The flow along the length l = 80 mm is confined by a free surface on one side and a plate on the opposite side. The bubbles are injected at the nozzle of the capillary channel via six capillary tubes of 100 μm in inner diameter. Different liquid and gas flow rates were tested leading to different liquid free surface shape and bubble size.     

An inclined wall jet: Mean flow characteristics and effects of acoustic excitation

 The mean velocity field of a 30° inclined wall jet has been investigated using both hot-wire and laser Doppler anemometry (LDA). Provided that the nozzle aspect ratio is greater than 30 and the inclined wall angle (β) is less than 50°, LDA measurements for various β show that the reattachment length is independent of the nozzle aspect ratio and the nozzle exit Reynolds number (in the range 6670–13,340). There is general agreement between the reattachment lengths determined by LDA and those determined using wall surface oil film visualisation technique. The role of coherent structures arising from initial instabilities of a 30° wall jet has been explored by hot-wire spectra measurements. Results indicate that the fundamental vortex roll-up frequency in both the inner and outer shear layer corresponds to a Strouhal number (based on nozzle exit momentum thickness and velocity) of 0.012. The spatial development of instabilities in the jet has been studied by introducing acoustic excitation at a frequency corresponding to the shear layer mode. The formation of the fundamental and its first subharmonic has been identified in the outer shear layer. However, the development of the first subharmonic in the inner shear layer has been severely suppressed. Distributions of mean velocities, turbulence intensities and Reynolds shear stress indicate that controlled acoustic excitation enhances the development of instabilities and promotes jet reattachment to the wall, resulting in a substantially reduced recirculation flow region. Received: 24 November 1998/Accepted: 24 August 1999     

Detached direct numerical simulations of turbulent two-phase bubbly channel flow

DNS simulations of two-phase turbulent bubbly channel flow at Re_τ = 180 (Reynolds number based on friction velocity and channel half-width) were performed using a stabilized finite element method (FEM) and a level set approach to track the air/water interfaces.     

悬臂输流管道在基础激励与脉动内流联合作用下的参激振动

首先建立了悬臂输流管道在基础激励与脉动内流联合作用下的运动方程;然后基于Galerkin法研究了该系统的非线性动力学行为,分析了系统运动状态随激励频率和相位差的变化,以及混沌百分比随频率比(基础激励频率与脉动频率之比)和相位差的变化。结果表明,无论以激励频率还是以相位差为分岔参数,系统都具有多种形式的动态响应,包括周期运动、概周期运动和混沌运动,但进入和脱离混沌的途径不同。相位差和频率比对系统的混沌百分比有重要影响:相位差为π/2时系统混沌百分比最大;频率比为1时系统混沌百分比最小,频率比较小或较大时系统混沌百分比与只有基础激励时接近。     

Damping of tubes due to internal two-phase flow

Two-phase internal flow is present in many piping system components. Although two-phase damping is known to be a significant constituent of the total damping, the energy dissipation mechanisms that govern two-phase damping are not well understood. In this paper, damping of three different clamped–clamped tubes subjected to two-phase air–water internal axial flow is investigated. Experimental data are reported, showing a strong dependence of two-phase damping on void fraction, flow velocity and flow regime. Data-points plotted on two-phase flow pattern maps indicate that damping is greater in a bubbly flow regime. The two-phase damping ratio reaches a maximum value at the highest void fraction before the transition to a churn flow regime. An analytical model that relates the two-phase damping ratio to the interface surface area is proposed. The model is based on rigid spherical bubbles in cubic elementary flow volumes. The analytical results are well correlated with the experiments.     

Sensitivity of plane Poiseuille flow to vibration of the channel walls

An accurate quantitative investigation of the disturbances induced in a Poiseuille flow by vibration of the walls is made on the basis of the Fourier transformation method.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 12–19, July–August, 1992.     

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.     

Void fraction and pressure fluctuations of bubbly flow in a vertical annular channel

In this investigation some hydrodynamic characteristics of two phase, two component, air water bubbly flow in a vertical annulus were studied. In particular, the void fraction profiles, and the pressure fluctuations were measured by the electrical resistivity probe and a capacitive type differential transducer respectively. These measurements were assessed under various system parameters, viz the air and water flux, the perforation ratio (Area of holes/channel cross sectional area) and the dimensionless axial distance. In addition, the pressure drop calculated from the void fraction measurements was in very good agreement with the corresponding one measured by the pressure transducers.List of symbols D _eq equivalent diameter of the annular channel (m) - j flux (discharge/channel cross sectional area) (m/s) - m mass flow rate (kg/s) - P pressure (Pa) - AP static pressure difference along the test section (Pa) - P pressure fluctuations (Pa) - P ^* dimensionless pressure (P _m/P _S.P. ) - P dimensionless pressure fluctuations (P _max /P _T.P. ) - r radius (m) - z axial distance (m) Greek symbols void fraction - dimensionless axial distance (Z/Dimeq) - perforation ratio (area of holes/channel cross sectional area) - density (kg/m3) - time (s) - dimensionless radial distance (r–r _i )/(r _o-r _i ) Suffix g gas - i inner - L liquid - m mean - Max Maximum - O outer - S.P. single-phase - T.P. two-phase     

Stability of a flexible insert in one wall of an inviscid channel flow

A hybrid of computational and theoretical methods is extended and used to investigate the instabilities of a flexible surface inserted into one wall of an otherwise rigid channel conveying an inviscid flow. The computational aspects of the modelling combine finite-difference and boundary-element methods for structural and fluid elements respectively. The resulting equations are coupled in state-space form to yield an eigenvalue problem for the fluid–structure system. In tandem, the governing equations are solved to yield an analytical solution applicable to inserts of infinite length as an approximation for modes of deformation that are very much shorter than the overall length of the insert. A comprehensive investigation of different types of inserts – elastic plate, damped flexible plate, tensioned membrane and spring-backed flexible plate – is conducted and the effect of the proximity of the upper channel wall on stability characteristics is quantified. Results show that the presence of the upper-channel wall does not significantly modify the solution morphology that characterises the corresponding open-flow configuration, i.e. in the absence of the rigid upper channel wall. However, decreasing the channel height is shown to have a very significant effect on instability-onset flow speeds and flutter frequencies, both of which are reduced. The channel height above which channel-confinement effects are negligible is shown to be of the order of the wavelength of the critical mode at instability onset. For spring-backed flexible plates the wavelength of the critical mode is much shorter than the insert length and we show very good agreement between the predictions of the analytical and the state-space solutions developed in this paper. The small discrepancies that do exist are shown to be caused by an amplitude modulation of the critical mode on an insert of finite length that is unaccounted for in the travelling-wave assumption of the analytical model. Overall, the key contribution of this paper is the quantification of the stability bounds of a fundamental fluid–structure interaction (FSI) system which has hitherto remained largely unexplored.     

Kinetic model of bubbly flow

A kinetic approach based on the approximate calculation of the fluid flow potential and formulation of Hamilton’s equations for generalized coordinates and momenta of bubbles is employed to describe processes of collective interaction of gas bubbles moving in an inviscid incompressible fluid. Kinetic equations governing the evolution of the distribution function of bubbles are derived. These equations are similar to Vlasov equations. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 5, pp. 130–138, September–October, 2000.     

Analytic solution to the micropolar-fluid flow through a semi-porous channel with an expanding or contracting wall

The flow of a micropolar fluid in a semi-porous channel with an expanding or contracting wall is investigated. The governing equations are reduced to ordinary ones by using similar transformations. To get the analytic solution to the problem, the homotopy analysis method （HAM） is employed to obtain the expressions for velocity fields. Graphs are sketched and discussed for various parameters, especially the effect of the expansion ratio on velocity and micro-rotation fields.     

Boundary-layer flow of a micropolar fluid due to a stretching wall

Summary  A Theoretical analysis is carried out to study the boundary-layer flow over a continuously moving surface through an otherwise quiescent micropolar fluid. The transformed boundary-layer equations are solved numerically for a power-law surface velocity using the Keller-box method. The effects of the micropolar K and exponent m parameters on the velocity and microrotation field as well as on the skin-friction group are discussed in a detailed manner. It is shown that there is a near-similarity solution of this problem. The accuracy of the present solution is also discussed. Accepted for publication 1 April 1996     

Effect of internal flow on vortex-induced vibration of risers

Although there are many studies dedicated to the problem of vortex-induced vibration (VIV) of marine risers, VIV experiments with internally flowing fluid have not been carried out before. In order to investigate this area, the present experiment with an internally flowing fluid and external current was designed. The riser was towed in the water flume with varying internal flow speeds. The tests in still water and in a current were conducted successfully. Various measurements were obtained including the frequency responses and the time-domain tracing of in-line and cross-flow responses. The experimental results exhibit several valuable features. First, with an increase in internal flow speed, the response amplitude increases while the vibration frequency decreases. Secondly, internally flowing fluid lessens the correlation of the vibration between different sections. In addition, by plotting both in-line strain and cross-flow strain simultaneously, a figure-of-eight for bending strain is also observed, and the trajectories in different cycles are more concordant with the increase of internal flow speed.     

Turbulent shear stress profiles in a bubbly channel flow assessed by particle tracking velocimetry

Particle tracking velocimetry (PTV) is applied to a bubbly two-phase turbulent flow in a horizontal channel at Re = 2 × 10⁴ to investigate the turbulent shear stress profile which had been altered by the presence of bubbles. Streamwise and vertical velocity components of liquid phase are obtained using a shallow focus imaging method under backlight photography. The size of bubbles injected through a porous plate in the channel ranged from 0.3 to 1.5 mm diameter, and the bubbles show a significant backward slip velocity relative to liquid flow. After bubbles and tracer particles are identified by binarizing the image, velocity of each phase and void fraction are profiled in a downstream region. The turbulent shear stress, which consists of three components in the bubbly two-phase flow, is computed by analysis of PTV data. The result shows that the fluctuation correlation between local void fraction and vertical liquid velocity provides a negative shear stress component which promotes frictional drag reduction in the bubbly two-phase layer. The paper also deals with the source of the negative shear stress considering bubble’s relative motion to liquid.     

Enhancement of heat transfer due to bubbles passing through a narrow vertical rectangular channel

An experimental study has been made of enhancement of heat transfer for water at a subcooling of 70, 40, 20, 10 K and for ethanol at a subcooling of 60, 30 K in a narrow vertical rectangular channel (0.02 m wide, 0.002 m space and 0.2 m long) at atmospheric pressure, in which a vertical heated surface (0.02 m wide and 0.01 m long) is located on one side at a position of 0.15 m from its entrance and air bubbles are forcibly passed through it at a designated period from 0.1 to 10 s. The length of the bubble can be controlled to be equal to 0.03, 0.02, and 0.01 m in length. In addition, a theoretical analysis on the basis of a transient thermal conduction model has been made to elucidate the mechanism of the enhancement. The experiment shows that the heat transfer coefficients at the subcooling of 40 and 70 K are determined only by the period of the passing bubble within the experimental range. The experimental result is a similar tendency to a theoretical one. The difference between both results arises from the effect of the flow which is not taken into account in the analysis.Es wurde eine experimentelle Untersuchung der Verbesserung der Wärmeübertragung von Wasser bei einer Unterkühlung von 70, 40, 20, 10 K und von Ethanol bei einer Unterkühlung von 60 und 30 K in einem schmalen, senkrechten, rechtwinkeligen Kanal (0,02 m breit, 0,002 m Zwischenraum und 0.2 m lang) bei Atmosphären-Druck durchgeführt. In dem Kanal ist eine senkrecht beheizte Oberfläche (0,02 m breit und 0,01 m lang), 0,15 m vom Kanaleingang entfernt, an einer Seite plaziert. Luftblasen strömen in einer Periode von 0,1 bis 10 s an der Oberfläche vorbei. Die Größe der Luftblasen liegt im Bereich von 0,03, 0,02 und 0,01 m. Um die Wirkungsweise der Verbesserung in der Wärmeübertragung erklären zu können, wurde eine theoretische Untersuchung auf der Basis eines transienten thermischen Leitungsmodells durchgeführt. Das Experiment zeigt, daß die Wärmeübertragungskoeffizienten bei einer Unterkühlung von 40 und 70 K nur von der Periode der vorbeiströmenden Blasen abhängen. Das experimentelle Ergebnis zeigt ähnliche Tendenzen wie das Theoretische. Der Unterschied zwischen beiden Ergebnissen entsteht durch den Effekt, den die Strömung auf die Untersuchung ausübt und der bei der Analyse nicht in Betracht gezogen wurde.     

Enhancement of heat transfer due to bubbles passing through a narrow vertical rectangular channel (Change in heat transfer along flow)

The objective of this paper is to make clear how heat transfer coefficient changes along the flow with the passing bubbles through a narrow vertical rectangular channel (20 mm wide, 2 mm deep and 450 mm long). The experiments were done using subcooled water of 80, 60, and 40 K at atmospheric pressure in which the air bubbles were injected into the channel at a designated period from 0.125 to 1.0 s and their length was controlled to be equal to 0.03, 0.02, and 0.01 m. The experiment shows that the heat transfer coefficients decrease along the flow and then reach a constant value beyond a certain distance from the leading edge of the heated surface where the flow becomes fully developed in both the velocity and the thermal conditions. Under the fully developed conditions, the heat transfer coefficients are predicted well by the existing theoretical analysis in which both the convective term and evaporation on the interface are ignored.     

Effect of wall conductances on convective magnetohydrodynamic channel flow

Summary The problem of convective magnetohydrodynamic channel flow in a vertical channel subjected simultaneously to an axial temperature gradient and a pressure gradient is examined when the thermal and the electrical conductance of the channel walls are arbitrary. The effects of wall conductances on the flow rate and heat transfer are found and discussed. When the vertical temperature gradient is negative, which is the case of heating from below, there exists a critical Rayleigh number at which the fluid becomes unstable. The critical Rayleigh number is also found as a function of the wall conductances.On leave from the State University of New York at Buffalo.     

简谐激励作用下轴向运动梁横向振动特性分析

为分析简谐激励作用下轴向运动梁的横向振动问题,采用单元数目及长度固定不变、节点参数在不同时间步下无缝传递的节点生死方法,建立了时变系统的动力学有限元模型,通过已有实例验证了模型的准确性和有效性。在此基础上,分析了架设速度、激励力频率和幅值对某型平推式军用桥梁架设过程横向动力响应的影响规律。结果表明,在架设过程中,当桥梁的时变固有频率与激励力频率接近时,桥梁位移动态响应呈共振的特点,据此提出了减小某型平推式军用桥梁架设过程动力响应的措施。     

Development of a robust image processing technique for bubbly flow measurement in a narrow rectangular channel

This paper presents a robust image processing technique for bubbly flow measurement over a wide range of void fractions. The proposed algorithm combines geometrical, optical and topological information recorded with high speed cameras to separate and reconstruct the overlapping bubbles. The common difficulties such as overlapping, irregular bubble shape, surface deformation and large clustering in digital image processing are solved by combining different information based on a preset decision table and flow chart. Test with synthetic bubble images is performed to evaluate the reliability of the algorithm and quantify the uncertainty of the data. The result shows that the proposed algorithm can accurately measure bubbly flows with void fraction up to 18% for large bubbles. Four runs of bubbly flow images in a 30 mm  ×  10 mm rectangular channel are then recorded by three high speed cameras. The area-averaged void fraction of these test runs range from 2.4% to 9.1%. The axial and lateral distributions of bubble number density are obtained by the present algorithm for studying the characteristics of these flows.