Air Entrainment in Dynamic Wetting: Knudsen Effects and the Influence of Ambient Air Pressure

Recent experiments on coating flows and liquid drop impact both demonstrate that wetting failures caused by air entrainment can be suppressed by reducing the ambient gas pressure. Here, it is shown that non-equilibrium effects in the gas can account for this behaviour, with ambient pressure reductions increasing the gas' mean free path and hence the Knudsen number $Kn$. These effects first manifest themselves through Maxwell slip at the gas' boundaries so that for sufficiently small $Kn$ they can be incorporated into a continuum model for dynamic wetting flows. The resulting mathematical model contains flow structures on the nano-, micro- and milli-metre scales and is implemented into a computational platform developed specifically for such multiscale phenomena. The coating flow geometry is used to show that for a fixed gas-liquid-solid system (a) the increased Maxwell slip at reduced pressures can substantially delay air entrainment, i.e. increase the `maximum speed of wetting', (b) unbounded maximum speeds are obtained as the pressure is reduced only when slip at the gas-liquid interface is allowed for and (c) the observed behaviour can be rationalised by studying the dynamics of the gas film in front of the moving contact line. A direct comparison to experimental results obtained in the dip-coating process shows that the model recovers most trends but does not accurately predict some of the high viscosity data at reduced pressures. This discrepancy occurs because the gas flow enters the `transition regime', so that more complex descriptions of its non-equilibrium nature are required. Finally, by collapsing onto a master curve experimental data obtained for drop impact in a reduced pressure gas, it is shown that the same physical mechanisms are also likely to govern splash suppression phenomena.Comment: Accepted for publication in the Journal of Fluid Mechanic

The Dynamics of Liquid Drops Coalescing in the Inertial Regime

We examine the dynamics of two coalescing liquid drops in the `inertial regime', where the effects of viscosity are negligible and the propagation of the bridge front connecting the drops can be considered as `local'. The solution fully computed in the framework of classical fluid-mechanics allows this regime to be identified and the accuracy of the approximating scaling laws proposed to describe the propagation of the bridge to be established. It is shown that the scaling law known for this regime has a very limited region of accuracy and, as a result, in describing experimental data it has frequently been applied outside its limits of applicability. The origin of the scaling law's shortcoming appears to be the fact that it accounts for the capillary pressure due only to the longitudinal curvature of the free surface as the driving force for the process. To address this deficiency, the scaling law is extended to account for both the longitudinal and azimuthal curvatures at the bridge front which, fortuitously, still results in an explicit analytic expression for the front's propagation speed. This new expression is then shown to offer an excellent approximation for both the fully-computed solution and for experimental data from a range of flow configurations for a remarkably large proportion of the coalescence process. The derived formula allows one to predict the speed at which drops coalesce for the duration of the inertial regime which should be useful for the analysis of experimental data.Comment: Accepted for publication in Physical Review

Kinetic effects in dynamic wetting

The maximum speed at which a liquid can wet a solid is limited by the need to displace gas lubrication films in front of the moving contact line. The characteristic height of these films is often comparable to the mean free path in the gas so that hydrodynamic models do not adequately describe the flow physics. This Letter develops a model which incorporates kinetic effects in the gas, via the Boltzmann equation, and can predict experimentally-observed increases in the maximum speed of wetting when (a) the liquid’s viscosity is varied, (b) the ambient gas pressure is reduced or (c) the meniscus is confined

Dynamics of liquid nano-threads : fluctuation-driven instability and rupture

