Directional spreading of a viscous droplet on a conical fibre

If a droplet smaller than the capillary length is placed on a substrate with a conical shape, it spreads by itself in the direction of growing fibre radius. We describe this capillary spreading dynamics by developing a lubrication flow approximation on a cone and by using the perturbation method of matched asymptotic expansions. Our results show that the droplet appears to adopt a quasi-static shape and the predictions of the droplet shape and the spreading velocity from the two mathematical models are in excellent agreement. At the contact line regions, a large pressure gradient is generated by the mismatch between the equilibrium contact angle and the apparent contact angle that maintains the viscous flow. It is the conical shape of the substrate that breaks the front/rear droplet symmetry in terms of the apparent contact angle, which is larger at the thicker part of the cone than at its thinner part. Consequently, the droplet is predicted to move from the cone tip to its base, consistent with experimental observations.


Introduction
A spherical droplet that comes in contact with a solid substrate will change its shape in order to minimize its total surface energy by generating a spreading motion. Droplet spreading on flat substrates has been widely studied and is quite well understood (Tanner 1979;Hocking 1983;Cox 1986;Chen & Wada 1989;Brenner & Bertozzi 1993;Bonn et al. 2009;Carlson, Bellani & Amberg 2012). If the flat substrate has a constant equilibrium contact angle, the centre of mass of the droplet will not change its position along the substrate. One way to generate a directional droplet spreading is to manipulate the chemical composition or the microscale/nanoscale structure of the substrate to make the equilibrium contact angle vary on the substrate (Brochard 1989;Chaudhury & Whitesides 1992;Dos Santos & Ondarçuhu 1995;Sumino et al. 2005;Moosavi & Mohammadi 2011;Li et al. 2017), where the contact angle difference at the front and back of the droplet generates its motion. An alternative to surface coatings is to instead change the macroscopic shape of the substrate to a geometric structure that breaks the front/rear symmetry, e.g. a cone-like structure. A droplet placed on a fibre with the shape of a cone spontaneously starts to move in the direction of a growing cone radius when the flow is dominated by capillary forces (Bico & Quéré 2002;Lorenceau & Quéré 2004). In fact, the principle of capillary induced self-propelled droplets by tuning the macroscopic geometry of the substrate has been widely exploited by living creatures, where plants  and animals (Zheng et al. 2010;Duprat et al. 2012;Wang et al. 2015) have evolved thin structures that generate droplet motion. Cacti that reside in arid regions have developed conical spines for water collection from humid air, which also transport the water droplets from the spine tip to its base for adsorption . Zheng et al. (2010) showed that a similar principle of directional water collection appears on wetted spider silk, where small water droplets condense at the thinner part (the joint) and move to the centre of the thicker part (knots) where multiple droplets in time coalescence to create a large droplet. It was recently shown by Chen et al. (2018b) that the plant Sarracenia trichome has developed a solution for droplet transport that generates velocities three orders of magnitude larger than that found on the spines of cacti by careful design of its macroscopic geometry and its microscopic structure. The biological transport solution in Sarracenia trichome was mimicked in microchannels as a solution for rapid droplet transport. Insects, on the other hand, are in general interested in getting rid of the unwanted weight of water droplets. Water striders have legs covered with tilted conical setae, which are elastic and hydrophobic, and aid the removal of water droplets (Wang et al. 2015). Controlled droplet motion has a broad industrial relevance for manipulation of chemical reactions, fabrication of materials that can maintain a dry or a wet state or new materials for water harvesting, where recent advances have found inspiration in nature (Park et al. 2016;Chen et al. 2018a,b).
When the droplet size, V 1/3 with V the droplet volume, is smaller than the capillary length ≡(γ /ρg) 1/2 with γ the surface tension coefficient, g the gravitational acceleration and ρ the liquid density, its shape and directional movement on a conical fibre is expected to be generated by the capillary forces. Lorenceau & Quéré (2004) studied this droplet dynamics experimentally, where the motion was rationalized by the aid of a theoretical model. The authors first consider the pressure p cy inside an equilibrium barrel-shaped droplet on a cylindrical fibre of radius R, which was derived by Carroll (1976), and expressed as p cy = (2γ /(R + H)) + p o , where H is the maximal thickness of the droplet and p o is the surrounding pressure in the air. It is then assumed that the pressure p co (x) inside a moving droplet on a conical fibre can be written in the same functional form, but with a replacement of R and H by the cone radius R(x) and the droplet thickness H(x), respectively. The substitution of these constants to variables gives rise to a pressure gradient dp co /dx = −(2γ /([R(x) + H(x)] 2 ))((dR/dx) + (dH/dx)) inside the droplet along the fibre's central axis, x. Despite this pressure gradient model being widely adopted to explain the motion of droplets on conical fibres (Zheng et al. 2010;Li et al. 2013;Li, Wu & Wang 2016;Chen et al. 2018b), it is conceptually not justified. For dynamical situations, the pressure gradient is determined through the coupling with the fluid flow inside the droplet and the interfacial curvature κ from the Laplace equation, i.e. p co − p o = γ κ. In this study, we provide a different physical picture to this phenomenon building on earlier seminal works on moving contact lines (Tanner 1979;Hocking 1983;Cox 1986;Eggers 2004Eggers , 2005b; when a droplet is placed in contact with a conical fibre, it quickly adopts to a quasi-static shape of uniform . Sketch of a droplet on a conical fibre with an angle α. The droplet shape is described by the height h(r, t) from the substrate to the free surface as a function of the radial distance from the vertex along the substrate r and time t. At the contact line positions, i.e. r = r r and r = r a , the free surface intersects with the substrate with an equilibrium (microscopic) contact angle θ e . The free surface deforms significantly in the vicinity of the contact line due to a large viscous stress. At the droplet scale, the contact angle appears as the apparent contact angles, θ r at the thinner part of the cone (receding contact line) and θ a at thicker part of the cone (advancing contact line).
pressure at the droplet scale. The conical geometry breaks the front/rear symmetry of the droplet shape. At the thicker part of the fibre, the apparent contact angle is larger than the equilibrium (microscopic) contact angle, thus a flow is generated and the contact line advances. The fluid recedes at the thinner part of the fibre, where the apparent contact angle is smaller than the equilibrium contact angle. Hence the droplet is expected to move from the tip to the base of a conical fibre, see figure 1. The pressure gradient in the bulk of the droplet only acts as a correction term for determining the droplet shape and plays a minor role on the dynamics. We start our study by considering a viscous droplet of dynamic viscosity η that moves on a fibre when inertia can be neglected, i.e. Reynolds number, Re ≡ ((ρUV 1/3 )/η) 1 with the droplet size V 1/3 (γ /ρg) 1/2 . The motion is dominated by the capillary force but hindered by viscous friction, i.e. the capillary number Ca ≡ (ηU/γ ) 1. A characteristic feature of slowly spreading viscous droplets is that it maintains a quasi-static shape during wetting (Tanner 1979;Bonn et al. 2009), where viscous stresses are predominantly located in the vicinity of the contact line and balanced by the capillary stress through interface deformations. Together, this makes the problem well suited to be studied by the perturbation method of matched asymptotic expansions, which has been widely used to describe viscous spreading of droplets on flat substrates (Tanner 1979;Hocking & Rivers 1982;Wilson 1982;Hocking 1983;Cox 1986;Eggers 2004Eggers , 2005bPismen & Thiele 2006;Savva & Kalliadasis 2009;Snoeijer & Eggers 2010;Chan, Gueudré & Snoeijer 2011). We develop here a similar approach, where we consider the viscous spreading of an axisymmetric droplet on a conical fibre by combining the lubrication theory and a perturbation method of asymptotic matching, to describe the directional droplet motion.

