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

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

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

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

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

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

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

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

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

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

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

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

液滴在梯度微结构表面上的铺展动力学分析

本文通过改变肋柱宽度和间距, 构造了二级和多级梯度微结构表面, 采用格子-Boltzmann方法对液滴在两种梯度表面上的铺展过程进行了研究, 探析液滴运动的机理和调控方法. 结果表明, 在改变肋柱间距的二级梯度表面上, 当液滴处于Cassie态时, 接触角滞后大小与粗糙度梯度成正比关系; 当液滴从Cassie态转换为Wenzel态或介于两者之间的不稳定态时, 这一正比关系不再遵循. 在改变肋柱宽度的二级梯度表面上, 接触角滞后大小与粗糙度梯度始终成正比关系. 在多级梯度表面上, 随液滴初始半径增大, 接触角滞后减小, 但液滴平衡位置相较于初始位置偏离增大. 对梯度微结构表面上液滴运动和接触角滞后的定量分析, 可为实现梯度微结构表面液滴运动调控提供理论依据.     

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

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

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

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

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

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

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

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

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

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

格子Boltzmann方法模拟气泡生长脱离

格子Boltzmann方法具有微介观的特性,能够自动捕捉和追踪界面,在多相流领域有广阔的应用前景。对于伴有气液相变的传热问题,格子Boltzmann方法已有模型较少,其中有一部分是通过质量变化间接反映相变过程,不能真实反映温度场的演化和影响。本文采用了模拟多相流动和分离的自由能模型,结合基于Stefan边界的相变传热模型,考虑壁面自由能对流体与壁面间润湿性的影响,模拟了静态接触角、过热壁面的气泡生长和脱离,并进行了更接近真实沸腾过程的随机生成气泡的模拟,以及流动沸腾的气泡生长和脱离的模拟。     

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

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

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

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

辐射传输的格子Boltzmann模拟

辐射动力学理论是描述辐射传输是重要手段,基于此,本文建立了辐射能和辐射动量守恒方程,并基于Chapman-Enskog多尺度展开方法实现了从辐射传输Boltzmann方程到宏观方程的推导,进而建立了适于一维辐射传输的2分量格子Boltzmann模型。数值结果与精确解吻合较好,表明本文提出的LBM方法具有很好的准确性和稳...     

三维双气泡融合的格子Boltzmann模拟

针对基于自由能模型的格子Boltzmann方法,推导D3Q15格子模型对应的平衡态分布函数;采用该模型模拟双气泡的融合过程.结果表明,气泡的融合不仅与它们之间的初始间距有关,还与表面张力有关.表面张力越大,气泡融合的临界距离也越大.此外,研究气泡融合速度与初始间距、表面张力及粘性系数的关系.     

疏水表面减阻的格子Boltzmann方法数值模拟

采用格子Boltzmann方法的多松弛模型和Shan-Chen多相流模型对雷诺数为100的疏水表面方柱绕流进行了数值模拟, 分析了疏水表面接触角和来流含气率对方柱绕流流场的影响. 研究结果表明: 疏水表面接触角一定时, 来流含气率在一定范围内, 疏水表面具有减阻的能力, 超出这一范围时会出现阻力系数、升力系数升高的现象, 同时在方柱近壁面处伴随涡的形成产生了气团脱落; 当来流含气率处于适当水平时, 接触角越大, 绕流物体近壁面处含气率越稳定, 减阻效果越明显. 分析发现疏水表面减阻的关键在于保证近壁面处气层的稳定性, 此时接触角越大, 减阻效果越明显. 本文从含气率角度出发分析疏水表面的减阻现象, 为进一步探索疏水表面减阻机理提出了新的思路.     

格子Boltzmann守恒型非平衡态的外推边界

针对非平衡态外推边界处理方法在某些条件下发生质量泄漏的问题,结合质量守恒定律对其进行修正,建立一种新的质量守恒型的边界处理格式.其基本思想是在边界处定义一个虚拟密度来修正平衡态分布函数,使得系统质量守恒.以重力作用下的密度场及磁流体模拟问题为例,利用该格式进行分析,讨论质量泄露问题的本质及各种因素的影响.通过对三种典型的流动:二维Poiseuille流、Couette流、Womersley流的模拟,验证该边界格式的计算精度和稳定性.     

磁流体流动与传热的格子Boltzmann模拟

考虑外加磁场下磁流体中纳米磁性粒子所受的各种作用力,建立了用于模拟磁流体流动与传热特性的两相格子Boltzmann模型,模拟了外加不同方向梯度磁场下平板间磁流体的流动与传热过程,计算了磁流体与平板间对流换热的Nusselt数,分析了磁场梯度方向、大小对Nusselt数的影响．