液滴在不同润湿性表面上蒸发时的动力学特性

基于润滑理论,采用滑移边界条件建立了二维液滴厚度的演化模型和移动接触线动力学模型,利用数值计算方法模拟了均匀加热基底上固着液滴蒸发时的动力学特性,分析了液-气、固-气和液-固界面张力温度敏感性对壁面润湿性和液滴动态特性的影响.结果表明,液滴的运动过程受毛细力、重力、热毛细力和蒸发的影响,重力对液滴铺展起促进作用,而毛细力、热毛细力则起抑制作用;通过改变界面张力温度敏感性系数,可使液滴蒸发过程中的接触线呈现处于钉扎或部分钉扎模式,且接触线钉扎模式下的液滴存续时间低于部分钉扎模式;提高液-气与液-固界面张力温度敏感系数均可改善壁面润湿性能,加快液滴铺展速率;而增大固-气界面张力温度敏感系数则导致壁面润湿性能恶化、延缓液滴铺展过程;通过改变固-气界面张力温度敏感系数更有利于调控处于蒸发状态下的液滴运动.     

液滴在气固交界面变形移动问题的光滑粒子流体动力学模拟

为准确模拟液滴在气固交界面变形移动问题,对基于连续表面张力模型的表面张力光滑粒子流体动力学方法进行了改进.改进方法采用新的边界处理方式和界面法向修正方法,即将固体边界虚粒子色函数值根据液面的位置进行相应设定以保证气-液-固三相交界处流体粒子的界面法向沿接触线法线方向,引入Brackbill提出的壁面附着力边界条件处理方法,对在气-液-固三相交界处的流体粒子及部分固体边界虚粒子的界面法向进行修正,修正前后保持法向模值不变,得到了含壁面附着力边界条件的表面张力算法.模拟了受壁面附着力影响的水槽中液面的变化过程、液滴润湿壁面过程和剪切气流驱动液滴在固体表面变形脱落过程,并与流体体积函数方法进行了对比.结果表明,该方法在处理壁面附着力问题时精度较高,稳定性较好,适合处理工程中液滴在气固交界面变形移动问题.     

前进接触线薄液膜结构与运移机制

气-液-固三相接触线薄液膜微观结构认识不清使得动态接触角预测问题始终存在争议.本文通过大规模分子动力学模拟研究液滴动态润湿过程,获得了完全浸润液滴前进接触线区域的前驱液膜结构、部分浸润液滴前进接触线区域的纳米级凸起结构.首次在模拟中获得了微观接触角与宏观接触角随时间的动态演变规律.分子轨迹追踪揭示了接触线薄液膜区域流体在固体表面摩擦作用下,随速度增大从滑动为主向滚动为主的模式转变,进而在底层出现气体裹挟导致滑移、发生溅射的机制.本研究为进一步构建普适的动态接触角预测模型提供了关键物理依据.     

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

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

悬垂液滴研究及表面张力和润湿角测定

对光滑固体表面下悬挂的液滴进行了理论分析,建立了悬垂液滴特征尺寸R和H与液滴表面张力σL和固液界面润湿角θ之间的关系式,计算发现对于特定的ρ,σL,θ值,液滴质量m与固液界面润湿半径R、液滴高度H满足特殊的曲线关系.利用此关系可以同时测量液滴的表面张力σL以及固液界面间的润湿角θ.     

神奇的表面浸润

<正>1.自然界中的浸润现象当把一种液体沉积于固体表面上时,由于液体、气体以及固体分子间的作用力,液体最后会处于平衡状态,这时气、液、固三相交点处气、液界面的切线会与固液界面线成一定角度的夹角,这种现象称为浸润,形成的角度称为接触角。对于水这种常见液体,一般人们把水滴接触角小于90°的表面称为亲水表面,接触角大于90°的表     

固着液滴蒸发薄液膜区传质传热特性数值研究

