The moving contact line: the slip boundary condition

1976 ◽  
Vol 77 (4) ◽  
pp. 665-684 ◽  
Author(s):  
E. B. Dussan V.

The singularity at the contact line which is present when the usual fluidmechanical modelling assumptions are made is removed by permitting the fluid to slip along the wall. The aim of this study is to assess the sensitivity of the overall flow field to the form of the slip boundary condition. Explicit solutions are obtained for three different slip boundary conditions. Two length scales emerge: the slip length scale and the meniscus length scale. It is found that on the slip length scale the flow fields are quite different; however, when viewed on the meniscus length scale, i.e. the length scale on which almost all fluidmechanical measurements are made, all of the flow fields appear the same. It is found that the characteristic of the slip boundary condition which affects the overall flow field is the magnitude of the slip length.

1988 ◽  
Vol 197 ◽  
pp. 157-169 ◽  
Author(s):  
P. A. Durbin

It has previously been shown that the no-slip boundary condition leads to a singularity at a moving contact line and that this forces one to admit some form of slip. Present considerations on the energetics of slip due to shear stress lead to a yield stress boundary condition. A model for the distortion of the liquid state near solid boundaries gives a physical basis for this boundary condition. The yield stress condition is illustrated by an analysis of a slender drop rolling down an incline. That analysis provides a formula for the frictional drag resisting the drop movement. With the present boundary condition the length of the slip region becomes a property of the fluid flow.


1997 ◽  
Vol 337 ◽  
pp. 49-66 ◽  
Author(s):  
Q. CHEN ◽  
E. RAMÉ ◽  
S. GAROFF

The dynamics of a spreading liquid body are dictated by the interface shape and flow field very near the moving contact line. The interface shape and flow field have been described by asymptotic models in the limit of small capillary number, Ca. Previous work established the validity and limitations of these models of the interface shape (Chen et al. 1995). Here, we study the flow field near the moving contact line. Using videomicroscopy, particle image velocimetry, and digital image analysis, we simultaneously make quantitative measurements of both the interface shape and flow field from 30 μm to a few hundred microns from the contact line. We compare our data to the modulated-wedge solution for the velocity field near a moving contact line (Cox 1986). The measured flow fields demonstrate quantitative agreement with predictions for Ca[les ]0.1, but deviations of ∼5% of the spreading velocity at Ca≈0.4. We observe that the interface shapes and flow fields become geometry independent near the contact line. Our experimental technique provides a way of measuring the interface shape and velocity field to be used as boundary conditions for numerical calculations of the macroscopic spreading dynamics.


2013 ◽  
Vol 715 ◽  
pp. 283-313 ◽  
Author(s):  
Yi Sui ◽  
Peter D. M. Spelt

AbstractUsing a slip-length-based level-set approach with adaptive mesh refinement, we have simulated axisymmetric droplet spreading for a dimensionless slip length down to $O(1{0}^{\ensuremath{-} 4} )$. The main purpose is to validate, and where necessary improve, the asymptotic analysis of Cox (J. Fluid Mech., vol. 357, 1998, pp. 249–278) for rapid droplet spreading/dewetting, in terms of the detailed interface shape in various regions close to the moving contact line and the relation between the apparent angle and the capillary number based on the instantaneous contact-line speed, $\mathit{Ca}$. Before presenting results for inertial spreading, simulation results are compared in detail with the theory of Hocking & Rivers (J. Fluid Mech., vol. 121, 1982, pp. 425–442) for slow spreading, showing that these agree very well (and in detail) for such small slip-length values, although limitations in the theoretically predicted interface shape are identified; a simple extension of the theory to viscous exterior fluids is also proposed and shown to yield similar excellent agreement. For rapid droplet spreading, it is found that, in principle, the theory of Cox can predict accurately the interface shapes in the intermediate viscous sublayer, although the inviscid sublayer can only be well presented when capillary-type waves are outside the contact-line region. However, $O(1)$ parameters taken to be unity by Cox must be specified and terms be corrected to ${\mathit{Ca}}^{+ 1} $ in order to achieve good agreement between the theory and the simulation, both of which are undertaken here. We also find that the apparent angle from numerical simulation, obtained by extrapolating the interface shape from the macro region to the contact line, agrees reasonably well with the modified theory of Cox. A simplified version of the inertial theory is proposed in the limit of negligible viscosity of the external fluid. Building on these results, weinvestigate the flow structure near the contact line, the shear stress and pressure along the wall, and the use of the analysis for droplet impact and rapid dewetting. Finally, we compare the modified theory of Cox with a recent experiment for rapid droplet spreading, the results of which suggest a spreading-velocity-dependent dynamic contact angle in the experiments. The paper is closed with a discussion of the outlook regarding the potential of using the present results in large-scale simulations wherein the contact-line region is not resolved down to the slip length, especially for inertial spreading.


