倾斜壁面浸润性梯度驱动液滴运动过程的格子Boltzmann模型研究

基于介观模型的多组分伪势格子Boltzmann方法,模拟了倾斜壁面浸润性梯度驱动液滴的运动过程,研究了壁面浸润性梯度、壁面倾斜角度对液滴运动过程的影响．结果表明,对于一定倾斜角度的壁面,当壁面上浸润性梯度足够大时,液滴能够克服重力的作用实现“爬坡”;液滴在运动过程中,其前进及后退接触角与当地静态接触角间存在差值;增大壁面浸润性梯度时,液滴能够获得更快的加速,并且前进及后退接触角与当地静态接触角之间的差值也随之增大;增大壁面倾斜角度时,液滴的运动受到阻碍,前进及后退接触角与当地静态接触角的差值小幅减小．     

Effect of the contact-line dynamics on the oscillations of a compressed droplet

This paper studies the natural and forced oscillations of a deformed droplet of an inviscid liquid surrounded by a different liquid and bounded in the axial direction by solid planes. In equilibrium, the droplet is a figure of revolution and the ratio of its radius to height is significant. The equilibrium contact angle between the side surface of the droplet and the solid surface is different from a right angle. The motion of the contact line is taken into account by setting an effective boundary condition. It is shown that three characteristic ranges of natural frequencies exist.     

基于SPH方法的湿润性固壁边界条件模拟研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。     

基于LBM的铝微滴斜柱沉积水平偏移研究

金属微滴沉积制造技术采用逐点堆砌方式成型, 为斜柱沉积提供无支撑制造方式, 具有高度灵活性. 本文针对铝液滴斜柱连续沉积过程, 建立格子玻尔兹曼模型进行数值模拟, 研究液滴在凝固表面上的水平偏移运动. 根据表面能充放过程, 沉积运动被划分为下落、快速扩张、慢速扩张、回弹4个阶段, 其受力状态由表面能、重力势能、动能和黏性耗散趋势得到. 液滴内部流动在扩张阶段中表现为滑动状态, 而在回弹阶段中表现为滚动状态. 液滴偏移运动的加速阶段主要发生在扩张阶段, 而偏移距离则在回弹阶段中产生. 扩张阶段的受力状态表明偏移运动的主要推动力为重力和毛细力. 随着液滴轴线距离的增大, 扩张阶段中的加速段时间缩短、速度峰值提高, 使水平偏移距离呈先增大后减小的趋势, 这种阶段化特征源于加速段时长和速度极大值的竞争关系. 不同沉积高度和固液浸润度下, 偏移距离均保持相同的演化趋势. 在相同的轴线距离下, 偏移距离随固液浸润度的增大而减小, 随沉积高度的增大而减小. 通过拟合水平偏移距离演化规律、优化扫描步距, 能够实现斜柱的均匀沉积, 并使倾角与理论结果一致.      

基于LS-DYNA的热成型钢断裂失效预测研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。     

STATIC SPREADING OF DROPLET IMPACT ON SOLID SURFACE：INFLUENCE FACTOR

通过建立液滴撞击固体平壁的静态铺展力学平衡的数学模型,从理论上得到了静态铺展半径与液滴物性参数、以及液滴与固体壁面接触角之间关系的数学表达式,将理论结果与数值模拟的结果进行了比较,两者吻合较好.比较了不同条件下液滴的静态铺展半径的变化规律,分别得到了液滴密度、体积、表面张力和接触角等因素对液滴静态铺展半径的影响规律.     

平壁液滴静态铺展影响因素的研究

通过建立液滴撞击固体平壁的静态铺展力学平衡的数学模型,从理论上得到了静态铺展半径与液滴物性参数、以及液滴与固体壁面接触角之间关系的数学表达式,将理论结果与数值模拟的结果进行了比较,两者吻合较好.比较了不同条件下液滴的静态铺展半径的变化规律,分别得到了液滴密度、体积、表面张力和接触角等因素对液滴静态铺展半径的影响规律.      

纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.      

A numerical study of droplet motion/departure on condensation of mixture vapor using lattice Boltzmann method

Droplet motion/departure, which is governed by external force acceleration coefficient, droplet radius and surface wettability on solid surfaces under external forces such as gravitational force, play a significant role in characterizing condensation heat transfer, especially when high fractional non-condensable gases (NCG) present. However, due to the challenge in visualizing the vapor/steam velocity field imposed by droplet motion/departure, the detailed mechanism of droplet motion/departure on condensing surfaces has not been completely investigated experimentally. In this study, droplet motion/departures on solid surfaces under external forces and their interactions with steam flow are simulated using two dimensional (2D) multiphase lattice Boltzmann method (LBM). Large external force acceleration coefficient, droplet radius and contact angle, lead to large droplet deformation and high motion/departure velocity, which significantly shortens the droplet residual time on the solid surface. Our simulation shows that steam vortices (lateral velocity) induced by droplet motion/departure can greatly disturb the vapor flow and would be intensified by increasing external force acceleration coefficient, droplet radius, and contact angle. In addition, the location of vortex center shifts in the ascending direction with increase of these factors. The average lateral velocities induced by droplet motion/departure at various conditions are obtained. The mass transfer resistance is substantially reduced owing to the droplet motion/departure, leading to an enhanced heat flux. The experimental results are compared to validate the influence of droplet motion/departure on condensation heat transfer performance, especially for steam–air mixture with the presence of high fractional NCG.     

Experimental investigation of inclined liquid water jet flow onto vertically located superhydrophobic surfaces

In this study, the behaviour of an inclined water jet, which is impinged onto hydrophobic and superhydrophobic surfaces, has been investigated experimentally. Water jet was impinged with different inclination angles (15°–45°) onto five different hydrophobic surfaces made of rough polymer, which were held vertically. The water contact angles on these surfaces were measured as 102°, 112°, 123°, 145° and 167° showing that the last surface was superhydrophobic. Two different nozzles with 1.75 and 4 mm in diameters were used to create the water jet. Water jet velocity was within the range of 0.5–5 m/s, thus the Weber number varied from 5 to 650 and Reynolds number from 500 to 8,000 during the experiments. Hydrophobic surfaces reflected the liquid jet depending on the surface contact angle, jet inclination angle and the Weber number. The variation of the reflection angle with the Weber number showed a maximum value for a constant jet angle. The maximum value of the reflection angle was nearly equal to half of the jet angle. It was determined that the viscous drag decreases as the contact angle of the hydrophobic surface increases. The drag force on the wall is reduced dramatically with superhydrophobic surfaces. The amount of reduction of the average shear stress on the wall was about 40%, when the contact angle of the surface was increased from 145° to 167°. The area of the spreading water layer decreased as the contact angle of the surface increased and as the jet inclination angle, Weber number and Reynolds number decreased.     

Experiments and analysis on self-motion behaviors of liquid droplets on gradient surfaces

A surface with surface energy gradient was fabricated by using chemical vapor deposition technology with dodecyltrichlorosilane (C₁₂H₂₅Cl₃Si), and its property was characterized by sessile drop method and Atomic Force Microscope scanning. Visualization experiments were carried out to investigate the motion behaviors of water and ethylene glycol droplets on horizontal and inclined gradient surfaces. And system free energy transition was analyzed to understand the mechanics of the droplet self-motion. The results show that the height and density of the silane molecules groups determined surface energy distribution on the surface. The liquid droplets were self-propelled to move horizontally or uphill from hydrophobic zone to hydrophilic zone on horizontal and inclined gradient surface. The motion process of the droplet experienced an accelerating stage and a creeping decelerating stage; the velocity and the displacement as well as the creeping frequency were proportional to the droplet size. The velocity of 2 ml water droplet reached 42 mm/s on the horizontal surface and 18 mm/s on the inclined surface, while that for ethylene glycol droplet reached 7 mm/s on the horizontal surface. The droplet motion was resulted from the energy transition among interfacial energy, kinetic energy, gravitational potential energy, and viscous dissipation energy. The interfacial energy released from deformation of the droplet is the main source for the motion.     

