大液滴撞击不同润湿性小球的三维格子Boltzmann方法模拟

采用三维多松弛(MRT)Shan-Chen的格子玻尔兹曼方法(LBM)模拟液滴撞击不同润湿性固体小球(液滴直径与小球直径比大于1)的动态过程,在多种方法验证计算模型准确性的基础上,探究液滴韦伯数、小球表面润湿性对液滴撞击小球动态过程的影响,以及液滴撞击小球后形成片层的几何特征与球面润湿性的关系。结果表明：液滴撞击小球后将形成钟罩形的片层,撞击亲水小球后的片层则会出现凹痕。当小球表面润湿程度一定时,片层的最大底径随韦伯数增大而增加。撞击亲水小球时,片层形成后会完全闭合,韦伯数越大,且片层完全闭合的时刻越早,并在闭合过程中片层容易发生破裂。而撞击疏水小球时,韦伯数越大,则片层发生破裂的时刻越早;在片层打开过程中,破裂后的片层成环条状下落,即片层的几何特征受韦伯数和小球润湿程度的共同影响。     

含活性剂液膜去润湿演化的稳定性特征

针对含非溶性活性剂的液膜在固体基底上的去润湿过程,基于润滑理论建立了基态和扰动态下液膜厚度和表面活性剂浓度的演化模型,应用非模态理论分析了演化过程的稳定性特征,探讨了分子间力对液膜去润湿过程的影响. 研究表明,微扰动波的引入（k=1）有利于液膜去润湿过程的稳定进行,扰动能量逐渐衰减,然而,该效果随着扰动波数的增加而显著改变,k ≥ 2时,液膜演化的稳定性反而恶化,扰动能量被逐步放大,演化呈现出非稳定特征. 增大初始液膜厚度可以有效改善液膜流动的稳定性. 范德华力放大了液膜表面的微扰动,使得液膜演化的稳定性下降;相反,Born斥力和静电斥力具有增强去润湿稳定性的作用. 关键词： 活性剂 去润湿 分离压 稳定性     

界面剪切力作用下波状液膜流的水动力稳定性

液膜流的水动力稳定性作为保障其高效传热传质性能的重要因素之一,受多种因素的制约和影响. 当气液界面处存在因气流流动而产生剪切力作用时,剪切力将通过改变界面处的边界条件,从而影响液膜流动的稳定性. 基于边界层理论,采用积分法建立了剪切力作用下降液膜表面波演化方程,分析了界面剪切力对水动力稳定性的影响. 研究表明,正向剪切力为不稳定性因素,反向剪切力在较小雷诺数时为不稳定因素,在大雷诺数时为稳定性因素;正向剪切力使临界波数和临界波速增大,反向剪切力使其减小;剪切力对临界波速的影响在不同雷诺数下也有所不同.      

连续凹槽基底对含非溶性活性剂薄液膜流动特性的影响

针对连续凹槽基底上含非溶性活性剂液膜的流动过程, 采用润滑理论建立了液膜厚度和活性剂浓度演化模型, 利用PDECOL程序数值模拟得到了液膜流动的动力学特征及基底结构参数的影响. 研究表明: 活性剂液膜流经连续凹槽时, 负向台阶处形成凹陷, 正向台阶处形成隆起, 且随时间逐渐增大; 与平整基底相比, 连续凹槽下的活性剂液膜铺展速度加快; 基底凹槽的高度增加或斜度减小可加速液膜破断的可能性; 增大凹槽宽度可促进液膜流动; 减小斜度会使液膜进入第一凹槽前形成隆起特征; 重力在液膜的爬坡和下坡过程中具有相反的作用, 但均使得流动稳定性变差; 分子间作用力形成的结合压可加速液膜流动, 进而引发去湿润现象, 分离压则与之相反. 关键词： 活性剂液膜 非平整基底 铺展 分子间作用力     

薄液滴在润湿性受限轨道上的热毛细迁移特性

通过固体表面改性可对液滴热毛细迁移过程进行调控.基于润滑理论和滑移模型建立了均匀温度梯度作用下液滴在润湿性受限轨道上运动的数学模型,通过将基底划分成亲水区域和疏水区域构建了润湿性受限轨道.结合接触线动力学提出了三维液滴在不同方向上接触线移动速度的计算方法,得到了液滴热毛细迁移的演化历程,分析了轨道宽度和润湿性对液滴迁移特性的影响.研究表明:液滴主体受温度梯度作用由高温区向低温区迁移,液滴后缘在移动过程中与主体部分间形成一层薄液膜,即薄液膜拖尾.液滴在垂直于轨道方向上的铺展受到抑制,收缩到轨道边缘后保持定扎状态.前进接触线移动速度开始时迅速减小,后缓慢降低趋于平稳;前进接触线移动速度与轨道宽度呈负相关.垂直于轨道方向上的壁面润湿性限制导致的排挤作用,在初始的短暂时刻对液滴在轨道上的热毛细迁移具有加速作用.液滴前进接触线移动速度与轨道润湿性呈正相关.增强轨道润湿性使得后退接触线移动速度的初始值增大,但对其稳定值影响不大.相比于改变轨道润湿性,改变轨道宽度更易于调控液滴热毛细迁移过程.     

