Do liquid drops roll or slide on inclined surfaces?

Self-Accelerating Drops on Silicone-Based Super Liquid-Repellent Surfaces

Design of super liquid-repellent surfaces has relied on an interplay between surface topography and surface energy. Perfluoroalkylated materials are often used, but they are environmentally unsustainable and notorious for building up static charge. Therefore, there is a need for understanding the performance of sustainable low surface energy materials with antistatic properties. Here, we explore drop interactions with perfluoroalkyl- and silicone-based surfaces, focusing on three modes of drop-to-surface interactions. The behavior of drops rolling under gravity is compared to those subjected to lateral and normal forces under constant slide (i.e., friction) and detachment (i.e., adhesion) velocities. We demonstrate that a drop's characteristic and dynamic mobility depends on surface chemistry, with sequential drop interactions being particularly affected. By utilizing force-and-charge instruments, we show how rolling drops are primarily governed by adhesion and its associated electrostatic effects, instead of friction. Perfluoroalkylated surfaces continuously accumulate charges, while silicone surfaces rapidly saturate. Consequently, sequentially contacting drops accumulate significant charges on the former while rapidly diminishing on the latter. The drop charge suppressing behavior of silicones enhances drop mobility despite their higher surface energy compared to perfluoroalkyls. Quantum mechanical density functional theory calculations show significant differences in surface charge distributions at the atomic level. Simulations suggest that variations in the lifetimes of surface hydroxyl ions likely drive the markedly different drop charging behaviors. Our findings demonstrate the critical role of surface chemistry and its coupled electrostatics in drop mobility, providing valuable insights for designing environmentally friendly, antistatic, super liquid-repellent surfaces.

Dynamics of Droplets Moving on Lubricated Polymer Brushes

Understanding the dynamics of drops on polymer-coated surfaces is crucial for optimizing applications such as self-cleaning materials or microfluidic devices. While the static and dynamic properties of deposited drops have been well characterized, a microscopic understanding of the underlying dynamics is missing. In particular, it is unclear how drop dynamics depends on the amount of uncross-linked chains in the brush, because experimental techniques fail to quantify those. Here we use coarse-grained simulations to study droplets moving on a lubricated polymer brush substrate under the influence of an external body force. The simulation model is based on the many body dissipative particle dynamics (MDPD) method and designed to mimic a system of water droplets on poly(dimethylsiloxane) (PDMS) brushes with chemically identical PDMS lubricant. In agreement with experiments, we find a sublinear power law dependence between the external force F and the droplet velocity v, F ∝ vα with α < 1; however, the exponents differ (α ∼ 0.6-0.7 in simulations versus α ∼ 0.25 in experiments). With increasing velocity, the droplets elongate and the receding contact angle decreases, whereas the advancing contact angle remains roughly constant. Analyzing the flow profiles inside the droplet reveals that the droplets do not slide but roll, with vanishing slip at the substrate surface. Surprisingly, adding lubricant has very little effect on the effective friction force between the droplet and the substrate, even though it has a pronounced effect on the size and structure of the wetting ridge, especially above the cloaking transition.

Sliding and rolling of yield stress fluid droplets on highly slippery lubricated surfaces

HYPOTHESIS

Droplets of yield stress fluids (YSFs), i.e. fluids that can flow only if they are subjected to a stress above a critical value and otherwise deform like solids, hardly move on solid surfaces due to their high viscosity. The use of highly slippery lubricated surfaces can shed light on the mobility of YSF droplets, which include everyday soft materials, such as toothpaste or mayonnaise, and biological fluids, such as mucus.

EXPERIMENTS

The spreading and mobility of droplets of aqueous solutions of swollen Carbopol microgels were studied on lubricant infused surfaces. These solutions represent a model system of YSFs. Dynamical phase diagrams were established by varying the concentration of the solutions and the inclination angle of the surfaces.

FINDINGS

Carbopol droplets deposited on lubricated surfaces could move even at low inclination angles. The droplets were found to slide because of the slip of the flowing oil that covered the solid substrate. However, as the descending speed increased, the droplets rolled down. Rolling was favored at high inclinations and low concentrations. A simple criterion based on the ratio between the yield stress of the Carbopol suspensions and the gravitational stress acting on the Carbopol droplets was found to nicely identify the transition between the two regimes.

Efficient Water Collection from Biodesigned and Natural Inclined Surfaces: Influence of Inclination Angle on Atmospheric Water Collection

