含气泡油滴撞击油膜壁面时气泡的变形与破裂

油--气润滑过程中润滑油液滴受高速气流扰动易形成含气泡油滴,微气泡将对油滴撞击壁面时的运动过程以及壁面油膜 层的形成质量产生重要影响. 基于耦合的水平集--体积分数 方法,对含气泡油滴撞击油膜壁面行为进行数值模拟研究, 考察含气泡油滴撞击油膜壁面时气泡的变形运动过程,探讨气泡破裂的动力学机制,分析气泡大小、碰撞速度和液体黏度等因素对含气 泡油滴撞壁过程中气泡变形特征参数的影响规律. 研究表明:含气泡油滴撞击油膜壁面后气泡会发生变形,并破裂形成膜液滴;气泡随同 液滴运动过程中,气泡内外压力和速度梯度变化是使气泡发生破裂的主要诱因. 气泡大小对气泡破裂方式影响较大,气泡较小时发生单 点破裂,而气泡较大时更容易发生多处破裂. 不同大小气泡受力差异较大,气泡大小与破裂发生时刻没有明显相关性. 碰撞速度和液体 黏度对气泡的变形、破裂和破裂发生时刻都具有一定的影响. 碰撞速度越大,油滴动能越大,更容易产生气泡变形和破裂现象. 液体黏 度增大,在油滴撞壁运动前期促进气泡变形,而在运动后期可以阻延气泡破裂行为发生.     

浮力气泡对水平壁面的回弹动力学特性

黏性液体中的气泡浮升运动有趣而又复杂,而气泡与固壁边界的相互作用更是广泛存在于实际工程中.基于轴对称数值计算,模拟了浮力驱动下气泡在黏性液体中上升并与顶部水平固壁面碰撞、回弹的过程.采用考虑表面张力的不可压、变密度Navier-Stokes方程来描述气液两相流流动,并通过基于分级八叉树的有限体积法进行数值求解.为准确捕捉气泡在回弹过程中局部而迅速的拓扑变化,采用了动态自适应网格技术耦合流体体积法(volume of fluid,VOF)来重构气泡的形状. 从气泡对壁面的碰撞和回弹的基本现象入手,研究了伽利略数 Ga和接触速度$U_{a}$对气泡回弹动力学特性的影响, 分析了气泡碰撞过程中涡结构的变化.用回弹高度$H$、回弹周期$T$、长宽比{$A_{r}$}、浮升速度$U$、轴向位置$z$和回复系数$C_{r}$等参数来表征不同条件时气泡的运动和形状特性. 研究结果表明,气泡的回弹运动特性对 Ga十分敏感. Ga的增大可加剧气泡形变, 促进气泡的回弹运动, 增多回弹次数,增大回弹参数($T$和$H)$, 提升回复系数. 然而,接触速度并非决定气泡回弹动力学的控制参数, $U_{a}$的改变并不会改变回复系数.      

浮力气泡对水平壁面的回弹动力学特性

黏性液体中的气泡浮升运动有趣而又复杂,而气泡与固壁边界的相互作用更是广泛存在于实际工程中.基于轴对称数值计算,模拟了浮力驱动下气泡在黏性液体中上升并与顶部水平固壁面碰撞、回弹的过程.采用考虑表面张力的不可压、变密度Navier-Stokes方程来描述气液两相流流动,并通过基于分级八叉树的有限体积法进行数值求解.为准确捕捉气泡在回弹过程中局部而迅速的拓扑变化,采用了动态自适应网格技术耦合流体体积法(volume of fluid,VOF)来重构气泡的形状. 从气泡对壁面的碰撞和回弹的基本现象入手,研究了伽利略数 Ga和接触速度$U_{a}$对气泡回弹动力学特性的影响, 分析了气泡碰撞过程中涡结构的变化.用回弹高度$H$、回弹周期$T$、长宽比{$A_{r}$}、浮升速度$U$、轴向位置$z$和回复系数$C_{r}$等参数来表征不同条件时气泡的运动和形状特性. 研究结果表明,气泡的回弹运动特性对 Ga十分敏感. Ga的增大可加剧气泡形变, 促进气泡的回弹运动, 增多回弹次数,增大回弹参数($T$和$H)$, 提升回复系数. 然而,接触速度并非决定气泡回弹动力学的控制参数, $U_{a}$的改变并不会改变回复系数.     