Mathematical formulation
We consider a droplet in contact with a solid conical fibre with an angle α as shown in figure 1. We assume the shape of the droplet is symmetric around the central axis of the cone. The droplet shape is described by the height h(r, t) from the substrate to the free surface as a function of the radial distance from the vertex along the substrate r and time t. As the droplet spontaneously starts to move on the cone, an incompressible viscous flow is generated. Since Re 1, the flow inside the droplet is described by Stokes equations (2.1) and the continuity equation reads where u is the velocity and p is the pressure. Moreover, we only consider small droplets, with a shape unaffected by gravity, i.e. the Bond number Bo ≡ ρgV 2/3 /γ 1. We neglect any influence of the air surrounding the droplet as its viscosity is orders of magnitude smaller than the liquid viscosity.
To describe the fluid flow and the droplet motion, equations (2.1)-(2.2) need to be accompanied by several boundary conditions. At the free surface, the tangential stress is zero as we neglect the viscous effects in the air. The normal stress σ f n is described by the Young-Laplace law where κ is the curvature of the interface. At the wetted substrate, the normal velocity is zero and we assume a tangential velocity u s t described by the Navier slip condition (Lauga, Brenner & Stone 2008) where σ s t is the shear stress parallel to the substrate and λ is the slip length. Slippage of fluid along the substrate is well known to regularize the viscous stress singularity at the moving contact line where the liquid, air and solid phase intersect. The slip length has been measured to be in a range of a few nanometres or less for simple fluids (Lauga et al. 2008). Other models, such as the diffuse interface model (Qian, Wang & Sheng 2006; Carlson, Do-Quang & Amberg 2011), a precursor film (Eggers 2005a) and the molecular-kinetic theory (Blake 2006) have been proposed to tackle the hydrodynamic singularity at the moving contact line. Nevertheless, like the slip length, these models introduce a characteristic length that is typically at the nanoscale, several orders of magnitude smaller than a droplet size of a few millimetres.
We also need to specify the slope of the free surface at the contact line. We assume that molecules at the contact lines quickly redistribute so that an equilibrium angle θ e is achieved and given by Young's law cos θ e = (γ SL − γ SV )/γ , and is independent of the contact line velocity, where γ SL and γ SV are, respectively, the liquid/solid and solid/air surface tension coefficients, as used in previous studies (Cox 1986;Eggers 2005b;Duez et al. 2007). The justification of this assumption is beyond the extent of our hydrodynamic model, and could be solved by other modelling approaches such as the diffuse interface model (Qian et al. 2006;Carlson et al. 2011) and molecular dynamics simulations (Johansson, Carlson & Hess 2015). Moreover, the surface of the fibre is assumed to be chemically homogeneous and smooth, allowing us to neglect any contact angle hysteresis.
Directional spreading of a viscous droplet on a cone 894 A26-5 2.1. Lubrication approximation on a cone For polar angles θ 1 and an equilibrium contact angle θ e 1, the flow inside the droplet is primarily in the radial direction and the droplet is fairly flat. By using these approximations, the Stokes equations (2.1) simplifies to the lubrication equations here given in spherical coordinates (a detailed derivation is given in A.1), where u = u(r, θ ) is the radial velocity. The boundary conditions are Equation (2.6) implies that p is independent of θ . Solving (2.5) and (2.6) with the boundary conditions (2.7) gives us the velocity The dynamics of the droplet's interface, i.e. θ = α + h(r, t)/r, is obtained by imposing mass conservation of the liquid (2.9) The characteristic velocity scale U in the radial direction is much smaller than the capillary velocity γ /η as Ca 1, see A.2 for a further description. Hence the free surface relaxes much faster than the motion of the droplet. We then expect that the spreading will be quasi-steady, and the entire droplet moves at a contact line velocity u cl . In the frame of the moving droplet, the droplet shape is stationary for a small time increment, i.e. ∂h/∂t = 0, and the velocity inside the liquid is u − u cl , hence we can reduce the time-dependent lubrication equation (2.9) to a stationary form In addition we have imposed a zero flux condition at the contact line. To evaluate (2.10), we substitute (2.3) and (2.8) with the normal stress σ f n = −p into (2.10) and get ∂κ ∂r = Ca F(h, r, α, λ) , (2.11) with Ca = ηu cl /γ and T. S. Chan, F. Yang and A. Carlson whereφ = h/(αr) andλ = λ/(αr). In the limit of the film thickness being much smaller than the cone radius, i.e.φ 1, F ≈ h(h + 3λ)/3. Hence (2.11) is reduced to the standard two-dimensional steady lubrication equation ∂κ/∂r = 3Ca/(h(h + 3λ)). The curvature of the free surface κ is expressed as with () ≡ ∂()/∂r. We note the second term of the curvature is derived by using a rotation matrix with an angle α 1, see the derivation in A.3. We use (2.13) as a description of the curvature as we will see in the following that the droplet shape away from the contact line is determined by κ = constant.
The boundary conditions for h(r) at the receding contact line r = r r read (2.14b) and at the advancing contact line r = r a read (2.15b) In the following, all the lengths are rescaled by V 1/3 with the volume V given by (2.16) For simplicity of notation, we keep the same symbols for all rescaled quantities, i.e. h, r, λ and κ. Note that (2.11)-(2.15) have the same forms after rescaling and the model parameters that dictate the dynamics are the cone angle α, the equilibrium contact angle θ e and the slip length λ. The droplet profile h(r) and the capillary number Ca are determined by solving (2.11) with the boundary conditions (2.14) and (2.15) by using the shooting method (Press et al. 2007).

