Substrate elastic deformation due to vertical component of liquid-vapor interfacial tension

Young’s equation is a fundamental equation in capillarity and wetting,which reflects the balance of the horizontal components of the three interfacial tensions with the contact angle(CA).However,it does not consider the vertical component of the liquid-vapor interfacial tension(VCLVIT).It is now well understood that the VCLVIT causes the elastic deformation of the solid substrate,which plays a significant role in the fabrication of the microfluidic devices because of the wide use of the soft materials.In this paper,the theoretical,experimental,and numerical aspects of the problem are reviewed.The effects of the VCLVIT-induced surface deformation on the wetting and spreading,the deflection of the microcantilever,and the elasto-capillarity and electroelasto-capillarity are discussed.Besides a brief review on the historical development and the recent advances,some suggestions on the future research are also provided.     

Equilibrium shapes of strained islands with non-zero contact angles

The equilibrium morphology of a strained island on an elastic substrate is determined. The island is assumed to partially wet the substrate (Volmer-Weber growth) and thus makes a non-zero contact angle with the surface. Both isotropic and anisotropic misfit strain are allowed. Two- and three-dimensional equilibrium island shapes are determined by using expressions for the elastic strain energy in the small-slope approximation. In this limit, the problem can be reduced to a singular integral-differential equation for the island thickness. We find that when there is a non-zero contact angle, all island shapes, for a given ratio of the elastic stress to surface energy, attain a form that is independent of the specific contact angle under an appropriate scaling. We show that for islands with non-zero contact angles, as the island volume increases, the shape approaches the geometry of a completely wetting island. But when the volume decreases, these islands approach a point while islands with a zero contact angle, approach a finite length line segment of zero volume. Multiple-hump equilibrium shapes are found. Single-humped islands are shown to have a lower chemical potential than multiple-humped islands, implying that they are the most stable. This conclusion is shown to be consistent with a stability analysis of the two-dimensional case. The effects of a tetragonal misfit strain on the three-dimensional island shape is investigated.     

A finite element method for free surface flows of incompressible fluids in three dimensions. Part II. Dynamic wetting lines

To date, few researchers have solved three‐dimensional free surface problems with dynamic wetting lines. This paper extends the free surface finite element method (FEM) described in a companion paper [Cairncross RA, Schunk PR, Baer TA, Sackinger PA, Rao RR. A finite element method for free surface flows of incompressible fluid in three dimensions. Part I. Boundary fitted mesh motion. International Journal for Numerical Methods in Fluids 2000; 33 : 375–403] to handle dynamic wetting. A generalization of the technique used in two‐dimensional modeling to circumvent double‐valued velocities at the wetting line, the so‐called kinematic paradox, is presented for a wetting line in three dimensions. This approach requires the fluid velocity normal to the contact line to be zero, the fluid velocity tangent to the contact line to be equal to the tangential component of web velocity, and the fluid velocity into the web to be zero. In addition, slip is allowed in a narrow strip along the substrate surface near the dynamic contact line. For realistic wetting line motion, a contact angle that varies with wetting speed is required because contact lines in three dimensions typically advance or recede at different rates depending upon location and/or have both advancing and receding portions. The theory is applied to capillary rise of static fluid in a corner, the initial motion of a Newtonian droplet down an inclined plane, and extrusion of a Newtonian fluid from a nozzle onto a moving substrate. The extrusion results are compared with experimental visualization. Copyright © 2000 John Wiley & Sons, Ltd.     

微液滴在PDMS软基体表面的动态摩擦行为研究

通过固液界面摩擦力测试装置研究了微液滴在PDMS软基体表面运动时的动态摩擦学行为,并对微液滴体积、滑动速度及软基体力学性能对固液界面动态摩擦行为的影响进行了分析. 结果表明:微液滴在软基体表面运动时表现出最大静摩擦力和动态摩擦力. 最大静摩擦力与微液滴黏度和速度梯度呈正比,动态摩擦力与微液滴体积、滑动速度和基体力学性能有关. 随着微液滴体积的增加,三相接触线长度增加,动态摩擦力增加;随着相对滑动速度增加,三相接触线长度及接触角滞后增加,动态摩擦力增加;随着软基体弹性模量降低,固液界面黏附力增加,固液界面运动能量耗散增加,动态摩擦力增加. 研究结果可为PDMS软基体表面微液滴的精确驱动和运动参数优化提供理论指导,也可进一步丰富固液界面摩擦理论.      

