疏水表面滑移流动及减阻特性的格子Boltzmann方法模拟

采用格子Boltzmann方法研究了固体壁面对流体的作用强度与其润湿性的关系,在此基础上进一步模拟了疏水表面微通道内的流体流动,获得了润湿性对疏水表面滑移流动及减阻特性的影响规律,证实了疏水表面表观滑移的存在性并揭示了其产生机理.结果表明,疏水性作用在疏水表面的近壁区诱导了一个低密度层,而表观滑移则发生在低密度层上.表观滑移是疏水表面具有减阻作用的直接原因,减阻效果随滑移长度的增大而增大.对于特定的流体系统,滑移长度是疏水表面的固有属性,仅是壁面润湿性的单一函数,而与流动本身的性质无关.     

Effect of the wall roughness on slip and rheological properties of hexadecane in molecular dynamics simulation of couette shear flow between two sinusoidal walls

Molecularly thin liquid films of alkanes in extreme conditions in a boundary lubrication regime have been investigated. The wall is modeled as a rough atomic sinusoidal wall. The effect on the boundary condition of the roughness characteristics, given by the period and amplitude of the sinusoidal wall, is studied here. The effect of the molecular length of the lubricating fluid is also examined here. The results show that the relative size of the fluid molecules and wall roughness determines the slip or nonslip boundary conditions. The effect of wall roughness characteristics on the rheological properties of the lubrication film is also studied.     

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

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

Behavior of adhesive boundary lubricated surfaces under shear: Effect of grafted diblock copolymers

The shear behavior and the normal interaction between mica surfaces covered by surfactant or surfactant-polymer mixtures were studied with a Surface Forces Apparatus (SFA) nanotribometer. If the surfaces are compressed while fully immersed in an aqueous surfactant solution that adsorbs in the form of flat bilayers, hemifusion can be induced. When the hemifused surfaces are subject to shear, at least five different dynamic regimes can be recognized. The general behavior may be described by a model based on the kinetics of formation and rupture of adhesive bonds between the shearing surfaces, with an additional viscous term. Once the adsorbed surfactant layer is decorated with physigrafted copolymers, the number of sliding regimes may be reduced to only one, in which the shear stress increases sublinearly with the driving velocity. The adhesion energy and the resistance to hemifusion of the adsorbed surfactant-polymer layers are also strongly modified as the grafting density increases.Received: 4 March 2004, Published online: 13 October 2004PACS: 46.55. + d Tribology and mechanical contacts - 81.40.Pq Friction, lubrication, and wear     

Structured surfaces for a giant liquid slip

We study experimentally how two key geometric parameters (pitch and gas fraction) of textured hydrophobic surfaces affect liquid slip. The two are independently controlled on precisely fabricated microstructures of posts and grates, and the slip length of water on each sample is measured using a rheometer system. The slip length increases linearly with the pitch but dramatically with the gas fraction above 90%, the latter trend being more pronounced on posts than on grates. Once the surfaces are designed for very large slips (>20 microm), however, further increase is not obtained in regular practice because the meniscus loses its stability. By developing near-perfect samples that delay the transition from a dewetted (Cassie) to a wetted (Wenzel) state until near the theoretical limit, we achieve giant slip lengths, as large as 185 microm.     

Temperature dependence of the slip length in polymer melts at attractive surfaces

Using Couette and Poiseuille flows, we extract the temperature dependence of the slip length, delta, from molecular dynamics simulations of a coarse-grained polymer model in contact with an attractive surface. delta is dictated by the ratio of bulk viscosity and surface mobility. At weakly attractive surfaces, lubrication layers form; delta is large and increases upon cooling. Close to the glass transition temperature Tg, very large slip lengths are observed. At a more attractive surface, a sticky surface layer is built up, giving rise to small slip lengths. Upon cooling, delta decreases at high temperatures, passes through a minimum, and grows for T-->Tg. At strongly attractive surfaces, the Navier-slip condition fails to describe Couette and Poiseuille flows simultaneously. The simulations are corroborated by a schematic, two-layer model suggesting that the observations do not depend on details of the computational model.     

Statistical thermodynamics of isolated rigid rod-like molecules near a surface

