Spreading on viscoelastic solids: are contact angles selected by Neumann's law?

The spreading of liquid drops on soft substrates is extremely slow, owing to strong viscoelastic dissipation inside the solid. A detailed understanding of the spreading dynamics has remained elusive, partly owing to the difficulty in quantifying the strong viscoelastic deformations below the contact line that determine the shape of moving wetting ridges. Here we present direct experimental visualisations of the dynamic wetting ridge using shadowgraphic imaging, complemented with measurements of the liquid contact angle. It is observed that the wetting ridge exhibits a rotation that follows exactly the dynamic liquid contact angle - as was previously hypothesized [Karpitschka et al., Nat. Commun., 2015, 6, 7891]. This experimentally proves that, despite the contact line motion, the wetting ridge is still governed by Neumann's law. Furthermore, our experiments suggest that moving contact lines lead to a variable surface tension of the substrate. We therefore set up a new theory that incorporates the influence of surface strain, for the first time including the so-called Shuttleworth effect into the dynamical theory for soft wetting. It includes a detailed analysis of the boundary conditions at the contact line, complemented by a dissipation analysis, which shows, again, the validity of Neumann's balance.

Dynamic Solid Surface Tension Causes Droplet Pinning and Depinning.

The contact line of a liquid drop on a solid exerts a nanometrically sharp surface traction. This provides an unprecedented tool to study highly localized and dynamic surface deformations of soft polymer networks. One of the outstanding problems in this context is the stick-slip instability, observed above a critical velocity, during which the contact line periodically depins from its own wetting ridge. Time-resolved measurements of the solid deformation are challenging, and the mechanism of dynamical depinning has remained elusive. Here we present direct visualisations of the dynamic wetting ridge formed by water spreading on a PDMS gel. Unexpectedly, it is found that the opening angle of the wetting ridge increases with speed, which cannot be attributed to bulk rheology, but points to a dynamical increase of the solid's surface tensions. From this we derive the criterion for depinning that is confirmed experimentally. Our findings reveal a deep connection between stick-slip processes and newly identified dynamical surface effects.

Thermal fluctuations of an interface near a contact line.

The effect of thermal fluctuations near a contact line of a liquid interface partially wetting an impenetrable substrate is studied analytically and numerically. Promoting both the interface profile and the contact line position to random variables, we explore the equilibrium properties of the corresponding fluctuating contact line problem based on an interfacial Hamiltonian involving a "contact" binding potential. To facilitate an analytical treatment, we consider the case of a one-dimensional interface. The effective boundary condition at the contact line is determined by a dimensionless parameter that encodes the relative importance of thermal energy and substrate energy at the microscopic scale. We find that this parameter controls the transition from a partial wetting to a pseudopartial wetting state, the latter being characterized by a thin prewetting film of fixed thickness. In the partial wetting regime, instead, the profile typically approaches the substrate via an exponentially thinning prewetting film. We show that, independently of the physics at the microscopic scale, Young's angle is recovered sufficiently far from the substrate. The fluctuations of the interface and of the contact line give rise to an effective disjoining pressure, exponentially decreasing with height. Fluctuations therefore provide a regularization of the singular contact forces occurring in the corresponding deterministic problem.

Physical processes causing the formation of penitentes.

Snow penitentes form in sublimation conditions by differential ablation. Here we investigate the physical processes at the initial stage of penitente growth and perform the linear stability analysis of a flat surface submitted to the solar heat flux. We show that these patterns do not simply result from the self-illumination of the surface-a scale-free process-but are primarily controlled by vapor diffusion and heat conduction. The wavelength at which snow penitentes emerge is derived and discussed. We found that it is controlled by aerodynamic mixing of vapor above the ice surface.

Droplets move over viscoelastic substrates by surfing a ridge.

