# [

## Abstract

The spreading of an incompressible viscous liquid over an isotropic homogeneous unsaturated porous substrate is considered. It is shown that, unlike the dynamic wetting of an impermeable solid substrate, where the dynamic contact angle has to be specified as a boundary condition in terms of the wetting velocity and other flow characteristics, the ‘effective’ dynamic contact angle on an unsaturated porous substrate is completely determined by the requirement of existence of a solution, i.e. the absence of a nonintegrable singularity in the spreading fluid’s pressure at the ‘effective’ contact line. The obtained velocity dependence of the ‘effective’ contact angle determines the critical point at which a transition to a different flow regime takes place, where the fluid above the substrate stops spreading whereas the wetting front inside it continues to propagate.

Dynamic contact angle on an unsaturated wettable porous substrate]
Dynamic contact angle of a liquid spreading on an unsaturated
wettable porous substrate
Y.D. Shikhmurzaev and J.E. Sprittles
]
Y\lsU\lsL\lsI\lsI\nsD.\nsS\lsH\lsI\lsK\lsH\lsM\lsU\lsR\lsZ\lsA\lsE\lsV^{1}^{2}

## 1 Introduction

The modelling of the spreading of liquids over porous substrates in the framework of continuum mechanics requires, and is based on, the separation of scales between the ‘macroscopic’ (or ‘Darcy-scale’) and ‘microscopic’ (or ‘pore-scale’) processes (Barenblatt et al., 1990). In the continuum approximation, the description of the spreading phenomenon brings in the notion of an ‘effective’ smooth penetrable solid substrate, which is how the actual porous medium is represented, together with the notions of ‘effective’ contact lines and ‘effective’ contact angles that the free surface of the pure fluid above and the wetting front inside the substrate form with it. These notions, being averages in the sense of mechanics of multiphase media (Whitaker, 1999), are fundamentally different from the concepts of ‘contact line’ and ‘contact angle’ used in the modelling of dynamic wetting on the actual, as opposed to ‘effective’, impermeable solid surfaces. In particular, the notion of a ‘static’ (or ‘equilibrium’) contact angle, central to the modelling of dynamic wetting where it is a measure of wettability of the solid (Ralston et al., 2008), becomes meaningless for the Darcy-scale description of the liquid spreading over a porous medium: as experiments show (e.g. Clarke et al., 2002; Starov et al., 2003; Markicevic et al., 2010), if, say, a drop of liquid is deposited onto an unsaturated wettable porous substrate, the eventual equilibrium state will be the drop completely imbibed into the porous medium, with no liquid left above it and hence no ‘effective’ static (equilibrium) contact angle between the, now non-existent, free surface of the pure liquid and the effective surface of the substrate.

The main theoretical implication of this absence of a meaningful effective static (equilibrium) contact angle is that the effective dynamic contact angle that the free surface of a spreading liquid forms with a porous substrate on the Darcy-scale cannot be regarded as essentially a perturbation of the static contact angle, which is what one invariably finds in all models dealing with the dynamic wetting of impermeable solid surfaces (see, Dussan V., 1979; de Gennes, 1985; Blake, 2006; Shikhmurzaev, 2011, for reviews). Therefore, it becomes important to consider the problem from first principles, without implying a-priori that the case of a porous substrate can be described by adjusting the concepts borrowed from the modelling of dynamic wetting of impermeable solid surfaces, such as the equilibrium contact angle.

In the present paper, we show that, unlike the situation one has in
dynamic wetting, where the ‘microscopic’ dynamic contact angle has
to be specified as an additional boundary condition^{3}

## 2 Problem formulation

Consider an incompressible Newtonian fluid of density and viscosity spreading at a speed over an isotropic homogeneous unsaturated porous substrate characterized by an effective pore size . The gas displaced by the fluid from the substrate and from the inside of the porous matrix is assumed to be ideal and at a constant pressure with respect to which the pressure in the pure fluid and in the fluid inside the porous medium will be measured. In order to be able to model the process in the framework of continuum mechanics, we need a separation of scales between the macroscopic and the pore-scale phenomena, i.e. we have to consider the continuum limit

(1) |