The configurational behaviour and thermodynamic properties of a dilute gas of rigid rod-like molecules in the vicinity of a macroscopic planar adsorption surface are investigated using statistical mechanics. The interaction energy between the surface and a rod-like molecule is determined as a function of both its molecular centre of mass separation R and its orientation relative to the surface. In calculating this interaction energy, each rod segment and molecule comprising the surface is assumed to interact through a Lennard-Jones pair potential. The average molecular order parameter is then determined as a function of R. We find that an isolated rod-like molecule tends to align nearly parallel to the surface for small separations. However, as R increases the order parameter first passes through a maximum then decays to zero as R ^-5 for large R. The configurational behaviour of an isolated rod-like molecule located between two parallel adsorption surfaces is also considered briefly. The surface spreading pressure, excess surface energy and entropy are also obtained for a dilute gas of rod-like molecules near a surface. We find that the extent of surface binding increases nearly exponentially with molecular length at constant temperature and surface density, and that the excess surface energy and entropy are essentially proportional to the molecular length.     

An updated overview of diamond-like carbon coating in tribology

During the last two decades, the industry (including scientists) has focused on diamond-like carbon (DLC) coating because of its wide range of application in various fields. This material has numerous applications in mechanical, electrical, tribological, biomedical, and optical fields. Severe friction and wear in some machine parts consumes high amount of energy, which makes the process energy inefficient. Thus, DLC coating can be an effective means to lower the friction and wear rate. Some important process variables that affect the tribological characteristics of DLC coating are adhesion promoter intermediate layer, substrate surface roughness, hydrogen incorporation or hydrogen non involvement, and coating deposition parameters (e.g., bias voltage, etching, current, precursor gas, time, and substrate temperature). Working condition of DLC-coated parts also affects the tribological characteristics, such as temperature, sliding speed and load, relative humidity, counter surface, and lubrication media (DLC additive interaction). Different types of lubricated oils and additives are used in engine parts to minimize friction and wear. DLC can be coated to the respective engine parts; however, DLC does not behave accordingly after coating because of lubricant oil and additive interaction with DLC. Some additive interacts positively and some behave negatively because of the tribochemical reactions between DLC coating and additives. Numerous conflicting views have been presented by several researchers regarding this coating additive interaction, resulting in unclear determination of true mechanism of such interaction. However, lubricant additive has been established to be more inert to DLC coating compared with uncoated metal surface because the additive is fabricated in such a way that it can react with metal surfaces. In this article, the tribological characteristics of different types of DLC coating in dry and lubricated conditions will be presented, and their behavior will be discussed in relation to working condition and processing parameters.     

Nature of mechanical instabilities and their effect on kinetic friction

It has long been recognized that kinetic friction F(k) between two solids must be due to instabilities, sudden "pops" of certain degrees of freedom. Here, such pops are studied with a focus on boundary lubrication. The pops' characteristics and consequently the friction-velocity relationship depend qualitatively on dimensionality, commensurability, and details of the lubricant wall interaction. It is found that F(k) should be small between commensurate surfaces. F(k) is large for incommensurate surfaces, unless the lubricant's motion is confined to 1D. The effects of thermal noise are discussed and computer simulations are employed to show the relevance of the predictions to less idealized models.     

Tunable-slip boundaries for coarse-grained simulations of fluid flow

On the micro- and nanoscale, classical hydrodynamic boundary conditions such as the no-slip condition no longer apply. Instead, the flow profiles exhibit "slip" at the surface, which is characterized by a finite slip length (partial slip). We present a new, systematic way of implementing partial-slip boundary conditions with arbitrary slip length in coarse-grained computer simulations. The main idea is to represent the complex microscopic interface structure by a spatially varying effective viscous force. An analytical equation for the resulting slip length can be derived for planar and for curved surfaces. The comparison with computer simulations of a DPD (dissipative particle dynamics) fluid shows that this expression is valid from full slip to no slip.     

A New Probe: AFM Measurements for Random Disorder Systems

We study the quenched random disorder(QRD) effects created by aerosil dispersion in the octylcyanobiphenyl(8CB) liquid crystal(LC) using atomic force microscopy technique. Gelation process in the 8CB+aerosil gels yields a QRD network which also changes the surface topography. By increasing the aerosil concentration, the original smooth pattern of LC sample surfaces is suppressed by the emergence of a fractal aerosil surface effect and these surfaces become more porous, rougher and they have more and larger crevices. The dispersed aerosil also serves as pinning centers for the liquid crystal molecules. It is observed that via the diffusion-limitedaggregation process, aerosil nano-particles yield a fractal-like surface pattern for the less disordered samples. As the aerosil dispersion increases, the surface can be described by more aggregated regions, which also introduces more roughness. Using this fact, we show that there is a net correlation between the short-range ordered x-ray peak widths(the results of previous x-ray diffraction experiments) and the calculated surface roughness. In other words, we show that these QRD gels can also be characterized by their surface roughness values.     

Association of polyethylene friction and thermal unfolding of interfacial albumin molecules

