skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Thermodynamically consistent phase-field modelling of contact angle hysteresis
In the phase-field description of moving contact line problems, the two-phase system can be described by free energies, and the constitutive relations can be derived based on the assumption of energy dissipation. In this work we propose a novel boundary condition for contact angle hysteresis by exploring wall energy relaxation, which allows the system to be in non-equilibrium at the contact line. Our method captures pinning, advancing and receding automatically without the explicit knowledge of contact line velocity and contact angle. The microscopic dynamic contact angle is computed as part of the solution instead of being imposed. Furthermore, the formulation satisfies a dissipative energy law, where the dissipation terms all have their physical origin. Based on the energy law, we develop an implicit finite element method that is second order in time. The numerical scheme is proven to be unconditionally energy stable for matched density and zero contact angle hysteresis, and is numerically verified to be energy dissipative for a broader range of parameters. We benchmark our method by computing pinned drops and moving interfaces in the plane Poiseuille flow. When the contact line moves, its dynamics agrees with the Cox theory. In the test case of oscillating drops, the contact line transitions smoothly between pinning, advancing and receding. Our method can be directly applied to three-dimensional problems as demonstrated by the test case of sliding drops on an inclined wall.  more » « less
Award ID(s):
1522604
PAR ID:
10351804
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
899
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The energetics of drop deposition are considered in the capillary-ballistic regime characterized by high Reynolds number and moderate Weber number. Experiments are performed impacting water/glycol drops onto substrates with varying wettability and contact-angle hysteresis. The impacting event is decomposed into three regimes: (i) pre-impact, (ii) inertial spreading and (iii) post contact-line (CL) pinning, conveniently framed using the theory of Dussan & Davis ( J. Fluid Mech. , vol. 173, 1986, pp. 115–130). During fast-time-scale inertial spreading, the only form of dissipation is CL dissipation ( $$\mathcal {D}_{CL}$$ ). High-speed imaging is used to resolve the stick-slip dynamics of the CL with $$\mathcal {D}_{CL}$$ measured directly from experiment using the $$\Delta \alpha$$ - $$R$$ cyclic diagram of Xia & Steen ( J. Fluid Mech. , vol. 841, 2018, pp. 767–783), representing the contact-angle deviation against the CL radius. Energy loss occurs on slip legs, and this observation is used to derive a closed-form expression for the kinetic K and interfacial $$\mathcal{A}$$ post-pinning energy $$\{K+\mathcal {A}\}_p/\mathcal {A}_o$$ independent of viscosity, only depending on the rest angle $$\alpha _p$$ , equilibrium angle $$\bar {\alpha }$$ and hysteresis $$\Delta \alpha$$ , which agrees well with experimental observation over a large range of parameters, and can be used to evaluate contact-line dissipation during inertial spreading. The post-pinning energy is found to be independent of the pre-impact energy, and it is broken into modal components with corresponding energy partitioning approximately constant for low-hysteresis surfaces with fixed pinning angle $$\alpha _p$$ . During slow-time-scale post-pinning, the liquid/gas ( $lg$ ) interface is found to vibrate with the frequencies and mode shapes predicted by Bostwick & Steen ( J. Fluid Mech. , vol. 760, 2014, pp. 5–38), irrespective of the pre-impact energy. Resonant mode decay rates are determined experimentally from fast Fourier transforms of the interface dynamics. 
    more » « less
  2. We studied the evolution of capillary bridges between nominally flat plates undergoing multiple cycles of compression and stretching in experiments and simulations. We varied the distance between the plates in small increments to study the full evolution of the bridge shape. Experiments show that contact angle hysteresis determines the shape of the bridge. In sliding drops, hysteresis can be modeled using a contact angle-dependent resistive force F̃R applied at the contact line. We developed a model that accurately captures the evolution of the bridge shape by combining F̃R and constrained energy minimization. Unlike previous work, this allows for both complete and partial contact line pinning. We also explored the effect of using nonparallel plates. The asymmetry in the bridge shape causes the movement of the center of mass of the bridge and can be explained by contact angle hysteresis. We find that even a slight misalignment between the flat plates can have a measurable effect. 
    more » « less
  3. Abstract Inkjet printing of electronic materials is of interest for digital printing of flexible electronics and sensors, but the width of the inkjet-printed lines is still large, limiting device size and performance. Decreasing the drop volume, increasing the drop spacing, and increasing the ink-substrate contact angle are all approaches by which the line width can be lowered, however these approaches are limited by the nozzle geometry, ink coalescence and bead instabilities, and contact angle hysteresis, respectively. Here we demonstrate a novel approach for stable inkjet printing of very narrow lines on ink-substrate combinations with a high contact angle, utilizing the de-wetting of the ink due to the decreased contact angle hysteresis. After printing and drying an initial layer of disconnected seed drops of silver nanoparticle ink, we print an additional layer of bridging drops of the same ink in between the dried seed drops. The bridging drops expand to touch the dried seed drops and then retract into a line, due to the pinning of the wet ink on the dried seed ink but not on the substrate, forming a continuous silver trace. The trace width is decreased from 60μm with a traditional printing approach down to 12.6μm with this seed-bridge approach. The electrical conductivity of the silver trace is similar to that of a conventionally printed trace. Due to poor adhesion on the print substrate, the trace was transferred to a separate polymer substrate with a simple hot-pressing procedure, which preserves the electrical conductivity of the trace. 
    more » « less
  4. Contact line dynamics is crucial in determining the deposition patterns of evaporating colloidal droplets. Using high-speed interferometry, we directly observe the stick-slip motion of the contact line in situ and are able to resolve the instantaneous shape of the inkjet-printed, evaporating pico-liter drops containing nanoparticles of varying wettability. Integrated with post-mortem optical profilometry of the deposition patterns, the instantaneous particle volume fraction and hence the particle deposition rate can be determined. The results show that the stick-slip motion of the contact line is a strong function of the particle wettability. While the stick-slip motion is observed for nanoparticles that are less hydrophilic ( i.e. , particle contact angle θ ≈ 74° at the water–air interface), which results in a multiring deposition, a continuous receding of the contact line is observed for more hydrophilic nanoparticles ( i.e. , θ ≈ 34°), which leaves a single-ring pattern. A model is developed to predict the number of particles required to pin the contact line based on the force balance of the hydrodynamic drag, interparticle interactions, and surface tension acting on the particles near the contact line with varying particle wettability. A three-fold increase in the number of particles required for pinning is predicted when the particle wettability increases from the wetting angle of θ ≈ 74° to θ ≈ 34°. This finding explains why particles with greater wettability form a single-ring pattern and those with lower wettability form a multi-ring pattern. In addition, the particle deposition rate is found to depend on the particle wettability and vary with time. 
    more » « less
  5. Abstract The synthesis of nanostructured surfaces via block copolymer (BCP) self-assembly enables a precise control of the surface feature shape within a range of dimensions of the order of tens of nanometers. This work studies how to exploit this ability to control the wetting hysteresis and liquid adhesion forces as the substrate undergoes chemical aging and changes in its intrinsic wettability. Via BCP self-assembly we fabricate nanostructured surfaces on silicon substrates with a hexagonal array of regular conical pillars having a fixed period (52 nm) and two different heights (60 and 200 nm), which results in substantially different lateral and top surface areas of the nanostructure. The wetting hysteresis of the fabricated surfaces is characterized using force–displacement measurements under quasistaic conditions and over sufficiently long periods of time for which the substrate chemistry and surface energy, characterized by the Young contact angle, varies significantly. The experimental results and theoretical analysis indicate that controlling the lateral and top area of the nanostructure not only controls the degree of wetting hysteresis but can also make the advancing and receding contact angles less susceptible to chemical aging. These results can help rationalize the design of nanostructured surfaces for different applications such as self-cleaning, enhanced heat transfer, and drag reduction in micro/nanofluidic devices. 
    more » « less