三维气泡羽流建立区流动形态的观察

本文对三维气泡羽流(Bubble plume)建立区流动形态进行了实验研究,给出实验参数的确定、实验方法与装置以及所取得的结果。     

壁面约束对裙带气泡动力学的影响

基于流体体积法(volume of fluid,VOF),数值模拟了装满黏性液体的圆柱形汽缸中的裙带气泡的浮升运动,研究了侧壁面约束对裙带气泡浮升动力学的影响.用雷诺数(Re)、韦伯数(We)、长宽比(χ)、裙带厚度(T/d)和裙带长度(L/d)等参数来表征不同约束比条件下(1.1≤Cr≤10)裙带气泡的运动和变形特性,分别在全局参考系和局部参考系下分析了壁面对气泡内外流场的影响.模拟结果显示,当Cr≥8时,裙带气泡的行为特性与在无界流域条件下的情况相当,可视作壁面无关的.当Cr8时,壁面对裙带气泡的浮升速度和形状演化有显著影响.随着壁面的靠近,裙带气泡受到的阻力增大,造成浮升速度下降.约束比降低使裙带厚度增厚而长度变短直至裙带消失,裙带气泡受挤压而被拉长并逐渐变为椭圆球帽形最后到子弹形.相反,约束比增大时,裙带气泡尾流效应增强,气泡边缘处流场产生明显的循环流动(涡环),促使裙带的形成.研究表明壁面会加剧裙带气泡产生破碎,印证了前人的推断.模拟结果与已有的经验公式吻合良好,分析了前人公式的适用性.     

表面粗糙度对微细管内气体流动特性的影响

采用了表面粗糙度粘性系数模型对微细管内的气体流动进行数值模拟，以研究微管内壁表面粗糙度对微管内气体流动的影响。运用本文改进的表面粗糙度粘性系数模型，数值模拟与实验数据十分吻合。计算结果表明，进出口压力一定时，表面粗糙度对流场的压力、密度及温度分布的影响不大，但是对速度场的影响十分显著，表面粗糙度使气体流动速度减小，并使壁面附近的速度梯度减小，从而使通过管道的气体质量流量减小，在微管内的气体流动中，表面粗糙度的影响是不能被忽略的。     

多孔介质中气泡的生成机理及流动特性研究

根据组成多孔介质的颗粒直径、颗粒排列方式、孔喉比、束缚水饱和度,将多孔介质简化为一簇变截面毛管组成的毛管束.利用质量守恒原理和动量定理,推导出单个孔隙单元内液相压力梯度和壁面剪切应力的计算公式,得到泡沫稳定流动状态时多孔介质中的压力分布以及泡沫当量直径的计算方法.通过对喉道处气泡的受力分析,推导了截断再生气泡的当量直径计算公式.利用长U型填砂管对多孔介质中稳定泡沫的流动特性进行了实验研究.研究结果表明:稳定泡沫流动时多孔介质中的压力分布呈线性下降,孔喉结构和泡沫干度是影响泡沫封堵能力的主要因素.     

滑动油膜轴承动态特性系数的迭代识别法

本文主要研究油膜轴承动态特性系数的现场识别方法,针对现场识别无法施加人工激振、测点较少,转子不平衡量未知,测点偏离轴承中心等实际限制条件,文中提出一种迭代识别法,直接利用转子不平衡响应(位移时间历程)来识别实际工况下油膜轴承的动态特性系数,最后通过计算机仿真算例分析,取得了满意的结果.     

喷管对水下爆轰气泡形态和压力特征影响的实验研究