Liquid drops on soft solids generate strong deformations below the contact line, resulting from a balance of capillary and elastic forces. The movement of these drops may cause strong, potentially singular dissipation in the soft solid. Here we show that a drop on a soft substrate moves by surfing a ridge: the initially flat solid surface is deformed into a sharp ridge whose orientation angle depends on the contact line velocity. We measure this angle for water on a silicone gel and develop a theory based on the substrate rheology. We quantitatively recover the dynamic contact angle and provide a mechanism for stick-slip motion when a drop is forced strongly: the contact line depins and slides down the wetting ridge, forming a new one after a transient. We anticipate that our theory will have implications in problems such as self-organization of cell tissues or the design of capillarity-based microrheometers.

Nonlocal effects in sand flows on an inclined plane.

The flow of sand on a rough inclined plane is investigated experimentally. We directly show that a jammed layer of grains spontaneously forms below the avalanche. Its properties and its relation with the rheology of the flowing layer of grains are presented and discussed. In a second part, we study the dynamics of erosion and deposition solitary waves in the domain where they are transversally stable. We characterize their shapes and velocity profiles. We relate their translational velocity to the stopping height and to the mass trapped in the avalanche. Finally, we use the velocity profile to get insight into the rheology very close to the jamming limit.

Moving contact line of a volatile fluid.

Interfacial flows close to a moving contact line are inherently multiscale. The shape of the interface and the flow at meso- and macroscopic scales inherit an apparent interface slope and a regularization length, both named after Voinov, from the microscopic inner region. Here, we solve the inner problem associated with the contact line motion for a volatile fluid at equilibrium with its vapor. The evaporation or condensation flux is then controlled by the dependence of the saturation temperature on interface curvature-the so-called Kelvin effect. We derive the dependencies of the Voinov angle and of the Voinov length as functions of the parameters of the problem. We then identify the conditions under which the Kelvin effect is indeed the mechanism regularizing the contact line motion.

Friction-induced amplification of acoustic waves in a low Mach number granular flow.

We show experimentally that, contrary to ordinary fluids, low Mach number granular pipe flows are linearly unstable toward the emission of acoustic waves. Exponential amplification of the waves propagating in the direction opposite to the flow is directly demonstrated. We relate the observed instability to the friction of the grains on the pipe, which pumps energy from the mean flow to coherent elastic waves. We show that the most amplified wavelength is proportional to the ratio of the tube radius to the friction coefficient.

Booming dune instability.

Sand avalanches flowing down the leeward face of some desert dunes spontaneously produce a loud sound with a characteristic vibrato around a well-defined frequency, a phenomenon called the "song of dunes." Here, we show through theory that a homogenous granular surface flow is linearly unstable towards growing elastic waves when a localized shear band forms at the interface between the avalanche and the static part of the dune. We unravel the nature of the acoustic amplifying mechanism at the origin of this booming instability. The dispersion relation and the shape of the most unstable modes are computed and compared to field measurements.

Evidence of rayleigh-hertz surface waves and shear stiffness anomaly in granular media.

Using the nonlinear dependence of sound propagation speed with pressure, we evidence the anomalous elastic softness of a granular packing in the vicinity of the jamming transition. Under gravity and close to a free surface, the acoustic propagation is only possible through surface modes guided by the stiffness gradient. These Rayleigh-Hertz modes are evidenced in a controlled laboratory experiment. The shape and the dispersion relation of both transverse and sagittal modes are compared to the prediction of nonlinear elasticity including finite size effects. These results allow one to access the elastic properties of the packing under vanishing confining pressure.

Thick films of viscous fluid coating a plate withdrawn from a liquid reservoir.

We consider the deposition of a film of viscous liquid on a flat plate being withdrawn from a bath, experimentally and theoretically. For any plate speed U, there is a range of "thick" film solutions whose thickness scales like U{1/2} for small U. These solutions are realized for a partially wetting liquid, while for a perfectly wetting liquid the classical Landau-Levich-Derjaguin film is observed, whose thickness scales like U{2/3}. The thick film is distinguished from the Landau-Levich-Derjaguin film by a dip in its spatial profile at the transition to the bath. We calculate the phase diagram for the existence of stationary film solutions as well as the film profiles and find excellent agreement with experiment.

