移动接触线的物理力学研究

移动接触线,指两种互不相溶的流体在固体表面形成移动的三相接触区域.移动接触区域跨越多个尺度,其中三相物质之间的相互作用影响着整个流场的动力学特征.由于在能源、航天、生物等领域中的重要应用和迅速发展,移动接触线在新的应用背景下发展了新的难题.标度分析是度量接触线自相似扩展的重要手段.本文以移动接触线的标度关系为主线,介绍了"力–电–热–化学"多场耦合环境下,亲水内角、微柱阵列、可溶解固体、水力压裂滞后区等复杂几何结构的刚性/柔性固体表面,采用物理力学方法对于移动接触线动力学属性研究的进展.通过跨尺度实验研究、大规模分子动力学模拟和分子动理论/水动力学理论相结合的方法,发现了类固体前驱膜、单分子前驱水链、锯齿形接触线等新现象.从原子尺度的界面结构到连续尺度的流动特性,讨论了移动接触线自相似扩展的标度关系,以及其驱动来源、能量耗散、边界条件等物理机制和规律,为多物理场中的"Huh-Scriven佯谬"探索了解答,为移动接触线的前景和应用提出了展望.     

接触力学的发展概况

接触问题的解法,最先在弹性理论中得到发展且比较完善;以后,考虑到粘弹性、塑性、各向异性、动力学和运动学特性,而大大地发展了各种接触问题,远远超过了弹性理论的范围.本文对接触问题进行了分类,介绍了经典、非经典的弹性接触理论的发展概况.最后,比较详细地分析了国内外当前的研究结果和发展动向.      

连接结构接触界面非线性力学建模研究

连接界面上存在的跨尺度、多物理场和非线性行为是引起结构复杂非线性动力学的主要原因。由于连接界面力学行为的复杂性,以及对连接界面进行直接试验观测的困难,连接界面的力学建模一直是非常具有挑战性的科学问题。本文首先从分析结合面的跨尺度物理机理入手,将名义的光滑平面视作凹凸不平的粗糙面,考虑单个微凸体的黏滑摩擦行为,建立接触载荷与变形的非线性关系,然后采用GW(Greenwood和Williamson, GW)模型数理统计方法建立整个粗糙界面的跨尺度力学模型,并与公开文献中试验结果进行对比。然后考虑连接界面典型非线性特征,提出一种改进的Iwan唯象模型,利用精细有限元方法获得非线性特征结果,采用系统辨识理论建立连接结构的降阶力学模型,并利用有限元结果进行模型验证。结果表明,本文提出的粗糙界面跨尺度模型在法向载荷较小时与试验结果吻合较好,改进的Iwan模型能够较好描述连接界面非线性特征,并与有限元结果吻合较好。      

发展新材料的一些物理力学问题：材料的力学性质和材料设计

本文论述了发展新材料必须进行材料设计，在设计当中要进行各个学科层次的分工与合作。结构材料设计的主要问题是材料的力学性质，这包括高强度、优良的加工性能和高的阻尼本领，其关键环节是加强物理力学和力学物理的研究。     

低维碳和氮化硼材料的物理力学研究进展

碳纳米管、石墨烯和六方氮化硼等低维材料具有优异的力学和电学性质,已经引起广泛的科学兴趣。然而由电荷、分子轨道、电子结构和自旋态构成的低维材料的局域场与力学变形、机械运动和物理化学环境等外场间往往存在强烈耦合,这导致低维材料会呈现出新颖独特的物理力学性能。论文对近年来碳纳米管、石墨烯和六方氮化硼等低维材料的力学性能、力电耦合与器件原理、表面和界面结构性能调控、层间相互作用、能量耗散和摩擦等物理力学方面的研究进展进行了简要综述,并讨论了利用低维材料多场耦合特性和结构性能关联发展新型功能器件的方法和途径,以及纳米力学和纳尺度物理力学的前沿和发展趋势。     

芶清泉与物理力学在中国的发展

芶清泉从20 世纪60 年代初致力于将物理力学学科研究建立在原子物理和高压物理的基础之上,系统地开展了原子结构、相互作用及原子碰撞等方面的科学研究工作;提出了物理力学发展微观化的思路和观点;修订并规范了物理力学研究的内容、方法和目的;组织并召开了4 届全国物理力学学术会议,为物理力学的有组织发展、科研队伍的建立与壮大以及物理力学的科学研究等做出了重大贡献.     

粗糙表面接触力学问题的重新分析

