Scaling relations in the diffusive infiltration in fractals

Abstract

In a recent work on fluid infiltration in a Hele-Shaw cell with the pore-block geometry of Sierpinski carpets (SCs), the area filled by the invading fluid was shown to scale as F ∼ t n , with n < 1/2, thus providing a macroscopic realization of anomalous diffusion [Filipovitch et al, Water Resour. Res. 52 5167 (2016)]. The results agree with simulations of a diffusion equation with constant pressure at one of the borders of those fractals, but the exponent n is very different from the anomalous exponent ν = 1/D W of single particle diffusion in the same fractals (D W is the random walk dimension). Here we use a scaling approach to show that those exponents are related as n = ν (D F − D B ), where D F and D B are the fractal dimensions of the bulk and of the border from which diffusing particles come, respectively. This relation is supported by accurate numerical estimates in two SCs and in two generalized Menger sponges (MSs), in which we performed simulations of single particle random walks (RWs) with a rigid impermeable border and of a diffusive infiltration model in which that border is permanently filled with diffusing particles. This study includes one MS whose external border is also fractal. The exponent relation is also consistent with the recent simulational and experimental results on fluid infiltration in SCs, and explains the approximate quadratic dependence of n on D F in these fractals. We also show that the mean-square displacement of single particle RWs has log-periodic oscillations, whose periods are similar for fractals with the same scaling factor in the generator (even with different embedding dimensions), which is consistent with the discrete scale invariance scenario. The roughness of a diffusion front defined in the infiltration problem also shows this type of oscillation, which is enhanced in fractals with narrow channels between large lacunas.

Models of infiltration into homogeneous and fractal porous media with localized sources

“…Subsequently, it motivated the proposal of the gradient percolation problem [23,24]. In a recent work, RWI in deterministic fractals was studied, also with the sources at flat boundaries [15], which helped to understand scaling properties of infiltration models.…”

“…6(a) shows I as a function of T ; the linear fit confirms that it has the same scaling as I DE in Eq. (15). The characteristic length of the infiltration [Eq.…”

“…In Ref. [15], a scaling approach was used to relate the infiltration exponent n [Eq. (1)] and the random walk dimension D W in cases where the source was one of the outer boundaries of the fractal.…”

“…1(b). A model for this process in a two-dimensional lattice was proposed by Sapoval, Rosso, and Gouyet [14], and it was recently extended to fractal lattices [15]: the solute is represented by particles that execute random walks in the lattice and that have a constant concentration at an external boundary. Here, this is termed random walk infiltration (RWI).…”

“…Ref. [15] used a scaling approach and numerical simulations to show the relation between the exponent n, the fractal dimensions of the infiltrated medium (D F ) and of the infiltration boundary (D B ), and the random walk dimension in that medium (D W ).…”