Simona Tonini

Single droplet impacts onto thin wall-films are a common phenomenon in many applications. For sufficiently high impact velocities, the droplet impact process consists of three phases, i.e., initial contact stage, droplet deformation with radial momentum transfer inducing an upward rising lamella, and crown propagation. Here, we present the results of a combined numerical and experimental study focusing on the early dynamics of the impact process. Specifically, the effects of the initial droplet shape, wall-film thickness, and contact line motion are analyzed. Prior to impact, an oblate spheroidal droplet shape was observed. Using direct numerical simulation, we show that the droplet shape affects the impact dynamics only during the first two phases, as it is one of the key parameter influencing the correct prediction of the impact zone. The contact line propagation is described by a square-root-time dependence R¯CL=ατ for both, dry and wetted surfaces. On dry surfaces, the advancement of the contact line is determined by the rolling motion of the truncated droplet. On wetted surfaces, the value of the α-parameter is controlled by two concurrent effects, namely, rolling motion and wall-film inertia. For impact onto thin films, the rolling motion prevails. With increasing wall-film height, the droplet penetrates into the soft substrates and wall-film inertia becomes the controlling factor. These insights into the early impact dynamics on wetted surface are important for the formulation of a unified modeling approach.

Pendant and sessile drops form a spherical cap only in the absence of gravity. The effect of gravity on drop shape is often neglected on the basis of the assumption that the drop size is smaller than the capillary length [Lc=(σ/gρ)1/2], although the deformation may not be fully negligible even in those cases. This paper focuses on evaluation of the effect that deformation due to gravity has on the evaporation characteristics of pendant and sessile drops. The drop shape is described by the Bashforth–Adams equation, a non-linear second order ordinary differential equation, which is solved numerically using a Runge–Kutta method with variable time steps. Under quasi-steady approximation, the species and energy conservation equations in the gas phase have analytical solutions, even for temperature-dependent gas thermophysical properties, once the solution of a basic Laplace problem is known. The Laplace equation is solved in axial symmetric geometry by using COMSOL Multiphysics®, for a wide range of drop sizes and contact angles, yielding vapor distribution, vapor fluxes, and evaporation rates. Comparison with the results from drops of same size in microgravity (i.e., having a spherical cap shape) shows that the effect is also perceptible for drops with a size smaller than the capillary length and that it can become quite important for those with larger sizes. Complementary results are found for sessile and pendant drops with respect to wall wettability, suggesting that the phenomenon can be analyzed using a unitary approach.