二维微柱阵列壁面对活性剂液滴铺展的影响

针对二维微柱阵列壁面上含不溶性活性剂液滴的铺展过程,采用润滑理论建立了液膜厚度和浓度演化模型,采用数值计算方法得到了液滴的铺展特征及相关参数的影响. 研究表明:活性剂液滴在微柱阵列壁面上铺展时,在壁面凸起处衍生出隆起结构,壁面凹槽处衍生出凹陷结构,随时间持续,隆起和凹陷均向两侧移动,且数量不断增加. 活性剂液膜流经凸起时,隆起高度呈驼峰形变化. 增大预置液膜厚度或活性剂初始浓度,铺展区域隆起和凹陷数量增多,液滴铺展速度加快. 增加凹槽深度或减小斜度会使毛细力作用增强,液膜破断可能性加大;增大凹槽宽度可加速活性剂液滴的铺展,加剧液膜表面波动幅度.      

STABILITY OF ULTRATHIN LIQUID FILM EVOLUTION WITH SURFACTANT

针对固体基底上厚度小于100 nm的含活性剂超薄液膜演化过程, 基于润滑理论推导出包含分离压影响的液膜厚度和活性剂浓度的演化方程, 采用正则模态法导出了描述液膜线性稳定性的特征方程, 分析了多个特征参数对线性稳定性的影响, 数值模拟了液膜厚度和活性剂浓度演化历程, 对比了模拟所得非线性结果与线性分析预测结果的一致性.结果表明:范德华力具有促进扰动增长的作用, 较强的玻恩斥力促使扰动衰减, 使液膜趋于稳定;较小的毛细力数易使液膜凹陷处发生二次失稳, 并最终导致去润湿现象发生;液膜厚度和溶于液膜内部的活性剂浓度初值越大, 液膜稳定性越强, 液膜表面活性剂浓度影响则相反;增大吸附系数不利于液膜稳定性.     

CHARACTERISTIC OF DROPLET OSCILLATION ON THE SURFACE OF RECTANGULAR HYDROPHOBIC GROOVES

液滴振荡行为是液滴运动中的重要伴随现象,具有重要科研价值.由于液滴撞击疏水沟槽板时运动行为与光滑表面明显不同,可以推测疏水沟槽表面液滴振荡特性也将会呈现与众不同的行为特点.采用高速摄像技术,研究了矩形疏水沟槽表面上水滴高度和接触线振荡行为随沟槽尺寸和撞击速度的变化规律.结果发现,矩形疏水沟槽造成的各向润湿异性使得振荡过程中水滴在平行沟槽方向上的接触线长度大于垂直方向,但并不影响水滴高度方向上衰减振荡的周期,即水滴振荡周期与沟槽间距无关;同时由于疏水沟槽表面上存在能垒束缚效应,致使水滴振荡过程中接触线的铺展和回缩运动不服从典型阻尼振荡规律,而呈现振荡数次后直接趋稳的特点.如水滴以0.61 m/s撞击时,接触线经历2次振荡后即维持稳定,但此时水滴仍在持续振荡中.另外,还初步分析了水滴振荡周期与沟槽间距无关的原因.     

Effects of surfactant on liquid film thickness in the Bretherton problem

The effects of insoluble and soluble surfactant on the motion of a long bubble propagating through a capillary tube are investigated computationally using a finite-difference/front-tracking method. Emphasis is placed on the effects of surfactant on the liquid film thickness between the bubble and the tube wall. The numerical method is designed to solve the evolution equations of the interfacial and bulk surfactant concentrations coupled with the incompressible Navier–Stokes equations. A non-linear equation of state is used to relate surface tension coefficient to surfactant concentration at the interface. Computations are first performed for soluble cases and then repeated for the corresponding clean and insoluble cases for a wide range of governing non-dimensional parameters in order to investigate the effects of surfactant and surfactant solubility. The computed film thickness for the clean case is found to be in a good agreement with Taylor’s law indicating the accuracy of the numerical method. We found that both the insoluble and soluble surfactant generally have a thickening effect on the film thickness, which is especially pronounced at low capillary numbers. This thickening effect strengthens with increasing sensitivity of surface tension to interfacial surfactant coverage mainly due to the enhanced Marangoni stresses along the liquid film. It is also observed that film thickening shows a non-monotonic behavior for variations in Peclet number. The validity of insoluble surfactant assumption is assessed for various non-dimensional numbers and it is demonstrated that insoluble assumption is valid only when capillary number is very low, i.e., Ca ≪ 1 and when surface tension is highly sensitive to interfacial surfactant coverage, i.e., the elasticity number is large.     

