三相流体的轴对称格子Boltzmann模型及其在Rayleigh-Plateau不稳定性的应用

基于多组分相场理论提出了一类模拟三相流体流动的轴对称格子Boltzmann模型.该模型利用两个粒子分布函数来捕捉三种不同流体之间的相界面,另一个粒子分布函数来求解流体动力学方程以获得流场信息.为了刻画坐标变换引起的轴对称效应,巧妙地设计了演化方程中平衡态分布函数和外力项分布函数,从理论上保证本文模型可以正确恢复三相流体系统的宏观控制方程,并且轴对称效应产生的源项中不包含任何复杂的梯度项,从而比现有的轴对称格子Boltzmann模型更加简单高效.首先通过模拟一系列轴对称多相流的基准算例,包括静态的双液滴、液体透镜的扩展和二元流体Rayleigh-Plateau不稳定性,来验证本文模型的有效性与正确性.接下来,利用该模型研究了三相流体的Rayleigh-Plateau不稳定性的增长过程,定量分析了波数和液柱半径比对复合液体线程破裂过程中界面动力学行为、界面破裂时间以及生成子液滴尺寸的影响.可以发现复合的液体线程在波数较大时破裂生成一个复合主液滴和卫星液滴,而在波数较小时可以生成更多数量的卫星液滴,这导致复合主液滴和卫星液滴的尺寸随着波数的增加呈现先增大而后减少的趋势.另外,我们发现内部流体...     

等直径微液滴碰撞过程的改进光滑粒子动力学模拟

运用一种改进光滑粒子动力学(SPH)方法模拟了相溶和不相溶两种情况下的等直径微液滴碰撞动力学过程. 为提高传统SPH方法的数值精度和稳定性, 采用一种不涉及核导数计算的核梯度改进形式; 为处理液滴界面张力采用修正的van der Waals表面张力模型. 通过模拟牛顿液滴碰撞聚并变形过程并与相关文献或试验结果进行对比, 验证了改进SPH 方法模拟微液滴碰撞过程的可靠性. 随后, 研究了基于van der Waals模型相溶聚合物微液滴碰撞聚并变形过程及不相溶微液滴碰撞后的反弹、分离过程, 讨论了碰撞过程中碰撞速度、碰撞角度、密度比等参数对碰撞变形过程的影响, 分析了流体桥、旋转角度等因素的变化情况. 关键词： 光滑粒子动力学 微液滴 聚合物液滴 碰撞     

气粒两相流传热问题的光滑离散颗粒流体动力学方法数值模拟

在气粒两相流动问题中,颗粒间以及气体与颗粒间的传热问题不可忽略.光滑离散颗粒流体动力学(SDPH)模型作为一种新的求解气粒两相流动问题的方法,已经成功应用于模拟风沙运动等问题.在此基础上,提出了SDPH方法的热传导模型,模拟了气粒两相流动问题中的热传导过程以及颗粒蒸发过程.首先引入各相的能量方程,利用有限差分与光滑粒子流体动力学一阶导数相结合的方法,处理各相内部热传导项中的二阶导数问题,基于气粒两相间温度差及对流换热系数计算颗粒与气体间的热传导量,推导得到了含热传导模型的气粒两相流SDPH计算方程组,模拟计算了圆盘形颗粒团算例及鼓泡流化床内部热传导算例,并与双流体模型计算结果进行对比,结果基本符合;其次利用离散液滴模型中的颗粒蒸发传质传热定律计算颗粒的蒸发过程,数值模拟了颗粒射流蒸发过程,并与离散颗粒模型结果进行对比,两者符合得较好,验证了该方法的准确性及实用性.     

基于摩擦-动力应力模型模拟喷动床气固流动

建立同时考虑颗粒和颗粒之间瞬间碰撞作用(瞬间接触)和非流化下颗粒滑动和滚动运动产生的颗粒半接触作用的颗粒碰撞-摩擦应力模型。推导具有普适性的固相动量守恒方程,建立高颗粒浓度气固两相流动模型,对喷动床内气体颗粒流动过程进行了数值模拟。模拟计算预测喷动床内喷射区、环隙区和喷泉区颗粒流动特性。预测喷动床内颗粒浓度和速度分布与他人实验结果相吻合。     

流动聚焦微通道中双重乳液乳化行为研究

本文基于VOF(Volume of Fluid)相界面追踪方法,建立了流动聚焦微通道中双乳液液滴形成的模型,研究了互不相溶的三种流体在微流控系统中形成双重液滴的典型物理过程,给出了双重液滴形成的典型工况,分析了流动参数和工质物性对液滴形成过程的影响。结果表明:双重乳液液滴形成过程可分为滴式与喷式两种模式;外流体流量增大时液滴形成模式由滴式转变为喷式,所形成的双重液滴尺寸减小;中间流体流量变化对液滴形成模式、所形成的双重乳液内液滴尺寸影响不明显,外液滴随中间流体流量增大而增大。     

