Slipping moving contact lines: critical roles of de Gennes’s ‘foot’ in dynamic wetting

2019 ◽  
Vol 873 ◽  
pp. 110-150
Author(s):  
Hsien-Hung Wei ◽  
Heng-Kwong Tsao ◽  
Kang-Ching Chu

In the context of dynamic wetting, wall slip is often treated as a microscopic effect for removing viscous stress singularity at a moving contact line. In most drop spreading experiments, however, a considerable amount of slip may occur due to the use of polymer liquids such as silicone oils, which may cause significant deviations from the classical Tanner–de Gennes theory. Here we show that many classical results for complete wetting fluids may no longer hold due to wall slip, depending crucially on the extent of de Gennes’s slipping ‘foot’ to the relevant length scales at both the macroscopic and microscopic levels. At the macroscopic level, we find that for given liquid height $h$ and slip length $\unicode[STIX]{x1D706}$, the apparent dynamic contact angle $\unicode[STIX]{x1D703}_{d}$ can change from Tanner’s law $\unicode[STIX]{x1D703}_{d}\sim Ca^{1/3}$ for $h\gg \unicode[STIX]{x1D706}$ to the strong-slip law $\unicode[STIX]{x1D703}_{d}\sim Ca^{1/2}\,(L/\unicode[STIX]{x1D706})^{1/2}$ for $h\ll \unicode[STIX]{x1D706}$, where $Ca$ is the capillary number and $L$ is the macroscopic length scale. Such a no-slip-to-slip transition occurs at the critical capillary number $Ca^{\ast }\sim (\unicode[STIX]{x1D706}/L)^{3}$, accompanied by the switch of the ‘foot’ of size $\ell _{F}\sim \unicode[STIX]{x1D706}Ca^{-1/3}$ from the inner scale to the outer scale with respect to $L$. A more generalized dynamic contact angle relationship is also derived, capable of unifying Tanner’s law and the strong-slip law under $\unicode[STIX]{x1D706}\ll L/\unicode[STIX]{x1D703}_{d}$. We not only confirm the two distinct wetting laws using many-body dissipative particle dynamics simulations, but also provide a rational account for anomalous departures from Tanner’s law seen in experiments (Chen, J. Colloid Interface Sci., vol. 122, 1988, pp. 60–72; Albrecht et al., Phys. Rev. Lett., vol. 68, 1992, pp. 3192–3195). We also show that even for a common spreading drop with small macroscopic slip, slip effects can still be microscopically strong enough to change the microstructure of the contact line. The structure is identified to consist of a strongly slipping precursor film of length $\ell \sim (a\unicode[STIX]{x1D706})^{1/2}Ca^{-1/2}$ followed by a mesoscopic ‘foot’ of width $\ell _{F}\sim \unicode[STIX]{x1D706}Ca^{-1/3}$ ahead of the macroscopic wedge, where $a$ is the molecular length. It thus turns out that it is the ‘foot’, rather than the film, contributing to the microscopic length in Tanner’s law, in accordance with the experimental data reported by Kavehpour et al. (Phys. Rev. Lett., vol. 91, 2003, 196104) and Ueno et al. (Trans. ASME J. Heat Transfer, vol. 134, 2012, 051008). The advancement of the microscopic contact line is still led by the film whose length can grow as the $1/3$ power of time due to $\ell$, as supported by the experiments of Ueno et al. and Mate (Langmuir, vol. 28, 2012, pp. 16821–16827). The present work demonstrates that the behaviour of a moving contact line can be strongly influenced by wall slip. Such slip-mediated dynamic wetting might also provide an alternative means for probing slippery surfaces.

2020 ◽  
Vol 229 (10) ◽  
pp. 1945-1977 ◽  
Author(s):  
Yulii D. Shikhmurzaev

Abstract After a brief overview of the ‘moving contact-line problem’ as it emerged and evolved as a research topic, a ‘litmus test’ allowing one to assess adequacy of the mathematical models proposed as solutions to the problem is described. Its essence is in comparing the contact angle, an element inherent in every model, with what follows from a qualitative analysis of some simple flows. It is shown that, contrary to a widely held view, the dynamic contact angle is not a function of the contact-line speed as for different spontaneous spreading flows one has different paths in the contact angle-versus-speed plane. In particular, the dynamic contact angle can decrease as the contact-line speed increases. This completely undermines the search for the ‘right’ velocity-dependence of the dynamic contact angle, actual or apparent, as a direction of research. With a reference to an earlier publication, it is shown that, to date, the only mathematical model passing the ‘litmus test’ is the model of dynamic wetting as an interface formation process. The model, which was originated back in 1993, inscribes dynamic wetting into the general physical context as a particular case in a wide class of flows, which also includes coalescence, capillary breakup, free-surface cusping and some other flows, all sharing the same underlying physics. New challenges in the field of dynamic wetting are discussed.