Asymptotic analysis
We now turn to a description based on the method of matched asymptotic expansions. The droplet size and the slip length differ by several orders of magnitude and we expect that the governing forces are different at these two length scales. In the vicinity of the contact line, denoted as the inner region, the flow is maintained by the balance of capillarity and the viscous force. Away from the contact line, the droplet is considered to be quasi-static with a shape only determined by the capillary force, denoted as the outer region. In the following, we summarize our analysis, which builds on the work by Eggers (2005b) and this matching procedure has also been justified by a detailed analysis of the behaviour of the solutions in the framework of the lubrication approximation (Sibley, Nold & Kalliadasis 2015).

Inner solutions
In the inner regions, the characteristic length is the slip length λ, which is assumed to be much smaller than the local radius of curvature of the cone, i.e. h ∼ λ rα, hence (2.17) . (2.18) We approximate the curvature of the free surface as κ ≈ h ∼ θ 2 e /λ > 1/(rα), equation (2.18) is then consistent with the lubrication equation for a two-dimensional flow (Batchelor 1967;Oron, Davis & Bankoff 1997). The behaviour of the solution of (2.18) has been discussed in detail by Eggers (2005b). To match the solution in the outer region, only the asymptotic behaviour when h is much larger than λ is required.
At the droplet's trailing edge, i.e. the receding inner region, we define the interfacial profile as h r = h r (x r ), where x r ≡ r − r r is the distance along the substrate from the contact line position, i.e. the profile is determined by the lubrication equation . (2.19) The prime symbol represents the derivative with respect to the independent variable. Equation (2.19) is complemented by the boundary conditions at the substrate where the height is h r (x r = 0) = 0, (2.20) and the profile slope is given by the equilibrium angle (2.21) By following the analysis by Eggers (2005b), the asymptotic behaviour of h r for θ e x r /λ 1 is where Here, Ai is the Airy function and s 1 needs to be determined from the asymptotic matching.
At the advancing droplet front, x a ≡ r a − r is defined as the distance from the contact line along the substrate, and has a profile h a = h a (x a ) described by and complemented with the boundary conditions at the substrate where the height is h a (x a = 0) = 0, (2.25) and the slope is given by the equilibrium angle h a (x a = 0) = θ e . (2.26) The asymptotic behaviour for θ e x a /λ 1 is in a functional form of the Cox-Voinov relation (Eggers 2005b) h a (x a ) 3 = θ 3 e + 9Ca ln(eθ e x a /3λ),