为了克服基于统计学参数的接触模型的尺度依赖性以及现有接触分形模型推导过程中初始轮廓表征受控于接触面积或取样长度的不足,基于粗糙表面轮廓分形维数$D$、尺度系数$G$ 和最大微凸体轮廓基底尺寸$l$,建立了新的粗糙表面接触分形模型,探讨了微凸体变形机制、粗糙表面的真实接触面积和接触载荷的关系,揭示了接触界面的孔隙率和真实接触面积随端面形貌、表面接触压力等参数变化的规律,给出了不同形貌界面被压实的最大变形量. 结果表明:微凸体变形从弹性变形开始,并随着平均接触压力$p_{\rm m}$ 的增大逐步向弹塑性变形和完全塑性变形转变;接触界面的初始孔隙率$\phi_{0}$ 随$D$ 的增大而增大,压实孔隙所需要的最大变形量$\delta $ 也随之增大;接触压力$p_{\rm c}$ 增大,孔隙率$\phi$ 减小,并随着$D$ 的增大和$G$ 减小,$\phi$ 快速减小,直至填实,变为零;$D$ 较小时,$G$ 的增大对真实接触面积的增大影响较小;$D$ 较大时,$G$ 的增大对真实接触面积的增大作用明显. 研究成果为端面摩擦副的润滑与密封设计提供了理论基础.      

针尖的化学物理力学研究

扫描探针显微镜的发明, 使人们了解纳米、分子和原子尺度的超微结构, 探测原子、分子间的力、电、磁及其复杂的物理和化学性质, 以及进行单原子、单分子操纵 成为现实. 本文从探针技术与表面化学物理力学耦合的角度出发, 首先对针尖的化学物理力 学研究领域进行了概述, 接着介绍了针尖的几种重要的化学修饰方法, 包括金属薄膜、自组 装单分子膜、胶体粒子及碳纳米管修饰; 然后以理论与实验结合的方式介绍了几个研究活跃 的领域: 表面力及分子间力(主要包括Van der Waals力, 双电层力, 憎水亲合力, Casimir力和单键力等)的理论描述与测量, 针尖的化 学物理力学在化学力滴定和表面化学识别等研究中的应用. 讨论了针尖与基底材料对测量力 的影响. 而这些复杂的原子、分子相互作用和物理、化学、力学及生物特性的实现均发生于 小小针尖上, 由此我们提出了``针尖力学'的概念. 并且指出多场(如电场、磁场、超声、 微波等)作用下, 针尖的化学物理力学研究将成为力学交叉学科研究的重点和热点.     

点、线接触真实粗糙表面的弹流润滑研究

本文给出了点、线接触的真实粗糙表面的微弹流数值解。在给定随机粗糙表面样本后,求解大小不同的载荷和粗糙的弹流问题。从计算结果可以看出,由于Reynolds方程中速度项的作用,在表面对应粗糙的位置处引起了压力变化,从而因其产生的弹性变形使粗糙变得平滑。对此光滑表面解可以看出,因粗糙引起的压力和膜厚的变化在光滑解附近波动。载荷较大时,压力分节接近固体接触情况。     

计算力学的有限元线法

本文对新近提出的有限元线法作一简介,对其正在取得的进展及发展前景作了评述及展望。      

Physical mechanics investigations of moving contact lines

移动接触线,指两种互不相溶的流体在固体表面形成移动的三相接触区域.移动接触区域跨越多个尺度,其中三相物质之间的相互作用影响着整个流场的动力学特征.由于在能源、航天、生物等领域中的重要应用和迅速发展,移动接触线在新的应用背景下发展了新的难题.标度分析是度量接触线自相似扩展的重要手段.本文以移动接触线的标度关系为主线,介绍了"力-电-热-化学"多场耦合环境下,亲水内角、微柱阵列、可溶解固体、水力压裂滞后区等复杂几何结构的刚性/柔性固体表面,采用物理力学方法对于移动接触线动力学属性研究的进展.通过跨尺度实验研究、大规模分子动力学模拟和分子动理论/水动力学理论相结合的方法,发现了类固体前驱膜、单分子前驱水链、锯齿形接触线等新现象.从原子尺度的界面结构到连续尺度的流动特性,讨论了移动接触线自相似扩展的标度关系,以及其驱动来源、能量耗散、边界条件等物理机制和规律,为多物理场中的"Huh-Scriven佯谬"探索了解答,为移动接触线的前景和应用提出了展望.     

Finite element methods for a class of continuum models for immiscible flows with moving contact lines

In this paper, we present a finite element method for two‐phase incompressible flows with moving contact lines. We use a sharp interface Navier–Stokes model for the bulk phase fluid dynamics. Surface tension forces, including Marangoni forces and viscous interfacial effects, are modeled. For describing the moving contact lines, we consider a class of continuum models that contains several special cases known from the literature. For the whole model, describing bulk fluid dynamics, surface tension forces, and contact line forces, we derive a variational formulation and a corresponding energy estimate. For handling the evolving interface numerically, the level‐set technique is applied. The discontinuous pressure is accurately approximated by using a stabilized extended finite element space. We apply a Nitsche technique to weakly impose the Navier slip conditions on the solid wall. A unified approach for discretization of the (different types of) surface tension forces and contact line forces is introduced. Results of numerical experiments are presented, which illustrate the performance of the solver. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical simulation of the flow with contact lines

