Hydrodynamic theory of liquid slippage on a solid substrate near a moving contact line

In this Letter a hydrodynamic theory of liquid slippage on a solid substrate near a moving contact line is proposed. A family of spatially varying slip lengths in the Navier slip law recovers the results of past formulations for slip in continuum theories and molecular dynamics simulations and is consistent with well-established experimental observations of complete wetting. This formulation gives a general approach for continuum hydrodynamic theories. New fluid flow behaviors are also predicted yet to be seen in experiment. © 2013 American Physical Society.

Moffatt vortices induced by the motion of a contact line

A 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 α and capillary number Ca. The limiting case α→O 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.

An In-Situ Observation On Initial Aggregation Process Of Colloidal Particles Near Three-Phase Contact Line Of Air, Water And Vertical Substrate

The self-assembling process near the three-phase contact line of air, water and vertical substrate is widely used to produce various kinds of nanostructured materials and devices. We perform an in-situ observation on the self-assembling process in the vicinity of the three phase contact line. Three kinds of aggregations, i.e. particle-particle aggregation, particle-chain aggregation and chain-chain aggregation, in the initial stage of vertical deposition process are revealed by our experiments. It is found that the particle particle aggregation and the particle-chain aggregation can be qualitatively explained by the theory of the capillary immersion force and mirror image force, while the chain-chain aggregation leaves an opening question for the further studies. The present study may provide more deep insight into the self-assembling process of colloidal particles.

Contact zone dynamics and the evolution of reproductive isolation in a North American treefrog, the spring peeper (Pseudacris crucifer)

Despite over seven decades of speciation research and 25 years of phylogeographic studies, a comprehensive understanding of mechanisms that generate biological species remains elusive. In temperate zones, the pervasiveness of range fragmentation and subsequent range expansions suggests that secondary contact between diverging lineages may be important in the evolution of species. Thus, such contact zones provide compelling opportunities to investigate evolutionary processes, particularly the roles of geographical isolation in initiating, and indirect selection against hybrids in completing (reinforcement), the evolution of reproductive isolation and speciation. The spring peeper (Pseudacris crucifer) has six well-supported mitochondrial lineages many of which are now in secondary contact. Here I investigate the evolutionary consequences of secondary contact of two such lineages (Eastern and Interior) in Southwestern Ontario using genetic, morphological, acoustical, experimental, and behavioural evidence to show accentuated divergence of the mate recognition system in sympatry. Mitochondrial and microsatellite data distinguish these two lineages but also show ongoing hybridization. Bayesian assignment tests and cline analysis imply asymmetrical introgression of Eastern lineage nuclear markers into Interior populations. Male calls are divergent between Eastern and Interior allopatric populations and show asymmetrical reproductive character displacement in sympatry. Female preference of pure lineage individuals is also exaggerated in sympatry, with hybrids showing intermediate traits and preference. I suggest that these patterns are most consistent with secondary reinforcement. I assessed levels of post-zygotic isolation between the Eastern and Interior lineages using a laboratory hybridization experiment. Hybrid tadpoles showed equal to or greater fitness than their pure lineage counterparts, but this may be countered through competition. More deformities and developmental anomalies in hybrid tadpoles further suggest post-zygotic isolation. Despite evidence for pre-mating isolation between the two lineages, isolation appears incomplete (i.e. hybridization is ongoing). I hypothesize that potentially less attractive hybrids may circumvent female choice by adopting satellite behaviour. Although mating tactics are related to body size, genetic status may play a role. I show that pure Eastern males almost always engage in calling, while hybrids adopt a satellite tactic. An absence of assortative mating, despite evidence of female preference, suggests successful satellite interception possibly facilitating introgression.

INVESTIGATING THE ROLE OF THE CONTACT LINE IN HETEROGENEOUS NUCLEATION WITH HIGH SPEED IMAGING

While nucleation of solids in supercooled liquids is ubiquitous [15, 65, 66], surface crystallization, the tendency for freezing to begin preferentially at the liquid-gas interface, has remained puzzling [74, 18, 68, 69, 51, 64, 72, 16]. Here we employ high-speed imaging of supercooled water drops to study the phenomenon of heterogeneous surface crystallization. Our geometry avoids the "point-like contact" of prior experiments by providing a simple, symmetric contact line (triple line defined by the substrate-liquid-air interface) for a drop resting on a homogeneous silicon substrate. We examine three possible mechanisms that might explain these laboratory observations: (i) Line Tension at the triple line, (ii) Thermal Gradients within the droplets and (iii) Surface Texture. In our first study we record nearly perfect spatial uniformity in the immersed (liquid-substrate) region and, thereby, no preference for nucleation at the triple line. In our second study, no influence of thermal gradients on the preference for freezing at the triple line was observed. Motivated by the conjectured importance of line tension (τ) [1, 66] for heterogeneous nucleation, we also searched for evidence of a transition to surface crystallization at length scales on the order of δ ∼ τ/σ, where σ is the surface tension [14]; poorly constrained τ [49] leads to δ ranging from microns to nanometers. We demonstrate that nano-scale texture causes a shift in the nucleation to the three-phase contact line, while micro-scale texture does not. The possibility of a critical length scale has implications for the effectiveness of nucleation catalysts, including formation of ice in atmospheric clouds [7].