Outer solution
At the length scale of the droplet size, i.e. the outer region, the dominant force is capillarity. Viscous effects appear only as a small correction in the higher-order terms of Ca as we have Ca 1. We define the outer solution ash(r) and expand it in series of Cah (r) =h 0 (r) + Cah 1 (r) + O(Ca 2 ). (2.28) We solve for the leading-order termh 0 (r) from the condition of uniform interfacial curvature κ 0 , which can be expressed by the relation (2.29) Note that this expression is the same as (2.13) with a replacement of h(r) byh 0 (r).
The value of the curvature κ 0 depends on the position of the droplet and is determined together with the relation between the volume V 0 and the profileh 0 of the droplet, that is V 0 = π r a r rh 0 (h 0 + 2αr) dr.
( 2.30) The contact angle of the static profileh 0 (r) with the substrate at the receding part is used to define the receding apparent contact angle θ r , while the advancing droplet edge is defined by the advancing apparent contact angle θ a . Two boundary conditions are required to solve the second-order ordinary differential equation (2.29), where we impose the following conditions at the receding contact line: h 0 (r = r r ) = 0, (2.31) h 0 (r = r r ) = θ r . (2.32) The position of the advancing contact line and the advancing apparent contact angle are given by the conditionsh 0 (r = r a ) = 0, (2.33) h 0 (r = r a ) = −θ a . (2.34) To match to the inner solution at the receding contact line region, only the leadingorder termh 0 (r), i.e. the static profile, is required. The asymptotic behaviour ofh 0 (r) near the contact line is determined by a Taylor expansion: At the advancing droplet front, the first-order correctionh 1 (r) is required in order to match to the logarithmic behaviour of the inner solution in (2.27) (Eggers 2005b). However, there is no analytical solution for the correction term for a cone geometry. Instead of computing the exact expression as done by Eggers (2005b) for a flat substrate, we assume that the outer solution near the advancing contact line has the functional formh where c a is an adjustable parameter but appears in the logarithm and has only a weak effect on the results. The excellent agreement between the results from matching and the numerical results from solving the full lubrication equation suggests that the assumption we have made here is justified.

