Self-affine fronts in self-affine fractures: large and small-scale structure.

The evolution and spatial structure of displacement fronts in fractures with self-affine rough walls are studied by numerical simulations. The fractures are open and the two faces are identical but shifted along their mean plane, either parallel or perpendicular to the flow. An initially flat front advected by the flow is progressively distorted into a self-affine front with a Hurst exponent equal to that of the fracture walls. The lower cutoff of the self-affine regime depends only on the aperture, while the upper cutoff grows with the lateral shift and linearly with the width of the front.

Tracer dispersion in two-dimensional rough fractures.

Tracer diffusion and hydrodynamic dispersion in two-dimensional fractures with self-affine roughness are studied by analytic and numerical methods. Numerical simulations were performed via the lattice-Boltzmann approach, using a boundary condition for tracer particles that improves the accuracy of the method. The reduction in the diffusive transport, due to the fractal geometry of the fracture surfaces, is analyzed for different fracture apertures. In the limit of small aperture fluctuations we derive the correction to the diffusive coefficient in terms of the tortuosity, which accounts for the irregular geometry of the fractures. Dispersion is studied when the two fracture surfaces are simply displaced normally to the mean fracture plane and when there is a lateral shift as well. Numerical results are analyzed using the Lambda parameter, related to convective transport within the fracture, and simple arguments based on lubrication approximation. At very low Péclet number, in the case where fracture surfaces are laterally shifted, we show using several different methods that convective transport reduces dispersion.

Anomalous diffusion with absorption: exact time-dependent solutions

Recently, analytical solutions of a nonlinear Fokker-Planck equation describing anomalous diffusion with an external linear force were found using a nonextensive thermostatistical Ansatz. We have extended these solutions to the case when an homogeneous absorption process is also present. Some peculiar aspects of the interrelation between the deterministic force, the nonlinear diffusion, and the absorption process are discussed.

Permeability of self-affine rough fractures

The permeability of two-dimensional fractures with self-affine fractal roughness is studied via analytic arguments and numerical simulations. The limit where the roughness amplitude is small compared with average fracture aperture is analyzed by a perturbation method, while in the opposite case of narrow aperture, we use heuristic arguments based on lubrication theory. Numerical simulations, using the lattice Boltzmann method, are used to examine the complete range of aperture sizes, and confirm the analytic arguments.

Exact expression for the diffusion propagator in a family of time-dependent anharmonic potentials.

We have obtained the exact expression of the diffusion propagator in the time-dependent anharmonic potential V(x,t)=1 / 2a(t)x(2)+b ln x. The underlying Euclidean metric of the problem allows us to obtain analytical solutions for a whole family of the elastic parameter a(t), exploiting the relation between the path integral representation of the short time propagator and the modified Bessel functions. We have also analyzed the conditions for the appearance of a nonzero flow of particles through the infinite barrier located at the origin (b<0).

Experimental evidence of power-law trapping-time distributions in porous media.

We present experimental results of solute transport in porous samples made of packings of activated carbon porous grains. Exchange experiments, where the tagged solution initially saturating the medium is replaced with the same solution without tracer, are accurately described by macroscopic transport equations. On the other hand, in desorption experiments, where the tagged solution is replaced by water, the solute concentration exhibits a power-law decay for long times, which requires a more detailed, mesoscopic description. We reproduce this behavior within a continuous-time random-walk approach, where the waiting time distribution is related to the desorption isotherm. Results are compatible with a power-law trapping time distribution with divergent first moment, characteristic of anomalous (sub)diffusion.