Water is one of the most important and crucial indicators of sustainable development goals (SDGs) for humans and other living organisms. Water demand has outstripped supply, resulting in shortage on a worldwide scale, particularly in arid regions. This water scarcity has impeded agricultural productivity and other developmental projects with the ongoing global warming and other anthropogenic activities making it more complicated. To address the worldwide water crisis, it is worthwhile to convert atmospheric air to drinking water. Sequel to that, a hydrophobic surface was designed using facile lithography to compare its water harvesting efficiency with a hydrophilic surface at different orientation angles. For the research, the hydrophobic designed surface is called biodesigned material, while the hydrophilic natural surface is a Mangifera indica leaf (MIL). It is against this background that we seek to investigate the most suitable orientation angle good for efficient water harvesting. To that end, a 60° inclination angle is the most efficacious orientation for water collection as it outperformed the 30 and 45° orientation angles. To minimize re-evaporation, absorption, suction, and other environmental challenges that impede efficient collection, atmospheric moisture should be collected immediately from functional surfaces.

Droplet Rolling Dynamics over a Hydrophobic Surface with a Minute Width Channel

Unidirectional and stabilize droplet rolling over hydrophobic surfaces is critical for self-cleaning applications of large areas. Introducing minute size channels on hydrophobic surfaces in the droplet rolling direction can minimize droplet wobbling and enables unidirectional rolling. The droplet rolling behavior over an inclined hydrophobic surface having a minute size channel is investigated. The flow field developed inside the droplet fluid is numerically simulated in a three-dimensional domain pertinent to experimental conditions. Experiments are carried out using a high-speed facility to monitor and evaluate droplet motion over channeled and flat hydrophobic surfaces. The findings revealed that predictions of the droplet translational velocity and those obtained from the experiments are in good agreement. The presence of a minute channel on the hydrophobic surface gives rise to droplet fluid inflection into the minute channel, which in turn modifies the center of mass of the droplet during rolling. This lowers the droplet wobbling height and enables the droplet to roll unidirectionally along the channel length. Enlarging the channel width on the hydrophobic surface increases droplet kinetic energy dissipation while reducing the droplet rolling speed. The complex flow structures formed in the droplet fluid modifies the pressure along the droplet centerline; however, the droplet fluid pressure remains almost the same order as the Laplace pressure in the upper region of a rolling droplet over the channeled hydrophobic surface.

Mapping outcomes of liquid marble collisions

Liquid marbles (LMs) have many promising roles in the ongoing development of microfluidics, microreactors, bioreactors, and unconventional computing. In many of these applications, the coalescence of two LMs is either required or actively discouraged, therefore it is important to study liquid marble collisions and establish parameters which enable the desired collision outcome. Recent reports on LM coalescence have focused on either two mobile LMs colliding, or an accelerating LM hitting a sessile LM with a backstop. A further possible scenario is the impact of a mobile LM against a non-supported static LM. This paper investigates such a collision, using high-speed videography for single-frame analysis. Multiple collisions were undertaken whilst varying the modified Weber number (We*) and offset ratios (X*). Parameter ranges of 1.0 < We* < 1.4 and 0.0 < X* < 0.1, resulted in a coalescence rate of approximately 50%. Whereas, parameter ranges X* > 0.25, and We* < 0.95 or We* > 1.55 resulted in 100% non-coalescence. Additionally, observations of LMs moving above a threshold velocity of 0.6 m s-1 have revealed a new and unusual deformation. Comparisons of the outcome of collisions whilst varying both the LM volume and the powder grain size have also been made, revealing a strong link. The results of this work provide a deeper understanding of LM coalescence, allowing improved control when designing future collision experiments.

Motion of a Droplet on an Anisotropic Microgrooved Surface

