A discrete network method based on previous asymptotic analysis for computing fluid flows between confined rough surfaces is proposed. This random heterogeneous geodesic network method could be either applied to surfaces described by a continuous random field or finely discretized on a regular grid. This method tackles the difficult problem of fluid transport between rough surfaces in close contact. We describe the principle of the method as well as detail its numerical implementation and performances. Macroscopic conductances are computed and analyzed far from the geometrical percolation threshold. Numerical results are successfully compared with the effective medium approximation, the application of which is also studied analytically.

In this paper, we present a network method for computing two-phase flows between two rough surfaces with significant contact areas. Low-capillary number drainage is investigated here since one-phase flows have been previously investigated in other contributions. An invasion percolation algorithm is presented for modeling slow displacement of a wetting fluid by a non wetting one between two rough surfaces. Short-correlated Gaussian process is used to model random rough surfaces.The algorithm is based on a network description of the fracture aperture field. The network is constructed from the identification of critical points (saddles and maxima) of the aperture field. The invasion potential is determined from examining drainage process in a flat mini-channel. A direct comparison between numerical prediction and experimental visualizations on an identical geometry has been performed for one realization of an artificial fracture with a moderate fractional contact area of about 0.3. A good agreement is found between predictions and observations.

We study the electrical and hydraulic conductivity percolation exponents in a Gaussian fracture using the method proposed in Plouraboué [Phys. Rev. E 73, 036305 (2006)]. Nonuniversal conductivity percolation exponents are found: they differ from the theoretical predictions for infinite system size for frozen power-law distributions of local conductivities, as with their finite size corrections. In the hydraulic case, we also find that the probability density function of the conductivity follows a power-law distribution near the percolation threshold.