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

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

Electrowetting-induced oil film entrapment and instability

We investigate the spreading at variable rate of a water drop on a smooth hydrophobic substrate in an ambient oil bath driven by electrowetting. We find that a thin film of oil is entrapped under the drop. Its thickness is described by an extension of the Landau-Levich law of dip coating that includes the electrostatic pressure contribution. Once trapped, the thin film becomes unstable under the competing effects of the electrostatic pressure and surface tension and dewets into microscopic droplets, in agreement with a linear stability analysis. Our results recommend electrowetting as an efficient experimental approach to the fundamental problem of dynamic wetting in the presence of a tunable substrate-liquid interaction.     

非晶态合金表面的水润湿动力学

采用分子动力学模拟方法研究了改进的Simple pointcharge模型SPC/E水滴在Cu_(50)Zr_(50)非晶薄膜上的润湿行为和铺展过程.通过与CuZr(110)和(100)晶面对比研究发现,水滴在Cu_(50)Zr_(50)非晶薄膜表面上表现出较高的铺展速度.水滴在非晶合金表面的铺展过程中形成了明显的吸附层;而在晶态表面,水滴铺展前沿呈脚状形态.分析结果表明非晶表面的水分子在吸附层内呈现完全无序的单层排列方式,而在晶态表面,特别是(100)晶面,吸附层水分子呈双层有序排列.这种吸附层结构的差异导致了吸附层内水分子方向的差异:非晶表面吸附层内水分子方向倾向平行于表面,而晶态基底上吸附层内的水分子倾向于指向液滴内部.由此造成了非晶表面吸附层中的水分子与液滴内部以相对较弱的氢键相互作用,这使得上层水分子能够较容易扩散至吸附层前沿,促进液滴铺展.     

Spreading of diblock copolymer droplets: A probe of polymer micro-rheology

