David N. Sibley

Slip or not slip? A methodical examination of the interface formation model using two-dimensional droplet spreading on a horizontal planar substrate as a prototype system

We consider the spreading of a thin two-dimensional droplet on a planar substrate as a prototype system to compare the contemporary model for contact line motion based on interface formation of Shikhmurzaev [Int. J. Multiphase Flow 19, 589–610 (1993)]10.1016/0301-9322(93)90090-H, to the more commonly used continuum fluid dynamical equations augmented with the Navier-slip condition. Considering quasistatic droplet evolution and using the method of matched asymptotics, we find that the evolution of the droplet radius using the interface formation model reduces to an equivalent expression for a slip model, where the prescribed microscopic dynamic contact angle has a velocity dependent correction to its static value. This result is found for both the original interface formation model formulation and for a more recent version, where mass transfer from bulk to surface layers is accounted for through the boundary conditions. Various features of the model, such as the pressure behaviour and rolling motion at the...

Fluid structure in the immediate vicinity of an equilibrium three-phase contact line and assessment of disjoining pressure models using density functional theory

We examine the nanoscale behavior of an equilibrium three-phase contact line in the presence of long-ranged intermolecular forces by employing a statistical mechanics of fluids approach, namely density functional theory (DFT) together with fundamental measure theory (FMT). This enables us to evaluate the predictive quality of effective Hamiltonian models in the vicinity of the contact line. In particular, we compare the results for mean field effective Hamiltonians with disjoining pressures defined through (I) the adsorption isotherm for a planar liquid film, and (II) the normal force balance at the contact line. We find that the height profile obtained using (I) shows good agreement with the adsorption film thickness of the DFT-FMT equilibrium density profile in terms of maximal curvature and the behavior at large film heights. In contrast, we observe that while the height profile obtained by using (II) satisfies basic sum rules, it shows little agreement with the adsorption film thickness of the DFT results. The results are verified for contact angles of 20, 40 and 60 degrees.

The contact line behaviour of solid-liquid-gas diffuse-interface models

A solid-liquid-gas moving contact line is considered through a diffuse-interface model with the classical boundary condition of no-slip at the solid surface. Examination of the asymptotic behaviour as the contact line is approached shows that the relaxation of the classical model of a sharp liquid-gas interface, whilst retaining the no-slip condition, resolves the stress, and pressure singularities associated with the moving contact line problem while the fluid velocity is well defined (not multi-valued). The moving contact line behaviour is analysed for a general problem relevant for any density dependent dynamic viscosity and volume viscosity, and for general microscopic contact angle and double well free-energy forms. Away from the contact line, analysis of the diffuse-interface model shows that the Navier–Stokes equations and classical interfacial boundary conditions are obtained at leading order in the sharp-interface limit, justifying the creeping flow problem imposed in an intermediate region in the seminal work of Seppecher [Int. J. Eng. Sci. 34, 977–992 (1996)]. Corrections to Seppechers work are given, as an incorrect solution form was originally used.

A comparison of slip, disjoining pressure, and interface formation models for contact line motion through asymptotic analysis of thin two-dimensional droplet spreading

The motion of a contact line is examined, and comparisons drawn, for a variety of models proposed in the literature. Pressure and stress behaviours at the contact line are examined in the prototype system of quasistatic spreading of a thin two-dimensional droplet on a planar substrate. The models analysed include three disjoining pressure models based on van der Waals interactions, a model introduced for polar fluids, and a liquid–gas diffuse-interface model; Navier-slip and two non-linear slip models are investigated, with three microscopic contact angle boundary conditions imposed (two of these contact angle conditions having a contact line velocity dependence); and the interface formation model is also considered. In certain parameter regimes it is shown that all of the models predict the same quasistatic droplet spreading behaviour.

On the moving contact line singularity: Asymptotics of a diffuse-interface model

The behaviour of a solid-liquid-gas system near the three-phase contact line is considered using a diffuse-interface model with no-slip at the solid and where the fluid phase is specified by a continuous density field. Relaxation of the classical approach of a sharp liquid-gas interface and careful examination of the asymptotic behaviour as the contact line is approached is shown to resolve the stress and pressure singularities associated with the moving contact line problem. Various features of the model are scrutinised, alongside extensions to incorporate slip, finite-time relaxation of the chemical potential, or a precursor film at the wall.Graphical abstract

Nanoscale fluid structure of liquid-solid-vapour contact lines for a wide range of contact angles