1997 ◽  
Vol 334 ◽  
pp. 211-249 ◽  
Author(s):  
YULII D. SHIKHMURZAEV

A general mathematical model which describes the motion of an interface between immiscible viscous fluids along a smooth homogeneous solid surface is examined in the case of small capillary and Reynolds numbers. The model stems from a conclusion that the Young equation, σ1 cos θ = σ2 − σ3, which expresses the balance of tangential projection of the forces acting on the three-phase contact line in terms of the surface tensions σi and the contact angle θ, together with the well-established experimental fact that the dynamic contact angle deviates from the static one, imply that the surface tensions of contacting interfaces in the immediate vicinity of the contact line deviate from their equilibrium values when the contact line is moving. The same conclusion also follows from the experimentally observed kinematics of the flow, which indicates that liquid particles belonging to interfaces traverse the three-phase interaction zone (i.e. the ‘contact line’) in a finite time and become elements of another interface – hence their surface properties have to relax to new equilibrium values giving rise to the surface tension gradients in the neighbourhood of the moving contact line. The kinematic picture of the flow also suggests that the contact-line motion is only a particular case of a more general phenomenon – the process of interface formation or disappearance – and the corresponding mathematical model should be derived from first principles for this general process and then applied to wetting as well as to other relevant flows. In the present paper, the simplest theory which uses this approach is formulated and applied to the moving contact-line problem. The model describes the true kinematics of the flow so that it allows for the ‘splitting’ of the free surface at the contact line, the appearance of the surface tension gradients near the contact line and their influence upon the contact angle and the flow field. An analytical expression for the dependence of the dynamic contact angle on the contact-line speed and parameters characterizing properties of contacting media is derived and examined. The role of a ‘thin’ microscopic residual film formed by adsorbed molecules of the receding fluid is considered. The flow field in the vicinity of the contact line is analysed. The results are compared with experimental data obtained for different fluid/liquid/solid systems.


2014 ◽  
Vol 746 ◽  
Author(s):  
E. Kirkinis ◽  
S. H. Davis

AbstractA recent hydrodynamic theory of liquid slippage on a solid substrate (Kirkinis & Davis, Phys. Rev. Lett., vol. 110, 2013, 234503) gives rise to a sequence of eddies (Moffatt vortices) that co-move with a moving contact line (CL) in a liquid wedge. The presence of these vortices is established through secular equations that depend on the dynamic contact angle $\alpha $ and capillary number Ca. The limiting case $\alpha \rightarrow 0$ is associated with the appearance of such vortices in a channel. The vortices are generated by the relative motion of the interfaces, which in turn is due to the motion of the CL. This effect has yet to be observed in experiment.


2018 ◽  
Vol 841 ◽  
pp. 767-783 ◽  
Author(s):  
Yi Xia ◽  
Paul H. Steen

Contact-line mobility characterizes how fast a liquid can wet or unwet a solid support by relating the contact angle $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ to the contact-line speed $U_{CL}$. The contact angle changes dynamically with contact-line speeds during rapid movement of liquid across a solid. Speeds beyond the region of stick–slip are the focus of this experimental paper. For these speeds, liquid inertia and surface tension compete while damping is weak. The mobility parameter $M$ is defined empirically as the proportionality, when it exists, between $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ and $U_{CL}$, $M\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}=U_{CL}$. We discover that $M$ exists and measure it. The experimental approach is to drive the contact line of a sessile drop by a plane-normal oscillation of the drop’s support. Contact angles, displacements and speeds of the contact line are measured. To unmask the mobility away from stick–slip, the diagram of $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ against $U_{CL}$, the traditional diagram, is remapped to a new diagram by rescaling with displacement. This new diagram reveals a regime where $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ is proportional to $U_{CL}$ and the slope yields the mobility $M$. The experimental approach reported introduces the cyclically dynamic contact angle goniometer. The concept and method of the goniometer are illustrated with data mappings for water on a low-hysteresis non-wetting substrate.


Author(s):  
Ndivhuwo Musehane ◽  
Rhameez Herbst

Computational Fluid Dynamics (CFD) is used to study the spreading process of a water droplet with a radius of 0.00275mm impacting a wax surface at a velocity of 1.18ms−1 . This type of flow is considered to be Multiphase, incompressible, laminar, surface tension dominated and is governed by the Navier stokes and continuity equations. To accurately model the spreading process 3 different contact angle models are investigated, two of which take into account the moving contact line. The governing equations are solved using the open source C++ library OpenFOAM, which uses a Finite Volume Method (FVM) of discretization and a Volume Of Fluid (VOF) interface capturing method. The VOF method is known to produce unphysical velocities when high pressure gradients exist between the two phases, thus a numerical improvement is implemented to reduce the magnitudes of the unphysical velocities. The improvement reduces the magnitudes of the unphysical velocities and as shown in literature their magnitudes increase with an increase in surface tension dominance. The improvement is implemented together with different contact angle models and results obtained show that contact angle models that take into account the moving contact line gives a good correlation of the spreading diameter obtained numerically with the one obtained experimentally.