针对具有不同类型喷管的爆轰管在水下爆轰中形成的燃气射流问题,搭建了爆轰实验平台,研究了单次爆轰过程中尾部喷管对水下气泡形态与压力特征的影响。采用数字粒子图像测速技术对高速摄影机拍摄得到的气泡脉动图片进行流场可视化分析,得到各喷管工况下的气泡速度场。为了确认爆轰管内是否形成稳定爆轰波,并观察爆轰波在气液两相界面上的透反射特性,爆轰管尾部安装有2个动态压力传感器,同时在距离喷管一定距离处设置一个水下爆炸传感器,以监测水中传播的压力波。结果表明：扩张喷管工况下的气泡脉动过程与直喷管工况基本一致,但扩张喷管提高了燃气射流速度,气泡膨胀体积更大;因为燃气射流的持续性,收敛喷管工况下的气泡脉动过程具有明显差异,气泡膨胀体积较小,但气泡二次脉动时长相较于一次脉动时长衰减更小;扩张喷管提高了气泡脉动强度,扩张喷管工况下的气泡脉动压力与透射冲击波压力远大于直喷管工况下的气泡脉动压力与透射冲击波压力;收敛喷管工况下的气泡脉动压力与透射冲击波压力都较小,但收敛喷管燃气射流的持续性减缓了气泡脉动压力的衰减速度。相比于直喷管,扩张喷管工况下的气泡脉动时间、气泡脉动压力与透射冲击波压力都更大。收敛喷管工况下的气泡...     

穿孔壁面通道流动阻力特性研究

利用计算流体力学方法对穿孔壁面通道内部流动进行了数值模拟,分析了其流动沿程摩擦阻力特性。计算中,假设流体不通过小孔进出穿孔壁面背后的吸声材料。经实验对比,采用Realizable k--\" 湍流模型结合增强型壁面函数,可获得比较理想的计算结果,数值模拟与实验测量的压损值相对偏差在10% 以内。对于穿孔壁面通道流动阻力计算,提出了等效粗糙度概念,结合达西公式、克罗布鲁克公式计算获得了等效粗糙度数值。研究了穿孔板小孔孔径、穿孔率对等效粗糙度的影响,发现等效粗糙度与孔径成二次方关系、与穿孔率成线性关系。     

A model for droplet entrainment in heated annular flow

The ability to accurately predict droplet entrainment in annular two-phase flow is required to effectively calculate the interfacial mass, momentum, and energy transfer, which characterizes nuclear reactor safety, system design, analysis, and performance. Most annular flow entrainment models in the open literature are formulated in terms of dimensionless groups, which do not directly account for interfacial instabilities. However, many researchers agree that there is a clear presence of interfacial instability phenomena having a direct impact on droplet entrainment. The present study proposes a model for droplet entrainment, based on the underlying physics of droplet entrainment from upward co-current annular film flow that is characteristic to light water reactor safety analysis. The model is developed based on a force balance and stability analysis that can be implemented into a transient three-field (continuous liquid, droplet, and vapor) two-phase heat transfer and fluid flow systems analysis computer code.     

The characteristics of a vapour bubble moving in a non-uniform flow field

The present paper describes the dynamic process of a vapour bubble moving in a non-uniform flow field. The coupling between the bubble moving as a whole and the deformation of the bubble surface is considered. The effect of the pressure gradient on the bubble movement is analysed. For a given flow field the numerical calculation is carried out until the vapour bubble is split by a micro-jet.     

Numerical study on effects of check valve closure flow conditions on pressure surges in pumping station with air entrainment

This paper proposes a numerical procedure to better compute the characteristics of pressure surges when check valves close under different flow conditions in a pumping station. Studies have shown that the effects of check valve closure on the pressure transients are predominantly dependent on the magnitude and gradient of the flow velocities immediately downstream of the check valve at the instant of valve closure. Through the present study, it was noted that the transient flow velocities near the check valve of a fluid system are also dependent on the characteristics of the air entrained into the fluid system. An improved numerical computational procedure for the fluid system with air entrainment under different transient conditions downstream of the check valve is also proposed in this paper. With a fluid system operating within the critical range of air entrainment values, the present analysis showed that there is a possibility of ‘high pressure surges’ when the check valves were closed at flow rates other than the positive flow conditions. This phenomenon was confirmed through field observations. This study thus concludes that a detailed numerical transient analysis of the fluid system, with various assumed amounts of entrained air, is necessary whenever there is the possibility of air entrainment into the fluid system, and that the flow conditions at the instant of check valve closure need to be modelled. Copyright © 2001 John Wiley & Sons, Ltd.     