液滴碰撞对液膜内汽泡运动与换热影响

采用数值模拟的方法,研究了沸腾雾化喷射过程中热壁面薄液膜层受到液滴碰撞扰动时液膜层内汽泡运动,相界面变化和由此引起的壁面换热特性.模拟结果显示汽泡生长初期相界面变化与液膜层内二次核化特征与文献结果吻合良好,汽泡生长后期相界面变化存在滞后.讨论了液滴下降速度.液滴直径与初始位置,多液滴碰撞对液膜层内流动与壁面换热的影响.     

氯气/BHP液滴化学反应流动一维数值模拟

为研究横向均匀液滴发生器的机理,对Cl₂/He气体横向流过BHP液滴场的物理化学行为建立一维化学反应流动模型.其中气流场按不可压缩假设近似为一维均匀流动;发生于液滴内的化学反应由四种物质的反应扩散方程描述;气相中各种物质浓度由对流扩散方程描述.液滴是这些物质的源或汇,两个模型方程通过此源项相联系.对模型方程作数值求解.结果显示,氯气利用率、单态氧产率及效率与文献所提供的实验结果较为接近,反映了所建模型具备模拟混合气体与BHP液滴反应的基本能力.     

油水两相分散流的相界面浓度输运方程研究

通过对连续相湍流导致的油滴聚合和破裂现象的合理模化,本文建立了适用于垂直上升管油水两相分散流的相界面浓度输运方程.其中对连续相湍流驱动油滴随机碰撞引起的聚合现象模化时,考虑了两油滴间连续相液膜排出过程中的凹面膜效应.而在描述油滴与湍流涡随机碰撞导致的破裂现象时,考虑了破裂过程所产生油滴的大小不等的因素.利用所建立的相界面浓度输运方程对沿实验段轴向的相界面浓度变化进行了预测,结果与实验数据比较表明吻合很好.     

偏晶合金液-液相分离的格子玻尔兹曼方法模拟

本文建立了一个模拟在弥散相液滴的扩散长大、碰撞凝并和Ostwald熟化等因素的作用下偏晶合金液-液相分离过程的二维格子玻尔兹曼方法 (lattice Boltzmann method, LBM) 模型.该模型结合了Shan-Chen的两相流模型和Qin的介观粒子相互作用势模型的优点,并在LB演化方程中引入了反映相变的源项.应用该模型模拟研究了偏晶合金液-液相分离过程中单液滴的生长、两液滴的合并和多液滴的生长规律.结果表明在两液相区中第二相单个液滴的生长是一个通过扩散从非平衡态到平衡态过渡的过程.两液滴合并 关键词： 偏晶合金 液-液相分离 格子玻尔兹曼方法     

基于大密度差多相流SPH方法的二维液滴碰撞数值模拟

该文结合了Ott提出的修正连续性方程和Adami改进的动量方程, 对空气中的液滴碰撞问题进行了二维数值模拟. 为有效提高计算精度, 推导了适用于大密度差多相流的人工黏性和人工应力方程. 通过表面张力作用下方形液滴自然变化和空气中两液滴互溶的算例, 验证了算法的有效性; 对不同韦伯数 (8.8, 19.8)、不同碰撞参数 (0, 0.5)下的液滴碰撞过程进行了数值模拟, 并与VOF方法对比,取得了较为一致的结果; 进一步计算多个韦伯数、多个碰撞参数下的液滴碰撞, 得到了空气中二维液滴碰撞结果分布图,与实验结果相符合. 结果表明, 该算法对于求解涉及大密度差多相流的液滴碰撞破碎问题十分有效,而且该方法容易拓展到三维, 从而为进一步模拟火箭发动机的二次雾化过程奠定了基础. 关键词： 光滑粒子流体动力学 大密度差 多相流 液滴碰撞     

倾斜均质表面上等径水滴聚合过程的可视化实验

对倾斜均匀表面上等直径水滴的聚合过程及特性进行了可视化实验研究,获得了水滴直径和表面倾角等参数对液滴聚合过程中液滴液桥半径、接触角和接触线变化特性的影响,分析了水滴聚合对其运动的影响.实验结果表明:表面倾角越大,下滑的临界半径越小;液滴的直径越大,液滴聚合后越容易下滑;液滴聚合可以加快液滴的运动,使下滑临界半径减小.     

倾斜均质表面上非等径液滴的聚合特性

对倾斜均匀表面上非等径液滴的聚合过程及特性进行了可视化实验研究,获得了液滴半径和表面倾角等参数对液滴聚合过程中液滴液桥半径、接触角和接触线变化特性的影响规律,进一步说明了倾斜表面上液滴聚合可以加快液滴的运动.     

Droplet transport and coalescence kinetics in emulsions subjected to acoustic fields