The paper proposes a physical model for the motion of the contact line and the gas-liquid interface. The local motion of the contact line at the solid wall is assumed and the interface between gas and liquid is traced by a level function. The numerically. The motion of the water column in a vertical pipe is computed and the results are in good agreement with experimental data.     

Recent progress in the moving contact line problem: a review

As pointed out long ago by Laplace, viscosity may become a large perturbation to capillary phenomena, especially close to solid surfaces where molecules may stick. A spectacular consequence of this is the impossibility for a triple line to move on a solid if the liquid/vapor interface is considered as a material surface and if the usual no slip boundary condition is enforced. As shown recently this specific phenomenon of contact line motion can be described by coupled van der Waals and fluid equations, yielding a rational theory that is divergence free and consistent with the equilibrium results. Far from the triple line, the equations of fluid mechanics are recovered in their usual form. In this approach, the contact line move close to the solid by evaporation or condensation, which requires (for evaporation) the molecules to jump above a high potential barrier on their way from the liquid to the vapor. An Arrhenius factor makes this process intrinsically slow, compared to molecular speeds. For (realistic) very small Arrhenius factors, the motion of the triple line induces a dynamical change of the functions in the van der Waals equations. This may lead to dynamical wetting and dewetting transitions, that is, to a change of the contact angle from a finite to a zero value or conversely. The dynamical wetting transition has been observed in liquids flowing down a plate (see Blake and Ruschak, Nature 282 (1979) 489–491) cusps on the contact line appear when it recedes faster than the speed of transition. Similar ideas account well also for the known sensitivity of contact line mobility to vapor pressure. To cite this article: Y. Pomeau, C. R. Mecanique 330 (2002) 207–222.     

Finite element framework for describing dynamic wetting phenomena

The finite element simulation of dynamic wetting phenomena, requiring the computation of flow in a domain confined by intersecting a liquid–fluid free surface and a liquid–solid interface, with the three‐phase contact line moving across the solid, is considered. For this class of flows, different finite element method (FEM) implementations have been used in the literature, and in some cases, these produced apparently contradictory results. In the present paper, a robust framework for the FEM simulation of dynamic wetting flows is developed, which, by consistently adhering to the FEM methodology, leaves no room for ad hoc ‘optional’ variations in the numerical handling of these flows. The developed approach makes it possible to conduct a convergence study, assess the spatial resolution required to achieve a preset accuracy and provide the corresponding benchmark calculations. This analysis allows one to identify numerical artefacts, which had previously been interpreted as physical effects, and demonstrates that suppressing numerical errors using a ‘strong’ implementation of a boundary condition creates bigger and less detectable errors elsewhere in the computational domain. We provide practical recommendations on the spatial resolution required by a numerical scheme for a given set of non‐dimensional similarity parameters and give a user‐friendly step‐by‐step guide specifying the entire implementation, which allows the reader to easily reproduce all presented results including the benchmark calculations. It is also shown how the developed framework accommodates generalizations of the mathematical model accounting for additional physical effects, such as gradients in surface tensions. Copyright © 2011 John Wiley & Sons, Ltd.     

Elastic-plastic contact mechanics of indentations accounting for phase transformations

A contact mechanics model is developed which takes into account possible phase transformations in materials induced by hydrostatic and shear stresses associated with indentation. The proposed model allows prediction of the average thickness and approximate shape of the phase transformation zone in semiconductors and ceramics under various types of diamond indenters. The results of theoretical calculation are in good agreement with the available experimental data.     

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

An efficient and energy stable scheme for a phase‐field model for the moving contact line problem

In this paper, we propose for the first time a linearly coupled, energy stable scheme for the Navier–Stokes–Cahn–Hilliard system with generalized Navier boundary condition. We rigorously prove the unconditional energy stability for the proposed time discretization as well as for a fully discrete finite element scheme. Using numerical tests, we verify the accuracy, confirm the decreasing property of the discrete energy, and demonstrate the effectiveness of our method through numerical simulations in both 2‐D and 3‐D. Copyright © 2015 John Wiley & Sons, Ltd.     

Boundary element method for moving and rolling contact of 2D elastic bodies with defects

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes. The project supported by the National Natural Science Foundation of China (19772025)     

一类功能梯度材料表面接触力学研究进展

针对一类杨氏模量梯度变化的功能梯度材料, 考虑摩擦、微动磨损和黏附等因素, 综述了梯度材料有限尺寸、梯度变化规律、梯度涂层厚度、界面摩擦热、压头形状等对表面接触力学性能的影响; 根据不同接触模型中接触界面应力场分布, 分析多种因素影响下功能梯度材料表面抵抗磨损的能力; 最后给出了功能梯度材料接触力学研究中仍存在的主要科学问题及进一步研究展望.