分离压对波状基底上活性剂液滴铺展过程的影响

针对波状基底上含不溶性活性剂液滴的铺展过程,引入受活性剂浓度影响的分离压模型,应用润滑理论建立了液滴高度和活性剂浓度演化方程组,通过数值计算方法得到了分离压作用下含活性剂液滴过程的演化特征. 研究表明:分离压作用下的液滴演化时间显著缩短,铺展速率加快,铺展前沿处衍生出的子波结构明显减少,铺展更加稳定;分离压对液滴铺展稳定性的影响与活性剂关联强度密切相关,减小引力强度系数α₁有利于促进液滴的铺展,而减小斥力强度系数α₂则起抑制作用,且放大了液滴的演化扰动能量,致使液滴铺展呈现不稳定特征;增加基底高度D或波数k均使液滴铺展速率减慢.      

MOVEMENT AND SPLASHING OF A DROPLET IMPACTING ON A WET WALL

壁面液体层的存在对液滴撞击壁面的运动具有重要的影响。采用气液两相流动相界面追踪的水平集和流体体积复合方法和壁面润湿模型，实现了液滴撞击湿润壁面运动的数值求解；在此基础上，开展了液滴撞击湿润壁面运动的研究。研究结果表明：液滴以不同速度撞击湿润壁面时，会呈现出黏附铺展、波动运动、皇冠几何体运动以及飞溅运动等几种不同的运动形态，液滴撞击湿润壁面后的压力分布是不同运动形态形成的主要原因；飞溅运动是一定条件下皇冠几何体运动的一种特殊形态，液滴从皇冠几何体侧壁顶端的飞溅分离满足毛细破碎理论；撞击速度对分离液滴的运动方向影响较小，而对壁面液体层厚度的影响则较大；撞击速度和壁面液体层厚度对分离液滴形态、飞溅分离位置、飞溅速度以及飞溅发生时刻等都具有一定的影响。     

乏油对纯滑动粗糙表面点接触热弹流润滑的影响

本文求解了纯滑动点接触热弹性流体动力润滑问题.分析中假设运动表面为光滑表面,静止表面在接触中心有1个垂直于卷吸速度方向的横向突起.将计算域划为乏油区、压力区和气穴区3个子区,乏油区和气穴区的压力计算算法一致,但是温度的计算方法不同.研究结果显示:在乏油工况下,油膜中压力和温度均在入口处新月形液面处开始建立,入口区乏油程度的增加会导致油膜中压力分布趋近于干接触状态,油膜中平均温度得以提高,最大压力和温度均出现在赫兹接触圆两侧的宏观马蹄形区域.同全膜润滑相比,油膜厚度被明显降低,并且随着乏油程度的增加,静止表面的突起逐渐被压平.     

STATIC SPREADING OF DROPLET IMPACT ON SOLID SURFACE：INFLUENCE FACTOR

通过建立液滴撞击固体平壁的静态铺展力学平衡的数学模型,从理论上得到了静态铺展半径与液滴物性参数、以及液滴与固体壁面接触角之间关系的数学表达式,将理论结果与数值模拟的结果进行了比较,两者吻合较好.比较了不同条件下液滴的静态铺展半径的变化规律,分别得到了液滴密度、体积、表面张力和接触角等因素对液滴静态铺展半径的影响规律.     

Impact of liquid drops on a rough surface comprising microgrooves

Impact of water drops on a stainless steel surface comprising rectangular shaped parallel grooves is studied experimentally. Geometric parameters of the surface groove structure such as groove depth, groove width and solid pillar width separating any two successive grooves were kept at 7.5, 136 and 66 μm, respectively. The study was confined to the impact of drops in inertia dominated flow regime with Weber number in the range 15–257. Experimental results of drop impact process obtained for the grooved surface were compared with those obtained for a smooth surface to elucidate the influence of surface grooves on the impact process. The grooves definitely influence both spreading and receding processes of impacting liquid drops. A more striking observation from this study is that the receding process of impacting liquid drops is dramatically changed by the groove structure for all droplet Weber number.     