液滴撞击超疏水壁面反弹及破碎行为研究

液滴撞击壁面时,壁面亲水性对液滴撞击壁面后的变化历程具有重要的影响。利用相界面追踪的复合Level Set-VOF方法对液滴撞击超疏水壁面的运动进行了研究。研究结果表明,撞击速度较小时,液滴撞壁后发生反弹;撞击速度较大时,液滴撞壁后会发生破碎现象;初始粒径的增大和表面张力的减小,有利于液滴撞壁后产生铺展破碎现象;撞击角度对撞壁后的液滴行为具有较大的影响。通过数值模拟,给出了一定条件下液滴垂直及倾斜撞击超疏水壁面反弹及破碎的临界条件。     

Axisymmetric oscillations of a cylindrical droplet with a moving contact line

Forced oscillations of a cylindrical droplet of an inviscid liquid surrounded by another liquid and bounded in the axial direction by rigid planes are investigated. The system is affected by vibrations whose force is directed parallel to the axis of symmetry of the droplet. The velocity of motion of the contact line is proportional to the deviation of the contact angle from the value at which the droplet is in equilibrium. Linear and nonlinear oscillations are considered. The conditions of the occurrence of resonance are determined.     

液滴撞击具有浸润性分布的倾斜固壁铺展行为的格子Boltzmann研究

采用基于单组分多相伪势模型的格子Boltzmann方法,模拟了三维液滴撞击左右两侧浸润性不同的倾斜固壁的铺展过程,获得了液滴在壁面两侧的铺展因子、相对铺展宽度、相对高度和液滴运动速度随时间的变化情况,研究了壁面浸润性分布和壁面倾斜角度对液滴铺展过程的影响．结果表明,液滴在倾斜壁面的铺展过程受到重力和表面力的综合作用,重力影响液滴的铺展和沿壁面向下的滑动,壁面浸润性分布影响液滴向壁面亲水侧横向移动．     

不同壁面条件下液滴撞击铺展特性的模拟研究

液滴撞击不同润湿性壁面的传热流动问题在自然界和工业生产中广泛存在。研究采用CLSVOF方法,引入描述壁面润湿特性的动态接触角,并考虑液滴物性参数随温度的变化,建立液滴撞壁模型,模拟研究液滴撞击流动行为,通过与实验对比验证,确定模型有效性。在此基础上,对传热作用下考虑壁面润湿性的液滴撞击问题展开研究,探讨壁面传热作用对液滴撞击铺展特性的影响。研究表明,在撞击过程中,液滴先铺展后逐渐收缩,与静态接触角模型相比,采用动态接触角模型所得的液滴流动特性与实验结果更加吻合;随着接触角增大,液滴在撞壁初期不易铺展,随后则易于收缩;虽然固液传热作用会影响液滴铺展直径,但不改变液滴的运动趋势。     

Drift of an inclined plate counter oncoming waves

The dependence of the velocity of the motion of a tow with an inclined plate mounted in a wave water channel on the wave parameters, the submergence depth, and the angle of inclination and dimensions of the plate is experimentally investigated. The effect of tow motion counter to the waves is detected and theoretically justified. The free surface profiles for periodic waves above an inclined plate obtained using the elolutionary system of the Boussinesq approximation equations correspond to the measured ones. The pulse generated as a result of wave breakup is estimated.     

旋转射弹高速倾斜入水多相流场与弹道数值模拟

