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

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

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

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

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

本文采用格子Boltzmann方法对液滴撞击液膜过程进行了研究, 主要考察了雷诺数(Re)、韦伯数(We)、相对液膜厚度 (h) 以及表面张力 (σ) 等物理参数对界面运动过程的影响. 首先, 随着Re数和We数的增加, 可以明显观察到液滴撞击液膜过程中形成的皇冠状水花以及卷吸现象; 当Re数较大时, 液体会发生飞溅, 由液体飞溅形成的小液滴则会继续下落, 并与液膜再次发生碰撞. 其次, 当相对液膜厚度较小时, 液滴撞击液膜并最终导致液膜断裂; 然而随着相对液膜厚度的增大, 尽管撞击过程溅起的液体会越来越多, 但是液膜并不会发生断裂. 再次, 随着表面张力的增大, 界面变形阻力增大, 撞击过程中产生的界面形变也逐渐减弱. 最后还发现皇冠(由液滴溅起形成)半径r 随时间满足r/(2R) ≈ α√Ut/(2R), 这一结果与已有结论是一致的.     

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

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

基于MRT-LB伪势模型的液滴撞击液膜数值研究

运用改进的格子玻尔兹曼（LB）伪势多松弛多相模型,研究单/双液滴撞击液膜的流动特性.考察单液滴在不同气液相密度比时撞击液膜的发展.随着密度比的减小,冠状水花顶端开始向内弯曲,且底部半径显著减小.在大密度比情况下研究双液滴撞击液膜.结果表明：双液滴撞击液膜有中心射流的形成;液滴水平间距的增大,延缓中心射流的出现,并降低初始中心射流的高度;随Re数的增加,中心射流的高度明显增大.     

液滴撞击不同湿润性固壁面上静止液滴的格子Boltzmann研究

采用单组分多相的伪势格子Boltzmann方法,在大小液滴粒径比为1.5的情况下,对大液滴竖直撞击壁面上静止小液滴的过程进行模拟,研究亲水与超疏水壁面上大液滴竖直碰撞小液滴的过程,得到液滴铺展因子和相对高度随时间的变化.结果表明:增大We数会使液滴的铺展因子增大,铺展直径变大,相对高度减小;并且随着We数的增加,在超疏...     

蒸汽剪切驱动液滴行为的格子Boltzmann模拟

滴状冷凝过程中,存在蒸汽流动对液滴的吹扫作用,液滴在蒸汽剪切作用下克服壁面黏附变形和运动,液滴运动速度越大,冷凝传热性能越高。但是液滴在蒸汽作用下变形和运动的细节还不清晰,蒸汽速度对液滴变形和运动的影响机理还不明确。本文采用自由能格子Boltzmann方法研究了在不同蒸汽速度剪切作用下,液滴在具有不同润湿性固体表面上的变形和运动过程,分析了蒸汽速度和接触角对液滴变形和运动的影响机制,结果显示随着蒸汽速度的增加,液滴变形越大,液滴在固体表面的运动速度越大,停留时间越短,有利于液滴的移除和表面更新,相同蒸汽速度的作用下,液滴在接触角大的固体表面上变形和运动速度越大,也有利于液滴的移除和表面更新。从而定性或半定量地揭示了蒸汽速度影响蒸汽滴状冷凝传热的物理机制。     

基于格子Boltzmann方法的液滴沿固壁铺展动态过程模拟

矿井喷雾降尘是利用水雾使粉尘润湿沉降的过程,考虑到固体与液体间分子作用力,本文采用格子Boltzmann方法对液滴沿固壁铺展的动力学行为进行了数值模拟,结果发现铺展直径及动态接触角随时间呈指数规律,确定了液滴表面张力与铺展最大直径间的关系,固壁润湿性对铺展最大速度值影响较大,这些与物理试验及文献结果符合良好. 进一步考察了疏水性强的固壁,发现当液滴表面张力足够小时,铺展接触角可以在90°以下,与理论公式符合. 研究发现铺展过程中伴随着振荡,且铺展到最大时液膜有回缩趋势. 关键词： 液滴 格子Boltzmann方法 铺展 数值模拟     

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

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

用格子Boltzmann方法数值模拟三维空化现象