Surface elastic waves in granular media under gravity and their relation to booming avalanches.

Due to the nonlinearity of Hertzian contacts, the speed of sound c in granular matter is expected to increase with pressure as P(1/6). A static layer of grains under gravity is thus stratified so that the bulk waves are refracted toward the surface. The reflection at the surface being total, there is a discrete number of modes (both in the sagittal plane and transverse to it) localized close to the free surface. The shape of these modes and the corresponding dispersion relation are investigated in the framework of an elastic description taking into account the main features of granular matter: Nonlinearity between stress and strain and the existence of a yield transition. We show in this context that the surface modes localized at the free surface exhibit a waveguide effect related to the nonlinear Hertz contact. Recent results about the song of dunes are reinterpreted in light of the theoretical results. The predicted propagation speed is compared with measurements performed in the field. Taking into account the finite depth effects, we show that the booming instability threshold can be explained quantitatively by a waveguide cutoff frequency below which no sound can propagate. Therefore, we propose another look at a recent controversy, confirming that the song of dunes can well originate from a coupling between avalanching grains and surface elastic waves once the specificity of surface waves (we baptized Rayleigh-Hertz) is correctly taken into account.

Comment on "Minimal size of a barchan dune".

It is now an accepted fact that the size at which dunes form from a flat sand bed as well as their "minimal size" scales on the flux saturation length. This length is by definition the relaxation length of the slowest mode toward equilibrium transport. The model presented by Parteli, Durán, and Herrmann [Phys. Rev. E 75, 011301 (2007)] predicts that the saturation length decreases to zero as the inverse of the wind shear stress far from the threshold. We first show that their model is not self-consistent: even under large wind, the relaxation rate is limited by grain inertia and thus cannot decrease to zero. A key argument presented by these authors comes from the discussion of the typical dune wavelength on Mars (650 m) on the basis of which they refute the scaling of the dune size with the drag length evidenced by Claudin and Andreotti [Earth Planet. Sci. Lett. 252, 30 (2006)]. They instead propose that Martian dunes, composed of large grains (500 microm), were formed in the past under very strong winds. We emphasize that this saltating grain size, estimated from thermal diffusion measurements, is far from straightforward. Moreover, the microscopic photographs taken by the rovers on Martian Aeolian bedforms show a grain size of 87+/-25 microm together with hematite spherules at millimeter scale. As those so-called "blueberries" cannot be entrained more frequently than a few hours per century, we conclude that the saltating grains on Mars are the small ones, which gives a second strong argument against the model of Parteli.

Flow rule, self-channelization, and levees in unconfined granular flows.

Unconfined granular flows along an inclined plane are investigated experimentally. During a long transient, the flow gets confined by quasistatic banks but still spreads laterally towards a well-defined asymptotic state following a nontrivial process. Far enough from the banks a scaling for the depth averaged velocity is obtained, which extends the one obtained for homogeneous steady flows. Close to jamming it exhibits a crossover towards a nonlocal rheology. We show that the levees, commonly observed along the sides of the deposit upon interruption of the flow, disappear for long flow durations. We demonstrate that the morphology of the deposit builds up during the flow, in the form of an underlying static layer, which can be deduced from surface velocity profiles, by imposing the same flow rule everywhere in the flow.

Hierarchical crack pattern as formed by successive domain divisions. II. From disordered to deterministic behavior.

Hierarchical crack patterns, such as those formed in the glaze of ceramics or in desiccated layers of mud or gel, can be understood as a successive division of two-dimensional domains. We present an experimental study of the division of a single rectangular domain in drying starch and show that the dividing fracture essentially depends on the domain size, rescaled by the thickness of the cracking layer e. Utilizing basic assumptions regarding the conditions of crack nucleation, we show that the experimental results can be directly inferred from the equations of linear elasticity. Finally, we discuss the impact of these results on hierarchical crack patterns, and in particular the existence of a transition from disordered cracks at large scales--the first ones--to a deterministic behavior at small scales--the last cracks.