EXPERIMENTAL INVESTIGATION ON THE EVOLUTION OF INTERNAL SOLITARY WAVE OVER A RIDGE AND ITS ACTING FORCES ON A SUBMERGED SLENDER BODY

内孤立波沿山脊地形传播过程中的结构变化使得影响海洋结构物及水下航行器安全的不确定因素增加,用大型重力式分层流试验水槽系统研究内孤立波过山脊地形及其对细长潜体作用力特性,可有效提高对其危害性机理的认识.为此,利用电导率探头阵列,结合染色标识方法,测量内孤立波演变特性,同时利用三分量传感器测量水下细长体模型受力特性.研究结果表明：山脊地形显著改变下凹型内孤立波结构,表现为坡前波幅增大,坡顶背风波面抬升,坡后波长变长;内波结构变化直接影响细长体受力特性,表现为坡前向下作用力明显增强,坡顶向上作用力突显,坡后作用力持续时间变长;细长体受力特性变化影响其运动趋势,坡前下沉运动增强、朝山脊方向纵荡运动显著.     

镓铟锡液滴撞击泡沫金属表面的运动学特性研究

常温下为液态的镓铟锡合金以其优异的导热性能在具有特殊要求的传热领域有着重要的应用价值,与传统流动介质相比较大的表面张力使得其产生的流动现象必有所区别.本文研究镓铟锡所形成的液滴撞击泡沫金属表面后所产生的铺展、回缩及回弹现象.采用高速相机拍摄液滴投影轮廓随液滴运动的变化过程,并通过图像处理获得不同撞击速度、底板表面孔径下的液滴铺展系数、中心位置轮廓高度以及液滴回弹后在空中的振动特性.研究结果表明:具有较高表面张力的镓铟锡液滴的铺展系数随无量纲时间的变化在铺展初始阶段仍满足常规流体的1/2次幂关系,只在铺展后期与底板的无量纲孔径有关系;液滴的最大铺展系数在较小无量纲孔径底板大于在光滑镍板,且随底板无量纲孔径增大而逐渐减小;在回弹过程,由于底板孔隙结构的存在使得液滴回弹后在空中的振动呈现3种形态:规则的横向和纵向振动、带旋转的横向和纵向振动以及旋转振动;最后,通过对振动频率的拟合和分析,进一步拓展了传统振动频率理论公式在非规则振动过程预测中的应用.      

离散型织构表面液滴的铺展及其接触线的力学特性分析

针对离散型织构表面上液滴的铺展过程,采用数值模拟和润湿性实验相结合的方法,引入织构润湿因子θ*,得到了不同类型的离散型织构对固体表面润湿性的影响,在此基础之上分析了液滴铺展过程中接触线的力学特性,以期从微观界面力学的角度解释微织构对液滴铺展过程的促进作用.研究表明:离散型织构增大了液滴铺展过程中的固-液接触面积,位于铺展前沿的液体分子部分浸润织构内部,导致液面曲率和液滴内部的拉普拉斯压力增大,相邻离散型织构间的液体获得了额外的驱动力和能量,铺展速度加快,平衡接触角减小;槽状离散型织构对表面润湿性的影响程度最大,液滴在其上铺展过程具有各向异性特性.另外数值仿真分析表明,接触线的钉扎效应与固体表面粗糙度的大小和微织构类型密切相关,表面粗糙度越大,钉扎效应越明显,其中槽状织构对接触线的钉扎作用还具有方向性.      

Surfactant spreading on a thin weakly viscoelastic film