本文基于扩展Young-Laplace方程和动力学理论,以水为工质,建立固着液滴蒸发薄液膜区、宏观弯月面区、固体基底内部及气液、液固、气固界面的耦合传热模型,并探究固体基底厚度和热物性对蒸发薄液膜区传热特性的影响。结果表明,由于液滴蒸发过程中的冷却效应,三相接触线附近的基底温度低于周围固体温度;蒸发速度越快,温度差越大。蒸发薄液膜的长度、最大质量流量和传热量与基底厚度与液滴半径比值h_r有关;当h_r=0.4时,薄液膜区传热对总体传热贡献最大;基底的热导率越小,蒸发液滴具有更长的蒸发薄液膜长度及更强的界面传热,薄液膜区传热对总体传热影响越显著.     

底板属性对液滴蒸发过程的影响

液滴蒸发机制目前引起人们广泛重视。通过研究纳米颗粒质量溶度为0.05%的Al_2O_3液滴在铁表面和玻璃表面的蒸发过程,发现液滴在铁表面蒸发时呈现4种蒸发形态;定接触线形态、定接触角形态、黏滑形态和不稳定形态;而在玻璃表面的蒸发只有不稳定形态。主要原因是液滴在铁表面蒸发时,固-液表面张力更大,形成定接触线蒸发形态;随着蒸发进行,液滴温度下降,气-液表面张力变大,因而产生黏滑形态。     

液气共存的耗散粒子动力学模拟

传统的耗散粒子动力学方法(DPD)由于采用了纯排斥的守恒力相互作用,从而不能适应液气共存或者带有自由面流体的模拟.这里研究了DPD方法中新近提出的一种短程排斥、长程吸引相互作用,探索了这种改进势能对于DPD方法模拟液气共存的能力.模拟了这种新势能所形成的液气过渡界面,计算了过渡界面区的应力分布,发现应力分布与多体DPD方法所得结果一致.进一步对表面张力进行了研究,验证了这种势能所形成的界面满足Laplace定律,而通过理论公式与Laplace定律分别所得到的表面张力也彼此相符。     

凝固前沿和气泡相互作用的大密度比格子玻尔兹曼方法模拟

建立了二维双组分两相流的大密度比格子玻尔兹曼方法 (lattice Boltzmann method, LBM)模型. 该模型基于改进的Shan-Chen伪势多相流LBM模型, 结合采用不同时间步长的方法, 实现密度比达800以上的气液两相流模拟. 为了对模型进行验证, 模拟了在不同气液相互作用系数和密度比条件下气泡内外压力差与其半径之间的关系, 其结果满足Laplace定律. 将所建立的大密度比LBM与介观尺度的元胞自动机(cellular automaton, CA)和有限差分法(FDM)相耦合, 用LBM模拟气液两相流, 用CA方法模拟固相生长, 用有限差分法模拟温度场, 采用LBM-CA-FDM耦合模型对定向凝固过程中凝固前沿的气泡与液-固界面之间的相互作用进行模拟研究. 结果表明, 绝热气泡的存在影响了温度场分布, 使得凝固前沿接近气泡时, 液-固界面凸起, 在不同的固相生长速度条件下, 出现凝固前沿淹没气泡或气泡脱离凝固前沿的不同情况, 模拟结果与实验结果符合良好. 关键词： 格子玻尔兹曼方法 元胞自动机 凝固 气泡     

Nonlinear shear wave interaction at a frictional interface: energy dissipation and generation of harmonics

Analytical and numerical modeling of the nonlinear interaction of shear wave with a frictional interface is presented. The system studied is composed of two homogeneous and isotropic elastic solids, brought into frictional contact by remote normal compression. A shear wave, either time harmonic or a narrow band pulse, is incident normal to the interface and propagates through the contact. Two friction laws are considered and the influence on interface behavior is investigated: Coulomb's law with a constant friction coefficient and a slip-weakening friction law which involves static and dynamic friction coefficients. The relationship between the nonlinear harmonics and the dissipated energy, and the dependence on the contact dynamics (friction law, sliding, and tangential stress) and on the normal contact stress are examined in detail. The analytical and numerical results indicate universal type laws for the amplitude of the higher harmonics and for the dissipated energy, properly non-dimensionalized in terms of the pre-stress, the friction coefficient and the incident amplitude. The results suggest that measurements of higher harmonics can be used to quantify friction and dissipation effects of a sliding interface.     