Comparative study of two lattice Boltzmann multiphase models for simulating wetting phenomena: implementing static contact angles based on the geometric formulation

Wetting phenomena are widespread in nature and industrial applications. In general, systems concerning wetting phenomena are typical multicomponent/multiphase complex fluid systems. Simulating the behavior of such systems is important to both scientific research and practical applications. It is challenging due to the complexity of the phenomena and difficulties in choosing an appropriate numerical method. To provide some detailed guidelines for selecting a suitable multiphase lattice Boltzmann model, two kinds of lattice Boltzmann multiphase models, the modified S-C model and the H-C-Z model, are used in this paper to investigate the static contact angle on solid surfaces with different wettability combined with the geometric formulation (Ding, H. and Spelt, P. D. M. Wetting condition in diffuse interface simulations of contact line motion. Physical Review E, 75(4), 046708 (2007)). The specific characteristics and computational performance of these two lattice Boltzmann method (LBM) multiphase models are analyzed including relationship between surface tension and the control parameters, the achievable range of the static contact angle, the maximum magnitude of the spurious currents (MMSC), and most importantly, the convergence rate of the two models on simulating the static contact angle. The results show that a wide range of static contact angles from wetting to non-wetting can be realized for both models. MMSC mainly depends on the surface tension. With the numerical parameters used in this work, the maximum magnitudes of the spurious currents of the two models are on the same order of magnitude. MMSC of the S-C model is universally larger than that of the H-C-Z model. The convergence rate of the S-C model is much faster than that of the H-C-Z model. The major foci in this work are the frequently-omitted important details in simulating wetting phenomena. Thus, the major findings in this work can provide suggestions for simulating wetting phenomena with LBM multiphase models along with the geometric formulation.     

The wetting problem of fluids on solid surfaces. Part 2: the contact angle hysteresis

In part 1 (Gouin, [13]), we proposed a model of dynamics of wetting for slow movements near a contact line formed at the interface of two immiscible fluids and a solid when viscous dissipation remains bounded. The contact line is not a material line and a Young-Dupré equation for the apparent dynamic contact angle taking into account the line celerity was proposed. In this paper we consider a form of the interfacial energy of a solid surface in which many small oscillations are superposed on a slowly varying function. For a capillary tube, a scaling analysis of the microscopic law associated with the Young-Dupré dynamic equation yields a macroscopic equation for the motion of the contact line. The value of the deduced apparent dynamic contact angle yields for the average response of the line motion a phenomenon akin to the stick-slip motion of the contact line on the solid wall. The contact angle hysteresis phenomenon and the modelling of experimentally well-known results expressing the dependence of the apparent dynamic contact angle on the celerity of the line are obtained. Furthermore, a qualitative explanation of the maximum speed of wetting (and dewetting) can be given.Received: 5 June 2001, Accepted: 24 May 2003, Published online: 29 July 2003PACS: 02.90, 47.50, 66.20, 68.03, 68.08     

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

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

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

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

Motion of a three-phase contact line along a wet surface

The motion of the contact line in gas-liquid-solid systems is theoretically investigated for small values of the capillary number and Reynolds number. The possible existence on the solid substrate of a residual microscopic film formed by adsorbed liquid molecules is taken into account and the spreading characteristics of the liquid on dry and wet substrates are compared. It is shown that, in accordance with the experimental data, in the model employed the motion of the liquid during wetting is rolling motion, and that the increase in the dynamic contact angle is slower for a wet than for a dry substrate. The maximum dynamic contact angle is much less than 180°. The flow structure in the neighborhood of the moving contact line is analyzed and it is shown that under certain conditions regions with closed streamlines may be formed. The reason for this is the self-induced Marangoni effect — the reaction of the surface tension gradient on the liquid-solid boundary caused by the liquid flow on the flow that caused it.Based on a paper read at the Seventh Congress on Theoretical and Applied Mechanics, Moscow, August 1991. Presented by R. I. Nigmatulin.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 81–89, November–December, 1992.     

Modelling the wetting of a solid occlusion by a liquid film