Boundary conditions in the vicinity of a dynamic contact line: experimental investigation of viscous drops sliding down an inclined plane.

To probe the microscopic balance of forces close to a moving contact line, the boundary conditions around viscous drops sliding down an inclined plane are investigated. At first, the variation of the contact angle as a function of the scale of analysis is discussed. The dynamic contact angle is measured at a scale of 6 mum all around sliding drops for different volumes and speeds. We show that it depends only on the capillary number based on the local liquid velocity, measured by particle tracking. This velocity turns out to be normal to the contact line everywhere. It indirectly proves that, in comparison with the divergence involved in the normal direction, the viscous stress is not balanced by intermolecular forces in the direction tangential to the contact line, so that any motion in this last direction gets damped.

The song of dunes as a wave-particle mode locking.

Singing dunes, which emit a loud sound as they avalanche, constitute a striking and poorly understood natural phenomenon. We show that, on the one hand, avalanches excite elastic waves at the surface of the dune, whose vibration produces the coherent acoustic emission in the air. The amplitude of the sound (approximately 105 dB) saturates exactly when the vibration makes the grains take off the flowing layer. On the other hand, we show that the sound frequency (approximately 100 Hz) is controlled by the shear rate inside the sand avalanche, which for granular matter is equivalent to the mean rate at which grains make collisions. This proves the existence of a feedback of elastic waves on particle motion, leading to a partial synchronization of the avalanching sand grains. It suggests that the song of dunes results from a wave-particle mode locking.

Corridors of barchan dunes: Stability and size selection.

Barchans are crescentic dunes propagating on a solid ground. They form dune fields in the shape of elongated corridors in which the size and spacing between dunes are rather well selected. We show that even very realistic models for solitary dunes do not reproduce these corridors. Instead, two instabilities take place. First, barchans receive a sand flux at their back proportional to their width while the sand escapes only from their horns. Large dunes proportionally capture more sand than they lose, while the situation is reversed for small ones: therefore, solitary dunes cannot remain in a steady state. Second, the propagation speed of dunes decreases with the size of the dune: this leads, through the collision process, to a coarsening of barchan fields. We show that these phenomena are not specific to the model, but result from general and robust mechanisms. The length scales needed for these instabilities to develop are derived and discussed. They turn out to be much smaller than the dune field length. As a conclusion, there should exist further, yet unknown, mechanisms regulating and selecting the size of dunes.

Constitutive property of the local organization of leaf venation networks.

The leaf venation of dicotyledons forms complex patterns. In spite of their large variety of morphologies these patterns have common features. They are formed of a hierarchy of structures, which are connected to form a reticulum. Excellent images of these patterns can be obtained from leaves from which the soft tissues have been removed. A numerical image processing has been developed, specially designed for a quantitative analysis of this type of network. It provides a precise characterization of its geometry. The resulting data reveals a surprising property of reticula's nodes: the angles between vein segments are very well defined and it is shown that they are directly related by the radii of the segments. The relation between radii and angles can be expressed very simply using a phenomenological analogy to mechanics. This local organization principle is universal; all leaf venation patterns studied show the same behavior. The results are compared with physical networks such as fracture arrays or soap froth in terms of hierarchy and reorganization.

Selection of velocity profile and flow depth in granular flows.

The dynamics of a two-dimensional pile constituted by spherical grains organized in parallel layers is investigated theoretically. Only three effects are taken into account in the model: driving by gravity, nonlocal dissipation due to shocks, and trapping of grains by the bumps of the underneath layer. This is sufficient to recover the basic properties of granular avalanches: the transition between static and flowing state is hysteretic; the pile does not flow on the whole height but only in a layer at the surface; the velocity profile inside the flowing layer is approximately linear and is followed by a creep motion in the (quasi) static part. The flow height increases as a function of the pile angle and tends to infinity for a critical angle straight phi(infinity). The dependence of this critical angle with the static angle straight phi(s), the restitution coefficient rho, and the moment of inertia J, is investigated.