垂直上升气液两相弹状流模型

基于等效弹单元思想,改进了预测垂直上升管中充分发展气液弹状流流动特性的模型。 模型中考虑了界面切应力对液膜运动的影响;并在液弹空隙度预测中引入临界气体夹带速度的概念,以此来描述弹状流中大气泡尾部的混合特性。本文提出的模型还考虑了管径对液弹空隙度的影响。弹状流模型的计算结果得到本文及其他作者实验数据的验证。     

Void fraction in vertical gas-liquid slug flow: Influence of liquid slug content

In this study we develop a model for computing the mean void fraction and the liquid slug void fraction in vertical upward gas-liquid intermittent flow. A new model for the rate of gas entrained from the Taylor bubble to the liquid slug is formulated. It uses the work done by the pressure force at the rear of the Taylor bubble. Then an iterative approach is employed for equating the gas entrainment flux and the gas flux obtained via conservation equations. Model predictions are compared with experimental data. The developed iterative method is found to provide reasonable quantitative predictions of the entrainment flux and of the void fraction at low and moderate liquid slug void fraction conditions. However, with an increased liquid slug void fraction experimental data indicate that the flow in the liquid slug transits to churn-heterogeneous bubbly flow thus gas entrainment flux tends to zero. Considering this effect in the iterative model significantly improved the predictions for large liquid slug void fraction conditions.     

矩形截面螺旋通道内气液两相流局部含气率分布实验研究

应用电导探针测量技术,对矩形截面螺旋通道内气液两相流局部含气率进行实验研究。在不同的气相折算速度下,应用电导探针测量了弹状流弹单元的长度,并与可视化方法进行对比,验证了电导探针的可靠性,并为信号处理选择合适的阈值。分别在泡状流、弹状流及环状流三种流型的条件下,分析了气相与液相折算速度对局部含气率分布的影响。实验结果发现,螺旋通道气液两相局部含气率呈非对称的抛物线形分布,这种非对称性受流型和液相折算速度的影响。     

用统计方法研究未充分发展弹状流动

在一内径19 mm、长2 m的垂直有机玻璃管内,采用自制的电导探针对未充分发展的气-液二相弹状流中的弹状气泡上升速度、液塞上升速度、弹状气泡长度和液塞长度进行了测量。得到了各自随表观气速或表观液速的变化规律。结果表明:在未充分发展的弹状流状态下,弹状气泡的上升速度略高于液塞的上升速度:弹状气泡长度随表观气速的增大而增大,随表观液速的增大而减小。文章对弹状气泡长度进行了统计分析。未充分发展弹状流中弹状气泡长度符合正态分布律。     

Influence of hydrodynamics on liquid mixing during Taylor flow in a microchannel

The miscible liquid‐liquid two phases based on Taylor flow in microchannels was investigated by high‐speed imaging techniques and Villermaux/Dushman reaction. The mixing based on Taylor flow was much better compared with that without introducing gas in microchannels, even the ideal micromixing performance could be obtained under optimized superficial gas and liquid velocities. In the mixing process based on Taylor flow, the superficial gas and liquid velocities affected the lengths and the velocities of Taylor bubble and liquid slug, and finally the micromixing performance. The formation process of Taylor flow in the inlets, the initial uniform distribution of reactants and the internal circulations in the liquid slug, and the thin liquid films all improved the mixing performance. Furthermore, a modified Peclet number that represented the relative importance of diffusion and convection in the mixing process was proposed for explaining and anticipating micromixing efficiency. © 2011 American Institute of Chemical Engineers AIChE J, 58: 1660–1670, 2012     

Experimental and numerical investigation of liquid circulation induced by a bubble plume in a baffled tank