Under the articulation of artificial joints, ultra-high molecular weight polyethylene (UHMWPE) acts as a bearing surface under the lubrication of synovial fluid containing various proteins. Albumin is the most abundant composition and acts as the interfacial molecule in the boundary lubrication regime. The dissipated energy including thermal energy from the tribological process may lead to the conformational change of albumin molecules.In this study, a series of experiments were designed and carried out to investigate the association of thermal unfolding albumin and the frictional characteristics of highly-crosslinked UHMWPE (x-UHMWPE). An accelerated oxidation experiment was used to prepare x-UHMWPE with an oxidized surface. Analysis of the albumin protein by circular dichroism (CD) spectroscopy was performed to detect the conformational changes during a thermal process. In addition, a molecular simulation was performed to understand the structural change of albumin at various temperatures and the exposed hydrophobic contact areas. Linear reciprocating frictional tests were carried out to obtain the start-up friction coefficients. The results indicate that a decrease of α-helix content and an unfolding of the secondary structure of albumin were observed with increasing temperatures which may come from the frictional heat of joint articulation process. The conformational change of albumin differentiates the frictional characteristics for x-UHMWPE with different oxidation levels. A model, describing that the properties of the lubricating molecules and articulating surfaces may affect the adsorption of the boundary lubrication thin film which is critical to the tribological behavior, is proposed.     

Determination of Slip Length in Couette Flow Based on an Analytical Simulation Incorporating Surface Interaction

An analytical simulation based on a new model incorporating surface interaction is conducted to study the slip phenomenon in the Couette flow at different scales.The velocity profile is calculated by taking account of the micro-force between molecules and macro-force from the viscous shearing effect,as they contribute to the achievement of the slip length.The calculated results are compared with those obtained from the molecular dynamics simulation,showing an excellent agreement.Further,the effect of the shear rate on the slip is investigated.The results can well predict the fluid flow behaviors on a solid substrate,but has to be proved by experiment.     

金属纳米结构表面吸附的CO分子在外电场中的相互作用

建立金属纳米颗粒在外电场中的排列结构模型，用经典理论分析纳米结构金属表面上吸附的CO分子在外电场中的相互作用能，包括有效偶极子间的相互作用和与局域电场的相互作用，并讨论和计算了纳米颗粒表面附近的局域电场. 用Monte-Carlo方法进行数值计算和模拟，具体给出纳米颗粒表面CO分子的分布和相互作用能，表明金属表面纳米结构使CO产生凝聚，并使分子相互作用能增加，为解释异常红外吸收效应提供依据. 关键词： 纳米结构金属 吸附分子 相互作用 局域电场     

气体-表面相互作用的分子动力学模拟研究

采用分子动力学模拟方法研究了气体分子Ar在光滑和粗糙Pt表面上的散射规律.提出了一种速度抽样方法,计算了不同温度条件下气体分子对光滑和粗糙表面的切向动量适应系数和吸附概率.结果显示：光滑表面条件下,气体分子的切向动量系数和吸附概率都随着温度的升高而降低;粗糙度对气体分子切向动量与表面的适应具有极大的促进作用,当粗糙度足够大时,切向动量适应系数的大小趋近于1.0,对温度的敏感性也逐渐降低.采用粒子束方法对气体分子在光滑和粗糙表面上的散射规律进行了定量分析.总结了散射过程中气体分子的典型轨迹和动量变化规律,将气体分子在光滑表面的散射分为两种类型：单次碰撞后散射和多次碰撞后散射.单次碰撞后散射的气体分子平均切向动量有所减小,而经过多次碰撞后散射的气体分子则倾向于保持原有的平均切向动量.对于粗糙表面,粗糙度的存在使气体分子与表面间的动量和能量适应更加充分,导致气体分子在较粗糙表面上散射后的平均切向动量大幅减小并接近于0,且气体分子在表面上经历的碰撞次数越多,其散射后的能量损失越严重.     

微通道疏水表面滑移的耗散粒子动力学研究

了解疏水表面的滑移规律对其在流动减阻方面的应用至关重要.利用耗散粒子动力学(dissipative particle dynamics, DPD)方法研究了微通道疏水表面的滑移现象.采用固定住的粒子并配合修正的向前反弹机制,构建了DPD固体壁面边界模型,利用该边界模型模拟了平板间的Couette流动.研究结果表明,通过调整壁面与流体间排斥作用强度,壁面能实现从无滑移到滑移的转变,壁面与流体间排斥作用越强,即疏水性越强,壁面滑移越明显,并且滑移长度与接触角之间存在近似的二次函数关系.无滑移时壁面附近密度分布均匀,有滑移时壁面附近存在低密度区域,低密度区域阻碍了动量传递,致使壁面产生滑移.     