We investigate in this work how the presence of an occlusion affects the dynamics of the wetting front of a liquid film draining down a vertical surface. This numerical study is developed in the context of the lubrication approximation. Through a parametric study, we show that depending on the asymptotic film thickness and the fluid properties, there exists a critical substrate contact angle below which separation of the contact line from the occlusion wall is observed which results in the appearance of a dry zone in the wake of the occlusion. In analogy with external aerodynamics, we also show that a sharp corner in the occlusion can induce this contact line separation. Our numerical results also highlight the importance of the occlusion wettability on the morphology of the wetting front suggesting a possible mechanism to control and mitigate the often undesirable fingering instability.     

Effects of oxidation and surface roughness on contact angle

Contact angle is known to be a parameter that effects boiling. This study was undertaken to measure contact angle of high and low surface tension fluids on copper and aluminum surfaces.Data were taken for polished, oxidized, and rough surfaces. A simple, yet fairly accurate method of measuring the static equilibrium contact angle of a solid/liquid interface is presented. The principles of a line light source and tilting plate were modified and then combined in the design of this apparatus. The angles obtained and their variation with the solid surface properties were in good agreement with previously published data. The contact angle of distilled water o of the organic fluids and refrigerants tested were in the range of 2–5°. Roughness and oxidation reduce the contact angle. If the depth of the roughness is less than 0.5 μm contact angle. The apparatus is fairly simple in construction, is inexpensive, and has good reproductibity. The measured angles were then compared to those measured with the sessile drop method.     

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.     

Green’s formula and singularity at a triple contact line. Example of finite-displacement solution

The various equations at the surfaces and triple contact lines of a deformable body are obtained from a variational condition, by applying Green’s formula in the whole space and on the Riemannian surfaces. The surface equations are similar to the Cauchy’s equations for the volume, but involve a special definition of the ‘divergence’ (tensorial product of the covariant derivatives on the surface and the whole space). The normal component of the divergence equation generalizes the Laplace’s equation for a fluid–fluid interface. Assuming that Green’s formula remains valid at the contact line (despite the singularity), two equations are obtained at this line. The first one expresses that the fluid–fluid surface tension is equilibrated by the two surface stresses (and not by the volume stresses of the body) and suggests a finite displacement at this line (contrary to the infinite-displacement solution of classical elasticity, in which the surface properties are not taken into account). The second equation represents a strong modification of Young’s capillary equation. The validity of Green’s formula and the existence of a finite-displacement solution are justified with an explicit example of finite-displacement solution in the simple case of a half-space elastic solid bounded by a plane. The solution satisfies the contact line equations and its elastic energy is finite (whereas it is infinite for the classical elastic solution). The strain tensor components generally have different limits when approaching the contact line under different directions. Although Green’s formula cannot be directly applied, because the stress tensor components do not belong to the Sobolev space H¹(V)$H^1(\text{V})$, it is shown that this formula remains valid. As a consequence, there is no contribution of the volume stresses at the contact line. The validity of Green’s formula plays a central role in the theory.     

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.     

Three-Phase Effective Contact Angle in a Model Pore

We consider a model pore (2D) in which a sharp interface between two fluids contact a third fluid which wets the solid boundary. If the configuration is capillary dominated, the geometry can be determined analytically in terms of the effective contact angle. This angle depends not only on the interfacial tensions, but also on the capillary pressures. However, if the height of the cusp formed by the wetting fluid is much smaller than the pore width, the effective contact angle is a simple function of the interfacial tensions. It turns out to be the same function as in the case of an undeformed wetting layer of molecular thickness. The analytical expression for the effective contact angle has been confirmed by a numerical technique, known as the lattice-Boltzmann method. This method, in turn, has been validated with Neumann's law for the three-phase contact angles.     

The wetting problem of fluids on solid surfaces. Part 1: the dynamics of contact lines