Bubble induced liquid circulation is important in applications such as bubble columns and air-lift reactors. In this work, we describe an experimental and numerical investigation of liquid circulation induced by a bubble plume in a tank partitioned by a baffle. The baffle divides the tank into two compartments. Liquid can flow from one compartment to the other through openings at the top and the bottom of the baffle. Gas (air) was injected in the riser section in the form of bubbles at one corner of the tank. The temporal and spatial variation of velocity field in the liquid as a function of the gas flow rate was measured using particle image velocimetry (PIV). At a constant gas flow rate, the liquid flow field is unsteady due to the interaction with the bubbles. The time scales associated with the velocity-time series and the bubble plume thickness variation were calculated. The time averaged-velocity field was used to quantify the variation of the liquid circulation rate with gas flow rate. The turbulence in the liquid was measured in terms of turbulent intensities. These were calculated from the experimental data and were observed to be less than 3 cm/s. A 2-d Euler-Euler two-fluid model with buoyancy and drag as the interaction terms was used to simulate the flow. The parameters chosen for the simulations were selected from literature. It is shown that inclusion of turbulence model such as k-ε is necessary to capture the overall flow behavior. Good agreement was observed between experimentally obtained velocity profiles and the recirculation rates with the simulation results.     

ELECTROOSMOTIC HELICAL FLOW PRODUCED BY COMBINED USE OF LONGITUDINAL AND TRANSVERSAL ELECTRIC FIELDS IN A RECTANGULAR MICROCHANNEL

We theoretically devise and simulate a microelectrode system that produces electroosmotic helical flow in a straight rectangular microchannel. In addition to a pair of primary electrodes that generate a longitudinal electric field, sets of secondary electrodes are installed to produce a transversal electric field. The secondary electrode system consists of point electrodes embedded along two edges of the bottom surface of the channel. The transversal electric field developed across the bottom surface causes the electroosmotic motion of fluid at the bottom of the channel along the transversal direction. As a combined effect of the primary electrode system that produces unidirectional longitudinal flow along the channel and the secondary electrode system that produces transversal flow across the bottom surface, the flow inside a rectangular microchannel follows a helical pattern. After simulating the electroosmotic helical flows developed in microchannels we analyze the mixing of sample liquid in such flow fields by calculating the trajectories of fluid particles.     

ELECTROOSMOTIC HELICAL FLOW PRODUCED BY COMBINED USE OF LONGITUDINAL AND TRANSVERSAL ELECTRIC FIELDS IN A RECTANGULAR MICROCHANNEL

We theoretically devise and simulate a microelectrode system that produces electroosmotic helical flow in a straight rectangular microchannel. In addition to a pair of primary electrodes that generate a longitudinal electric field, sets of secondary electrodes are installed to produce a transversal electric field. The secondary electrode system consists of point electrodes embedded along two edges of the bottom surface of the channel. The transversal electric field developed across the bottom surface causes the electroosmotic motion of fluid at the bottom of the channel along the transversal direction. As a combined effect of the primary electrode system that produces unidirectional longitudinal flow along the channel and the secondary electrode system that produces transversal flow across the bottom surface, the flow inside a rectangular microchannel follows a helical pattern. After simulating the electroosmotic helical flows developed in microchannels we analyze the mixing of sample liquid in such flow fields by calculating the trajectories of fluid particles.     

Toward Understanding of Two‐Phase Eccentric Helical Reactor Performance

Mass transfer investigations in a two‐phase gas‐liquid Couette‐Taylor flow (CTF) reactor and a numerical flow simulation are reported. The CTF reactor is characterized by high values of the mass transfer parameters. Previous mass transfer investigations have yielded high values of the volumetric mass transfer coefficients (of the order of 10^–1 s^–1) and the specific interfacial area, compared to those obtained in a stirred tank (10³ m² m^–3). In order to intensify mass transfer in the CTF reactor, an eccentric rotor (rotating inner cylinder) was used. In the eccentric annulus with rotating inner cylinder, due to frequent variation of the hydrodynamic flow field parameters, nonlinear hydrodynamic conditions occurred. These conditions can influence the rate of mass transfer. The experimental results of benzaldehyde oxidation in an eccentric CTF reactor confirmed an increase in mass transfer, as against a concentric CTF reactor. Numerical simulation of the Couette‐Taylor (helical) flow was performed in a concentric and in an eccentric annulus. Calculation of parameters such as velocity, static pressure, kinetic energy and energy dissipation rate revealed a significant effect of gap eccentricity on the flow behavior.     

Electrohydrodynamic Spraying (EHDS) of Conductive Liquids and Mass Transfer in EHDS Distillation Columns

