Towards the Computation of Viscous Flow Resistance of a Liquid Bridge

Towards the Computation of Viscous Flow Resistance of a Liquid Bridge

C. Chen P. Duru M. Prat P. Joseph S. Geoffroy

Université de Toulouse, INPT, UPS, IMFT, France

Université de Toulouse, UPS, INSA, LMDC, France

CNRS, IMFT, Toulouse, France


Université de Toulouse, LAAS, France

| |
| | Citation



Flow within thick liquid films present owing to capillary effects in the pore space is of key importance in many multiphase flow applications in porous media, for example, drying or oil recovery processes. The viscous resistance to the flow is a key parameter for modelling fluid transport in such situations. It is well known for liquid films wetting the corners of tubes of polygonal cross-section. In this latter case, the liquid films shape is simple and can be readily obtained. The situation is much more involved when considering a realistic pore space, as in a packing of spherical particles, for example. In this case, X-ray tomography observations have shown that most of the liquid is confined around contact points between particles at intermediate liquid saturation. Nonetheless, a connectivity of all the liquid bridges throughout the particle packing can exist, allowing liquid transport across the porous medium. The ultimate goal of the present research is to provide the viscous flow resistance for such capillary liquid cluster of complex shape. As a first step in this direction, we present in this paper direct numerical simulation of the Stokes flow in liquid bridges obtained between two cylindrical pillars confined between two hori- zontal plates. The liquid bridge shape is obtained under conditions of hydrostatic equilibrium thanks to the Surface Evolver software. Then simulations of the viscous flow within the bridge are performed using Comsol Multiphysics® Creeping flow solver.


capillary flow, liquid bridge, porous media, viscous resistance


[1] Chauvet, F., Duru, P., Geoffroy, S. & Prat, M., Three periods of drying of a single square capillary tube. Phys. Rev. Lett. 103, pp. 124502, 2009. doi:

[2] Prat, M., On the influence of pore shape, contact angle and film flows on drying of capillary porous media. Int. J. Heat Mass Transfer 50, pp. 1455–1468, 2007. doi:

[3] Yiotis, A.G., Salin, D., Tajer, E.S. & Yortsos, Y.C., Drying in porous media with gravity stabilized fronts: experimental results. Phys. Rev. E 86, pp. 026310, 2012. doi:

[4] Scheel, M., Seemann, R., Brinkmann, M., Michiel, M.D., Sheppard, A. & Herminghaus, S., Liquid distribution and cohesion in wet granular assemblies beyond the capillary bridge regime. J. Phys. Condensed Matter 20, pp. 494236, 2008. doi:

[5] Brakke, K., The Surface Evolver. Exp. Math., 1(2), pp. 141–165, 1992. doi:

[6] Ransohoff, T.C. & Radke, C.J., Laminar flow of a wetting liquid along the corners of a predominantly gas-occupied noncircular pore. J. Colloid Interface Sci. 121(2), pp. 392–401, 1988. doi: