润湿性梯度驱动液滴运动的格子Boltzmann模拟

运用考虑了固体与液体间分子作用力的格子Boltzmann方法,数值研究了由于固液界面上表面张力梯度引起的Marangoni效应驱动的液滴运动.当表面张力梯度较小时,计算结果和前人的理论预测符合较好.而表面张力梯度较大时,由于液滴不变形和准平衡态等假设不再满足,理论预测的液滴运动速度高于数值模拟的结果.计算结果显示,在向亲水端运动过程中液滴内部出现旋涡结构,当润湿性梯度较大时,其前进速度和接触角随时间变化出现振荡. 关键词： 润湿性 格子Boltzmann方法 Marangoni效应 液滴     

润湿性梯度驱动液滴运动的格子Boltzmann模拟

运用考虑了固体与液体间分子作用力的格子Boltzmann方法,数值研究了由于固液界面上表面张力梯度引起的Marangoni效应驱动的液滴运动.当表面张力梯度较小时,计算结果和前人的理论预测符合较好.而表面张力梯度较大时,由于液滴不变形和准平衡态等假设不再满足,理论预测的液滴运动速度高于数值模拟的结果.计算结果显示,在向亲水端运动过程中液滴内部出现旋涡结构,当润湿性梯度较大时,其前进速度和接触角随时间变化出现振荡.     

池沸腾中气泡生长过程的格子Boltzmann方法模拟

在通过引入精确差分方法的单组分多相格子Boltzmann模型的基础上耦合能量方程,并考虑流体与固壁间的相互作用力来调节气泡与固壁间的接触角,从而建立了一种新的描述气液相变的格子Boltzmann理论模型. 为验证该模型的正确性,利用其对工质为水的相变过程进行了模拟,发现模拟结果与实验值符合良好;进而利用其验证Laplace定律,发现计算所得的水的表面张力与实验值甚为符合. 为考察该模型处理复杂相变问题的能力,利用其对工质为水的池沸腾中的气泡生长过程进行模拟,发现气泡脱离直径与g^{-0 关键词： 格子Boltzmann方法 池沸腾 气泡生长过程 接触角}     

楔形体上复合液滴润湿铺展行为的格子Boltzmann方法研究

固壁上液滴的润湿铺展行为是自然界中普遍存在的现象,针对楔形体上的复合液滴,采用基于相场理论的格子Boltzmann方法对其润湿铺展行为进行探究.通过理论分析和数值模拟,发现液滴润湿面积随接触角、楔形体顶角的减小而增大,液滴也越容易分裂.处于理论分裂临界状态附近的液滴,在一定密度比、黏度比条件下将沿楔形体壁面分裂.基于模拟结果生成以密度比、黏度比为坐标的液滴分裂状态相图,比较发现相同条件下初始状态为平衡态的复合液滴更不易发生分裂.另外模拟还表明非对称界面张力及非对称运动黏度比也是影响液滴分裂结果的重要因素.     

改性疏水固壁润湿性反转现象的格子Boltzmann方法模拟

基于疏水固壁改性会引起润湿性反转的特点,采用考虑固体与液体间分子力的格子Boltzmann方法,从壁面的线性和瞬时改性两方面对润湿性反转现象进行了数值模拟,并结合流体体积方法处理界面层质量.结果表明:壁面线性改性的过程中润湿性反转变化平稳,润湿所需时间大幅减少,所得到的接触角与固液吸引力系数的关系与其他文献结果一致;壁面瞬时改性幅度越大说明固壁对液滴作用力越强,表现为润湿性变化越明显,瞬时改性后接触角随时间呈指数规律变化,这与现有结论相符合.研究发现:在改性条件下液膜铺展过程中伴随着振荡变化,线性改性的振动峰值与改性幅度相关;瞬时改性的液膜速度会在某一时刻突然增大,这种现象与夹带空气有关.     

格子Boltzmann方法模拟双液滴撞击液膜的流动过程