where is the characteristic length scale on which the phenomenon is described. The resulting model will be applicable to experiments if , with the actual value of determining its accuracy. In the 0th approximation in the above limit, one has (see Fig. 1) (a) a macroscopic ‘wetting front’ as a sharp interface separating the saturated porous medium from the unsaturated matrix , (b) an effective ‘contact line’ at which the free surface , confining the domain occupied by the pure fluid, and the wetting front meet, and (c) well-defined ‘contact angles’ and that the free surface and the wetting front form with the ‘effective’ surface of the solid. The reference frame in which the problem will be considered and the directions of unit normals , and to, respectively, , and are shown in the figure. Importantly, in the scheme outlined above we have already made an assumption that the two contact lines, i.e. the contact line CL1 formed with the substrate by the free surface and the contact line CL2 formed with it by the wetting front , coincide. This is always what happens when the fluid is first brought in contact with the substrate, and we will examine what follows from this initial situation and, later, the conditions when the assumption that CL1 and CL2 coincide no longer holds.

We will consider the case of small Reynolds and capillary numbers for the flow of the pure fluid, i.e. the limit and

(2) |

where is the surface tension of the fluid-gas interface. In this limit, to leading order, inertial effects can be neglected, and, as in the case of the dynamic wetting of an impermeable solid surface (Huh & Scriven, 1971; Shikhmurzaev, 1993), from the normal-stress boundary condition on the free surface we have that near the contact line the free surface () is locally planar, so that locally is a wedge region. For simplicity we will also neglect gravity, though its inclusion would not change the main results of the analysis below.

Scaling the pressure and velocity in with and respectively, one has that in the dimensionless pressure and velocity obey the Stokes equations

(3) |

and on the free surface satisfy the standard kinematic and tangential-stress boundary conditions,

(4) |

where is the metric tensor.

In the porous medium, the flow is driven by the capillary pressure in the menisci that collectively form the wetting front, and, for the problem in question, the characteristic pressure and velocity are and , respectively. Using the notation and for the dimensionless pressure and velocity in , in a frame moving with the contact line we have an equation of motion in the form of Darcy’s law,

(5) |

where is the nondimensionalized permittivity of the porous matrix ( is the actual permittivity) and is a unit vector directed along the velocity of the porous substrate. Given that in porous media (Probstein, 1989), we have that and, for the problem to be nontrivial, it is assumed to be finite in the limits (1), (2). In the present context, it is convenient to define by setting .

Since the porosity in a homogeneous matrix is constant, the mass balance equation has the standard form , so that, after substituting (5) into it, one arrives at Laplace’s equation for in :

(6) |

On the wetting front, one has the kinematic condition that the front propagates with the velocity of the fluid, i.e. in the reference frame described above where the process is steady,

(7) |

and needs to specify the dynamic condition on the pressure . To formulate this condition in a general case, one has to consider the modes of motion the menisci go through on the pore scale as the macroscopic wetting front propagates (Shikhmurzaev & Sprittles, 2012; ?). Here we will be considering the simplest case involving the main, ‘wetting’, mode. In this case, for an ‘effective’ pore with a circular cross-section the dimensionless pressure at the front is given by

(8) |

where is the contact angle formed on the pore scale by the representative meniscus with the pore wall (Fig. 1). Unlike the simplified approach pioneered by Washburn (1921), where it is assumed that is constant and equal to the prescribed static contact angle , we need to take into account that, as demonstrated by numerous experiment (e.g. see Ch. 3 of Shikhmurzaev, 2007, for a review), depends on the wetting speed, i.e. to consider the dependence

(9) |

where is the appropriate velocity scale depending on the material parameters of the system (as all velocities above, we have it nondimensionalized using ). The function has to be determined theoretically or empirically. For example, one can apply derived using the theory of dynamic wetting as a process of interface formation (Shikhmurzaev, 2007), which has been shown to reliably describe experimental data, though, in the present context, any function representing the experimentally observed dependence of the form (9) could be used, such as, for instance, the one that comes from the molecular-kinetic theory of wetting (Blake & Haynes, 1969). For more information about the dynamic wetting modelling we refer the reader to a recent review (Shikhmurzaev, 2011).

From the theory of flows with forming interfaces one has

(10) |

where ,

and , , , , , are material constants characterizing the contacting media. Their values for some systems can be found elsewhere (Shikhmurzaev, 2007; Blake & Shikhmurzaev, 2002).

