Institut für Theoretische Physik I
Universität Erlangen-Nürnberg
Staudtstraße 7
91058 Erlangen

How to find us

Phone: +49-9131-85 28442
Fax: +49-9131-85 28444

Department Physik
Teilbibliothek Physik Web-Opac
UnivIS FAU Erlangen-Nürnberg

Physikalisches Kolloquium
Kolloquium der Theor. Physik
Gruppenseminar der Theorie 1


  • Markus Spanner
  • Anomalous Transport in Porous Media near the Percolation Transition
A realization of the 2D Lorentz model: Obstacles (blue) are randomly placed, trajectories of a point-like ballistic tracer (orange) are simulated in the remaining void space having periodic boundaries.

Transport properties of non-interacting fluid particles in porous materials are studied via molecular dynamics simulation in three-dimensional space. The focus is on systems where the transport is characterized by sub-diffusive behavior, that is the mean square displacement follows a power-law in time with exponent smaller than unity. We show that the mean square displacement of ballistic tracer particles (representing the non-interacting fluid particles) confined to the void space of quenched equilibrium hard sphere configurations follows δr2(t)~t2/z with z=6.25 when the obstacle radius is increased so that the void space is near its percolation transition. This exponent is the same as the exponent observed for the Lorentz model, an even simpler model for transport in porous materials, where a point-like tracer moves through an array of uncorrelated quenched spherical obstacles, which was studied previously in great detail. This shows that the exponent z is not immediately affected by structural correlations of the porous medium, at least for these two stochastic models.

Log-log plot of mean square displacement in cherry pit systems of different system sizes Lbox with hard cores packing fraction η=0.25 and obstacle radius <σc>. Hatched areas indicate the range of MSD values for obstacle sizes σ=<σcδσc. Lines for Lbox=150, 100 and 50 were multiplied by 10, 10² and 10³ respectively. Anomalous diffusion is compatible with t2/z, z=6.25 for all system sizes, larger systems stay longer in the sub-diffusive domain.

Furthermore we have performed simulations of particles confined to the `infinite' cluster, i.e. the fraction of void space that percolates through the system of obstacles. The computational challenge is to identify the percolating cluster by a suitable Voronoi tessellation before the trajectories are calculated by molecular dynamics simulations. We find that anomalous transport at the percolation threshold is still observed to obeying sub-diffusion, but with a different exponent dw=4.81 known as the walk dimension of the system in the context of random walks in lattice percolation. Again our result is independent of the obstacle correlations which is shown by simulations for hard sphere systems for a wide range of densities and for the Lorentz-model.