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

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

液滴碰撞Janus颗粒球表面的行为特征

为研究液滴碰撞Janus颗粒(双亲性)球表面的独特行为特征,以粒径为5.0 mm铜球为材料制备了Janus颗粒,用直径为2.0 mm的液滴,在韦伯数(We)为2.7,10,20,30的测试情况下对Janus颗粒球表面进行了碰撞实验.结果表明:液滴碰撞Janus颗粒球表面后的运动可分为铺展、回缩、振荡和回弹4个过程.在不同We下,液滴碰撞Janus颗粒后的运动状态主要与表面润湿性相关,在Janus颗粒亲水侧表现为铺展特性且铺展系数γ随着时间t的增大而逐渐增大并趋于稳定;但在疏水侧,表现为回弹现象,铺展系数γ会出现类似"抛物线"形状;当液滴碰撞Janus颗粒球表面亲-疏水分界线时,液滴铺展和回弹同时发生.基于能量平衡和受力分析发现,液滴动能和表面能的互相转化是液滴铺展的关键,液滴会在重力、惯性力、表面张力、黏性力、接触力等力的综合作用下展现其独特的行为特征并最终达到平衡状态.     

Simulation solution for micro droplet impingementon a flat dry surface

This paper presents a computational fluid dynamics approach for micro droplet impacting on a flat dry surface. A two-phase flow approach is employed using FLUENT VOF multiphase model to calculate the flow distributions upon impact. The contact line velocity is tracked to calculate the dynamic contact angle through user defined function program. The study showed that the treatment of contact line velocity is crucial for the accurate prediction of droplet impacting on poor wettability surfaces. On the other hand, it has much less influence on the simulation of droplet impacting on good wettability surfaces. Good fit between simulation results and experimental data is obtained using this model.     

Structure of the surface microrelief of a droplet evaporating from a rough substrate as a possible cause of contact angle hysteresis

Based on the refraction images of a droplet evaporating on a rough substrate, we simultaneously observed the dynamics of its surface microrelief, contact angle, and contact line deformations along the entire perimeter of the contact line. This has led us conclude that the microrelief structure is directly related to the phenomenon of contact angle hysteresis and the jump-like pattern of contact line deformation. We suggest a possible mechanism for the occurrence of contact angle hysteresis during droplet evaporation and derive the relations that specify the range of possible contact angles at known microrelief parameters.     

An advected-field method for deformable entities under flow

We study dynamics of a deformable entity (such as a vesicles under hydrodynamical constraints). We show how the problem can be solved by means of Green's functions associated with the Stokes equations. A gauge-field invariant formulation makes the study of dynamics efficient. However, this procedure has its short-coming. For example, if the fluids are not Newtonian, then no Green's function is available in general. We introduce a new approach, the advected field one, which opens a new avenue of applications. For example, non-Newtonian entities can be handled without additional deal. In addition problems like budding, droplet break-up in suspensions, can naturally be treated without additional complication. We exemplify the method on vesicles filled by a fluid having a viscosity contrast with the external fluid, and submitted to a shear flow. We show that beyond a viscosity contrast (the internal fluid being more viscous), the vesicle undergoes a tumbling bifurcation, which has a saddle-node nature. This bifurcation is known for blood cells. Indeed red cells either align in a shear flow or tumble according to whether haematocrit concentration is high or low. Received 19 December 2001 / Received in final form 31 May 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: chaouqi.misbah@ujf-grenoble.fr     

Continuous droplet rebound on heated surfaces and its effects on heat transfer property: A lattice Boltzmann study

A thermal multiphase lattice Boltzmann(LB) model is used to study the behavior of droplet impact on hot surface and the relevant heat transfer properties.After validating the correctness of the codes through the D~2 law,the simulations of intrinsic contact angle and the temperature-dependent surface tension are performed.The LB model is then used to simulate the droplet impact on smooth and micro-hole heated surface.On the smooth surface,the impinging droplet is reluctant to rebound,unless the intrinsic wettability of the solid surface is fairly good.On the micro-hole surface,however,the micro-holes provide favorable sites for generating a high-pressure vapor cushion underneath the impinging droplet,which thereby facilitates the continuous droplet rebound.For the continuously rebounding droplet.The time evolution of volume and temperature display obvious oscillations.The achievable height of the rebounding droplet increases as the intrinsic wettability of the solid surface becomes better,and the maximum transient heat flux is found to be directly proportional to the droplet rebounding height.Within a certain time interval,the continuous rebounding behavior of the droplet is favorable for enhancing the total heat quantity/heat transfer efficiency,and the influence of intrinsic wettability on the total heat during droplet impingement is greater than that of the superheat.The LB simulations not only present different states of droplets on hot surfaces,but also guide the design of the micro-hole surface with desirable heat transfer properties.     