Departures from Amonton's law of friction

The method of transition from light friction (rolling) to heavy friction (sliding) is used for analyzing the principal departure from Amonton's law (force of sliding friction T proportional to the normal pressure force N) for several metal — nonmetal pairs when the sliding time t approaches zero. It appears that the departure mode is, under these conditions, the same for all pairs which have been tested here. Experimental and theoretical evidence is given to prove that all other modes of departure from Amonton's law represent consequences of changes in the rubbing surfaces which occur during sliding and that, therefore, such departures will appear only after a long enough period of sliding.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 11, pp. 94–98, November, 1974.     

石墨烯层间纳米摩擦性质的第一性原理研究

基于密度泛函理论的第一性原理计算方法研究了纳米尺度下石墨烯层间摩擦现象, 探讨了对称和非对称两种情况下双层石墨烯层间沿不同方向的摩擦性质. 研究发现对于对称的双层石墨烯, 层间摩擦沿不同方向同性; 摩擦因数依赖于正压力, 随正压力增大, 摩擦因数的变化曲线分为三个阶段, 在较小以及较大压力下, 摩擦因数遵循Amonton法则不随压力变化而变化; 而在中间3-6 nN阶段, 摩擦因数随压力增加线性增加. 整个研究压力范围内摩擦因数在0.05-0.25之间. 对于非对称性双层石墨烯层间摩擦, 不同压力下摩擦因数在0.006上下波动, 摩擦因数较两层对称性石墨烯大大降低. 上述研究结果与实验一致.     

Fabrication of durable hydrophobic surfaces through surface texturing

Low surface energy polymer thin-films can be applied to surfaces to increase hydrophobicity and reduce friction for a variety of applications. However, wear of these thin films, resulting from repetitive rubbing against another surface, is of great concern. In this study, we show that highly hydrophobic surfaces with persistent abrasion resistance can be fabricated by depositing fluorinated carbon thin films on sandblasted glass surfaces. In our study, fluorinated carbon thin films were deposited on sandblasted and as-received smooth glass using deep reactive ion etching equipment by only activating the passivation step. The surfaces of the samples were then rubbed with FibrMet abrasive papers in a reciprocating motion using an automatic friction abrasion analyzer. During the rubbing, the static and kinetic friction forces were also measured. The surface wetting properties were then characterized using a video-based contact angle measuring system to determine the changes in water contact angle as a result of rubbing. Assessment of the wear properties of the thin films was based on the changes in the water contact angles of the coated surfaces after repetitive rubbing. It was found that, for sandblasted glass coated with fluorinated carbon film, the water contact angle remained constant throughout the entire rubbing process, contrary to the smooth glass coated with fluorinated carbon film which showed a drastic decrease in water contact angle with the increasing number of rubbing cycles. In addition, the static and kinetic friction coefficients of the sandblasted glass were also much lower than those of the smooth glass.     

Friction Properties of Bio-mimetic Nano-fibrillar Arrays

Nano-fibrillar arrays are fabricated using polystyrene materials. The average diameter of each fiber is about 300 nm. Experiments show that such a fibrillar surface possesses a relatively hydrophobie feature with a water contact angle of 142°. Nanoscale friction properties are mainly focused on. It is found that the friction force of polystyrene nano-fibrillar surfaces is obviously enhanced in contrast to polystyrene smooth surfaces. The apparent coefficient of friction increases with the applied load, but is independent of the scanning speed. An interesting observation is that the friction force increases almost linearly with the real contact area, which abides by the fundamental Bowden-Tabor law of nano-seale friction.     

Interface dynamics of lipid membrane spreading on solid surfaces