We study the nanoscale behaviour of the density of a simple fluid in the vicinity of an equilibrium contact line for a wide range of Young contact angles between 40 and 135 degrees. Cuts of the density profile at various positions along the contact line are presented, unravelling the apparent step-wise increase of the film height profile observed in contour plots of the density. The density profile is employed to compute the normal pressure acting on the substrate along the contact line. We observe that for the full range of contact angles, the maximal normal pressure cannot solely be predicted by the curvature of the adsorption film height, but is instead softened -- likely by the width of the liquid-vapour interface. Somewhat surprisingly however, the adsorption film height profile can be predicted to a very good accuracy by the Derjaguin-Frumkin disjoining pressure obtained from planar computations, as was first shown in [Nold et al., Phys. Fluids, 26, 072001, 2014] for contact angles less than 90 degrees, a result which here we show to be valid for the full range of contact angles. This suggests that while two-dimensional effects cannot be neglected for the computation of the normal pressure distribution along the substrate, one-dimensional planar computations of the Derjaguin-Frumkin disjoining pressure are sufficient to accurately predict the adsorption height profile.

Films, layers, and droplets: The effect of near-wall fluid structure on spreading dynamics

We present a study of the spreading of liquid droplets on a solid substrate at very small scales. We focus on the regime where effective wetting energy (binding potential) and surface tension effects significantly influence steady and spreading droplets. In particular, we focus on strong packing and layering effects in the liquid near the substrate due to underlying density oscillations in the fluid caused by attractive substrate-liquid interactions. We show that such phenomena can be described by a thin-film (or long-wave or lubrication) model including an oscillatory Derjaguin (or disjoining or conjoining) pressure and explore the effects it has on steady droplet shapes and the spreading dynamics of droplets on both an adsorption (or precursor) layer and completely dry substrates. At the molecular scale, commonly used two-term binding potentials with a single preferred minimum controlling the adsorption layer height are inadequate to capture the rich behavior caused by the near-wall layered molecular packing. The adsorption layer is often submonolayer in thickness, i.e., the dynamics along the layer consists of single-particle hopping, leading to a diffusive dynamics, rather than the collective hydrodynamic motion implicit in standard thin-film models. We therefore modify the model in such a way that for thicker films the standard hydrodynamic theory is realized, but for very thin layers a diffusion equation is recovered.

Nonequilibrium molecular dynamics simulations of nanoconfined fluids at solid-liquid interfaces

We investigate the hydrodynamic properties of a Lennard-Jones fluid confined to a nanochannel using molecular dynamics simulations. For channels of different widths and hydrophilic-hydrophobic surface wetting properties, profiles of the fluid density, stress, and viscosity across the channel are obtained and analysed. In particular, we propose a linear relationship between the density and viscosity in confined and strongly inhomogeneous nanofluidic flows. The range of validity of this relationship is explored in the context of coarse grained models such as dynamic density functional-theory.

The vicinity of an equilibrium three-phase contact line using density-functional theory: density profiles normal to the fluid interface

ABSTRACT The paper by Nold et al. [Phys. Fluids 26 (7), 072001 (2014)] examined density profiles and the micro-scale structure of an equilibrium three-phase (liquid–vapour–solid) contact line in the immediate vicinity of the wall using elements from the statistical mechanics of classical fluids, namely density-functional theory. The present research note, building on the above work, further contributes to our understanding of the nanoscale structure of a contact line by quantifying the strong dependence of the liquid–vapour density profile on the normal distance to the interface, when compared to the dependence on the vertical distance to the substrate. A recent study by Benet et al. [J. Phys. Chem. C 118 (38), 22079 (2014)] has shown that this could explain the emergence of a film-height-dependent surface tension close to the wall, with implications for the Frumkin–Derjaguin theory. GRAPHICAL ABSTRACT

The pressure tensor across a liquid-vapour interface

Inhomogeneous fluids exhibit physical properties that are neither uniform nor isotropic. The pressure tensor is a case in point, key to the mechanical description of the interfacial region. Kirkwood and Buff and, later, Irving and Kirkwood, obtained a formal treatment based on the analysis of the pressure across a planar surface [J. G. Kirkwood and F. P. Buff, J. Chem. Phys. 17(3), 338 (1949); J. H. Irving and J. G. Kirkwood, J. Chem. Phys. 18, 817 (1950)]. We propose a generalisation of Irving and Kirkwoods argument to fluctuating, non-planar surfaces and obtain an expression for the pressure tensor that is not smeared by thermal fluctuations at the molecular scale and corresponding capillary waves [F. P. Buff et al., Phys. Rev. Lett. 15, 621-623 (1965)]. We observe the emergence of surface tension, defined as an excess tangential stress, acting exactly across the dividing surface at the sharpest molecular resolution. The new statistical mechanical expressions extend current treatments to fluctuating inhomogeneous systems far from equilibrium.