A physicomathematical model for electrohydrodynamic spraying (EHDS) of a conductive liquid and the mass transfer in a parallel plate electrodes system of an EHDS distillation column is presented. The influence of the electric field on the formation and general characteristics of the finely dispersed aerosol flow of the conductive liquid in EHDS distillation columns has been determined from the model. It allows to explain in detail the mechanism of the influence of the electric field on the mass transfer in the gas‐liquid system. In addition, the maximum possible liquid flow rate in the mass transfer section of the EHDS distillation column has been estimated. The role of process parameters is discussed in detail.     

Liquid flow pattern around Taylor bubbles in an etched rectangular microchannel

Liquid flow around Taylor bubbles and the motion of bubble interface in a rectangular microchannel etched on a microfluidic chip were investigated using a three-dimensional particle tracking method. The Taylor bubbles were generated by releasing the dissolved air in working the liquid (water) through heating the microfluidic chip to 35–55 °C and had low velocities (15–1500 μm/s). Three-dimensional velocity distributions of liquid recirculation flows surrounding the Taylor bubble head and tail were obtained by tracking submicron fluorescent particles seeded in the working liquid and the motion of the bubble interface was analyzed by monitoring the motions of the particles attached on the bubble interface. The high velocity film flow through the microchannel corners acted as a liquid jet in front of bubble head and drainage into the corners behind the bubble tail to drive the liquid recirculation flows. The bubble interface near the microchannel corners was also moved by the strong liquid shear induced from the high velocity liquid flow in the microchannel corners. This high velocity liquid flow through the corners could be considered to be driven by the pressure drop over the Taylor bubble. The pressure drop resulted from the decrease of bubble surface mobility due to tracer surfactant in the gas–liquid interface.     

Hydrogenation of 2‐ethylanthraquinone under Taylor flow in single square channel monolith reactors

The hydrogenation of 2‐ethylanthraquinone (EAQ) to 2‐ethylanthrahydroquinone (EAHQ) was carried out under Taylor flow in single square channel monolith reactors. The two opening ends of opaque reaction channel were connected with two circular transparent quartz‐glass capillaries, where Taylor flow hydrodynamics parameters were measured and further used to obtain practical flow state of reactants in square reaction channels. A carefully designed gas‐liquid inlet mixer was used to supply steady gas bubbles and liquid slugs with desired length. The effects of various operating parameters, involving superficial gas velocity, superficial liquid velocity, gas bubble length, liquid slug length, two‐phase velocity and temperature, on EAQ conversion were systematically researched. Based on EAQ conversion, experimental overall volumetric mass transfer coefficients were calculated, and also studied as functions of various parameters as mentioned earlier. The film model, penetration model, and existing semi‐empirical formula were used to predict gas‐solid, gas‐liquid, and liquid‐solid volumetric mass transfer coefficients in Taylor flow, respectively. The predicted overall volumetric mass transfer coefficients agreed well with the experimental ones. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

Viscous co‐current downward Taylor flow in a square mini‐channel

This article presents a computational study of the co‐current downward Taylor flow of gas bubbles in a viscous liquid within a square channel of 1 mm hydraulic diameter. The three‐dimensional numerical simulations are performed with an in‐house computer code, which is based on the volume‐of‐fluid method with interface reconstruction. The computed (always axi‐symmetric) bubble shapes are validated by experimental flow visualizations for varying capillary number. The evaluation of the numerical results for a series of simulations reveals the dependence of the bubble diameter and the interfacial area per unit volume on the capillary number. Correlations between bubble velocity and total superficial velocity are also provided. The present results are useful to estimate the values of the bubble diameter, the liquid film thickness and the interfacial area per unit volume from given values of the gas and liquid superficial velocities. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

微小三角形槽道内电渗流动特性分析

提出了微小三角形槽道内电渗流动理论计算方法,通过Galerkin法计算并分析了其内部的电势及速度分布,获得了温度、槽道尺寸、外加电势的电场强度、ζ电势以及电解质浓度对微小三角形槽道内电渗流动的影响规律。计算结果表明：微小三角形槽道内液体的质量流量随ζ电势、电场强度、流体温度及电解质浓度的增加而增加,随微小三角形槽道尺寸的增加先增加后减小。     

Liquid wall friction in two-phase turbulent gas laminar liquid stratified pipe flow