So far, the flow in the pure fluid and in the porous medium were considered separately, and to link them one has to specify three boundary conditions at . One condition that we obviously have on this surface is the continuity of mass flux:

(11) |

For the velocity components parallel to , i.e. for and , a number of boundary conditions have been discussed in the literature (e.g. Saffman, 1971; Jones, 1973; Murdoch & Soliman, 1999; Nield, 2009; Auriault, 2010), following the experiments reported by Beavers & Joseph (1967) and an empirical condition these authors proposed. As noted, for example, by Auriault (2010), these conditions are aiming at capturing the effects of order , i.e. go beyond the classical, i.e. 0th-order, approximation in the continuum limit (1). In the 0th approximation, all these conditions reduce to no-slip for the pure fluid,

(12) |

and it is this condition that we will be using here.

The condition of continuity of pressure on in the dimensionless form yields that

and hence, to leading order in the limits (1) and (2), one has

(13) |

Importantly, unlike the case of an impermeable solid substrate, where one has to specify the dynamic contact angle (as we need to specify on the pore scale), here, for the effective contact angle , we will require only that the flow parameters in the porous medium remain regular at the contact line and that the solution in the pure fluid exists.

## 3 Dynamic contact angle

Consider the asymptotic behaviour of the solution to (3)–(13) in the case of a two-dimensional flow as the distance to the contact line . In the polar coordinates shown in Fig. 1, to leading order, for the pressure in the porous medium one has

where the last condition follows from (7) and (5). The separable solution to this problem is obviously given by

(14) |

so that, using the dynamic boundary condition (8), we obtain that

Then, from (9), where now , we have an equation determining :

(15) |

In order to consider the flow in the pure fluid, we introduce a stream function

(16) |

so that for the leading-order term of the asymptotic expansion of as we have a biharmonic equation

(17) |

together with conditions (4), i.e.

(18) |

on the free surface, and conditions (11) and (12), i.e.

(19) |

on the surface of the solid substrate. In writing down the first of conditions (19) we made use of (14) and integrated along .

The separable solution to the problem (17)–(19) has the form

(20) |

where the constants are given by

(21) |

(22) |

Using (16), we can write down the radial projection of the second equation in (3) in the form

and, after substituting the solution (20), arrive at

Thus, the leading term in the coordinate expansion of the stream function will not give rise to a nonintegrable singularity of pressure if and only if , or, using (21) and (22),

Given the expression (15) for and introducing the dimensionless contact-line speed , we can also write this equation down as

(23) |

This equation specifies in terms of the speed of the contact line with respect to the substrate, the porosity of the matrix and the material properties of the contacting media accumulated in and . The velocity-dependence of for different porosities is shown in Fig. 2. As one can see, this dependence is much steeper than the velocity-dependence of the dynamic contact angle for the dynamic wetting of an impermeable substrate (e.g. see Ch. 3 of Shikhmurzaev, 2007, for a review). In particular, as .

The dynamic contact angle given by (23) ensures integrability of the normal stress on the free surface and hence the existence of a solution. The next term in the asymptotic expansion of has the form and therefore can give rise, at most, to a logarithmic (i.e. integrable) singularity of at the contact line, which does not affect the existence of the solution.

The obtained result has a clear physical meaning. In the case we are considering, the dynamics of imbibition determines the normal to the substrate component of the pure fluid’s velocity independently of the pure fluid’s bulk flow, and, given that the tangential component of the fluid’s velocity on the surface of the substrate is prescribed, as it satisfies the no-slip condition with a known speed of the substrate, we have a moving contact-line problem where on the solid surface both components of velocity are set. In this situation, in the pure fluid the solution exists only if from the family of stream functions described by (20) we choose the one that corresponds to a uniform flow (, ) with the projections of velocity on the normal and tangential to the substrate directions equal to the speed of imbibition and the speed of the solid substrate, respectively.

The solution (23) for a steady spreading of a fluid over an unsaturated substrate ceases to exist when the dimensionless velocity becomes equal to and both and reach . For we have a different regime with the two contact lines CL1 and CL2 no longer coinciding as the wetting front moves ahead (dashed line in Fig. 1), and the pure fluid finds itself on a saturated substrate. Then, in the vicinity of CL1 the pressure satisfies Laplace’s equation together with the boundary conditions

and the local solution with the regular pressure gradient has the form