Matching the inner and the outer solutions
We are now in a position to perform the matching between the asymptotic behaviour of the inner and outer solutions to determine the unknown quantities, Ca, s 1 and θ r . To do this, we see that three matching conditions are required. By comparing the asymptotic behaviour of the inner solution (2.22) and the outer solution (2.35) at the receding region, we find the matching conditions . (2.38) At the advancing contact line region, we match the inner solution (2.27) to the outer solution (2.36). This procedure of matching the cubes of the free surface slope has been shown to be justified for these flows described by lubrication theory (Sibley et al. 2015). The advancing apparent contact angle becomes θ 3 a = θ 3 e + 9Ca ln(c/λ), (2.39) with c ≡ eθ e /3c a that is treated as an adjustable parameter. We fix c = 0.04 when θ e = 0.2 rad by fitting the results of asymptotic matching and lubrication approximation. We also find that c depends on θ e , and we use c = 0.027 when θ e = 0.3 rad and c = 0.007 when θ e = 0.4 rad (see § 3).

2.2.4.
Completely wetting substrate, θ e = 0 rad In the case of a completely wetting droplet, a Landau-Levich-Derjaguin film (Landau & Levich 1942;de Gennes 1985) will be deposited on the substrate at the receding contact line as the droplet moves. We can assume that the receding apparent contact angle θ r is zero for any value of Ca and matching at the receding region is not needed. The advancing apparent contact angle θ a is obtained by solving for the static outer solution governed by (2.29) with θ r = 0 rad. Once θ a is computed, the capillary number Ca is determined by the condition (2.39) at the advancing contact line Ca = θ 3 a 9 ln(c/λ) . (2.40)