采用单相自由面格子Boltzmann方法中流场计算,自由面条件和表面张力处理方法,研究双液滴撞击液膜的流动过程,实现对垂直间隔的两液滴相继撞击液膜流动过程的模拟.数值结果表明,液滴相继冲击液膜时,液滴间上下距离对液膜状态有很大的影响,液滴间距不同,液膜形状会出现很大差异.数值结果与实验结论定性一致.     

基于光滑粒子动力学方法的液滴冲击固壁面问题数值模拟

采用改进的光滑粒子动力学(SPH)方法对液滴冲击固壁面问题进行了数值模拟. 为了提高传统SPH方法的计算精度和数值稳定性, 在传统的SPH方法的基础上对粒子方法中的密度和核梯度进行了修正, 采用了考虑黎曼解法的SPH流体控制方程, 构造了一种新型的粒子间相互作用力(IIF)模型来模拟表面张力的影响. 应用改进的SPH方法对液滴冲击固壁面问题进行了数值模拟. 计算结果表明：新型的IIF 模型能够较好地模拟表面张力的影响, 改进的SPH方法能够精细地描述液滴与固壁面相互作用过程中液滴的内部压力场演变和自由面形态变化, 液滴的铺展因子随初始韦伯数的增大而增大, 数值模拟结果与实验得到的结果基本一致. 关键词： 液滴 固壁面 光滑粒子动力学 表面张力     

液滴在粗糙固体表面上的铺展特性

本文建立了粗糙固体表面上液滴铺展过程的理论模型并采用了格子Boltzmann方法进行数值模拟,研究了粗糙与光滑表面上液滴铺展过程的形貌变化。研究结果表明,无论表面光滑还是粗糙,液滴都沿着半径方向逐渐铺展,整体首先呈帽形,随着时间的推移,液滴在表面铺展成一平面。粗糙表面与光滑表面铺展特性间的差异主要在于粗糙表面液滴受到粗糙元干扰,其铺展形状随时间推移由圆形向八边形过渡.粗糙元间隙的增大导致液滴润湿半径增大,液滴厚度沿半径方向变化梯度降低。     

气泡成核过程的格子Boltzmann方法模拟

利用精确差分格子Boltzmann模型探讨水在特定温度下的亚稳态及不稳定平衡态, 获得等温相变过程中形成气泡和液滴的条件, 模型预测结果与理论解符合良好. 在该等温模型的基础上耦合能量方程, 通过调节流体-壁面相互作用力获得不同的气泡与固壁间接触角, 从而建立了一种新的描述气液相变的格子Boltzmann理论模型. 利用该新模型模拟不同流体-壁面相互作用力下凹坑气泡成核过程, 再现了气泡成核过程中的三阶段特性; 探讨了接触角、曲率半径及气泡体积随气泡成核过程的变化关系, 获得了与文献结果定性符合的曲率-气泡体积关系曲线. 关键词： 格子Boltzmann方法 气泡成核过程 气液相变 接触角     

液滴振荡模型及与数值模拟的对比

由于碰撞壁面后液滴内部流动的复杂性, 以及气-液-固三相间的相互作用, 对液滴碰撞壁面形态变化的数学理论研究有较大的难度, 因此所见者多为实验和数值模拟. 本文通过对液滴受力状态的分析, 得到了惯性力、黏性力和表面张力带经验系数的表达式, 并进一步建立了液滴碰撞壁面振荡模型, 得到了液滴铺展半径的振荡表达式, 以及表面张力、黏性系数等参数对液滴铺展的影响. 最后通过与液滴衰减振荡数值模拟结果的对比, 确定了液滴振荡模型中的修正系数, 验证了模型的可行性. 关键词： 液滴碰撞 振荡 铺展半径/高度 数值模拟     

剪切变稀液滴撞击不同浸润性壁面的数值模拟研究