A novel methodology for simulating contact-line behavior in capillary-driven flows

Despite the wide swath of applications where multiphase fluid contact lines exist, there is still no consensus on an accurate and general simulation methodology. Most prior numerical work has imposed one of the many dynamic contact-angle theories at solid walls. Such approaches are inherently limited by the theory accuracy. In fact, when inertial effects are important, the contact angle may be history dependent and, thus, any single mathematical function is inappropriate. Given these limitations, the present work has two primary goals: 1) create a numerical framework that allows the contact angle to evolve naturally with appropriate contact-line physics and 2) develop equations and numerical methods such that contact-line simulations may be performed on coarse computational meshes.

Fluid flows affected by contact lines are dominated by capillary stresses and require accurate curvature calculations. The level set method was chosen to track the fluid interfaces because it is easy to calculate interface curvature accurately. Unfortunately, the level set reinitialization suffers from an ill-posed mathematical problem at contact lines: a ``blind spot'' exists. Standard techniques to handle this deficiency are shown to introduce parasitic velocity currents that artificially deform freely floating (non-prescribed) contact angles. As an alternative, a new relaxation equation reinitialization is proposed to remove these spurious velocity currents and its concept is further explored with level-set extension velocities.

To capture contact-line physics, two classical boundary conditions, the Navier-slip velocity boundary condition and a fixed contact angle, are implemented in direct numerical simulations (DNS). DNS are found to converge only if the slip length is well resolved by the computational mesh. Unfortunately, since the slip length is often very small compared to fluid structures, these simulations are not computationally feasible for large systems. To address the second goal, a new methodology is proposed which relies on the volumetric-filtered Navier-Stokes equations. Two unclosed terms, an average curvature and a viscous shear VS, are proposed to represent the missing microscale physics on a coarse mesh.

All of these components are then combined into a single framework and tested for a water droplet impacting a partially-wetting substrate. Very good agreement is found for the evolution of the contact diameter in time between the experimental measurements and the numerical simulation. Such comparison would not be possible with prior methods, since the Reynolds number Re and capillary number Ca are large. Furthermore, the experimentally approximated slip length ratio is well outside of the range currently achievable by DNS. This framework is a promising first step towards simulating complex physics in capillary-dominated flows at a reasonable computational expense.

Deployment strategy for controlled morphologies in sessile, mixed colloidal droplets

Colloidal systems offer an effective medium to micro-engineer complex structures without involving sophisticated fabrication procedures. This article presents a deployment strategy of multiple droplets of different colloidal composition and utilizes the inherent capillary flow driven self assembly of nanoparticles to construct stacks of multiple materials on a given glass substrate. Here we used aqueous nano-crystalline titania and nano-amorphous silica solutions as the two materials. Initially, a pure nanotitania (nanosilica) droplet is deployed and allowed to dry partially. Subsequently, a second droplet of pure nanosilica (nanotitania) is deployed co-axially on the partially dried precipitate. The proposed deployment strategy allowed significant morphological differences when the deployment order of nanosilica and nanotitania were interchanged. Compositional analysis performed using EDX (Energy Dispersive X-ray spectroscopy) showed preferential deposition of nanosilica and nanotitania along the radial as well as the axial plane of the final deposit pattern. The underlying mechanism for such a phenomenon could be attributed to the contact line dynamics of a sessile double droplet. We also observe heteroaggregation of the nanosilica-nanotitania interaction along a narrow interface which resulted in nanotitania particles clustering into isolated islands embedded into a matrix of nanosilica particles. Overall, this work elucidates the evaporation driven dynamics of a mixed colloidal system which displays both macroscopic as well as microscopic phenomena. Such a system could be used to generate ordered arrays of functional materials with engineered micro to nano-scale properties.

Thickness measurement of dynamic thin liquid films generated by plug-annular flow in non-wetting microchannels

