The effect of sharp solid edges on the droplet wettability

Water emerging from micrometer-sized pores into millimeter-sized gas-flow channels forms drops. The drops grow until the force from the flowing gas is sufficient to detach the drops as either (1) slugs that completely occlude the cross section of the channel and move at the superficial gas velocity, (2) drops that partially occlude the channel and move at a velocity that is less than the gas velocity, or (3) films that flow continuously, occluding part of the channel. At steady state, small residual water droplets, ∼100 μm in diameter, left in corners and on surface defects from previous drops, are key in determining the shape of water drops at detachment. Slugs are formed at low-gas-phase Reynolds numbers (ReG) in both hydrophilic and hydrophobic channels. Drops are shed in Teflon-coated hydrophobic channels for ReG > 30. Films are formed in acrylic hydrophilic channels for ReG > 30. Slugs form when growing drops encounter residual water droplets that nucleate the drop to slug transition. Drops are shed when the force exerted by the flowing gas on growing drops exceeds the force needed to advance the gas/liquid/solid contact line before they grow to the critical size for the drop to slug transition. Drops grow by "stick-slip" of the solid-liquid-gas contact lines and with pinned contact lines until the force on the drops results in either the downstream contact angle becoming greater than the dynamic advancing contact angle or the upstream contact angle becoming less than the dynamic receding contact angle. The upstream contact line never detaches for hydrophilic channels, which is why films form. The shape of water drops and the detachment energies are shown to be well approximated by the force balance between the force needed to advance the drop's contact lines and the force that the flowing gas exerts on a stationary drop.

We investigate the cavity formation during the impact of spheres and cylinders into a liquid pool by using a combination of experiments, simulations and theoretical analysis, with particular interest in contact-line pinning and its relation with the subsequent cavity evolution. The flows are simulated by a Navier–Stokes diffuse-interface solver that allows for moving contact lines. On the basis of agreement on experimentally measured quantities such as the position of the pinned contact line and the interface shape, we investigate flow details that are not accessible experimentally, identify the interface regions in the cavity formation and examine the geometric effects of impact objects. We connect wettability, inertia, geometry of the impact object, interface bending and contact-line position with the contact-line pinning by analysing the force balance at a pinned meniscus, and the result compares favourably with those from simulations and experiments. In addition to adjusting the interface bending, the object geometry also has a significant effect on the magnitude of low pressure in the liquid and the occurrence of flow separation. As a result, it is easier for an object with sharp edges to generate a cavity than a smooth object. A theoretical model based on the Rayleigh–Besant equation is developed to provide a quantitative description of the radial expansion of the cavity after the pinning of the contact line. The accuracy of the solution is greatly affected by the geometrical information on the interface connected to the pinned meniscus, showing the dependence of the global cavity dynamics on the local flows around the pinned contact line. Vertical ripple propagation on the cavity wall is found to follow the dispersion relation for the perturbation evolution on a hollow jet.

Evaporation of fluid from a pinned drop placed on solid surface proceeds via constant contact radius (CCR) mode, with a continuous reduction in the contact angle. The reduction of contact angle leads to an imbalance of interfacial tensions at the three-phase contact line. When the unbalanced force is sufficiently strong, the drop slips from the pinned contact line and slides inward. Depinning of the drop alters the mode of evaporation to constant contact angle (CCA) mode till it repins onto the surface. The change in evaporation mode from CCR to CCA is usually achieved by tuning the pinning energy barrier by controlling the surface properties of the substrate. Here, we demonstrate that the evaporation mode can be controlled by solely tailoring the surface tension of the drop, which is achieved in microgel particle-laden sessile drops that show spontaneous adsorption of microgels to the air/water interface, leading to a decrease in the interfacial tension. We show that droplets containing a sufficient number of microgels evaporate predominantly in CCR mode even on a hydrophobic surface, and the contact line remains pinned throughout the evaporation of the drop. Interestingly, the contact line dynamics can be controlled by tuning the softness of the microgels and the particle concentration in the drops.

The liquid drop captured at the capillary end, which is observed in capillary valve and pendant drop technique, is investigated theoretically and experimentally. Because of contact line pinning of the lower meniscus, the lower contact angle is able to rise from the intrinsic contact angle (θ*) so that the external force acting on the drop can be balanced by the capillary force. In the absence of contact angle hysteresis (CAH), the upper contact angle remains at θ*. However, in the presence of CAH, the upper contact angle can descend to provide more capillary force. The coupling between the lower and upper contact angles determines the equilibrium shape of the captured drop. In a capillary valve, the pinned contact line can move across the edge as the pressure difference exceeds the valving pressure, which depends on the geometrical characteristic and wetting property of the valve opening. When CAH is considered, the valving pressure is elevated because the capillary force is enhanced by the receding contact angle. For a pendant drop under gravity, the maximal capillary force is achieved as the lower contact angle reaches 180° in the absence of CAH. However, in the presence of CAH, four regimes can be identified by three critical drop volumes. The lower contact angle can exceed 180°, and therefore the drop takes on the shape of a light bulb, which does not exist in the absence of CAH. The comparisons between Surface Evolver simulations and experiments are quite well.

