Variational approach to contact line dynamics for thin films

This paper investigates a variational approach to viscous flows with contact line dynamics based on energy-dissipation modeling. The corresponding model is reduced to a thin-film equation and its variational structure is also constructed and discussed. Feasibility of this modeling approach is shown by constructing a numerical scheme in 1D and by computing numerical solutions for the problem of gravity driven droplets. Some implications of the contact line model are highlighted in this setting

Numerics of thin-film free boundary problems for partial wetting

We present a novel framework to solve thin-film equations with an explicit non-zero contact angle, where the support of the solution is treated as an unknown. The algorithm uses a finite element method based on a gradient formulation of the thin-film equations coupled to an arbitrary Lagrangian-Eulerian method for the motion of the support. Features of this algorithm are its simplicity and robustness. We apply this algorithm in 1D and 2D to problems with surface tension, contact angles and with gravity

Variational approach to contact line dynamics for thin films

This paper investigates a variational approach to viscous flows with contact line dynamics based on energy-dissipation modeling. The corresponding model is reduced to a thin-film equation and its variational structure is also constructed and discussed. Feasibility of this modeling approach is shown by constructing a numerical scheme in 1D and by computing numerical solutions for the problem of gravity driven droplets. Some implications of the contact line model are highlighted in this setting

Two-phase flows for sedimentation of suspensions

We present a two-phase flow model that arises from energetic-variational arguments and study its implication for the sedimentation of buoyant particles in a viscous fluid inside a Hele--Shaw cell and also compare corresponding simulation results to experiments. Based on a minimal dissipation argument, we provide a simplified 1D model applicable to sedimentation and study its properties and the numerical discretization. We also explore different aspects of its numerical discretization in 2D. The focus is on different possible stabilization techniques and their impact on the qualitative behavior of solutions. We use experimental data to verify some first qualitative model predictions and discuss these experiments for different stages of batch sedimentation

Model hierarchies and higher-order discretisation of time-dependent thin-film free boundary problems with dynamic contact angle

We present a mathematical and numerical framework for the physical problem of thin-film fluid flows over planar surfaces including dynamic contact angles. In particular, we provide algorithmic details and an implementation of higher-order spatial and temporal discretisation of the underlying free boundary problem using the finite element method. The corresponding partial differential equation is based on a thermodynamic consistent energetic variational formulation of the problem using the free energy and viscous dissipation in the bulk, on the surface, and at the moving contact line. Model hierarchies for limits of strong and weak contact line dissipation are established, implemented and studied. We analyze the performance of the numerical algorithm and investigate the impact of the dynamic contact angle on the evolution of two benchmark problems: gravity-driven sliding droplets and the instability of a ridge

Two-phase flows for sedimentation of suspensions

We present a two-phase flow model that arises from energetic-variational arguments and study its implication for the sedimentation of buoyant particles in a viscous fluid inside a Hele--Shaw cell and also compare corresponding simulation results to experiments. Based on a minimal dissipation argument, we provide a simplified 1D model applicable to sedimentation and study its properties and the numerical discretization. We also explore different aspects of its numerical discretization in 2D. The focus is on different possible stabilization techniques and their impact on the qualitative behavior of solutions. We use experimental data to verify some first qualitative model predictions and discuss these experiments for different stages of batch sedimentation

Sharp-interface limits of Cahn--Hilliard models and mechanics with moving contact lines

We construct gradient structures for free boundary problems with moving capillary interfaces with nonlinear (hyper)elasticity and study the impact of moving contact lines. In this context, we numerically analyze how phase-field models converge to certain sharp-interface models when the interface thickness tends to zero. In particular, we study the scaling of the Cahn--Hilliard mobility with certain powers of the interfacial thickness. In the presence of interfaces, it is known that the intended sharp-interface limit holds only for a particular range of powers However, in the presence of moving contact lines we show that some scalings that are valid for interfaces produce significant errors and the effective range of valid powers of the interfacial thickness in the mobility reduces

Droplets on liquids and their long way into equilibrium

The morphological paths towards equilibrium droplets during the late stages of the dewetting process of a liquid film from a liquid substrate is investigated experimentally and theoretically. As liquids, short chained polystyrene (PS) and polymethyl-methacrylate (PMMA) are used, which can be considered as Newontian liquids well above their glass transition temperatures. Careful imaging of the PS/air interface of the droplets during equilibration by \emph{in situ} scanning force microscopy and the PS/PMMA interface after removal of the PS droplets reveal a surprisingly deep penetration of the PS droplets into the PMMA layer. Droplets of sufficiently small volumes develop the typical lens shape and were used to extract the ratio of the PS/air and PS/PMMA surface tensions and the contact angles by comparison to theoretical exact equilibrium solutions of the liquid/liquid system. Using these results in our dynamical thin-film model we find that before the droplets reach their equilibrium they undergo several intermediate stages each with a well-defined signature in shape. Moreover, the intermediate droplet shapes are independent of the details of the initial configuration, while the time scale they are reached depend strongly on the droplet volume. This is shown by the numerical solutions of the thin-film model and demonstrated by quantitative comparison to experimental results

Impact of energy dissipation on interface shapes and on rates for dewetting from liquid substrates

We revisit the fundamental problem of liquid-liquid dewetting and perform a detailed comparison of theoretical predictions based on thin-film models with experimental measurements obtained by atomic force microscopy (AFM). Specifically, we consider the dewetting of a liquid polystyrene (PS) layer from a liquid polymethyl methacrylate (PMMA) layer, where the thicknesses and the viscosities of PS and PMMA layers are similar. The excellent agreement of experiment and theory reveals that dewetting rates for such systems follow no universal power law, in contrast to dewetting scenarios on solid substrates. Our new energetic approach allows to assess the physical importance of different contributions to the energy-dissipation mechanism, for which we analyze the local flow fields and the local dissipation rates.Comment: 15 pages, 5 figure

Stationary solutions of liquid two-layer thin film models

We investigate stationary solutions of a thin-film model for liquid two-layer flows in an energetic formulation that is motivated by its gradient flow structure. The goal is to achieve a rigorous understanding of the contact-angle conditions for such two-layer systems. We pursue this by investigating a corresponding energy that favors the upper liquid to dewet from the lower liquid substrate, leaving behind a layer of thickness $h_*$. After proving existence of stationary solutions for the resulting system of thin-film equations we focus on the limit $h_*\to 0$ via matched asymptotic analysis. This yields a corresponding sharp-interface model and a matched asymptotic solution that includes logarithmic switch-back terms. We compare this with results obtained using $\Gamma$-convergence, where we establish existence and uniqueness of energetic minimizers in that limit