We experimentally characterize the sliding angle of water droplets (volume 3.1-22.2 μL) migrating on inclined microgrooved surfaces along the longitudinal and transverse directions of the grooves. The rectangular microgrooves are manufactured on silicon wafers using standard photolithography techniques. We tilt the surface gradually using a rotating stage mechanism until the incipience of the sliding. The droplet migration in the longitudinal and transverse directions to the grooves is recorded using a high-speed camera. For the droplets migrating downward in the transverse direction, the contact line exhibits a "stick-slip" type motion, that is, the advancing contact line is attached to the surface, whereas the receding contact line is detached from the surface. However, no significant change in the relative position of the advancing and receding contact lines is observed in the case of the longitudinal migration of the droplets. The sliding behavior of the droplet in the longitudinal direction is similar to that observed in the case of a smooth surface. The sliding angle in the longitudinal direction of motion is found to be smaller as compared to that in the transverse motion of the droplet. In both longitudinal and transverse migrations, increasing the pitch of the grooves increases the contact angle, which in turn decreases the sliding angle. As the droplet volume is increased, the component of the gravitational force in the direction of inclination increases, which acts to decrease the sliding angle. A theoretical analysis is also conducted to predict the sliding angle of a droplet on microgrooved surfaces. The model predictions agree with the trends observed in our experiments and thus validate the proposed sliding mechanisms in the longitudinal and transverse migrations of the droplet.

Drop Dynamics on Liquid-Infused Surfaces: The Role of the Lubricant Ridge

We employ a free-energy lattice-Boltzmann method to study the dynamics of a ternary fluid system consisting of a liquid drop driven by a body force across a regularly textured substrate, infused by a lubricating liquid. We focus on the case of partial wetting lubricants and observe a rich interplay between contact line pinning and viscous dissipation at the lubricant ridge, which become dominant at large and small apparent angles, respectively. Our numerical investigations further demonstrate that the relative importance of viscous dissipation at the lubricant ridge depends on the drop to lubricant viscosity ratio, as well as on the shape of the wetting ridge.

Shear-Induced Migration of Rigid Particles near an Interface between a Newtonian and a Viscoelastic Fluid

Simulations of rigid particles suspended in two-phase shear flow are presented, where one of the suspending fluids is viscoelastic, whereas the other is Newtonian. The Cahn-Hilliard diffuse-interface model is employed for the fluid-fluid interface, which can naturally describe the interactions between the particle and the interface (e.g., particle adsorption). It is shown that particles can migrate across streamlines of the base flow, which is due to the surface tension of the fluid-fluid interface and a difference in normal stresses between the two fluids. The particle is initially located in the viscoelastic fluid, and its migration is investigated in terms of the Weissenberg number Wi (shear rate times relaxation time) and capillary number Ca (viscous stress over capillary stress). Four regimes of particle migration are observed, which can roughly be described by migration away from the interface for Wi = 0, halted migration toward the interface for low Wi and low Ca, particle adsorption at the interface for high Wi and low Ca, and penetration into the Newtonian fluid for high Wi and high Ca. It is found that the angular velocity of the particle plays a large role in determining the final location of the particle, especially for high Wi. From morphology plots, it is deduced that the different dynamics can be described well by considering a balance in the first-normal stress difference and Laplace pressure. However, it is shown that other parameters, such as the equilibrium contact angle and diffusion of the fluid, are also important in determining the final location of the particle.

Dynamics of a water droplet on a hydrophobic inclined surface: influence of droplet size and surface inclination angle on droplet rolling

An understanding of the dynamic motion of a water droplet is critical to reduce the effort required to remove dust particles from such surfaces.

Lattice Boltzmann simulations of droplet formation in confined channels with thermocapillary flows

Based on mesoscale lattice Boltzmann simulations with the "Shan-Chen" model, we explore the influence of thermocapillarity on the breakup properties of fluid threads in a microfluidic T-junction, where a dispersed phase is injected perpendicularly into a main channel containing a continuous phase, and the latter induces periodic breakup of droplets due to the cross-flowing. Temperature effects are investigated by switching on-off both positive-negative temperature gradients along the main channel direction, thus promoting a different thread dynamics with anticipated-delayed breakup. Numerical simulations are performed at changing the flow rates of both the continuous and dispersed phases, as well as the relative importance of viscous forces, surface tension forces, and thermocapillary stresses. The range of parameters is broad enough to characterize the effects of thermocapillarity on different mechanisms of breakup in the confined T-junction, including the so-called "squeezing" and "dripping" regimes, previously identified in the literature. Some simple scaling arguments are proposed to rationalize the observed behavior, and to provide quantitative guidelines on how to predict the droplet size after breakup.

No-Loss Transportation of Water Droplets by Patterning a Desired Hydrophobic Path on a Superhydrophobic Surface