(24) |

where , are constants. Then, according to (5) and (11), one has as . As a result, if as , then, to leading order as , we will have the ‘classical’ moving contact-line problem, with no imbibition and the no-slip boundary condition on the solid. This problem, as is well known (see, e.g., Ch. 3 of Shikhmurzaev, 2007), has no solution. The only way out of this situation is to conclude that, as CL1 and CL2 separate, the contact line CL1 stops moving. Then, for the pure fluid one simply has a static contact line with the imbibition velocity near it and . It is easy to verify that the solution to this problem exists. As CL1 and CL2 separate and CL2 moves ahead, the imbibition process near CL1 slows down, leading to in (24) and hence , so that no longer has to be equal to and can go down as the imbibition continues.

## 4 Discussion

The described scenario is in agreement with the available experimental data. In experiments on the spreading drops, it has been observed that the regime where the drop base and the saturated area underneath it expand together is followed by the regime where the two contact lines, CL1 and CL2, separate, as drop’s base stops expanding whereas the saturated area continues to grow (Clarke et al., 2002; Starov et al., 2003; Keshav & Basu, 2007; Markicevic et al., 2010). This is usually attributed to the ‘competition’ between spreading and imbibition, and the above analysis shows what this actually means.

Markicevic et al. (2010) report that in their experiment the spreading droplet maintained the shape of a spherical cap throughout the process and that its base stopped expanding when “the droplet is a half of sphere”. In other words, it stopped expanding when the dynamic contact angle became equal to . This is exactly what follows from the result of our analysis.

The dependence of the effective contact angle on the contact-line speed for given by (23) and illustrated in Fig. 2 awaits its experimental verification. The issue here is that, to extract this dependence from experiments with unsteady flows, one has to deal with very short time intervals ( ms according to Markicevic et al. (2010)) with a well-controlled spatial resolution and conditions of the Darcy-scale description satisfied. The problem is made more complicated by the fact that, as one can see in Fig. 2, the velocity-dependence of is very steep, which brings in additional conditions on the temporal resolution of the experiments. On the other hand, experiments with steady flows of the kind commonly used to study dynamic wetting of impermeable solid surface require sizeable substrates with well-reproducible properties. To date, no systematic measurements of with controlled spatial resolution have been reported.

As we have shown, the requirement of the absence of a nonintegrable singularity of pressure in the pure fluid at the contact line, which is equivalent to the requirement of the existence of a solution, uniquely specifies the velocity-dependence of the ‘effective’ dynamic contact angle formed by the free surface and the ‘effective’ surface of the porous substrate. In practical computations, it may be convenient to set (23) as a boundary condition while considering the normal-stress condition on the free surface as the equation determining the free surface shape. One implication of the obtained result is that, if one specifies as an additional boundary condition different from (23), then, for a solution to exist, it will become necessary to use a slip boundary condition instead of the no-slip condition (12) employed here (Davis & Hocking, 2000). As mentioned earlier, this would mean bringing in effects of into the essentially classical, i.e. as , formulation. If the velocity-dependence of different from (23) is imposed together with the no-slip condition (12) for the pure fluid (Alleborn & Raszillier, 2004; Reis Jr. et al., 2004), then there will be no solution to the problem, although this fact can be hidden behind, and as a result masked by, simplifying assumptions made in the process of finding the solution (e.g. lubrication approximation) and the numerical implementation.

From the theoretical viewpoint, it would be interesting to consider the essentially unsteady process of separation of CL1 and CL2, which marks the transition between the two regimes described above. As with every finite-time transition, this is a challenging problem deserving a detailed investigation.

### Footnotes

- E-mail: Y.D.Shikhmurzaev@bham.ac.uk
- E-mail: sprittles@maths.ox.ac.uk
- Here we are not discussing the so-called ‘apparent’ contact angle resulting from the free-surface bending near the contact line. This angle is not part of the mathematical problem formulation; it is an auxiliary concept introduced in some works and in different ways (see, e.g., ?) to interpret experimental results when the measured angle differs from the angle imposed in the problem formulation.

### References