A mathematical model is presented for surfactant-driven thin weakly viscoelastic film flows on a flat, impermeable plane. The Oldroyd-B constitutive relation is used to model the viscoelastic fluid. Lubrication theory and a perturbation expansion in powers of the Weissenberg number (We) are employed, which give rise to non-linear coupled evolution equations governing the transport of insoluble surfactant and thin liquid film thickness. Spreading on a Newtonian film is recovered to leading order and corrections to viscoelasticity are obtained at order We. These equations are solved numerically over a wide range of viscosity ratio (ratio of solvent viscosity to the sum of solvent and polymeric viscosities), pre-existing surfactant level and Peclet number (Pe). The effect of viscoelasticity on surfactant transport and fluid flow is investigated and the mechanisms underlying this effect are explored. Shear stress, streamwise normal stress and the temporal rate of change of extra shear stress generated from gradients in surfactant concentration dominate thin viscoelastic film flows whereas only shear stresses play a role in Newtonian thin film flows. Our results also reveal that, for weak viscoelasticity, the influence of viscosity ratio on the evolution of surfactant concentration and film thickness can be significant and varies considerably, depending on the concentration of pre-existing surfactant and surfactant surface diffusivity.     

单个微油滴弹流润滑行为的试验研究

利用光干涉技术研究了微油滴通过弹流润滑接触区的润滑行为,考察了油滴大小、卷吸速度和载荷等因素的影响.结果表明微油滴在入口区域因挤压或毛细力效应发生表面积扩展,从而影响润滑膜的形成.油滴越大,挤压扩展直径越大,形成的膜厚越大.卷吸速度越高,入口处微油滴表面积扩展越不充分,仅接触区局部形成油膜,微油滴在接触表面挤压出凹坑穿过接触区.     

Convective and diffusive surfactant transfer in multiphase liquid systems

Laser interferometry was used to investigate diffusive and convective mass transfer in a multicomponent fluid mixture with a liquid–liquid or liquid–gas interface. For this purpose, an immobile gas bubble or insoluble fluid droplet, having the shape of a short cylinder with a free lateral surface, was inserted into a thin liquid layer. In the case of non-uniform distribution of the dissolved surfactant component, the Marangoni convection near the drop/bubble was initiated by the surface tension inhomogeneities, depending on the surfactant concentration. The applied experimental techniques allowed us to study the structure and evolution of the convective flows and concentration fields in a liquid layer, which due to its small thickness were nearly two-dimensional. Making use of both the vertical and horizontal orientation of the liquid layer, we investigated the mass transfer process at different levels of the interaction between gravity and capillary forces. During the experiments, we detected new solutocapillary phenomena, which were found to be caused by oscillatory regimes of solutal convection occurring around air bubbles and chlorobenzene drops in heterogeneous aqueous solutions of alcohol with a vertical surfactant concentration gradient. The role of the oscillatory instability in the processes of drop saturation by the surfactant from its water solution and an inverse process of surfactant extraction from the drop into the surrounding homogeneous fluid (water) was determined. A reasonable explanation for the driving mechanisms of the discovered effects has been proposed.     

超疏水沟槽表面通气减阻实验研究

减阻是解决航行体提速和增程的主要技术途径之一, 对缓解日益严峻的能源危机极为重要. 在重力式管道实验系统中, 测试给出了湍流状态下不同通气速率时减阻率随雷诺数及沟槽无量纲间距的变化规律和气膜铺展状态, 对比分析了单纯超疏水表面与超疏水沟槽表面上通气时减阻效果的差异.实验板材质为无色亚克力, 沟槽结构采用机械方法加工, 并在表面喷涂超疏水涂层. 结果表明, 持续通气能解决超疏水沟槽表面气膜层流失问题, 实现气膜层长时间稳定维持; 恒定雷诺数下, 随通气速率增大, 超疏水沟槽表面气膜铺展更趋均匀, 减阻率上升; 由于通气速率影响气膜横向扩展能力, 致使恒定通气速率下, 减阻率随雷诺数的变化呈现两种模式; 在固定雷诺数及通气速率时, 减阻率随沟槽尺寸的扩大先增后减, $S^{+}\approx 76$时减阻率最大. 分析其原因在于, 沟槽结构增大沾湿面积的同时, 显著提升了通气状态下超疏水表面气膜层的稳定性, 因而展示出与超疏水表面和沟槽表面均不相同的减阻规律, 且效果更佳.      

On the dynamic behavior of a liquid droplet impacting upon an inclined heated surface