微液滴在超疏水表面的受迫振动及其接触线的固着-移动转变

利用高速摄影技术对超疏水表面液滴振动的动态行为进行观测,研究液滴在不同频率下的振动特性. 实验发现,液滴的共振频率满足Rayleigh方程,微液滴在超疏水表面具有自由液滴的振动性质. 在80–200 Hz的驱动频率范围内,接触线出现了明显的固着-移动现象,液滴的振动频率是驱动频率的一半,液滴振动时的形变较大. 当驱动频率大于200 Hz时,接触线基本固着,液滴的振动频率近似等于驱动频率,液滴共振时的形态边缘始终有节点存在. 分析表明,液滴对外界驱动的不同响应与接触线的振荡行为和变形程度密切相关 关键词： 超疏水表面 受迫振动 共振 接触线     

带凹槽的微通道中液滴运动数值模拟

运用改进的耗散粒子动力学方法模拟了液滴在由凹槽所构成的粗糙表面微通道内的运动行为.改进的耗散粒子动力学方法采用新近提出的一种短程排斥、长程吸引相互作用势能函数,从而可以模拟带有自由面的流体,如液滴等.模拟了新势能函数下液滴与固体壁面的静态接触角,并用2次多项式拟合了"接触角-a_wf/a_f"变化曲线.研究了液滴在带凹槽的微通道中运动时,微通道壁面浸润性、外场力、液滴温度对液滴流动特性的影响.研究表明壁面浸润性和外场力对液滴流动特性的影响较大,液滴温度对液滴流动特性的影响较小.研究结果对运用耗散粒子动力学方法模拟并分析微流体在复杂微通道的流动有一定的参考价值.     

A comparative study of droplet impact dynamics on a dual-scaled superhydrophobic surface and lotus leaf

The impact dynamics of water droplets on an artificial dual-scaled superhydrophobic surface was studied and compared with that of a lotus leaf with impact velocity V up to 3 m/s. The lower critical impact velocity for the bouncing of droplets was about 0.08 m/s on both surfaces. At relatively low impact velocities, regular rebound of droplets and air bubble trapping and flow jetting on both surfaces were observed as V was increased. For intermediate V, partial pinning and rebound of droplets were found on the artificial dual-scaled surface due to the penetration of the droplets into the micro- and nano-scale roughness. On the lotus leaf, however, the droplets bounced off with intensive vibrations instead of being partially pinned on the surface because of the irregular distribution of microbumps on the leaf. As the impact velocity was sufficiently high, droplet splashing occurred on both surfaces. The contact time and restitution coefficient of the impinging droplets were also measured and discussed.     

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

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

微粗糙结构表面液滴浸润特性的多体耗散粒子动力学研究

超疏水表面因其优异的自洁作用一直是表面科学领域关注的重点.本文使用多体耗散粒子动力学(many-body dissipative particle dynamics, MDPD)方法模拟研究了不同粗糙结构下液滴的浸润特性, 并与Cassie-Baxter理论进行了对比. 研究使用了一种具有长吸短斥作用的流固作用函数来获得不同的液滴浸润性, 并利用一种简洁的数值方法来测量接触角. 模拟结果表明本研究方法能够稳定地捕捉到液滴在粗糙表面的静态和动态特性. 模拟了粗糙结构对三相接触线运动的黏滞作用, 与物理实验结果相符合, 表明Cassie-Baxter理论在实际应用中尚存在一定局限性. 研究分析了在动态铺展过程中的能量转化关系, 并指出在低值表面容易引起液滴反弹.     

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

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

Effect of suspended particles on spreading of volatile droplet on solid substrate

Abstract  

We have been interested in behaviors of suspended particles in a volatile droplet placed on a smooth substrate. It is known that the particles gather and deposit in the vicinity of the macroscopic contact line of the droplet, which is generally called ‘coffee stain problem’. A convective flow induced by non-uniform evaporation through the interface brings suspended particles toward the pinned contact line in the drying droplet, which forms a ring stain. We have focused on the dynamics of the droplet with/without suspended particles spreading on the solid substrate and on the behaviors of particles in the evaporating droplet. Spreading process of the droplet is significantly affected by the suspended particles. We indicate flow patterns in the droplet, in which the flow exhibits a modal structure with a mode number in the azimuthal direction, and indicate particles depositions after the dryout of the droplet. Three-dimensional particle tracking velocimetry is applied to reconstruct such unique flow patterns in the spreading process of the droplet. Resultant patterns of the particles depositing on the substrate are introduced.     