Control of the swelling, chemical functionalization, and adhesivity of hydrogels are finding new applications in a wide range of material systems. We investigate experimentally the effect of adsorbed particles on hydrogels on the depinning of contact lines. In our experiments, a water drop containing polystyrene microspheres is deposited on a swelling hydrogel, which leads to the drop absorption and particle deposition. Two regimes are observed: a decreasing drop height with a pinned contact line followed by a receding contact line. We show that increasing the particles concentration increases the duration of the first regime and significantly decreases the total absorption time. The adsorbed particles increase the pinning force at the contact line. Finally, we develop a method to measure the receding contact angle with the consideration of the hydrogel swelling.

In the present study, the effect of contact-line mobility on the natural oscillations of sessile drops is investigated through interface-resolved simulations. A hyperbolic tangent function is employed to model the variation of contact angle (CA) with contact-line (CL) velocity. The range of CL mobility considered is sufficient to recover the free contact-line (FCL) and pinned contact-line (PCL) limits. For finite CL mobility between these asymptotic limits, the CL can move, and the CA varies with CL velocity. The equilibrium CA and CL mobility are systematically varied to investigate their effects on oscillation frequency, energy dissipation, and oscillation decay rate. The results show that the oscillation frequency increases monotonically with CL mobility, transitioning from the FCL to the PCL limit. The CL dissipation is zero at the FCL and PCL limits and exhibits a non-monotonic dependence on CL mobility. The oscillation decay rate generally decreases with increasing equilibrium CA. However, due to the non-monotonic variation of CL dissipation with CL mobility, the decay rate initially increases and then decreases as CL mobility increases.

Using a continuum model capable of describing the one-component liquid-gas hydrodynamics down to the contact line scale, we carry out numerical simulation and physical analysis for the droplet motion driven by thermal singularity. For liquid droplets in one-component fluids on heated or cooled substrates, the liquid-gas interface is nearly isothermal. Consequently, a thermal singularity occurs at the contact line and the Marangoni effect due to temperature gradient is suppressed. Through evaporation or condensation in the vicinity of the contact line, the thermal singularity makes the contact angle increase with the increasing substrate temperature. This effect on the contact angle can be used to move the droplets on substrates with thermal gradients. Our numerical results for this kind of droplet motion are explained by a simple fluid dynamical model at the droplet length scale. Since the mechanism for droplet motion is based on the change of contact angle, a separation of length scales is exhibited through a comparison between the droplet motion induced by a wettability gradient and that by a thermal gradient. It is shown that the flow field at the droplet length scale is independent of the statics or dynamics at the contact line scale.

Nanoparticles deposition patterns in evaporating nanofluid droplets on smooth heated hydrophilic substrates: A 2D immersed boundary-lattice Boltzmann simulation

МОДЕЛИРОВАНИЕ ТУРБУЛЕНТНОЙ ЕСТЕСТВЕННОЙ КОНВЕКЦИИВ ЗАМКНУТЫХ ВЫТЯНУТЫХ ПО ВЫСОТЕ ОБЛАСТЯХ В.М.Головизнин 1 , И.А.Короткин 2 , С.А.Финогенов 2 1 МГУ им.М.В.Ломоносова, Москва, Российская Федерация 2 Институт проблем безопасного развития атомной энергетики РАН, Москва, Российская Федерация В предыдущей работе авторов (см.журнал «Вычислительная механика сплошных сред», 2015, Т. 8, № 1, С. 60-70) было показано, что использование вихреразрешающей, свободной от настроечных параметров схемы КАБАРЕ для решения как двумерной, так и трехмерной задачи Дэвиса приводит к удивительно хорошему совпадению результатов расчетов на очень грубых сетках (20×20 и 20×20×20) с результатами экспериментов и прецизионных расчетов для чисел Релея вплоть до значения 10 14 .Настоящая работа посвящена исследованию чувствительности этого феномена к изменению формы каверны (переходу каверны от кубической формы к прямоугольной).Рассмотрены области в форме параллелепипеда с аспектными отношениями 1:4, 1:10 и 1:28,6.Представлено сравнение данных, полученных методом КАБАРЕ, с экспериментально установленными (для отношения 1:28,6), а также с найденными в ходе прямого численного моделирования (для отношения 1:4) и вычисленными по эмпирической формуле (для отношения 1:10).Во всех случаях метод КАБАРЕ продемонстрировал хорошее согласование интегральных характеристик течения в каверне с данными других авторов.При этом в расчетах применялись показательно грубые сетки, сгущающиеся по мере приближения к стенкам области.Показано, что приемлемая точность расчетов на экстремально грубых сетках сохраняется вплоть до аспектного отношения 1:10.При более высоких

Evaluation of unusual change in contact angle between MgO and molten magnesium

A study of the structure and surface properties of nanostructured biocompatible coatings on Ti alloys

Pinning and depinning in imbibition beyond a sharp edge: A lattice Boltzmann study