The instability and rupture of nanoscale liquid threads is shown to strongly depend on thermal fluctuations. These fluctuations are naturally occurring within molecular dynamics (MD) simulations and can be incorporated via fluctuating hydrodynamics into a stochastic lubrication equation (SLE). A simple and robust numerical scheme is developed for the SLE that is validated against MD for both the initial (linear) instability and the nonlinear rupture process. Particular attention is paid to the rupture process and its statistics, where the `double-cone’ profile reported by Moseler & Landmann [Science, 2000, 289(5482): 1165-1169] is observed, as well as other distinct profile forms depending on the flow conditions. Comparison to the Eggers’ similarity solution [Physical Review Letters, 2002, 89(8): 084502], a power law of the minimum thread radius against time to rupture, shows agreement only at low surface tension; indicating that surface tension cannot generally be neglected when considering rupture dynamics

Capillary breakup of a liquid bridge : identifying regimes and transitions

Computations of the breakup of a liquid bridge are used to establish the limits of applicability of similarity solutions derived for different breakup regimes. These regimes are based on particular viscous-inertial balances, that is different limits of the Ohnesorge number Oh. To accurately establish the transitions between regimes, the minimum bridge radius is resolved through four orders of magnitude using a purpose-built multiscale finite element method. This allows us to construct a quantitative phase diagram for the breakup phenomenon which includes the appearance of a recently discovered low-Oh viscous regime. The method used to quantify the accuracy of the similarity solutions allows us to identify a number of previously unobserved features of the breakup, most notably an oscillatory convergence towards the viscous-inertial similarity solution. Finally, we discuss how the new findings open up a number of challenges for both theoretical and experimental analysis

The Formation of a Bubble from a Submerged Orifice

The formation of a single bubble from an orifice in a solid surface, submerged in an in- compressible, viscous Newtonian liquid, is simulated. The finite element method is used to capture the multiscale physics associated with the problem and to track the evolution of the free surface explicitly. The results are compared to a recent experimental analysis and then used to obtain the global characteristics of the process, the formation time and volume of the bubble, for a range of orifice radii; Ohnesorge numbers, which combine the material parameters of the liquid; and volumetric gas flow rates. These benchmark calculations, for the parameter space of interest, are then utilised to validate a selection of scaling laws found in the literature for two regimes of bubble formation, the regimes of low and high gas flow rates.Comment: Accepted for publication in the European Journal of Mechanics B/Fluid

Bouncing off the walls : the influence of gas-kinetic and van der Waals effects in drop impact

A model is developed for liquid drop impact on a solid surface that captures the thin film gas flow beneath the drop, even when the film’s thickness is below the mean free path in the gas so that gas kinetic effects (GKE) are important. Simulation results agree with experiments, with the impact speed threshold between bouncing and wetting reproduced to within 5 least 50 mapped and provides experimentally verifiable predictions. There are two principal modes of contact leading to wetting and both are associated with a van der Waals driven instability of the film

Revisiting the Rayleigh-Plateau instability for the nanoscale

The theoretical framework developed by Rayleigh and Plateau in the 19th century has been remarkably accurate in describing macroscale experiments of liquid cylinder instability. Here we re-evaluate and revise the Rayleigh-Plateau instability for the nanoscale, where molecular dynamics experiments demonstrate its inadequacy. A new framework based on the stochastic lubrication equation is developed that captures nanoscale flow features and highlights the critical role of thermal fluctuations at small scales. Remarkably, the model indicates that classically stable (i.e. ‘fat’) liquid cylinders can be broken at the nanoscale, and this is confirmed by molecular dynamics

Lifetime of a nanodroplet : kinetic effects and regime transitions

A transition from a d2 to a d law is observed in molecular dynamics (MD) simulations when the diameter (d) of an evaporating droplet reduces to the order of the vapor’s mean free path; this cannot be explained by classical theory. This Letter shows that the d law can be predicted within the Navier-Stokes-Fourier (NSF) paradigm if a temperature-jump boundary condition derived from kinetic theory is utilized. The results from this model agree with those from MD in terms of the total lifetime, droplet radius, and temperature, while the classical d2 law underpredicts the lifetime of the droplet by a factor of 2. Theories beyond NSF are also employed in order to investigate vapor rarefaction effects within the Knudsen layer adjacent to the interface

Molecular simulation of thin liquid films : thermal fluctuations and instability

The instability of a thin liquid film on a solid surface is studied both by molecular dynamics simulations (MD) and a stochastic thin-film equation (STF), which models thermal fluctuations with white noise. A linear stability analysis of the STF allows us to derive a power spectrum for the surface fluctuations, which is quantitatively validated against the spectrum observed in MD. Thermal fluctuations are shown to be critical to the dynamics of nanoscale films. Compared to the classical instability mechanism, which is driven by disjoining pressure, fluctuations (a) can massively amplify the instability, (b) cause the fluctuation wavelength that is dominant to evolve in time (a single fastest-growing mode does not exist), and (c) decrease the critical wavelength so that classically stable films can be ruptured