Sliding behavior of water droplets on line-patterned hydrophobic surfaces

We prepared line-patterned hydrophobic surfaces using fluoroalkylsilane (FAS) and octadecyltrimethoxysilane (ODS) then investigated the effect of line direction on sliding behavior of water droplets by direct observation of the actual droplet motion during sliding. Water droplets slide down with a periodic large deformation of the contact line and sliding velocity fluctuation that occurred when they crossed over the 500-μm ODS line regions in FAS regions on a Si surface tilted at 35°. These behaviors are less marked for motion on a 100-μm line surface, or on lines oriented parallel to the slope direction. Smaller droplets slide down with greater displacement in the line direction on 500-μm line patterning when the lines were rotated at 13° in-plane for the slope direction. This sliding behavior depended on the droplet size and rotation angle, and is accountable by the balance between gravitational and retentive forces.     

Model for reversible colloidal gelation

We report a numerical study, covering a wide range of packing fraction Phi and temperature T, for a system of particles interacting via a square well potential supplemented by an additional constraint on the maximum number n(max) of bonded interactions. We show that, when n(max)<6, the liquid-gas coexistence region shrinks, giving access to regions of low Phi where dynamics can be followed down to low T without an intervening phase separation. We characterize these arrested states at low densities (gel states) in terms of structure and dynamical slowing down, pointing out features which are very different from the standard glassy states observed at high Phi values.     

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

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

Interface dynamics under nonequilibrium conditions: From a self-propelled droplet to dynamic pattern evolution

In this article, we describe the instability of a contact line under nonequilibrium conditions mainly based on the results of our recent studies. Two experimental examples are presented: the self-propelled motion of a liquid droplet and spontaneous dynamic pattern formation. For the self-propelled motion of a droplet, we introduce an experiment in which a droplet of aniline sitting on an aqueous layer moves spontaneously at an air-water interface. The spontaneous symmetry breaking of Marangoni-driven spreading causes regular motion. In a circular Petri dish, the droplet exhibits either beeline motion or circular motion. On the other hand, we show the emergence of a dynamic labyrinthine pattern caused by dewetting of a metastable thin film from the air-water interface. The contact line between the organic phase and the aqueous phase forms a unique spatio-temporal pattern characterized as a dynamic labyrinth. Motion of the contact line is controlled by diffusion processes. We propose a theoretical model to interpret essential aspects of the observed dynamic behavior.     

Non-Newtonian thin films with normal stresses: dynamics and spreading

The dynamics of a thin film on a horizontal solid substrate is investigated in the case of non-Newtonian fluids exhibiting normal stress differences, the rheology of which is strongly non-linear. Two coupled equations of evolution for the thickness of the film and the shear rate are proposed within the lubrication approximation. This framework is applied to the motion of an advancing contact line. The apparent dynamic contact angle is found to depend logarithmically on a lengthscale determined solely by the rheological properties of the fluid and the velocity of the contact line.     

液滴撞击加热壁面传热实验研究

本文采用高速摄像仪对水滴和乙醇液滴撞击加热壁面后的蒸发过程进行了实验观测, 分析了液滴撞击加热壁面后的蒸发特性参数. 实验中, 两种液体初始温度均为20 ℃, 不锈钢壁面初始温度范围为68-126℃. 水滴初始直径为2.07 mm, 撞击壁面时Weber 数为2-44; 乙醇液滴初始直径为1.64 mm, Weber数为3-88. 结果表明, 液滴受到重力、表面张力及流动性的影响, 在蒸发过程的大部分时间内, 水滴高度持续降低而接触直径几乎不变; 蒸发后期, 液滴发生回缩, 水滴的接触直径、高度和接触角出现振荡现象. 乙醇液滴的接触角随时间的增加呈现先减小随后保持不变的趋势, 而接触直径和高度则持续减小, 直到液滴完全蒸发. 液滴蒸发总时长与液体物性和壁面温度有关, 随壁面温度的升高而减小, 与液滴撞击壁面时的Weber 数无关. 同时, 随着壁面温度的升高, 液滴显热部分占总换热量的比重增大, 显热部分能量不可忽略, 本文实验条件下得到水滴的平均热流密度为0.014-0.110 W·mm^-2.