We present here the capillary cohesion resulting from a liquid bridge between two particles. The bridge between two particles of different sizes takes a complex shape as illustrated in figure 1. R1 and R2 are the particle radii, ζ1 and ζ2 are the filling angles (corresponding to the wetted surface of the particles), θ is the wetting angle, and δn is the distance between particles. The axis x coincides with the axis of rotation of the liquid bridge, and the coordinate y describes the profile of the meridian of the bridge as a function of x. The radius of the bridge is denoted y0. The coordinates x of the three-phase contact line (also called triple line), i.e. the line defining the solid-liquid-gas interface, are denoted by xc1 and xc2 respectively for particles 1 and 2.
In the following, we assume that
The pressure difference Δp = pgaz -pliquide (also called capillary suction) through the liquid-gas interface is related to the curvature of the liquid bridge and the surface tension of the liquid σ by the Young-Laplace equation (Hotta et al. , Soulié et al. ):
The volume of the bridge is given by:
The interparticle distance δn can be expressed in the coordinates xc1 and xc2:
The capillary force can be calculated at the gorge (Hotta et al. , Lian et al. , Mikami et al. ):
The relationship between the capillary force and the shape of the liquid bridge is described by a system of nonlinear equations (1) to (4). However, in discrete numerical simulations an explicit expression of fn is needed as a function of interparticle distance δn and volume of liquid bridge V for any set of parameters R1, R2, σ and θ.
Numerical resolutions of the system of equations have been proposed by several authors (Soulie et al. , Soulié et al. , Scholtes et al. [2009b,a]). The major drawback of these approaches is that they do not provide a physically interpretable form of the capillary law. Moreover, they are not very efficient in terms of computation time. Richefeu it et al. (Richefeu et al. , Richefeu et al.) proposed a simple expression for the capillary force (figure 2b):
where κ = 2πσ cosθ (Willett et al. , Bocquet et al. , Herminghaus ), and λ = λ(V,R1,R2) expresses a length scale that governs the exponential decay of the capillary force as a function of the gap (Richefeu et al.). This explicit expression gives an excellent approximation of the capillary force obtained by the integration of the Young-Laplace equation; see figure 2b.
To complete this expression, a failure criterion of the liquid bridge is required. This failure is governed by energy considerations. It corresponds to a minimal liquid-gas surface interface (Pepin et al. ). Using a numerical solution of the Young-Laplace equation, Erle et al. (Erle et al. ) and De Bisschop and Rigole (Bisschop and Rigole ) have proposed two very close solutions for the failure distance. They proposed an empirical criterion based on the filling angle and the radius of the gorge to evaluate the distance of separation. From similar considerations, Lian et al. (Lian et al. ) proposed the following relationship between the debonding distance Drupt, the volume of liquid bridge V and the wetting angle θ:
This failure criterion has been used in various numerical studies (Richefeu et al. , Soulie et al. ). Notice that the debonding distance is not the same as the distance at which the liquid bridge reforms. This hysteresis phenomenon has been discussed in the literature (Pepin et al. , Soulié ) and can be taken into account. A very common solution is to neglect the volume of the adsorbed liquid to the particle surface so that rejoining occurs at contact.
The distribution of liquid bridges in a granular medium is poorly studied (Fournier et al. , Kohonen et al. ). It plays an important role in the force transmission (Richefeu et al. , Radjai and Richefeu ) and in hydro-texturing phenomena (Rondet et al. [2009b,a]). The study of the distribution of the liquid is all the more difficult that the samples are polydisperse and that the configuration is three-dimensional. In the case of polydisperse media, the volume of all liquid bridges cannot be the same. A common way to allocate the water to the grains is to attribute to each liquid bridge a volume of liquid corresponding to a fraction of the total volume of liquid. This distribution cannot be consistent with the thermodynamic equilibrium between gas, solid and liquid phases since it implies a constant pressure in the liquid phase (the gas phase is percolating). For this reason, it is necessary to take into account the capillary suction in the formulation of local laws of cohesion.
Finally, it is sometimes necessary to take into account the viscosity of the liquid. The normal component fnd of the viscous force can be obtained from (Adams and Perchard ):
where μ is the dynamic viscosity of the liquid and R⋆ is the reduced radius. For the tangential component ftd, it is possible to use the following approximation (Goldman et al. ) which is valid for small values of δn and for small filling angles ζ1 and ζ2:
F. R. E. De Bisschop and W. J. L. Rigole. A physical model for liquid capillary bridges between adsorptive solid spheres: the nodoid of plateau. Journal of Colloid and Interface Science, 88:117–128, 1982.
Z. Fournier, D. Gerimichalos, S. Herminghaus, M.M. Kohonen, F. Mugele, M. Scheel, M. Schulz, B. Schulz, C. Schier, R. Seemann, and A. Shudelny. Mechanical properties of wet granular materials. Applied Physics: Condensed Matter, 17:S477–S502, 2005.
X. Pepin, D. Rossetti, S. M. Iveson, and S. J. R. Simons. Modeling the evolution and rupture of pendular liquid bridges in the presence of large wetting hysteresis. Journal of Colloid and Interface Science, 232:289–297, 2000.
X. Pepin, S. J. R. Simons, S. Blanchon, D. Rossetti, and G. Couarraze. Hardness of moist agglomerates in relation to interparticle friction, granule liquid content and nature. Powder Technology, 117:123–138, 2001.
F. Radjai and V. Richefeu. Bond anisotropy and cohesion of wet granular materials. Philosophical Transactions of the Royal Society A-mathematical Physical and Engineering Sciences, 367(1909):5123–5138, 2009.
E. Rondet, M. Rundgsiyopas, T. Ruiz, M. Delalonde, and J. P. Desfours. Hydrotextural description of an unsaturated humid granular media: Application for kneading, packing and drying operations. Kona-powder and Particle, 27:174–185, 2009b.
L. Scholtes, P. Y. Hicher, F. Nicot, B. Chareyre, and F. Darve. On the capillary stress tensor in wet granular materials. International Journal For Numerical and Analytical Methods In Geomechanics, 33(10):1289–1313, 2009b.
F. Soulié, F. Cherblanc, M.S. El Youssoufi, and C. Saix. Influence of liquid bridges on the mechanical behaviour of polydisperse granular materials. International Journal for Numerical and Analytical Methods in Geomechanics, 30(3):213–228, 2006. URL http://dx.doi.org/10.1002/nag.476.