We present an experimental study of the spreading dynamics of symmetric diblock copolymer droplets above and below the order-disorder transition. Disordered diblock droplets are found to spread as a homopolymer and follow Tanner’s law (the radius grows as R ∼ t ^m , where t is time and m = 1/10 . However, droplets that are in the ordered phase are found to be frustrated by the imposed lamellar microstructure. This frustration is likely at the root of the observed deviation from Tanner’s law: droplet spreading has a much slower power law ( m ∼ 0.05±0.01 . We show that the different spreading dynamics can be reconciled with conventional theory if a strain-rate-dependent viscosity is taken into account.     

Topography driven spreading

Roughening a hydrophobic surface enhances its nonwetting properties into superhydrophobicity. For liquids other than water, roughness can induce a complete rollup of a droplet. However, topographic effects can also enhance partial wetting by a given liquid into complete wetting to create superwetting. In this work, a model system of spreading droplets of a nonvolatile liquid on surfaces having lithographically produced pillars is used to show that superwetting also modifies the dynamics of spreading. The edge speed-dynamic contact angle relation is shown to obey a simple power law, and such power laws are shown to apply to naturally occurring surfaces.     

微液滴在不同能量表面上润湿状态的分子动力学模拟

微小液滴在不同能量表面上的润湿状态对于准确预测非均相核化速率和揭示界面效应影响液滴增长微观机理具有重要意义. 通过分子动力学模拟, 研究了纳米级液滴在不同能量表面上的铺展过程和润湿形态. 结果表明, 固液界面自由能随固液作用强度增加而增加, 并呈现不同液滴铺展速率和润湿特性. 固液作用强度小于1.6的低能表面呈现疏水特征, 继续增强固液作用强度时表面变为亲水, 而固液作用强度大于3.5的高能表面上液体呈完全润湿特征. 受微尺度条件下非连续、非对称作用力影响, 微液滴气液界面存在明显波动, 呈现与宏观液滴不同的界面特征. 统计意义下, 微小液滴在不同能量表面上铺展后仍可以形成特定接触角, 该接触角随固液作用强度增加而线性减小, 模拟结果与经典润湿理论计算获得的结果呈现相似变化趋势. 模拟结果从分子尺度为核化理论中的毛细假设提供了理论支持, 揭示了液滴气液界面和接触角的波动现象, 为核化速率理论预测结果和实验测定结果之间的差异提供了定性解释.     

Wetting of phospholipid membranes on hydrophilic surfaces - Concepts towards self-healing membranes

We report on the wetting behavior of phospholipid membranes on solid surfaces immersed in aqueous solution. Using fluorescence microscopy, the spreading velocity of fluid bilayers advancing from a lipid source is investigated. The kinetic spreading coefficient was measured as a function of temperature for pure DMPC membranes and as a function of charge density and cholesterol content for binary membranes. A theoretical model for the membrane flow is presented, which takes into account the liquid crystalline bilayer architecture of the lipid membrane. The spreading power results from the membrane-solid VdW interaction and is dissipated in hydrodynamic shear flow as well as by inter-monolayer friction within the bilayer. The frictional drag causes a dynamic tension gradient in the spreading membrane, which is manifested by a single exponential decay of the fluorescence intensity profile along the spreading direction. Obstacles are shown to act as pinning centers deforming the advancing line interface. However, no depinning was observed, since the centers are circumflown without abrupt relaxation. Received 6 November 1998     

Pb液滴在Ni基底润湿铺展行为的分子动力学模拟

采用分子动力学方法研究了Pb液滴在Ni(100)晶面、(110)晶面和(111)晶面的铺展润湿行为. 结果表明: Pb液滴在Ni(100)及(111)基底上的润湿铺展现象呈各向同性, 而在Ni(110)基底上的润湿铺展现象呈明显的各向异性, 且这种各向异性源于Ni(110)晶面点阵结构上Pb原子沿不同晶向的扩散机制及速度的明显差异; Pb液滴在Ni(111)晶面上铺展时, 未发生表面合金化, 液滴铺展动力学描述近似满足 R² ∝ t, 而液滴在(100)晶面和(110)晶面上铺展时表面产生合金化现象, 铺展动力学关系近似满足 R⁴ ∝ t, 且液滴在(100)晶面上的铺展速度高于(110)晶面上的铺展速度. 关键词： 分子动力学 润湿各向异性 铺展膜 扩散机制     

A mesoscopic model for (de)wetting

We present a mesoscopic model for simulating the dynamics of a non-volatile liquid on a solid substrate. The wetting properties of the solid can be tuned from complete wetting to total non-wetting. This model opens the way to study the dynamics of drops and liquid thin films at mesoscopic length scales of the order of the nanometer. As particular applications, we analyze the kinetics of spreading of a liquid drop wetting a solid substrate and the dewetting of a liquid film on a hydrophobic substrate. In all these cases, very good agreement is found between simulations and theoretical predictions.     

Cross-Over Between First-Order and Critical Wetting at the Liquid-Vapour Interface of n-Alkane/Methanol Mixtures: Tricritical Wetting and Critical Prewetting

A simple mean-field theory is presented which describes the basic observations of recent experiments revealing rich wetting behaviour of n-alkane/methanol mixtures at the liquid-vapour interface. The theory, qualitative and in part heuristic, is based on a microscopic lattice-gas model from which a Cahn–Landau approach is distilled. Besides the physics associated with the short-range components of the intermolecular interactions, effects of the long-range tails of the net van der Waals forces between interfaces are also taken into account. Further, gravitational thinning of the wetting phase is incorporated. The calculation of the spreading coefficient S is extended to the experimentally relevant situation in which the bulk adsorbate is slightly away from two-phase coexistence due to gravity. Analysis of this novel approximation to S for systems with short-range forces leads to the conclusion that the surface specific heat exponents _s =1,1/2, and 0, for first-order wetting, tricritical wetting and critical wetting, respectively, are robust with respect to (weak) gravitational thinning, consistently with experiment. For three different systems the adsorption is calculated as a function of temperature and compared with the experimentally measured ellipticity. Including weak long-range forces which favour wetting in the theory does not visibly alter the critical wetting transition for the nonane/methanol mixture, in contrast with the generic expectation of first-order wetting for such systems, but in good agreement with experiment. For decane/methanol weak long-range forces bring the transition very close to the prewetting critical point, leading to an adsorption behaviour closely reminiscent of short-range tricritical wetting, observed experimentally for alkane chain length between 9.6 and 10. Finally, for undecane/methanol the transition is clearly of first order. First-order wetting is also seen in the experiment.     

Capillary filling with pseudo-potential binary Lattice-Boltzmann model

We present a systematic study of capillary filling for a binary fluid by using a mesoscopic lattice Boltzmann model for immiscible fluids describing a diffusive interface moving at a given contact angle with respect to the walls. The phenomenological way to impose a given contact angle is analysed. Particular attention is given to the case of complete wetting, that is contact angle equal to zero. Numerical results yield quantitative agreement with the theoretical Washburn’s law, provided that the correct ratio of the dynamic viscosities between the two fluids is used. Finally, the presence of precursor films is experienced and it is shown that these films advance in time with a square-root law but with a different prefactor with respect to the bulk interface.     

Gels at interfaces

We discuss the adsorption of polymer gels on flat surfaces. Even in cases of complete wetting where the spreading power S is positive and where an equivalent liquid would spread, the elastic stresses due to the gel deformation upon adsorption oppose the spreading. The competition between elasticity characterized by the bulk shear modulus G and capillarity characterized by the spreading power S defines a typical length scale ℓ = S/G for the deformation in the gel. For loose gels ℓ can be of the order of 1 μm. Macroscopic gels larger than ℓ deform only at their edges over a region of size ℓ. Microscopic gels smaller than ℓ show a finite deformation despite the elastic stresses. The elastic stresses limit the spreading of the polymer, but solvent can be sucked out of a swollen gel by wetting the surface. The thin solvent film can extend rather far from the gel edge and carry solvent. We calculate the kinetics of the solvent film formation and of the solvent transfer from a more swollen gel to a less swollen gel. Received 16 July 2001     

Complete wetting of pits and grooves

For one-component volatile fluids governed by dispersion forces an effective interface Hamiltonian, derived from a microscopic density functional theory, is used to study complete wetting of geometrically structured substrates. Also the long range of substrate potentials is explicitly taken into account. Four types of geometrical patterns are considered: i) one-dimensional periodic arrays of rectangular or parabolic grooves and ii) two-dimensional lattices of cylindrical or parabolic pits. We present numerical evidence that at the centers of the cavity regions the thicknesses of the adsorbed films obey precisely the same geometrical covariance relation, which has been recently reported for complete cone and wedge filling. However, this covariance does not hold for the laterally averaged wetting film thicknesses. For sufficiently deep cavities with vertical walls and close to liquid-gas phase coexistence in the bulk, the film thicknesses exhibit an effective planar scaling regime, which as a function of undersaturation is characterized by a power law with the common critical exponent -1/3 as for a flat substrate, but with the amplitude depending on the geometrical features.     