含活性剂液滴在倾斜粗糙壁面上的铺展稳定性

针对含非溶性活性剂液滴在倾斜粗糙壁面上的铺展过程, 应用润滑理论推导出基态和扰动态下液膜厚度和活性剂浓度的演化方程组, 基于非模态理论研究了液滴铺展的稳定性特征, 探讨了相关参数的影响及其内在机理. 研究表明: 液膜厚度和活性剂浓度扰动量均呈现双驼峰型变化, 且峰值位于液滴底部凹陷处; 随扰动波数k增加, 最易失稳区域由液滴底部右侧凹陷处移至左侧凹陷, 壁面结构的影响逐渐减弱, 液滴铺展历程趋于稳定; 增加壁面倾角θ 将导致液滴铺展不稳定性加剧, 增大壁面高度D和壁面波数k₀均导致液滴铺展稳定性先增强后减弱; 随毛细数C减小, 液滴铺展稳定性下降, 重力的影响逐渐突显, 扰动量最大值呈现先增大后减小的变化趋势.     

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

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

波纹基底上含不溶性活性剂液滴的铺展稳定性

针对波纹基底上含不溶性活性剂液滴的铺展历程, 采用润滑理论建立了液滴铺展数理模型, 推导出基态和扰动态下液膜厚度和活性剂浓度的演化方程组, 基于非模态稳定性理论分析了液滴铺展的稳定性及参数的影响规律. 研究表明: 扰动量在液滴中心及铺展前沿处很小, 在液膜最薄处达到最大值且活性剂浓度的负扰动现象比较明显; 扰动波数可增强液滴铺展稳定性, 但随扰动波数增加, 该稳定性逐渐下降甚至转变为不稳定. 增加 Marangoni数将导致液滴铺展不稳定性加剧; 增大基底高度具有增强液滴铺展稳定的作用, Peclet数和基底波数取适中值时有利于液滴铺展的稳定性. 关键词： 活性剂液滴 非平整基底 铺展 非模态稳定性     

平衡接触角对受热液滴在水平壁面上铺展特性的影响

壁面温度是影响壁面润湿性的重要外部条件. 为解决液滴铺展中三相接触线处应力集中问题, 已有研究多采用预置液膜假设, 但无法探究壁面温度对润湿性的影响. 本文针对受热液滴在固体壁面上的铺展过程, 基于润滑理论建立了演化模型, 通过数值模拟, 从平衡接触角角度分析了温度影响壁面润湿性及铺展过程的内部机理. 研究表明: 随温度梯度增大, 液滴所受Marangoni效应增强, 致使液滴向低温区的铺展速率加快; 铺展过程中, 位于高温区的接触线与液滴主体部分间形成一层薄液膜, 重力与热毛细力先后主导该区域的铺展; 当液-固或气-液界面张力对温度的敏感度高于另两个界面时, 低温区方向的平衡接触角不断增大, 使壁面润湿性恶化, 导致液滴铺展减慢; 而当气-固界面张力对温度的敏感度高于其他两个界面时, 低温区方向上的平衡接触角将减小, 由此改善壁面润湿性, 加快液滴铺展; 在温度影响壁面润湿性和液滴铺展过程中, 平衡接触角起关键作用.     

垂直排液过程不稳定性的数值模拟

为分析线框排液实验中液膜表面出现的不稳定现象及其成因,针对含有不溶性活性剂的线框液膜排液过程,模拟液膜底部的不稳定现象,分析Marangoni效应、膨胀黏性和扰动波数因素的影响。结果表明：底部扰动在排液开始比较剧烈,而后快速减弱,到排液后期又逐渐增强。排液开始的扰动是由初始扰动引起,而排液后期的扰动与活性剂分布有关。较弱的Marangoni效应可增强表面扰动,而较强的Marangoni效应则减弱底部扰动,使液膜呈刚性,发生表面逆流现象;较高的膨胀黏性减慢液膜排液进程,降低表面速度,且能抑制Marangoni效应引起的逆流现象;波数较大的扰动使液膜在排液初期的扰动变强,但对排液后期的稳定性不产生影响。