A method to aid the separation of the oil phase from aqueous emulsions using a low-intensity, resonant ultrasonic field has recently been developed. The density and compressibility difference between the dispersed and continuous phases within the emulsion results in a net force on the oil drops that pushes them toward the pressure antinodes of the standing-wave field, where coalescence subsequently occurs. A trajectory model is developed to predict the relative motion of drops subjected to the acoustic field. Such trajectories are sensitive to the physical properties and relative size of interacting drops, the initial configuration of the drops, and acoustic field parameters. Model predictions are validated by comparing experimentally observed trajectories with those predicted by the model. The modeling approach is then extended to determine the temporal evolution of the size of the region surrounding a target drop cleared by coalescence as a function of physical and acoustic field parameters. These results form the basis of a population balance model that attempts to track the size-evolution of a drop population coalescing under the influence of an acoustic field.     

水平均质表面上液滴聚合过程的可视化实验研究

对水平均匀表面上液滴的聚合过程及特性进行了可视化实验研究,获得了液滴半径和液滴物性等参数对液滴聚合过程中液滴液桥半径和接触角变化特性的影响规律。实验结果表明:液滴聚合中液桥半径和接触角都呈衰减振荡变化; 聚合前液滴半径越小,液桥半径振荡频率越大,振幅越小,振荡时间越短;液滴的粘度越大,液桥半径的振荡频率越小, 振幅越小,振荡时间越短;液滴聚合前的接触角明显大于聚合液滴静止后的接触角,其差值与固体界面状况和气、固、液物性相关。     

On the coalescence of sessile drops with miscible liquids

Sessile drops sitting on highly wettable solid substrates fuse in qualitatively different ways after contact, depending on the surface tension gradients between the mixing droplets. In early time evolution the drop coalescence can be fast or delayed (intermittent). In long time evolution a secondary drop formation can occur. We study numerically droplet dynamics during coalescence in two and three spatial dimensions, within a phase field approach. We discuss criteria to distinguish different coalescence regimes. A comparison with recent experiments will be done.     

Binary collision of a burning droplet and a non-burning droplet of xylene: Outcome regimes and flame dynamics

The binary collisions of a burning droplet and a non-burning droplet of xylene are experimentally investigated. The experimental parameters span an extensive range of Weber number and impact parameter, covering the collision outcome regimes of coalescence, reflexive separation, and stretching separation. A high-speed camera captures the temporal details of the collision process, involving flame spread, visible radiation, and flame distributions around droplets. For reflexive separation and stretching separation, the flame from the droplet spreads to the ligament, surrounding it during the interaction process, and then spreads around separated droplets and satellite droplets. Highly-interactive flames are formed in-between the droplets, with very sooty flames generated for most collisions. For the coalescence case, a swirling flame forms around the rotating coalesced droplet. For similar Weber numbers, visible flame radiation is compared for different collision regimes. The visible flame radiation changes more significantly for the reflexive and stretching separation cases than it does for the coalescence case. The change of the averaged visible flame radiation for reflexive separation and stretching separation is more than two times higher than that for coalescence. The map of three different collision regimes is plotted in the Weber number versus impact parameter domain and compared with available theoretical model predictions. Although the different outcomes of collision with the presence of flame can be well predicted by the model, using fluid properties determined by the averaged properties of the two droplets, the dynamics of the detailed processes involved in the collisions are very interesting and have strong implications on overall combustion behavior that go well beyond the mapped regimes.     

Viscous to inertial crossover in liquid drop coalescence

Using an electrical method and high-speed imaging, we probe drop coalescence down to 10?ns after the drops touch. By varying the liquid viscosity over two decades, we conclude that, at a sufficiently low approach velocity where deformation is not present, the drops coalesce with an unexpectedly late crossover time between a regime dominated by viscous and one dominated by inertial effects. We argue that the late crossover, not accounted for in the theory, can be explained by an appropriate choice of length scales present in the flow geometry.     

Propagation of drop coalescence in a two-dimensional emulsion: a route towards phase inversion

The phase inversion that undergoes an emulsion while being sheared is a sudden phenomenon that is still puzzling. In this Letter, we report an experimental investigation on propagative coalescence by using a microfluidic device where a calibrated two-dimensional emulsion is created and destabilized. The velocity of propagation and the probability of the coalescence are reported as a function of the size and the spatial distribution of the drops, respectively. We then discuss the efficiency of this novel scenario of phase inversion and suggest that inversion can be favored by the existence of a drop size distribution.     

Nucleation and growth on a superhydrophobic grooved surface

The growth dynamics of water drops condensed on a superhydrophobic geometrically patterned surface were studied. Drop size evolution at early and intermediate times is self-similar. Drop growth laws do not differ for a flat surface because of a reduction of both drop and substrate dimensionality. A striking observation is the instantaneous drying of the top surface of grooves at a point in time due to coalescence of the drops with a completely filled channel. At late times, only a few large drops grow connected to the channels, in a mixed Wenzel-penetration regime.