空化是一种微观、瞬时、随机、多相的复杂现象,其过程中所产生的极端条件以及伴随的一系列空化效应,将对液流系统产生破坏性和建设性两方面的作用.采用基于Shan-Chen模型的单组分多相流格子Boltzmann方法对水体中的三维空化现象进行了数值模拟,研究了低压下水体中气核半径与空化现象的相互关系,成功再现了低压下水体中微小气核发展成气泡的过程,并进一步研究了水体依次流经低压区、高压区时空化产生、发展、溃灭的全过程.数值模拟结果和理论预测结果符合良好. 关键词： 单组分多相流 格子Boltzmann方法 三维空化     

低温均匀相变的格子Boltzmann方法模拟

文中基于气液两相接触时的传质机理和两相格子Boltzmann理论模型,提出相界面的密度数处于气液相密度数之间,对低温液滴相变进行了模拟研究。在低温状态下,对单个液滴进行了均匀蒸发现象模拟,分析了不同过热度对蒸发速率的影响,并模拟运动液滴的相变。结果表明,所采用的格子Boltzmann方法能够有效分析模拟低温介质的相变。     

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

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

激波的格子玻尔兹曼方法模拟

格子玻尔兹曼方法已被广泛应用于模拟流体,但主要应用于模拟不可压流体。应用二维正方格阵上的3速9点模型模拟激波的传播。考察了激波关系式、激波的宽度和结构。模拟结果与由Chapman-Enskog展开导出的流体方程的解相符。这说明格子玻尔兹曼方法能够模拟可压缩流体。     

润湿梯度作用下液滴在倾斜表面向上运动的格子Boltzmann模拟

采用基于Shan-Chen伪势模型的格子Boltzmann方法,对液滴在存在润湿梯度的倾斜表面上克服重力、自下而上运动的过程进行模拟。探究润湿梯度、液滴尺寸、Bond数以及表面倾斜角度对液滴运动的影响。计算结果表明：液滴在运动过程中,内部会出现沿斜面向上的速度矢量,润湿梯度越大,液滴运动速度越快,润湿长度也越长,且动态接触角减小速率越快。液滴尺寸和Bond数对液滴运动的影响较小,但存在临界Bond数,超过该临界Bond数时,液滴将沿梯度润湿表面向下运动。表面倾角对液滴运动有显著影响,倾角增大,液滴运动速度和润湿长度都明显减小。     

格子Boltzmann方法中三维运动边界的统一模型

格子Boltzmann方法(LBM)中边界条件的处理很复杂,在现有的边界条件处理方法中,动力学格式能够精确满足宏观边界条件,但由于要解一个不定方程,必须引入附加假设确保方程非奇异.作为动力学格式和反弹格式的一种扩展,提出一种处理三维任意速度运动边界的统一模型,其中人口速度和固体壁面速度是该模型的特殊情形.给出用于三维15速度的表达式.为了检验该模型,模拟对角顶盖驱动三维空腔流,并将结果与有限差分法计算的结果进行比较,说明所提出的统一模型是合理可行的.     

基于格子Boltzmann方法的平板射流大涡模拟

应用多GPU技术,将格子Boltzmann方法与大涡模拟相结合(LBM-LES),使用1.12×10⁸网格,对雷诺数Re=4 000,倾斜角α=30°,吹风比M=0.5工况下的平板单孔射流进行了大规模高性能数值模拟研究.合理的定性与定量结果验证了LBM-LES模拟平板射流的有效性与可行性.使用上亿的计算网格捕捉了精细的湍流拟序结构,有利于主流与射流之间的掺混机理研究.此外,使用6个K20M GPU并行计算,模拟了71 680 LBM时间步长,仅耗时15 402秒,计算性能达到521.24MLUPS,即每秒更新5.212 4×10⁸个网格点的数据.     

纳米液滴撞击柱状固体表面动态行为的分子动力学模拟

液滴撞击固体表面是一种广泛存在于工农业生产中的现象.随着微纳技术的发展,纳米液滴撞击行为的定量描述有待完善.采用分子动力学模拟纳米水滴撞击柱状粗糙铜固体表面的动态行为.分别在液滴速度为2—15?/ps,五种方柱高度和六种固体表面特征能的情况下分析液滴的动态特征.结果表明,随着液滴初始速度V0的增加,其最终稳定状态先由C...     

液滴撞击等温固体平壁的数值模拟

建立了液滴撞击等温固体平壁的物理和数学模型,数值模拟了撞击后液滴的变形和反弹过程,将模拟结果与文献中的实验结果做了比较,两者吻合较好,表明所采用的模型可以模拟液滴在水平壁面上的运动形态.通过比较不同条件下铺展因素随时间的变化规律,分析了液滴碰撞速度、初始接触角对变形过程的影响.     

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

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