The directional transportation of droplets on solid surfaces is essential in a wide range of engineering applications. It is convenient to guide liquid droplets in a given direction by utilizing the gradient of wettability, by which the binding forces can be produced. In contrast to the mass-loss transportation of a droplet moving along hydrophilic paths on a horizontal superhydrophobic surface, we present no-loss transportation by fabricating a hydrophobic path on the same surface under tangential wind. In experimental exploration and theoretical analysis, the conditions of no-loss transportation of a droplet are mainly considered. We demonstrate that the lower (or upper) critical wind velocity, under which the droplet starts on the path (or is derailed from the path), is determined by the width of the path, the length of the contact area in the direction parallel to the path, the drift angle between the path and the wind direction, and the surface wettability of the pattern. Meanwhile, the no-loss transportation of water droplets along the desired path zigzagging on a superhydrophobic surface can be achieved steadily under appropriate conditions. We anticipate that such robust no-loss transportation will find an extensive range of applications.

Universal evolution of a viscous-capillary spreading drop

The rate of spreading or retraction of a drop on a flat substrate is determined through a balance of surface tension and hydrodynamic flow. While asymptotic regimes are known, no general rate equation has hitherto been available. Here, we revisit this classic problem, in a regime governed by capillary and viscous forces, by performing an exhaustive numerical study of drop evolution as a function of the contact angle with the substrate. Our study reveals a universal evolution of the drop radius parameterised only by the substrate wettability. Two limits of this evolution recover the familiar exponential and algebraic regimes. Our results show quantitative comparison with the evolution derived from lubrication theory, indicating that dissipation at the contact line is the key determinant in drop evolution. Our work, both numerical and theoretical, provides a foundation for studying the full temporal dynamics of droplet evolution under the influence of external fields and thermal fluctuations, which are of importance in nanofluidics.

Effects of viscoelasticity on droplet dynamics and break-up in microfluidic T-Junctions: a lattice Boltzmann study

The effects of viscoelasticity on the dynamics and break-up of fluid threads in microfluidic T-junctions are investigated using numerical simulations of dilute polymer solutions at changing the Capillary number (Ca), i.e. at changing the balance between the viscous forces and the surface tension at the interface, up to Ca ≈ 3×10(-2). A Navier-Stokes (NS) description of the solvent based on the lattice Boltzmann models (LBM) is here coupled to constitutive equations for finite extensible non-linear elastic dumbbells with the closure proposed by Peterlin (FENE-P model). We present the results of three-dimensional simulations in a range of Ca which is broad enough to characterize all the three characteristic mechanisms of break-up in the confined T-junction, i.e. squeezing, dripping and jetting regimes. The various model parameters of the FENE-P constitutive equations, including the polymer relaxation time τP and the finite extensibility parameter L2, are changed to provide quantitative details on how the dynamics and break-up properties are affected by viscoelasticity. We will analyze cases with Droplet Viscoelasticity (DV), where viscoelastic properties are confined in the dispersed (d) phase, as well as cases with Matrix Viscoelasticity (MV), where viscoelastic properties are confined in the continuous (c) phase. Moderate flow-rate ratios Q ≈ O(1) of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, as the flow driving the break-up process upstream of the emerging thread can be sensibly perturbed by the polymer stresses.

A lattice Boltzmann study of the effects of viscoelasticity on droplet formation in microfluidic cross-junctions

Based on mesoscale lattice Boltzmann (LB) numerical simulations, we investigate the effects of viscoelasticity on the break-up of liquid threads in microfluidic cross-junctions, where droplets are formed by focusing a liquid thread of a dispersed (d) phase into another co-flowing continuous (c) immiscible phase. Working at small Capillary numbers, we investigate the effects of non-Newtonian phases in the transition from droplet formation at the cross-junction (DCJ) to droplet formation downstream of the cross-junction (DC) (Liu and Zhang, Phys. Fluids. 23, 082101 (2011)). We will analyze cases with Droplet Viscoelasticity (DV), where viscoelastic properties are confined in the dispersed phase, as well as cases with Matrix Viscoelasticity (MV), where viscoelastic properties are confined in the continuous phase. Moderate flow-rate ratios Q≈O(1) of the two phases are considered in the present study. Overall, we find that the effects are more pronounced with MV, where viscoelasticity is found to influence the break-up point of the threads, which moves closer to the cross-junction and stabilizes. This is attributed to an increase of the polymer feedback stress forming in the corner flows, where the side channels of the device meet the main channel. Quantitative predictions on the break-up point of the threads are provided as a function of the Deborah number, i.e., the dimensionless number measuring the importance of viscoelasticity with respect to Capillary forces.