As a model system for two-dimensional interface dynamics, we study the wetting front of a lipid membrane moving over a solid substrate that is structured with regularly spaced pinning centers. By analyzing the contour of the front, we derive the normal growth rate and the relaxation coefficient. Both exhibit a 1/t(1/2) time dependence. Moreover, the friction coefficient and the line tension can be determined. Randomly distributed pinning centers cause a fractal contour line, whereas on surfaces that are artificially roughened, self-affine contour lines are observed. The latter exhibit an anomalous roughness exponent of zeta = 0.81+/-0.05.     

A model for the characterization of friction contacts in turbine blades

Stresses produced by the forced vibrations can lead to a significant reduction of the life of turbo engine blades. To predict the vibration amplitudes of this components an accurate dynamic analysis is necessary. The forced response calculation of these dynamic systems is strongly affected by the presence of the contact interfaces (i.e., underplatform dampers, shrouds, root joints). Different contact models are available in literature. These models make use of contact parameters, contact stiffness and friction coefficient to evaluate the damping and stiffness related to the contact interfaces. In this paper a model is proposed to characterize friction contact of non-spherical contact geometries obeying the Coulomb friction law with constant friction coefficient and constant normal load. The hysteresis curves of the oscillating tangential contact forces vs. relative tangential displacements and the dissipated energy at the contact are obtained for different contact geometries. The developed model is suitable to be implemented in numerical solvers for the calculation of the forced response of turbine blades with embedded friction contacts.     

Phase field modeling of hysteresis in sessile drops

We propose a novel approach to describe wetting of plane solid surfaces by liquid drops. A two-dimensional nonconserved phase field variable is employed to distinguish between wetted and nonwetted regions on the surface. The imbalance in the Young's force provides for the exchange of relative stability of the two phases. The three-phase contact line tension arises from the gradient energy and contact angle hysteresis from the kinetic coefficient. Using this theory, we discuss contact angle hysteresis on chemically heterogeneous surfaces. We show significant departure from the classical Cassie theory, which is attributed to defect pinning of the continuous triple line.     

Nonlinear dynamics of a rub-impact micro-rotor system with scale-dependent friction model

Scale effects in dry friction at microscale and the coefficients of friction due to adhesion force and two- and three-body deformations are considered. A rub-impact micro-rotor model with scaling nonlinear rub-impact force is presented and the nonlinear dynamic characteristics of the system in micro-electro-mechanical systems (MEMS) are investigated when the rotating speed, imbalance, damping coefficient, scale length and fractal dimension are regarded as the control parameters. Effects of scale length, fractal dimension, velocity-dependent impact factor and contact form on the coefficients of dry friction are investigated and discussed, and used to study the nonlinear behavior of rub-related vibrations with a large number of numerical simulations. The effects of rotating speed, imbalance, damping coefficient, and friction coefficient on the micro-rotor system responses are studied. It is indicated that the rub-impact micro-rotor system with the scale effects in friction alternates among the periodic, quasi-periodic and chaotic motions as the system parameters change. The results can be effectively used to diagnose the rub-impact fault, reduce the failure and improve the characteristics of a micro-rotor system, and optimize the design of micro-rotating machinery in MEMS.     

固体表面液滴铺展与润湿接触线的移动分析

液滴在固体表面上的铺展行为与润湿特性对许多工业生产过程的研究具有重要意义.根据液滴在光滑表面上的受力情况,建立了液滴平壁铺展的动力学模型.应用润滑近似方法和二维Navier-Stokes方程,建立了液滴沿理想表面铺展的动量和连续性方程.根据建立的方程,应用数值解法求解并详细分析了液滴在铺展过程中膜厚、接触线铺展半径以及铺展速度随时间的变化关系.研究结果表明:液滴的铺展过程可分为扩展和收缩两个阶段,铺展过程伴随着表面能、动能以及各种势能的相互转化,液滴最终的铺展半径大小由固体基面固有的润湿特性所决定;液滴在铺展过程中出现的"坍塌效应"与弯曲液面处的Laplace压力差有关;铺展半径随时间变化的标定律近似满足"1/7"次方标度律.