The main result of the present work is an analytic expression for the mean liquid wall shear stress in two-phase turbulent gas/laminar liquid stratified pipe flow. The Navier-Stokes equations are solved assuming a flat fluid-fluid interface subject to a constant interfacial shear — approximating the interfacial drag exerted by the gas. The effect of a pipe inclination is accounted for, thereby retaining the interesting two-phase phenomenon of backflow in upwardly inclined pipes. The corresponding expression for the wall shear stress distribution is obtained by formal differentiation. Its limiting behaviour in the triple points, where the fluid-fluid interface meets the pipe wall, is determined by residue calculus. The expression for the mean wall shear stress is given by integration. It is found to be a linear combination of two terms. The first term accounts for the free surface liquid flow in the absence of the gas. The corresponding approximate hydraulic diameter model is found to be in surprisingly good agreement with this term. The second term represents the shear flow associated with the interfacial drag exerted by the gas (not accounted for by the hydraulic diameter approximation). The shear flow increases the flow rate near the interface on behalf of the flow rate near the pipe wall, thus reducing the wall shear stress below the free surface flow value. Expedient evaluation of the expression for the mean wall shear stress, suitable for use in a 1-D multiphase pipe flow simulator, is facilitated by replacing certain one-parameter integrals with highly accurate rational approximations.     

IMPROVED HYDRODYNAMIC MODEL OF TWO-PHASE SLUG FLOW IN VERTICAL TUBES

Merits of the Fernandes model(Fernandes et al.1983)for two-phase slug flow in verticaltubes are reviewed in this paper.While predicting many macroscopic parameters of slug flow in verti-cal tubes,it fails to present correctly the trend that the average voidage in liquid slugs increases asthe rising velocity of Taylor bubbles is increased.It is also desirable to extend its application toelectrolyte systems, and to churn flow conditions.Based on the diagnostic analysis,the model equa-tion for gas entrainment by falling liquid film is reformulated and the influence of surface tension isalso accounted for.Development of the falling liquid film is recognized in the revised model in or-der to suit the case of short Taylor bubbles as well.The modified model predicts the variation of av-erage voidage in liquid slugs in good agreement with available experimental data.     

湿气管道积液临界气速预测的新模型

气田开发经常采用湿气集输方案。针对湿气输送管道出现的积液问题,基于分层流最小界面剪切应力准则,利用气液平界面分层流液膜区的速度分布规律,建立了求解积液临界气速的新机理模型。由分层流液膜区的流场描述和气相动量方程得到气液界面剪切应力的表达式;利用界面剪切应力函数曲线特性,可以通过界面剪切应力关于持液率求导获得临界气速。以不同文献中收集的临界气速实验数据,对新模型和其他具有代表性的湿气管道积液模型进行验证对比,表明新模型的预测精度要优于其他模型。     

壁面润湿性对微通道内Taylor流动特性的影响

壁面润湿性不仅影响着Taylor气泡的形状,同时对通道内流体流动、相变换热等有着关键的作用。采用VOF模型对T型微通道内气液两相Taylor流动进行三维数值模拟,重点研究了接触角改变对Taylor气泡流体动力学特性的影响。模拟结果与他人实验数据对比基本吻合,验证了模型的有效性。结果表明：随着接触角增大,气泡周围液含量逐渐降低,相界面也由外凸形变为内凹形。壁面越接近润湿（或疏水）状态,气液接触面的曲率就越大;当120°≤θ≤150°时Taylor气泡稳定性变差。当θ≥150°时“拖曳流态”出现,分析指出在大接触角下气体更易贴附壁面导致接触区内流场发生变化,形成的涡流减弱了水对气相的水平剪切作用,进而引起流型转变。接触角对通道内压力有着重要影响,通道中心轴向压力曲线以θ=90°为过渡,润湿状态下呈凸函数递减且p _G>p _L,疏水状态下气液进口处的压力分配改变,曲线趋势相反。     

倾斜上升弹状流中Taylor气泡运动速度研究

用高速动态分析仪对倾斜上升管中气液两相弹状流中Ｔａｙｌｏｒ气泡的运动速度进行了研究。获得了无干扰流场下Ｔａｙｌｏｒ气泡运动参数的测量结果，并分析了混合物流速及管倾角对气泡头部位置的影响，以及由此而引起的气泡漂移速度及液体速度影响系数的变化情况，在理论分析的基础上，推荐了计算气泡运动速度的实验关联式。并且与可利用的结果进行了比较，两者符合较好。