Generation of randomly distributed geometrically complex porous media samples. Hierarchical data structure and fast multigrid-like solver. Computation on high-performance computing clusters (more than 100.000 processes). Multi-scale approach based on coupling of Darcy and Navier-Stokes equations. Interactive data exploration during the runtime of the simulation.