Author(s):  
Derek C. Tretheway ◽  
Luoding Zhu ◽  
Linda Petzold ◽  
Carl D. Meinhart

This work examines the slip boundary condition by Lattice Boltzmann simulations, addresses the validity of the Navier’s hypothesis that the slip velocity is proportional to the shear rate and compares the Lattice Boltzmann simulations to the experimental results of Tretheway and Meinhart (Phys. of Fluids, 14, L9–L12). The numerical simulation models the boundary condition as the probability, P, of a particle to bounce-back relative to the probability of specular reflection, 1−P. For channel flow, the numerically calculated velocity profiles are consistent with the experimental profiles for both the no-slip and slip cases. No-slip is obtained for a probability of 100% bounce-back, while a probability of 0.03 is required to generate a slip length and slip velocity consistent with the experimental results of Tretheway and Meinhart for a hydrophobic surface. The simulations indicate that for microchannel flow the slip length is nearly constant along the channel walls, while the slip velocity varies with wall position as a results of variations in shear rate. Thus, the resulting velocity profile in a channel flow is more complex than a simple combination of the no-slip solution and slip velocity as is the case for flow between two infinite parallel plates.


2010 ◽  
Vol 76 (769) ◽  
pp. 1370-1378
Author(s):  
Takahiro ITO ◽  
Hayato FUJIWARA ◽  
Akira HIBI ◽  
Yuzuru NADA ◽  
Susumu NODA

2018 ◽  
Vol 849 ◽  
pp. 805-833 ◽  
Author(s):  
Xianmin Xu ◽  
Yana Di ◽  
Haijun Yu

The sharp-interface limits of a phase-field model with a generalized Navier slip boundary condition for binary fluids with moving contact lines are studied by asymptotic analysis and numerical simulations. The effects of the mobility number as well as a phenomenological relaxation parameter on the boundary condition are considered. In asymptotic analysis, we consider both the cases that the mobility number is proportional to the Cahn number and the square of the Cahn number, and derive the sharp-interface limits for several set-ups of the boundary relaxation parameter. It is shown that the sharp-interface limit of the phase-field model is the standard two-phase incompressible Navier–Stokes equations coupled with several different slip boundary conditions. Numerical results are consistent with the analysis results and also illustrate the different convergence rates of the sharp-interface limits for different scalings of the two parameters.


Author(s):  
Shi-Ming Li ◽  
Danesh K. Tafti

A mean-field free-energy lattice Boltzmann method (LBM) is applied to simulate moving contact line dynamics. It is found that the common bounceback boundary condition leads to an unphysical velocity at the solid wall in the presence of surface forces. The magnitude of the unphysical velocity is shown proportional to the local force term. The velocity-pressure boundary condition is generalized to solve the problem of the unphysical velocity. The simulation results are compared with three different theories for moving contact lines, including a hydrodynamic theory, a molecular kinetic theory, and a linear cosine law of moving contact angle versus capillary number. It is shown that the current LBM can be used to replace the three theories in handling moving contact line problems.


Author(s):  
Qi Zhou ◽  
Chiu-On Ng

The hydrodynamic dispersion of a neutral non-reacting solute due to steady electro-osmotic flow in a circular channel with longitudinal step changes of zeta potential and hydrodynamic slippage is analyzed in this study. The channel wall is periodically micro-patterned along the axial position with alternating slip-stick stripes of distinct zeta potentials. Existing studies on electrically driven hydrodynamic dispersion are based on flow subject to either the no-slip boundary condition on the capillary surface or the simplification of lubrication approximation. Taking wall slippage into account, a homogenization analysis is performed in this study to derive the hydrodynamic dispersion coefficient without subject to the long-wave constraint of the lubrication approximation, but for a general case where the length of one periodic unit of wall pattern is comparable with the channel radius. The flow and the hydrodynamic dispersion coefficient are calculated numerically, using the packages MATLAB and COMSOL, as functions of controlling parameters including the period length of the wall pattern, the area fraction of the slipping region (EOF-suppressing) in a periodic unit, the ratio of the two zeta potentials, the intrinsic hydrodynamic slip length, the Debye parameter, and the Péclet number. The dispersion coefficient is found to show notable, non-monotonic in certain situations, dependence on these controlling parameters. It is noteworthy that the introduction of hydrodynamic slippage will generate much richer behaviors of the hydrodynamic dispersion than the situation with no-slip boundary condition, as slippage interacts with zeta potentials in the EOF-suppressing and EOF-supporting regions (either likewise or oppositely charged).


Sign in / Sign up

Export Citation Format

Share Document