Prediction of the slug-to-churn flow transition in vertical two-phase flow

An assessment is made of the various viewpoints on the slug-to-churn flow transition in vertical upward flow in the light of recent experimental results obtained at Harwell Laboratory. It is found that the flooding model of McQuillan & Whalley and the bubble entrainment model of Barnea & Brauner give satisfactory results at low and high liquid flow rates, respectively. An improved model for flooding, which takes account of the effect of the falling film, has been proposed. It is shown that this new model is in good agreement with experimental results at both low and high liquid flow rates.     

Interpretation of atomization rates of the liquid film in gas-liquid annular flow

For vertical gas-liquid annular flow the fraction of the liquid in the gas is controlled by the rate of atomization of the liquid film flowing along the wall and the rate of deposition of droplets entrained in the gas. Measurements of the rate of atomization are interpreted by a Kelvin-Helmholtz mechanism. Small wavelets on the liquid film are visualized to be entrained when wave-induced variations in the gas pressure cannot be counterbalanced by surface tension effects.     

Long waves on a viscoelastic film flow down a wavy incline

Long waves on a viscoelastic film flow down a wavy inclined plane is investigated. The analysis is performed to see how long non-linear waves on viscoelastic film down an uneven inclined wall are deformed due to the non-uniformity of the basic flow. The results are then compared with those corresponding to Newtonian film down a wavy inclined wall as well as viscoelastic film down a plane inclined wall.     

Disengagement of the gas phase in bubble columns

The technique of gas disengagement is popularly used to assess the bubble size distribution in bubble columns. The technique involves the dynamic measurement of dispersion height when the gas supply is stopped. In this paper a mathematical model has been proposed for the process of dynamic gas disengagement. It has been shown that the initial faster disengagement is due to the presence of internal liquid circulation and not due to the presence of very large bubbles. Further, slower disengagement has been attributed to the transition from heterogeneous dispersion to homogeneous dispersion. The new model also explains the effects of superficial gas velocity, column diameter, column height and liquid phase physical properties on the gas disengagement.     

Study of gas-sheared liquid film in horizontal rectangular duct using high-speed LIF technique: Three-dimensional wavy structure and its relation to liquid entrainment

The flow of a liquid film sheared by high velocity gas stream in a horizontal rectangular duct was investigated using a high-speed laser-induced fluorescence technique. Measurements of local film thickness were resolved in both longitudinal and transverse coordinates with high spatial and temporal resolution. It was found that the generation of fast and slow ripples by the disturbance waves was qualitatively the same as it was observed earlier in completely different conditions. The transverse size and curvature of the disturbance waves and ripples were measured. A relationship between the three-dimensional structure of ripples on top of disturbance waves and the two mechanisms of liquid entrainment, known as ‘bag break-up’ and ‘ligament break-up’, is proposed.     

The drainage and rupture of a non-foaming liquid film formed upon bubble impact with a free surface

Experiments are performed to measure the thickness of a thin liquid film formed between a free surface and the apex of a bubble which is approaching at its terminal velocity. Measurements are made using bubbles of < 1 mm dia both in distilled water and alcoholic solutions (methanol, ethanol). The results of the experiments show that the models based on the lubrication approximation fail to predict the initial stage of drainage. A better agreement is obtained by modifying the model proposed by Chesters.     

Weakly nonlinear stability analysis of a falling film with countercurrent gas flow

Weakly nonlinear stability analysis of a falling film with countercurrent gas–liquid flow has been investigated. A normal mode approach and the method of multiple scales are employed to carry out the linear and nonlinear stability solutions for the film flow system. The results show that both supercritical stability and subcritical instability are possible for a film flow system when the gas flows in the countercurrent direction. The stability characteristics of the film flow system are strongly influenced by the effects of interfacial shear stress when the gas flows in the countercurrent direction. The effect of countercurrent gas flow in a falling film is to stabilize the film flow system.