Motility of active fluid drops on surfaces

Drops of active liquid crystal have recently shown the ability to self-propel, which was associated with topological defects in the orientation of active filaments [Sanchez et al., Nature 491, 431 (2013)]. Here, we study the onset and different aspects of motility of a three-dimensional drop of active fluid on a planar surface. We analyze theoretically how motility is affected by orientation profiles with defects of various types and locations, by the shape of the drop, and by surface friction at the substrate. In the scope of a thin drop approximation, we derive exact expressions for the flow in the drop that is generated by a given orientation profile. The flow has a natural decomposition into terms that depend entirely on the geometrical properties of the orientation profile, i.e., its bend and splay, and a term coupling the orientation to the shape of the drop. We find that asymmetric splay or bend generates a directed bulk flow and enables the drop to move, with maximal speeds achieved when the splay or bend is induced by a topological defect in the interior of the drop. In motile drops the direction and speed of self-propulsion is controlled by friction at the substrate.

Nanodrop on a smooth solid surface with hidden roughness. Density functional theory considerations

A nanodrop of a test fluid placed on a smooth surface of a solid material of nonuniform density which covers a rough solid surface (hidden roughness) is examined, on the basis of the density functional theory (DFT), in the presence of an external perturbative force parallel to the surface. The contact angles which the drop profile makes with the surface at the leading edges of the drop are determined as functions of drop size and perturbative external force. A critical sticking force, defined as the largest value of the perturbative force for which the drop remains at equilibrium, is determined and its dependence on the size of the drop is explained on the basis of the shape of the interaction potential generated by the solid in vicinity of the leading edges of the drop. For even larger values of the perturbative force no drop-like solution of the Euler-Lagrange equation of the DFT was found. The upper bound of the inclination angle of a surface containing a macroscopic drop is estimated on the basis of results obtained for nanodrops and some experimental results are interpreted. The main conclusion is that the hidden roughness has a similar effect on the drop features as the traditionally considered physical and chemical roughnesses.

Wettability effects on contact line dynamics of droplet motion in an inclined channel

The present work discusses the implications of wall wettability and inclination of the surface on droplet dynamics. In this work, mesoscopic illustration of droplet dynamics in a duct with different inclination angles, based on the two-phase lattice Boltzmann model, is reported. Temporal evolution of wetted length, wetted area, and maximum height of the droplet for surfaces with different inclination angles and wettabilities is furnished in detail in order to elucidate the droplet displacement dynamics. It has been observed that the effect of inclination of the surface on droplet dynamics is more pronounced on a hydrophobic surface as compared to a hydrophilic surface. The time evolution of height and contact line motion of the droplet shows that higher angle of inclination of substrate affects the dynamics strongly irrespective of wettability.

Simulating droplet motion on virtual leaf surfaces

A curvilinear thin film model is used to simulate the motion of droplets on a virtual leaf surface, with a view to better understand the retention of agricultural sprays on plants. The governing model, adapted from Roy et al. (2002 J. Fluid Mech. 454, 235-261 (doi:10.1017/S0022112001007133)) with the addition of a disjoining pressure term, describes the gravity- and curvature-driven flow of a small droplet on a complex substrate: a cotton leaf reconstructed from digitized scan data. Coalescence is the key mechanism behind spray coating of foliage, and our simulations demonstrate that various experimentally observed coalescence behaviours can be reproduced qualitatively. By varying the contact angle over the domain, we also demonstrate that the presence of a chemical defect can act as an obstacle to the droplet's path, causing break-up. In simulations on the virtual leaf, it is found that the movement of a typical spray size droplet is driven almost exclusively by substrate curvature gradients. It is not until droplet mass is sufficiently increased via coalescence that gravity becomes the dominating force.

Combing of genomic DNA from droplets containing picograms of material

