The dynamics and stability of thin liquid films have fascinated scientists over many decades: the observations of regular wave patterns in film flows down a windowpane or along guttering, the patterning of dewetting droplets, and the fingering of viscous flows down a slope are all examples that are familiar in daily life. Thin film flows occur over a wide range of length scales and are central to numerous areas of engineering, geophysics, and biophysics; these include nanofluidics and microfluidics, coating flows, intensive processing, lava flows, dynamics of continental ice sheets, tear-film rupture, and surfactant replacement therapy. These flows have attracted considerable attention in the literature, which have resulted in many significant developments in experimental, analytical, and numerical research in this area. These include advances in understanding dewetting, thermocapillary- and surfactant-driven films, falling films and films flowing over structured, compliant, and rapidly rotating substrates, and evaporating films as well as those manipulated via use of electric fields to produce nanoscale patterns. These developments are reviewed in this paper and open problems and exciting research avenues in this thriving area of fluid mechanics are also highlighted.

Dynamics and stability of thin liquid films

Computational methods have been extended recently to allow for the presence of moving contact lines in simulated two-phase flows. The predictive capability offered by these methods is potentially large, joining theoretical and experimental methods. Several challenges rather unique to this area need to be overcome, however, notably regarding the conditions near a moving contact line and the very large separation of length scales in these flows. We first summarize the main models for moving contact lines and follow with an overview of computational methods that includes direct continuum approaches and macroscale models that resolve only the large-scale flow by modeling the effects of the conditions near the contact line using theory. Results are presented for contact-line motion on ideal as well as patterned and grooved surfaces and for extensions to account for complexities such as thermocapillarity and phase change.

Numerical Simulations of Flows with Moving Contact Lines

Numerical simulations of 3D flows with moving contact lines

Self-propulsion of liquid marbles filled with aqueous alcohol solutions and placed on a water surface is reported. The characteristic velocity of the marbles is ca. 0.1 m/s. The phenomenon of self-propulsion is related to the Marangoni solutocapillary flow caused by the condensation of alcohol, evaporated from the liquid marble, on a water surface. The Marangoni flow in turn enhances the evaporation of alcohol from the marbles. Addition of alcohol to the water supporting the marbles suppresses the self-propulsion. The propulsion of liquid marbles is mainly stopped by water drag. The velocity of the center of mass of the marbles increases with an increase in the concentration of alcohol within the marble. The self-propulsion velocity of marbles is independent of their volume. Self-propulsion was not observed when marbles were placed on a glycerol surface.

Self-Propulsion of Liquid Marbles: Leidenfrost-like Levitation Driven by Marangoni Flow

We study the dynamics of a colloidal fluid including inertia and hydrodynamic interactions, two effects which strongly influence the nonequilibrium properties of the system. We derive a general dynamical density functional theory which shows very good agreement with full Langevin dynamics. In suitable limits, we recover existing dynamical density functional theories and a Navier-Stokes-like equation with additional nonlocal terms.

/pdf/general-dynamical-density-functional-theory-for-classical-2a9xjpljp8.pdf

General Dynamical Density Functional Theory for Classical Fluids

We investigate the interplay between a stable horizontal thin liquid film on a solid substrate and an excitable or bistable reactive mixture on its free surface. Their coupling is twofold. On the one hand, flow in the film transports the reacting surfactants convectively. On the other hand, gradients in the surfactant concentration exert Marangoni stresses on the free surface of the film. A reduced model is derived based on the long-wave approximation. We analyze the linear stability of the coupled system as well as the nonlinear behavior, including the propagation of solitary waves, fronts, and pulses. We show, for instance, that the coupling of thin film hydrodynamics and surfactant chemistry can either stabilize instabilities occurring in the pure chemical system, or in a regime where the pure hydrodynamic and chemical subsystems are both stable, the coupling can induce instabilities.

Dynamics of a horizontal thin liquid film in the presence of reactive surfactants

We investigate the dynamics of a thin liquid film on an inclined planar substrate in the presence of an insoluble surfactant on its free surface. We consider both the linear and nonlinear regimes. The linear regime is examined through the Orr-Sommerfeld eigenvalue problem of the full Navier-Stokes and concentration equations and wall and free-surface boundary conditions. The nonlinear regime is investigated through two different models. The first one is obtained from the classical long-wave expansion and the second one through an integral-boundary-layer approximation combined with a simple Galerkin projection. Although accurate close to the instability threshold, the first model fails to describe the dynamics of the system far from criticality. On the other hand, the second model not only captures accurately the behavior close to the instability threshold, but is also valid far from criticality. Analytical and numerical results on the role of the surfactant on the free-surface dynamics are presented. In the linear regime, the Marangoni stresses induced by the surfactant reduce the domain of instability for the base flow. In the nonlinear regime, the system evolves into solitary pulses for both the free surface and surfactant concentration. The amplitude and velocity of these pulses decrease as the Marangoni effect becomes stronger.

Dynamics of a falling film with solutal Marangoni effect.

We investigate the equilibrium of a fluid in contact with a solid boundary through a density-functional theory. Depending on the conditions, the fluid can be in one phase, gas or liquid, or two phases, while the wall induces an external field acting on the fluid particles. We first examine the case of a liquid film in contact with the wall. We construct bifurcation diagrams for the film thickness as a function of the chemical potential. At a specific value of the chemical potential, two equally stable films, a thin one and a thick one, can coexist. As saturation is approached, the thickness of the thick film tends to infinity. This allows the construction of a liquid–gas interface that forms a well-defined contact angle with the wall.

Equilibrium gas-liquid-solid contact angle from density-functional theory

We investigate the equilibrium of a fluid in contact with a solid boundary through a density-functional theory. Depending on the conditions, the fluid can be in one phase, gas or liquid, or two phases, while the wall induces an external field acting on the fluid particles. We first examine the case of a liquid film in contact with the wall. We construct bifurcation diagrams for the film thickness as a function of the chemical potential. At a specific value of the chemical potential, two equally stable films, a thin one and a thick one, can coexist. As saturation is approached, the thickness of the thick film tends to infinity. This allows the construction of a liquid-gas interface that forms a well defined contact angle with the wall.

We present a derivation of the time evolution equation of a physical quantity defined as an integral of a density over a moving surface. The equation can be readily used to obtain the transport equation for an interfacial entity such as surfactants on the surface of a single-valued interface as in thin films. By utilizing the concept of the extension field of the density we generalize the transport equation to the case of multi-valued interfaces.

Antonio Pereira

Papers

Dynamics of a horizontal thin liquid film in the presence of reactive surfactants

Dynamics of a falling film with solutal Marangoni effect.

Equilibrium gas-liquid-solid contact angle from density-functional theory

Equilibrium gas-liquid-solid contact angle from density-functional theory

On the transport equation for an interfacial quantity