Results
Next, we solve the two mathematical models, the lubrication approximation on a cone (LAC) and the asymptotic matching (AM), to predict the directional spreading dynamics of droplets on a conical fibre. There are several physical parameters that can influence the spreading phenomenon, the cone angle α, the slip length λ and the equilibrium contact angle θ e and we show how these parameters influence the dynamics, i.e. the droplet's capillary number Ca and shape. The LAC given by (2.11)-(2.13) is solved by using the shooting method with boundary conditions (2.14), (2.15) and the constant volume condition (2.16). The solution for the AM is obtained by solving the matching conditions (2.37)-(2.39) together with the uniform curvature condition (2.29) and the constant volume condition (2.30) for the static outer solution, which allows us to determine Ca, s 1 , θ r and θ a . The droplet centre of mass x cm is used to quantify its position on the cone and is given by 3.1. Droplet spreading velocity on a cone The droplet shape during directional spreading is shown in figure 2(a), where the droplet moves from left to right for a cone with radius R = αx, i.e. the droplet moves to the thicker part of the cone but decelerates along the way, qualitatively consistent with experimental measurements (Lorenceau & Quéré 2004). We start by determining the droplet velocity along the cone, illustrated by the capillary number Ca as a function of x cm and shown for different slip lengths in figure 2(b). The predictions from the two mathematical models (LAC and AM) are in excellent . The capillary number Ca plotted as a function of the droplet's centre of mass x cm for four different equilibrium contact angles θ e , predicted by the lubrication approximation on a cone and by the asymptotic matching. The three symbols in black indicate the critical capillary numbers above which the film deposition occurs. The slip length λ = 10 −6 and the cone angle α = 0.01 rad (0.57 • ). Note that the completely wetting θ e = 0 rad is a special case as a Landau-Levich-Derjaguin (known as LLD) film is deposited on the substrate surface. The values of c used in AM for θ e = 0, 0.2, 0.3 and 0.4 rad, respectively, are 0.07, 0.04, 0.027 and 0.007. agreement although the slip length is varied by two orders of magnitude. We see that the droplet velocity is not a linear function of the cone radius R and there exists a critical capillary number Ca c above which no solution is predicted by the LAC or the AM. To understand why no solution is found, we show the receding apparent contact angle θ r as predicted by asymptotic matching as a function of Ca in figure 2(c). The receding contact angle approaches zero as the capillary number gets closer to the critical value Ca → Ca c . A vanishing receding contact angle implies that a film is formed at the droplet tail (Cox 1986;Eggers 2004;Snoeijer et al. 2008). Also the slip length λ can influence the spreading, but its influence on the results is weaker than the other parameters in the system and primarily influence the results for larger Ca, in concordance with other moving contact line models (Cox 1986;Eggers 2005b;Snoeijer & Andreotti 2013) as λ appears in the logarithmic term in (2.39).
The droplet wettability, i.e. the equilibrium contact angle θ e also influences the spreading dynamics and we solve the droplet velocity for different θ e , see figure 3. For partially wetting droplets, the capillary number Ca increases with θ e when compared at the same droplet position. The completely wetting droplet θ e = 0 rad is a special case in figure 3, with θ r = 0 rad at any position x cm on the cone, as a Landau-Levich-Derjaguin film formed at the receding contact line (de Gennes 1985). We can then compare a completely wetting droplet with a partially wetting droplet at θ r = 0 rad and at the same position on the conical fibre. Here, the macroscopic droplet shapes would be the same, as well as their advancing apparent contact angle, but the droplet velocity is determined by (2.39) explaining why a completely wetting droplet moves faster than a partially wetting droplet.  FIGURE 4. The capillary number Ca as a function of the cone radius R cm evaluated at the droplet's centre of mass for three different cone angles α by using the approaches of lubrication approximation on a cone and the asymptotic matching. The slip length λ = 10 −6 and the equilibrium contact angle θ e = 0.2 rad.
To illustrate the dependence on the slope of the cone by varying α, we first define R cm as the cone radius at the droplet's centre of mass, (3.2) Equation (3.2) allows us to compare Ca for different α at the same cone thickness R cm , as shown in figure 4. As expected, the droplet moves faster when α is larger for the same R cm . Our results highlight how the droplet velocity can be tuned by the control of the macroscopic geometry of the conical substrate shape.
3.2. Mismatch between the apparent and equilibrium contact angle Since the motion of the droplet is generated by the capillary force, which depends on the interfacial curvature, we would like to know how the droplet shape maintains the viscous flow. In figure 5(a), the interface shape is plotted at a scale similar to the droplet size. We see that the prediction from the two theoretical approaches (LAC and AM) are in excellent agreement, illustrating that the droplet shape is quasi-static. In the contact line regions, it is expected that the viscous stress is balanced by the capillary stress. We show the droplet shape and the slopes of the droplet profile in the receding and the advancing contact line regions as a function of the distance from the contact line positions, see figures 5(b) and 5(c). In the advancing contact line region, we compare the results from the LCA and the AM in the form of the Cox-Voinov relation (Voinov 1976;Cox 1986;Eggers 2005b). The good agreement between LAC and AM suggests that the force balance assumption in the contact line regions is correct. The interface deforms significantly near the contact lines. Large Laplace pressure gradients are generated to maintain the flow inside the two contact line regions. This large interfacial deformation is clearly illustrated in figure 6, where the local slope of the interface is plotted as a function of the distance from the contact line. As we can see from the curve computed by the LAC in figure 6(a), the local interface slope decreases from the equilibrium value (θ e = 0.2 rad) at the receding contact line position to a local minimum within a very small distance. The variation of local angle is due to the difference between the equilibrium contact angle and the receding apparent contact angle (θ r = 0.12 rad) determined from the macroscopic shape, see figure 6(a). In the advancing region, since the apparent contact angle is larger than the equilibrium angle, the interface slope increases with the distance from the contact line position, see figure 6(b). The flow and the droplet motion is generated by the mismatch of the equilibrium and apparent contact angles at the receding and the advancing contact lines. Because of the asymmetry of the conical shape, the apparent contact angle at the thicker part of the cone is larger than that at the thinner part and the droplet moves spontaneously from the thinner part to the thicker part.