This work focuses on the numerical study of electrohydrodynamic multiphase fluid flow problems. Liquid bodies experience deformations caused by strong external electric fields. These deformations can be used in industrial applications to accurately control the motion of single droplets. In the limiting case, strong electric fields can force atomisation of fine microdroplets from much larger liquid bulks. Due to differences in the electrical properties of the different fluid phases, the flow affects in turn the electric field distribution. The analysis of such coupled flows requires numerical modelling of the interdependent hydrodynamic and electrodynamic problems. Since most liquids exhibit some conductivity due to intrinsic ionic species and dissolved impurities, the electrodynamic problem must be modelled by an electroquasistatic model taking into account capacitive, resistive and convective electrical currents. This electroquasistatic problem is coupled to an incompressible fluid flow problem described by the Navier-Stokes equations. Both problems are solved on the same computational grid, using the finite volume method. The fluid-fluid interface is modelled using the volume of fluid method. The resulting diffuse interface captures the motion of the phase boundaries while efficiently handling topology changes. The motion of contact lines is furthermore modelled using a dynamic contact angle model including hysteresis effects. Pinned contact lines and stick-slip contact line motion in transient problems can therefore be represented. The developed solver is readily applicable to a large range of electrohydrodynamic flow problems. This work investigates electrohydrodynamic flows occurring in three technical applications. First, the dynamics of sessile droplets subjected to an AC electric field on the surface of an insulator are considered. The dynamics of the oscillating droplets are compared with experimental data. Partial discharge inception fields are then estimated for similar configurations. Secondly, the detachment dynamics of two liquids in an on-demand droplet generator, where droplet detachment is enforced by electric fields, are considered and compared to experimental data. The different conductivities of the liquids are shown to lead to substantially different detachment dynamics. Relevant parameters in the detachment dynamics are extracted from the simulations. Finally, transient electrosprays in the cone-jet mode are characterised for a number of liquids with different electromechanical properties. The charge-radius correlations of the first ejections are found to obey the scaling laws reported in the literature. Moreover, additional scaling laws are found for the subsequent ejections.

We report direct atomic-force-microscope measurements of capillary force hysteresis (CFH) and relaxation of a circular moving contact line (CL) formed on a long micron-sized hydrophobic fiber intersecting a liquid-air interface. By using eight different liquid interfaces with varying solid-liquid molecular interactions, we find a universal behavior of the asymmetric speed dependence of CFH and CL relaxation. A unified model based on force-assisted barrier crossing is used to connect the mesoscopic measurements of CFH and CL relaxation with the energy barrier height E_{b} and size λ associated with the surface defects. The experiment demonstrates that the CL pinning (relaxation) and depinning dynamics are closely related and can be described by a common microscopic framework.

A comprehensive numerical and experimental investigation on micrometer-sized water droplet impact dynamics and evaporation on an unheated, flat, dry surface is conducted from the standpoint of spray-cooling technology. The axisymmetric time-dependent governing equations of continuity, momentum, energy, and species are solved. Surface tension, wall adhesion effect, gravitational body force, contact line dynamics, and evaporation are accounted for in the governing equations. The explicit volume of fluid (VOF) model with dynamic meshing and variable-time stepping in serial and parallel processors is used to capture the time-dependent liquid-gas interface motion throughout the computational domain. The numerical model includes temperature- and species-dependent thermodynamic and transport properties. The contact line dynamics and the evaporation rate are predicted using Blake's and Schrage's molecular kinetic models, respectively. An extensive grid independence study was conducted. Droplet impingement and evaporation data are acquired with a standard dispensing/imaging system and high-speed photography. The numerical results are compared with measurements reported in the literature for millimeter-size droplets and with current microdroplet experiments in terms of instantaneous droplet shape and temporal spread (R/D(0) or R/R(E)), flatness ratio (H/D(0)), and height (H/H(E)) profiles, as well as temporal volume (inverted A) profile. The Weber numbers (We) for impinging droplets vary from 1.4 to 35.2 at nearly constant Ohnesorge number (Oh) of approximately 0.025-0.029. Both numerical and experimental results show that there is air bubble entrapment due to impingement. Numerical results indicate that Blake's formulation provides better results than the static (SCA) and dynamic contact angle (DCA) approach in terms of temporal evolution of R/D(0) and H/D(0) (especially at the initial stages of spreading) and equilibrium flatness ratio (H(E)/D(0)). Blake's contact line dynamics is dependent on the wetting parameter (K(W)). Both numerical and experimental results suggest that at 4.5 < We < 11.0 the short-time dynamics of microdroplet impingement corresponds to a transition regime between two different spreading regimes (i.e., for We < or = 4.5, impingement is followed by spreading, then contact line pinning and then inertial oscillations, and for We > or = 11.0, impingement is followed by spreading, then recoiling, then contact line pinning and then inertial oscillations). Droplet evaporation can be satisfactorily modeled using the Schrage model, since it predicts both well-defined transient and quasi-steady evaporation stages. The model compares well with measurements in terms of flatness ratio (H/H(E)) before depinning occurs. Toroidal vortices are formed on the droplet surface in the gaseous phase due to buoyancy-induced Rayleigh-Taylor instability that enhances convection.