The understanding of the spreading of liquids on solid surfaces is an important challenge for contemporary physics. Today, the motion of the contact line formed at the intersection of two immiscible fluids and a solid is still subject to dispute. In this paper, a new picture of the dynamics of wetting is offered through an example of non-Newtonian slow liquid movements. The kinematics of liquids at the contact line and equations of motion are revisited. Adherence conditions are required except at the contact line. Consequently, for each fluid, the velocity field is multivalued at the contact line and generates an equivalent concept of line friction but stresses and viscous dissipation remain bounded. A Young-Dupré equation for the apparent dynamic contact angle between the interface and solid surface depending on the movements of the fluid near the contact line is proposed.Received: 5 June 2001, Accepted: 24 May 2003, Published online: 29 July 2003PACS: 47.17., 47.50, 66.20, 68.08, 83.50     

Some generic capillary-driven flows

This paper deals with numerical simulations of some capillary-driven flows. The focus is on the wetting phenomenon in sintering-like flows and in the imbibition of liquids into a porous medium. The wetting phenomenon is modeled using the coupled Cahn–Hilliard/Navier–Stokes system. The Cahn–Hilliard equation is treated as a system where the chemical potential is solved first followed by the composition. The equations are discretised in space using piecewise linear functions. Adaptive finite element method is implemented with an ad hoc error criterion that ensures mesh resolution along the vicinity of the interface. In the 3D case we use parallel adaptive finite element method. First, a basic wetting of a liquid drop on a solid surface is shown and is established the independence of the dynamic contact angle on the interface width. In addition, the dependence of the dynamic contact angle on the Capillary number is matched with experimental data. Next, some generic sintering-like flows with a fixed matrix is presented. Different geometries in 2D and 3D are considered. We observed rapid wetting, precursor films, coalescence, breakup of melt drops as well as pore migration and elimination that are all microstructural characteristics of a liquid phase sintering. Finally, the effect of equilibrium contact angles on imbibition of liquid into a porous medium is studied.     

ASYMPTOTIC ANALYSIS OF MODE Ⅱ STATIONARY GROWTH CRACK ON ELASTIC-ELASTIC POWER LAW CREEPING BIMATERIAL INTERFACE

A mechanical model was established for mode Ⅱ interfacial crack static growing along an elastic-elastic power law creeping bimaterial interface. For two kinds of boundary conditions on crack faces, traction free and frictional contact, asymptotic solutions of the stress and strain near tip-crack were given. Results derived indicate that the stress and strain have the same singularity, there is not the oscillatory singularity in the field; the creep power-hardening index n and the ratio of Young' s module notably influence the cracktip field in region of elastic power law creeping material and n only influences distribution of stresses and strains in region of elastic material. When n is bigger, the creeping deformation is dominant and stress fields become steady, which does not change with n.Poisson ' s ratio does not affect the distributing of the crack- tip field.     

The Effect of Isopropyl Alcohol and Non-Ionic Surfactant Mixtures on the Wetting of Porous Coated Paper

The influence of isopropyl alcohol and non-ionic surfactant solutions on aqueous droplet wetting behaviour on porous coated paper was determined. Paper coatings provide a micro- and nano-porous surface structure, which strictly speaking cannot be described in simple roughness terms as sub-surface lateral absorption directly impacts on the apparent contact angle. It is this very deviation from an idealised system that leads to novel wetting phenomena. Isopropyl alcohol and surfactant-based systems, both of which are commonly used in the printing industry, show differences in wetting behaviour, on both short and long timescales, with changes in the relative composition of the mixtures. Small variations of 0.1?wt% in surfactant concentration have a dramatic influence on the dynamic surface tension, and thus the wetting. It was observed that the wetting kinetics for isopropyl alcohol and surfactant solutions were different in terms of both wetting area and the penetration rate, even in cases where the dynamic surface tension of the solutions was kept the same. Different stages in the wetting and following drying processes could be observed with near infrared spectral imaging. In addition, the surfactant chemistries such as their degrees of hydrophilicity and molecular weights generated comparative differences in the wetting kinetics. The dominating factor affecting the wetting was, as expected, the solid?Cliquid interfacial energy defined on the practical porous substrate, which differed from the direct comparison with dynamic surface tension, thus exemplifying the deviation from idealised surface roughness behaviour when considering porous materials. An apparent ??equivalent?? surface roughness value for the porous material was determined, and it was seen that an increase in this equivalent parameter enhanced the rate of wetting behaviour with decreasing solution surface tension, and so also affected the wetting evolution. The wetting was enhanced by cavities in the coating layer, which were enlarged by the penetrating liquids.