1994 ◽  
Vol 116 (4) ◽  
pp. 938-945 ◽  
Author(s):  
P. C. Wayner

A new physical model for the spreading dynamics of fluids with an apparent finite contact angle on solid substrates is presented. The model is based on the premise that both interfacial intermolecular forces and temperature control change-of-phase heat transfer and (therefore) motion in the moving contact line region. Classical change-of-phase kinetics and interfacial concepts like the Kelvin–Clapeyron, Young–Dupre, and augmented Young–Laplace equations are used to compare the effects of stress (change in apparent dynamic contact angle) and temperature (superheat). Explicit equations are obtained for the velocity, heat flux, and superheat in the contact line region as a function of the change in the apparent contact angle. Comparisons with experimental data demonstrate that the resulting interfacial model of evaporation/condensation not only describes the “apparently isothermal” contact line movement in these systems at 20°C but also describes the substrate superheat at the critical heat flux.


Author(s):  
Yi Lu ◽  
Aritra Sur ◽  
Dong Liu ◽  
Carmen Pascente ◽  
Paul Ruchhoeft

Electrowetting has drawn significant interests due to the potential applications in electronic displays, lab-on-a-chip devices and electro-optical switches, etc. Current understanding of electrowetting-induced droplet dynamics is hindered by the inadequacy of available numerical and theoretical models in properly handling the dynamic contact angle at the moving contact line. A combined numerical and experimental approach was employed in this work to study the spatiotemporal responses of a droplet subject to EW with both direct current and alternating current actuating signals. The time evolution of the droplet shape was measured using high-speed photography. Computational fluid dynamics models were developed by using the Volume of Fluid-Continuous Surface Force method in conjunction with a selected dynamic contact angle model. It was found that the numerical models were able to accurately predict the key parameters of the electrowetting-induced droplet dynamics.


2016 ◽  
Vol 20 (5) ◽  
pp. 1163-1182 ◽  
Author(s):  
Shixin Xu ◽  
Weiqing Ren

AbstractThe level set method is one of the most successful methods for the simulation of multi-phase flows. To keep the level set function close the signed distance function, the level set function is constantly reinitialized by solving a Hamilton-Jacobi type of equation during the simulation. When the fluid interface intersects with a solid wall, a moving contact line forms and the reinitialization of the level set function requires a boundary condition in certain regions on the wall. In this work, we propose to use the dynamic contact angle, which is extended from the contact line, as the boundary condition for the reinitialization of the level set function. The reinitialization equation and the equation for the normal extension of the dynamic contact angle form a coupled system and are solved simultaneously. The extension equation is solved on the wall and it provides the boundary condition for the reinitialization equation; the level set function provides the directions along which the contact angle is extended from the contact line. The coupled system is solved using the 3rd order TVD Runge-Kutta method and the Godunov scheme. The Godunov scheme automatically identifies the regions where the angle condition needs to be imposed. The numerical method is illustrated by examples in three dimensions.


2007 ◽  
Vol 574 ◽  
pp. 343-364 ◽  
Author(s):  
OLIVIER DEVAUCHELLE ◽  
CHRISTOPHE JOSSERAND ◽  
STEPHANE ZALESKI

We study the dewetting of a porous plate withdrawn from a liquid bath. The contact angle is fixed to zero and the flow is assumed to be almost parallel to the plate (lubrication approximation). The ordinary differential equation involving the position of the water surface is analysed in phase space by means of numerical integration. We show the existence of a stationary moving contact line with zero contact angle below a critical value of the capillary number η U/γ. Above this value, no stationary contact line can exist. An analytical model, based on asymptotic matching is developed, which reproduces the dependence of the critical capillary number on the angle of the plate with respect to the horizontal (3/2 power law), provided the capillary length is much larger than the square root of the porous-medium permeability. In addition, it is shown that the classical lubrication equation leads not only to the well-known Landau–Levich–Derjaguin films, but also to a family of films for which thickness is not imposed by the problem parameters.


Author(s):  
Fangjun Hong ◽  
Ping Cheng ◽  
Zhen Sun ◽  
Huiying Wu

In this paper, the electrowetting dynamics of a droplet on a dielectric surface was investigated numerically by a mathematical model including dynamic contact angle and contact angle hysteresis. The fluid flow is described by laminar N-S equation, the free surface of the droplet is modeled by the Volume of Fluid (VOF) method, and the electrowetting force is incorporated by exerting an electrical force on the cells at the contact line. The Kilster’s model that can deal with both receding and advancing contact angle is adopted. Numerical results indicate that there is overshooting and oscillation of contact radius in droplet spreading process before it ceases the movement when the excitation voltage is high; while the overshooting is not observed for low voltage. The explanation for the contact line overshooting and some special characteristics of variation of contact radius with time were also conducted.


Sign in / Sign up

Export Citation Format

Share Document