Surface tension forces are significant at millimeter length-scales, causing profoundly different flow morphologies in microchannels than in macroscale flows. The existence and morphology of thin liquid films is particularly relevant for predicting performance and operational stability of devices containing microscale two phase flows. Analytical, computational, and experimental methods previously employed in the study of thin liquid films are discussed. Thicknesses before and after a novel film morphology, referred to as a `shock,' are measured with a novel film thickness measurement technique that uses confocal microscopy. Film thicknesses predicted by previous work are compared to experimental results. Methods for increasing the accuracy of the confocal film thickness measurement technique are discussed.

On the dynamic contact angle in simulation of impinging droplets with sharp interface methods

Effects of dynamic contact angle models on the flow dynamics of an impinging droplet in sharp interface simulations are presented in this article. In the considered finite element scheme, the free surface is tracked using the arbitrary Lagrangian-Eulerian approach. The contact angle is incorporated into the model by replacing the curvature with the Laplace-Beltrami operator and integration by parts. Further, the Navier-slip with friction boundary condition is used to avoid stress singularities at the contact line. Our study demonstrates that the contact angle models have almost no influence on the flow dynamics of the non-wetting droplets. In computations of the wetting and partially wetting droplets, different contact angle models induce different flow dynamics, especially during recoiling. It is shown that a large value for the slip number has to be used in computations of the wetting and partially wetting droplets in order to reduce the effects of the contact angle models. Among all models, the equilibrium model is simple and easy to implement. Further, the equilibrium model also incorporates the contact angle hysteresis. Thus, the equilibrium contact angle model is preferred in sharp interface numerical schemes.

Oscillation dynamics of sessile droplets subjected to substrate vibration

Sessile droplets on a vibrating substrate are investigated focusing on axisymmetric oscillations with pinned contact line. Proper orthogonal decomposition is employed to identify the different modes of droplet shape oscillation and quantitatively assess the droplet oscillation and spectral response. We offer the first experimental evidence for the analogy of an oscillating sessile droplet with a non-linear spring mass damper system. The qualitative and quantitative agreement of amplitude response and phase response curves and limit cycles of the model dynamical system with that observed experimentally suggest that the bulk oscillations in the fundamental mode of a sessile droplet can be very well modeled by a Duffing oscillator with a hard spring, especially near the resonance. The red shift of the resonance peak with an increase in the glycerol concentration is clearly evidenced by both the experimental and predicted amplitude response curves. The influence of various operational parameters such as excitation frequency and amplitude and fluid properties on the droplet oscillation characteristics is adequately captured by the model. (C) 2014 Elsevier Ltd. All rights reserved.

Agglomeration front dynamics: Drying in sessile nano-particle laden droplets

The drying of sessile, nano-silica laden water droplet is studied under ambient conditions, in the absence of any convection. The drying process can be divided into two distinct regimes. During regime 1, the outer edge of the droplet remains pinned and particles agglomerate at the droplet periphery similar to the traditional coffee ring. However in regime 2, with further evaporation, both the liquid contact line and the agglomeration front starts moving radially inwards from the initial contact edge. The contact between the liquid and the agglomerate is maintained throughout regime 2 and the vaporisation driven liquid edge recession essentially drives the inward growth of the particle deposition. Fast kinetics of particle aggregation results in rapid growth of this agglomeration front as seen from the experiments. A theoretical formulation involving a simplistic model of the agglomeration front growth based on particle mass balance has been proposed. (C) 2014 Elsevier Ltd. All rights reserved,

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

The effect of flow field & geometry on the dynamic contact angle

A number of recent experiments suggest that, at a given wetting speed, the dynamic contact angle formed by an advancing liquid-gas interface with a solid substrate depends on the flow field and geometry near the moving contact line. In the present work, this effect is investigated in the framework of an earlier developed theory that was based on the fact that dynamic wetting is, by its very name, a process of formation of a new liquid-solid interface (newly “wetted” solid surface) and hence should be considered not as a singular problem but as a particular case from a general class of flows with forming or/and disappearing interfaces. The results demonstrate that, in the flow configuration of curtain coating, where a liquid sheet (“curtain”) impinges onto a moving solid substrate, the actual dynamic contact angle indeed depends not only on the wetting speed and material constants of the contacting media, as in the so-called slip models, but also on the inlet velocity of the curtain, its height, and the angle between the falling curtain and the solid surface. In other words, for the same wetting speed the dynamic contact angle can be varied by manipulating the flow field and geometry near the moving contact line. The obtained results have important experimental implications: given that the dynamic contact angle is determined by the values of the surface tensions at the contact line and hence depends on the distributions of the surface parameters along the interfaces, which can be influenced by the flow field, one can use the overall flow conditions and the contact angle as a macroscopic multiparametric signal-response pair that probes the dynamics of the liquid-solid interface. This approach would allow one to investigate experimentally such properties of the interface as, for example, its equation of state and the rheological properties involved in the interface’s response to an external torque, and would help to measure its parameters, such as the coefficient of sliding friction, the surface-tension relaxation time, and so on.