Interaction and local magnetic moments of metal phthalocyanine and tetraphenylporphyrin molecules on noble metal surfaces

In order to understand the Kondo effect observed in molecular systems, first-principles calculations have been widely used to predict the ground state properties of molecules on metal substrates. In this work, the interaction and the local magnetic moments of magnetic molecules （3d-metal phthalocyanine and tetraphenylporphyrin molecules） on noble metal surfaces are investigated based on the density functional theory. The calculation results show that the dz2 orbital of the transition metal atom of the molecule plays a dominant role in the molecule-surface interaction and the adsorption energy exhibits a simple declining trend as the adsorption distance increases. In addition, the Au（111） surface generally has a weak interaction with the adsorbed molecule compared with the Cu（ll 1） surface and thus serves as a better candidate substrate for studying the Kondo effect. The relation between the local magnetic moment and the Coulomb interaction U is examined by carrying out the GGA＋U calculation according to Dudarev＇s scheme. We find that the Coulomb interaction is essential for estimating the local magnetic moment in molecule-surface systems, and we suggest that the reference values of parameter U are 2 eV for Fe and 2-3 eV for Co.     

A dynamic rheological model for thin-film lubrication

In this study, the effects of the non-Newtonian rheological properties of the lubricant in a thin-film lubrication regime between smooth surfaces were investigated. The thin-film lubrication regime typically appears in Stribeck curves with a clearly observable minimum coefficient of friction (COF) and a low-COF region, which is desired for its lower energy dissipation. A dynamic rheology of the lubricant from the hydrodynamic lubrication regime to the thin-film lubrication regime was proposed based on the convected Maxwell constitutive equation. This rheology model includes the increased relaxation time and the yield stress of the confined lubricant thin film, as well as their dependences on the lubricant film thickness. The Deborah number (De number) was adopted to describe the liquid-solid transition of the confined lubricant thin film under shearing. Then a series of Stribeck curves were calculated based on Tichy's extended lubrication equations with a perturbation of the De number. The results show that the minimum COF points in the Stribeck curve correspond to a critical De number of 1.0, indicating a liquid-to-solid transition of the confined lubricant film. Furthermore, the two proposed parameters in the dynamic rheological model, namely negative slipping length b (indicating the lubricant interfacial effect) and the characteristic relaxation time λ 0 , were found to determine the minimum COF and the width of the low-COF region, both of which were required to optimize the shape of the Stribeck curve. The developed dynamic rheological model interprets the correlation between the rheological and interfacial properties of lubricant and its lubrication behavior in the thin-film regime.     

Generalized stacking fault energy surfaces in the molecular crystal αRDX

The generalized stacking fault (GSF) energy surfaces in the organic energetic molecular crystal, hexahydro-1,3,5-trinitro-s-triazine (RDX), were studied through atomistic simulations. Using a fully flexible molecular potential in highly damped molecular dynamics simulations, we determined quenched 0?K GSF energy surfaces and structures for a set of planes in the α-polymorph RDX crystal and subsequently compare predictions of slip or cleavage with available experimental observations. To account for the steric contributions and elastic shearing due to the presence of flexible molecules, a modified calculation procedure for the GSF energy surface is proposed that enables the distinction of elastic shear energy from the energy associated with the interfacial displacement discontinuity at the slip plane. Comparisons of the unstable stacking fault energy with the surface energy are used to differentiate cleavage planes from likely slip planes, and the calculations are found to be largely in agreement with available experimental data.     

Effect of dipole forces on the structure of the liquid phases of DNA

The dipole-dipole contribution to the energy of the pair interaction between DNA molecules has been calculated and analyzed. Rigid fragments of DNA, i.e., of a length of about the persistent length, which have discrete dipole moments of base pairs, are considered. For a given distance between the centers of mass of molecules, the energy of the dipole-dipole interaction is a function of three angular variables; the angles ?₁ and ?₂ between the central dipoles of both molecules and the z axis passing through the centers of the molecules, as well as the angle ξ between long axes of the molecules, are taken as these variables. It is shown that the dipole energy has minima when the mutual orientation of the molecules satisfies one of the following conditions: (i) ?₁ = ?₂ = 0 or (ii) ?₁ = ?₂ = π. The cholesteric twist angle ξ can be both positive and negative in dependence on the type of the minimum. For realistic cholesteric dispersion parameters, the dipole energy is only slightly lower than the experimentally known binding energy of the molecules in dispersion. The results allow the assumption that the dipole forces significantly affect the structure and other properties of DNA suspensions; in particular, they can lead to nontrivial texture phenomena, biaxial correlation, and multistability.