Deposition of linear DNA molecules is a critical step in many single-molecule genomic approaches including DNA mapping, fiber-FISH, and several emerging sequencing technologies. In the ideal situation, the DNA that is deposited for these experiments is absolutely linear and uniformly stretched, thereby enabling accurate distance measurements. However, this is rarely the case, and furthermore, current approaches for the capture and linearization of DNA on a surface tend to require complex surface preparation and large amounts of starting material to achieve genomic-scale mapping. This makes them technically demanding and prevents their application in emerging fields of genomics, such as single-cell based analyses. Here we describe a simple and extremely efficient approach to the deposition and linearization of genomic DNA molecules. We employ droplets containing as little as tens of picograms of material and simply drag them, using a pipet tip, over a polymer-coated coverslip. In this report we highlight one particular polymer, Zeonex, which is remarkably efficient at capturing DNA. We characterize the method of DNA capture on the Zeonex surface and find that the use of droplets greatly facilitates the efficient deposition of DNA. This is the result of a circulating flow in the droplet that maintains a high DNA concentration at the interface of the surface/solution. Overall, our approach provides an accessible route to the study of genomic structural variation from samples containing no more than a handful of cells.

Morphological evolution of domains in spinodal decomposition

Domain growth in spinodal decomposition is usually described by a single time-evolving length scale. We show that the evolution of morphology of domains is nonmonotonic. The domains elongate rapidly at first and then, with the help of hydrodynamics, return to a more circular shape. The initial elongation phase does not alter with hydrodynamics. A small deviation from critical composition changes the morphology dramatically.

A finite-element method model for droplets moving down a hydrophobic surface

We set up a 2D computational Finite-Element Method (FEM) model describing the initial descent of a droplet down an inclined hydrophobic substrate. We solve the full Navier-Stokes equations inside the drop domain, and use the arbitrary Lagrangian-Eulerian method to keep track of the droplet surface. The contact angle is included by using the Frennet-Serret equations. We investigate the behaviour of the drop velocity as a function of the slip length and compare with experimental results. Furthermore, we quantify the energy associated with centre-of-mass translation and internal fluid motion, and we also compute the local dissipation of energy inside the drop. The model predicts trajectories for tracer particles deposited inside the drop, and satisfactorily describes the sliding motion of steadily accelerating droplets. The model can be used for determining a characteristic slip parameter, associated with slip lengths and drag reduction for hydrophobic surfaces.

Efficient sample preparation from complex biological samples using a sliding lid for immobilized droplet extractions

Sample preparation is a major bottleneck in many biological processes. Paramagnetic particles (PMPs) are a ubiquitous method for isolating analytes of interest from biological samples and are used for their ability to thoroughly sample a solution and be easily collected with a magnet. There are three main methods by which PMPs are used for sample preparation: (1) removal of fluid from the analyte-bound PMPs, (2) removal of analyte-bound PMPs from the solution, and (3) removal of the substrate (with immobilized analyte-bound PMPs). In this paper, we explore the third and least studied method for PMP-based sample preparation using a platform termed Sliding Lid for Immobilized Droplet Extractions (SLIDE). SLIDE leverages principles of surface tension and patterned hydrophobicity to create a simple-to-operate platform for sample isolation (cells, DNA, RNA, protein) and preparation (cell staining) without the need for time-intensive wash steps, use of immiscible fluids, or precise pinning geometries. Compared to other standard isolation protocols using PMPs, SLIDE is able to perform rapid sample preparation with low (0.6%) carryover of contaminants from the original sample. The natural recirculation occurring within the pinned droplets of SLIDE make possible the performance of multistep cell staining protocols within the SLIDE by simply resting the lid over the various sample droplets. SLIDE demonstrates a simple easy to use platform for sample preparation on a range of complex biological samples.

Sliding drops across alternating hydrophobic and hydrophilic stripes

We perform a joint numerical and experimental study to systematically characterize the motion of 30 μl drops of pure water and of ethanol in water solutions, sliding over a periodic array of alternating hydrophobic and hydrophilic stripes with a large wettability contrast and a typical width of hundreds of microns. The fraction of the hydrophobic areas has been varied from about 20% to 80%. The effects of the heterogeneous patterning can be described by a renormalized value of the critical Bond number, i.e., the critical dimensionless force needed to depin the drop before it starts to move. Close to the critical Bond number we observe a jerky motion characterized by an evident stick-slip dynamics. As a result, dissipation is strongly localized in time, and the mean velocity of the drops can easily decrease by an order of magnitude compared to the sliding on the homogeneous surface. Lattice Boltzmann numerical simulations are crucial for disclosing to what extent the sliding dynamics can be deduced from the computed balance of capillary, viscous, and body forces by varying the Bond number, the surface composition, and the liquid viscosity. Beyond the critical Bond number, we characterize both experimentally and numerically the dissipation inside the droplet by studying the relation between the average velocity and the applied volume forces.