Discussion and conclusions
We have shown that the motion of a fairly flat and viscous droplet on a conical fibre is driven by the difference between the equilibrium and the apparent contact angle, which is not the same at the advancing front and the receding tail. Moreover, our analysis shows that the capillary pressure gradient at the scale of the droplet size is very small and scales with Ca 1, which predicts a quasi-static droplet shape. Instead, large pressure gradients and strong interfacial deformations are found in the vicinity of the contact line, along with the dominant part of the viscous stresses. Our findings contrast the model proposed by Lorenceau & Quéré (2004) that is based on a capillary pressure gradient at the drop scale to generate the droplet motion, which for complete wetting is described by the relation where h m ≡ max (h) is the maximum height of the droplet. We plot our results for the completely wetting case θ e = 0 rad and compare the results with the prediction (4.1) for several different cone angles α in figure 7. The first observation we make, is that our predictions do not follow the relation (4.1). However, our results point to the possibility that if measurements are made over a limited range of Ca or h m /((r a − r r )R 3 cm ) the interpretation of the data can easily be misinterpreted to follow (4.1) in particular for larger cone angles.
A mismatch between the equilibrium contact angle and the apparent contact angle is found to dominate the flow and since this difference is not the same at the receding tail and advancing front a directional droplet spreading motion is generated. Although the equilibrium contact angle is constant on the conical fibre, the apparent angle changes due to its macroscopic asymmetric shape. Our model provides a general . The capillary number Ca is computed by using asymptotic matching with (2.40) for a completely wetting droplet as a function of the quantity h m /(r a − r r )R 3 cm . The slip length is λ = 10 −6 . The straight line has a linear slope. description of droplet motion on a conical substrate and can help to better understand droplet motion on slender geometries with varying shapes, which are found in an abundance of engineering systems as well as in nature.
A.3. Derivation of the curvature Here we give the derivation of the expression of curvature (2.13). We first describe the interface as y =ȟ(x), whereȟ is the distance of the free interface from the axis of rotation. The curvature can be expressed as the following standard form for axisymmetric shape: to the order of α.