Modelling of spreading process: effect from hydrogen bonds

Lubricant spreading on solid substrates has drawn considerable attention not only for the microscopic wetting theory but also for the dramatic application in head-disk interface of magnetic storage drive systems. Molecular dynamic simulation based on a coarse-grained bead-spring model has been used to study such a spreading process. The spreading profiles indicate that the hydrogen bonds among lubricant molecules and the hydrogen bonds between lubricant molecules and polar atoms of solid substrates will complicate the spreading process in a tremendous degree. The hydrogen bonds among lubricant molecules will strengthen the lubricant combination intensity, which may hinder most molecules from flowing down to the substrates and diffusing along the substrates. And the hydrogen bonds between lubricant molecules and polar atoms of solid substrates will confine the lubricant molecules around polar atoms, which may hinder the molecules from diffusing along the substrates and cause precursor film to vanish.     

The Role of Surface Chemistry on Spreading Kinetics of Molten Silicides on Silicon Carbide

The spreading time for millimeter-sized droplets of nonreactive molten silicides on silicon carbide in high vacuum is several orders of magnitude higher than typical spreading times observed in nonreactive metal/ceramic systems. To explain this paradox, two types of experiments were performed: (i) wetting experiments for various nonreactive CuSi alloys on -SiC single crystals using the sessile drop and dispensed drop techniques, with emphasis on determining the initial contact angle; and (ii) characterization of surface chemistry of SiC after different heat treatments in high-vacuum furnaces. It is shown that spreading kinetics in these systems are controlled by the kinetics of removing of wetting barriers present or developed in situ on SiC surface.     