The dynamic behavior of a water droplet impinging upon a heated surface was shown to be significantly different, depending on the normal momentum of the impinging droplet before impact. This experimental study focused mainly on the effects of the impinging angle of the droplet on impact dynamics and its dependence on surface temperature. At the surface temperature of the nucleate boiling regime, disintegration of the droplet did not occur, whereas the deforming droplet adhered to the surface. The liquid film was spread and contracted several times on the horizontal surface, but the expanded droplet merely slipped without noticeable contraction on the inclined surfaces. In the film boiling regime, the impinging droplet spread over the surface as a liquid film separated from the surface by the vapor produced. Depending on the magnitude of the normal momentum of droplet, disintegration into several irregular shapes of liquid elements occurred in the case of the horizontal and 30°-inclined surfaces. The impinging droplet in the case of the 60°-inclined surface did not break up and tended to recover its original spherical shape. Received: 16 February 1999/Accepted: 9 November 1999     

Study of the impacts of droplets deposited from the gas core onto a gas-sheared liquid film

The results of an experimental study on droplet impactions in the flow of a gas-sheared liquid film are presented. In contrast to most similar studies, the impacting droplets were entrained from film surface by the gas stream. The measurements provide film thickness data, resolved in both longitudinal and transverse coordinates and in time together with the images of droplets above the interface and images of gas bubbles entrapped by liquid film. The parameters of impacting droplets were measured together with the local liquid film thickness. Two main scenarios of droplet-film interaction, based on type of film perturbation, are identified; the parameter identifying which scenario occurs is identified as the angle of impingement. At large angles an asymmetric crater appears on film surface; at shallow angles a long, narrow furrow appears. The most significant difference between the two scenarios is related to possible impact outcome: craters may lead to creation secondary droplets, whereas furrows are accompanied by entrapment of gas bubbles into the liquid film. In addition, occurrence of partial survival of impacting droplet is reported.     

Droplet impact on wetted structured surfaces

In this study, we numerically investigate the droplet impact onto a thin liquid film deposited on a structured surface with square pillars and cavities. The time evolution of crown geometry is strongly affected by the surface structure. When the thickness of the liquid film is larger than the structure height, the expanding speed of the crown base radius is independent of the structure width. However, if the liquid film thickness is equal to the structure height, the crown base expands slower as...     

Effect of surfactants on liquid loading in vertical wells

Most gas wells produce some amount of liquid. The liquid is either condensate or water. At high rates, the gas is able to entrain liquid to the surface; however, as gas well depletes, the liquid drops back in a gas well (called liquid loading) creating a back pressure on the reservoir formation. Addition of surfactants to the well to remove liquid is one of the common methods used in gas wells. Liquid loading in vertical gas wells with and without surfactant application was investigated in this study. Anionic, two types of amphoteric (amphoteric I and amphoteric II), sulphonate and cationic surfactants were tested in 2-inch and 4-inch 40-feet vertical pipes. Pressure gradient and liquid holdup are measured. Visual observation with a high speed camera was used to gain insight into the direction of foam flow in intermittent flow and foam film flow under annular flow conditions.Liquid loading is initiated when the liquid film attached to the wall in annular flow starts flowing downwards. Introduction of foam causes the gas velocity at which film reversal occurs to decrease; this shift increases with increasing surfactant concentration and it is more pronounced in 2-inch pipe than in 4-inch pipe. That is, the benefit of surfactants is much more pronounced in 2-inch pipe than in 4-inch pipe. The reason for postponement of liquid loading is reduction in the liquid holdup at low gas velocities which reduces the liquid holdup in foam flow compared to air-water flow. However, at higher gas velocities, the pressure drop in 2-inch compared to 4-inch pipe increases rapidly as the surfactant concentration increases. The selection of optimum concentration of the surfactant is a balance between the reductions in the gas velocity at which liquid loading occurs compared to increase in the frictional loss as the concentration increases. We provide guidelines about the selection of the surfactant concentration.Visual observations using high speed camera show differences in the behavior under foam flow conditions. Unlike air-water flow, the liquid film attached to the wall is replaced by thick foam capturing the gas bubbles. The type of roll waves which carry the liquid in 2-inch pipe is different than what was observed in 4-inch pipe. Compared to 4-inch pipe, the roll waves in 2-inch pipe are much thicker. This partly explains the differences in 2-inch versus 4-inch pipe behavior.