- Alleborn, N. & Raszillier, H. 2004 Spreading and sorption of a droplet on a porous substrate. Chem. Eng. Sci. 59, 2071–2088.
- Auriault, J.-L. 2010 About the Beavers and Joseph boundary condition. Transp. Porous Med. 83, 257–266.
- Barenblatt, G. I., Entov, V. M. & Ryzhik, V. M. 1990 Theory of Fluid Flows through Natural Rocks. Kluwer Acad. Publ., London.
- Beavers, G. S. & Joseph, D. D. 1967 Boundary conditions at a naturally permeable wall. J. Fluid Mech. 30, 197–207.
- Blake, T. D. 2006 The physics of moving wetting lines. J. Colloid Inerf. Sci. 299, 1–13.
- Blake, T. D. & Haynes, J. M. 1969 Kinetics of liquid/liquid displacement. J. Colloid & Interf. Sci. 30, 421–423.
- Blake, T. D. & Shikhmurzaev, Y. D. 2002 Dynamic wetting by liquids of different viscosity. J. Colloid Interf. Sci. 253, 196–202.
- Clarke, A., Blake, T. D., Carruthers, K. & Woodward, A. 2002 Spreading and imbibition of liquid droplets on porous surfaces. Lagmuir 18, 2980–2984.
- Davis, S. H. & Hocking, L. M. 2000 Spreading and imbibition of viscous liquid on a porous base. ii. Phys. Fluids 12, 1646–1655.
- Dussan V., E. B. 1979 On the spreading of liquids on solid surfaces: Static and dynamic contact lines. Annu. Rev. Fluid Mech 11, 371.
- de Gennes, P.-G. 1985 Wetting: statics and dynamics. Rev. Mod. Phys. 57, 827–863.
- Huh, C. & Scriven, L. E. 1971 Hydrodynamic model of steady movement of a solid/liquid/fluid contact line. J. Colloid Interf. Sci. 35, 85–101.
- Jones, I. P. 1973 Low reynolds number flow past a porous spherical shell. Proc. Camb. Philos. Soc. 73, 231–238.
- Keshav, T. R. & Basu, S. 2007 Spreading of liquid droplets on proton exchange membrane of a direct alcohol fuel cell. Chem. Eng. Sci. 62, 7515–7522.
- Markicevic, B., D’Onofrio, T. G. & Navaz, H. K. 2010 On spread extent of sessile droplet into porous medium: Numerical solution and comparison with experiments. Phys. Fluids 22, 012103.
- Murdoch, A. I. & Soliman, A. 1999 On the slip boundary condition for liquid flow over planar porous boundaries. Proc. R. Soc. London A 455, 1315–1340.
- Nield, D. A. 2009 The Beavers-Joseph boundary condition and related matters: A historical and critical note. Transp. Porous Med. 78, 537–540.
- Probstein, R. F. 1989 Physicochemical Hydrodynamics: An Introduction. Butterworth, Stonehaven, MA.
- Ralston, J., Popescu, M. & Sedev, R. 2008 Dynamics of wetting from an experimental point of view. Annu. Rev. Mater. Res. 38, 23–43.
- Reis Jr., N. C., Griffiths, R. F. & Santos, J. M. 2004 Numerical simulation of the impact of liquid droplets on porous surfaces. J. Comput. Phys. 198, 747–770.
- Saffman, P. G. 1971 On the boundary condition at the surface of a porous medium. Stud. Appl. Math. 50, 93–101.
- Shikhmurzaev, Y. D. 1993 The moving contact line on a smooth solid surface. Intl J. Multiphase Flow 19, 589–610.
- Shikhmurzaev, Y. D. 2007 Capillary Flows with Forming Interfaces. Chapman & Hall/CRC, Boca Raton-London-New York.
- Shikhmurzaev, Y. D. 2011 Some dry facts about dynamic wetting. Eur. Phys. J. - Special Topics 197, 47–60.
- Shikhmurzaev, Y. D. & Sprittles, J. E. 2012 Wetting front dynamics in an isotropic porous medium. J. Fluid Mech. 694, 399–407.
- Starov, V. M., Zhdanov, S. A., Kosvintsev, S. R., Sobolev, V. D. & Velarde, M. G. 2003 Spreading of liquid drops over porous substrates. Adv. Colloid & Interf. Sci. 104, 123–158.
- Washburn, E. W. 1921 The dynamics of capillary flow. Phys. Rev. 17, 273–283.
- Whitaker, S. 1999 The Method of Volume Averaging. Kluwer.