Finished Theses

Institut für Theoretische Physik I

**Address:**

Institut für Theoretische Physik I

Universität Erlangen-Nürnberg

Staudtstraße 7

91058 Erlangen

Germany

**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

our partners at FAU | |

- Anomalous Transport in Porous Media near the Percolation Transition
- supervised by Gerd Schröder-Turk, Klaus Mecke and Thomas Franosch
- finished 2010-05-05

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
*δr*^{2}(*t*)~*t*^{2/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 *L*_{box}=150, 100 and 50 were multiplied by 10, 10² and 10³ respectively. Anomalous diffusion is compatible with *t*^{2/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 *d*_{w}=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.