基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程,结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响,然后对不同入水角下倾斜入水过程进行分析,得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性,但会降低弹体侧向稳定性,使射弹受到的阻力系数、俯仰力矩系数变小;入水角越小,形成的空泡越不对称,由射弹运动状态的改变引起的空泡形态变化越明显,在超空泡航行阶段,弹体运动较稳定,不同角度下流体动力系数差别很小,当弹体下表面刺破空泡壁沾湿时,弹体运动状态发生较大变化,流体动力系数迅速增大,此时入水角度过小,弹体容易失稳;弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。     

Computational study of bouncing and non-bouncing droplets impacting on superhydrophobic surfaces

We numerically investigate bouncing and non-bouncing of droplets during isothermal impact on superhydrophobic surfaces. An in-house, experimentally validated, finite element method-based computational model is employed to simulate the droplet impact dynamics and transient fluid flow within the droplet. The liquid–gas interface is tracked accurately in Lagrangian framework with dynamic wetting boundary condition at three-phase contact line. The interplay of kinetic, surface and gravitational energies is investigated via systematic variation of impact velocity and equilibrium contact angle. The numerical simulations demonstrate that the droplet bounces off the surface if the total droplet energy at the instance of maximum recoiling exceeds the initial surface and gravitational energy, otherwise not. The non-bouncing droplet is characterized by the oscillations on the free surface due to competition between the kinetic and surface energy. The droplet dimensions and shapes obtained at different times by the simulations are compared with the respective measurements available in the literature. Comparisons show good agreement of numerical data with measurements, and the computational model is able to reconstruct the bouncing and non-bouncing of the droplet as seen in the measurements. The simulated internal flow helps to understand the impact dynamics as well as the interplay of the associated energies during the bouncing and non-bouncing. A regime map is proposed to predict the bouncing and non-bouncing on a superhydrophobic surface with an equilibrium contact angle of 155°, using data of 86 simulations and the measurements available in the literature. We discuss the validity of the computational model for the wetting transition from Cassie to Wenzel state on micro- and nanostructured superhydrophobic surfaces. We demonstrate that the numerical simulation can serve as an important tool to quantify the internal flow, if the simulated droplet shapes match the respective measurements utilizing high-speed photography.     

A height function based momentum balance model to simulate contact angle dynamics with hysteresis

We report a combined method to deal with the contact angle dynamics with hysteresis. The momentum balance model is applied to obtain the transient contact angle by balancing the inertia and the capillary force where the curvatures are estimated by the height function at the contact line. This integrated approach provides a great convenience that no need for a prior knowledge of the contact angle, and possesses the good property of the sharp interface approximation. In order to facilitate the wider use of the present method, we incorporate a dynamic contact line model to estimate the contact angle when the contact line starts to move. This combination plus the hysteresis region will take the ability to solve most of problems related to wetting with more physical sense. This proposed model is finally validated by the droplet equilibrium, spreading, and sliding tests.     

Direct numerical simulations of thermocapillary migration of a droplet attached to a solid wall

This paper is designated to gain further insight into the physical mechanisms of thermal droplet actuation on a wall through direct numerical simulation. Classical theory states that free droplets in a nonuniform temperature field always move towards the hot side. However, when attaching a droplet to a wall with a nonuniform temperature gradient, lubrication theory explains how such a droplet moves towards the colder side. This paper aims at further investigating and clarifying the physical mechanisms and acting forces in the environment of a nonuniform temperature field and offers some explanations. For the numerical simulations of a droplet attached to a wall with a linear temperature gradient and larger contact angles, the full Navier–Stokes equations and energy equation are solved in a Volume of Fluid framework. The solver is extended with a dynamic contact angle treatment and thoroughly validated. The droplet motion is studied both in two and three dimensions, where a movement towards the cold and the warm side can be observed. The forces acting in such a setting are identified and interpreted. A decomposition of the jump conditions shows that the tangential stress due to the temperature dependent surface tension alone would lead to a motion towards the cold side, whereas the normal component alone would move the droplet to the opposite direction. The differences between two- and three-dimensional simulations show that the problem at hand is clearly three-dimensional.