基于有限元法,采用水平集方法捕捉相界面的移动,构建了液滴撞击固体壁面的数值模型.通过修正的幂律模型描述流体的非牛顿剪切变稀特性,探讨了剪切变稀特性对液滴撞击固体壁面后铺展行为的影响,分析了撞击不同浸润性壁面时剪切变稀特性对液滴撞击壁面行为的影响差异.研究结果表明:随着幂律指数m的减小,液滴撞击过程中的黏性耗散减小,液滴的形貌变化及无量纲参数变化更为显著.接触角为55°的情况下:当m降低至0.85时,液滴铺展过程中开始出现显著区别于牛顿流体液滴的振荡现象;当m降低至0.80时,液滴在回缩过程中会出现中心液膜断裂的情况.接触角为100°时,剪切变稀液滴均会出现振荡行为,振荡幅度随着m的减小而增大.接触角为160°时,牛顿流体液滴与剪切变稀液滴均会在回缩过程中弹起,但剪切变稀液滴的弹起速度更快.此外,基于数值计算结果,本文提出了接触角为55°情况下剪切变稀液滴撞击壁面后的最大无量纲铺展直径预测模型.     

Estimating critical surface tension from droplet spreading area

Critical surface tension (CST) is a measure of solid surface tension and is mainly determined by measuring the contact angle of a droplet on a target solid surface. The concept of CST makes it possible to determine solid surface tension without any unprovable assumptions such as the Fowkes hypothesis. However, it requires somewhat special devices and skills for measuring the contact angle. In this work, we propose a simple method to determine the CST of a solid by measuring the droplet spreading area. This method is developed by combining the conventional CST with a simple analytical droplet model. The difference in estimated CSTs between our method and the conventional one is within 3.0%. Our method enables a quick and simple evaluation of the solid surface tension without special devices for measuring the contact angle.     

液滴撞击圆柱内表面的数值研究

针对液滴撞击圆柱内表面的过程,利用基于相场的格子Boltzmann方法模拟液滴以不同初速度、从不同初始高度、撞击不同大小的圆柱内表面时液滴的形态变化,分析了液滴自身物性(如密度和黏性等)和圆柱内表面润湿性等因素对撞击现象的具体影响.研究发现:撞击韦伯数、密度比及动力黏性比、圆柱半径等对液滴撞击后沿圆柱内表面的铺展均有一定影响,较高的韦伯数下液滴可能会发生分裂;液滴初始高度对大密度比和动力黏性比的撞击影响较小;液滴反弹现象可能出现在接触角较大时;重力作用会抑制撞击后液滴的振荡.     

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

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

低温光滑壁面上水滴撞击结冰行为

采用高速摄像技术测试低温光滑壁面上水滴撞击结冰过程, 分析了撞击速度、壁面温度和材料热导率对水滴撞击铺展、振荡及结冰行为的影响规律. 结果表明, 低温壁面造成水滴最大铺展直径缩小, 且结冰时间随温度降低而缩短; 当撞击We数提高时, 水滴最大铺展直径增大, 而振荡和结冰时间减小; 同时材料热导率越高, 最大铺展直径越小, 结冰越迅速. 另外, 从热力学角度推导出水滴撞击结冰时间的理论公式, 预测误差<5.3%.     

基于范德瓦尔斯表面张力模式液滴撞击疏水壁面过程的研究

液滴撞击疏水壁面过程的研究在介观流体力学和微流体作用材料科学的研究中具有重要的理论意义和工程价值. 论文在SPH方法中引入范德瓦尔斯状态方程处理液滴表面张力, 考虑流体粒子之间远程吸引, 近程排斥的内部作用力, 提出了流体粒子与疏水壁面粒子间势能函数与表面张力相结合的作用模式. 通过模拟真空条件下两个静止的等体积液滴相互融合的过程, 验证了计算模式在模拟液滴的表面张力中的有效性. 采用该模式模拟的液滴撞击疏水壁面过程, 不仅能够有效地模拟液滴撞击壁面后的变形过程, 而且清晰地模拟出液滴的回弹、腾空以及二次撞壁现象的完整过程. 模拟结果与液滴撞击疏水壁面的实验结果以及VOF模拟结果符合较好, 表明本文所提出的表面张力和疏水壁面作用力处理模式对模拟液滴撞壁过程具有实际应用价值.