Effect of solid surface on the formation of thin confined lubricating film of water with micro-content of oil

The formation of confined film between two contacting surfaces is significant for evaluating the lubricating ability of liquid. A micro-content of oil in water was experimentally demonstrated to be significantly effective to the film formation of water, which was much thicker than predicted by elastohydrodynamic lubrication theory. The effect of solid surface characteristics on the liquid film confined in a nanogap has been investigated. The film forming performances of such films were presented. The work of adhesion between two different phases was calculated, and the competitive wetting behaviours of water and oil on different solid surfaces were employed to understand the film formation mechanism.     

Dependence of wetting on cavitation during the spreading of a filler droplet on the ultrasonically agitated Al substrate

The cavitation characteristics during the spreading of a pure Sn liquid droplet subjected to ultrasonication were studied for the first time through high-speed photography to reveal the wetting mechanism. Ultrasonic vibration realized the spreading of Sn droplet on the nonwetting pure Al substrate. However, the oxide layer of the substrate at the spreading front is difficult to remove. The high-speed photography result shows that a noncavitation region consistently appears at the spreading front, because the acoustic pressure is below the cavitation threshold of 1.26 MPa. In particular, the width of the noncavitation region gradually increases as the size of the spreading area increases. Such a result accounts for the condition wherein the oxide layer at the spreading front is difficult to remove. Furthermore, the bubble density during spreading gradually decreases due to the decreased acoustic pressure of the thinned liquid. Finally, the bubble dynamics were calculated to verify the wetting mechanism.     

Sonoprocessing of wetting of SiC by liquid Al: A thermodynamic and kinetic study

The sonoprocessing of droplet spreading during the wetting process of molten aluminum droplets on SiC ceramic substrates at 700 °C is investigated in this paper. When wetting is assisted by a 20 kHz frequency ultrasonic field, the wettability of liquid metal gets enhanced, which has been determined by the variations in thermodynamic energy and wetting kinetics. Wetting kinetic characteristics are divided into two stages according to pinning and depinning states of substrate/droplet contact lines. The droplet is static when the contact line is pinning, while it is forced to move when the contact line is depinning. When analyzing the pinning stage, high-speed photography reveals the evidence of oxide films being rapidly crushed outside the aluminum droplet. In this work, atomic models of spherical Al core being wrapped by alumina shell are tentatively built, whose dioxide microstructures are being transformed from face-centered cubic into liquid at the atomic scale. At the same time, the wetting experiment reveals that the oxide films show changes in the period of sonoprocessing from 3rd to 5th second.During the ultrasonic spreading behavior in the late stage, there is a trend of evident expansion of the base contact area. The entire ultrasonic process lasts for no longer than 10 s. With the aid of ultrasonic sinusoidal waves, the wettability of metal Al gets a rapid improvement. Both molecular dynamic (MD) investigations and the experiments results reveal that the precursor film phenomenon is never found unless wetting is assisted by ultrasonic treatments. However, the precursor film appears near the triple line after using ultrasonics in the droplet wetting process, whose formation is driven by ultrasonic oscillations. Due to the precursor film, the ultrasonic wetting contact angle is lower than the non-ultrasonic contact angle. In addition, the time-variant effective ultrasonic energy has been quantitatively evaluated. The numerical expressions of thermodynamic variables are well verified by former ultrasonic spreading test results, which altogether provide an intrinsic explanation of the fast-decreasing contact angle of Al/SiC.     

Short-time dynamics of partial wetting

When a liquid drop contacts a wettable surface, the liquid spreads over the solid to minimize the total surface energy. The first moments of spreading tend to be rapid. For example, a millimeter-sized water droplet will wet an area having the same diameter as the drop within a millisecond. For perfectly wetting systems, this spreading is inertially dominated. Here we identify that even in the presence of a contact line, the initial wetting is dominated by inertia rather than viscosity. We find that the spreading radius follows a power-law scaling in time where the exponent depends on the equilibrium contact angle. We propose a model, consistent with the experimental results, in which the surface spreading is regulated by the generation of capillary waves.